Peer Revenue Returns subject to peer review

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Thursday 23 July

Cosmology and Nongalactic Astrophysics

Web-Halo Model Peak-Background Split (WHM-PBS): halo bias as a distribution, not a number

We present the Web-Halo Model Peak-Background Split (WHM-PBS), an analytic theory in which the large-scale bias of a dark-matter halo is inherited from its cosmic-web environment. Building on the Web-Halo Model, we use the Shen et al. moving barriers for ellipsoidal collapse generating the web hierarchy in which every halo sits inside a host filament, itself inside a sheet. Combined with the peak-background split, this picture replaces the deterministic bias-mass relation $b(M_h)$ with the bias of the host environment, averaged over the conditional mass function. As a result, halo bias $b(M_h)$ is no longer a number but a strongly skewed distribution. In this work we make use of this distribution in three different ways: as (i) a physically motivated prior on bias relations, (ii) a prediction on halo stochasticity, and (iii) a framework for assembly bias models. Regarding (i) we find that the density bias relations $b_2(b_1)$ and $b_3(b_1)$ stay tight, while the tidal bias $b_{s^2}(b_1)$ shows significant scatter, as found in $N$-body simulations. Regarding (ii), once including halo exclusion, our model reproduces the super- to sub-Poisson shot-noise trend of Baldauf et al. which we convert into a prior band on the EFT stochasticity amplitude parameters. Finally, regarding (iii) in the density sector it explains the bias-concentration-correlation inversion of Paranjape et al. at the characteristic mass ($M_\mathrm{h}\simeq 1.7\times 10^{13}\,h^{-1} M_\odot$), with no parameter tuned to assembly bias. We apply the resulting priors on synthetic data, demonstrating that the WHM-PBS priors mitigate projection effects arising when all nuisance parameters are varied freely, while fixed priors may fail severely. The code to reproduce our results is publicly available at https://github.com/SamuelBrieden/whmpbs.git.
Read the full abstract

High Energy Astrophysical Phenomena

Spirited Away: Advective Loss of Cosmic Rays into the Milky Way's Circumgalactic Medium Explains the Large Magellanic Cloud's Low Gamma-ray Luminosity

Models of galactic-scale cosmic ray production and transport have successfully reproduced the radio and $\gamma$-ray spectra of many galaxies; however, one notable exception is the Large Magellanic Cloud (LMC), where models that successfully fit other galaxies consistently overestimate its $\gamma$-ray flux. Here we investigate this discrepancy by applying the CRIPTIC cosmic ray transport code to recent magnetohydrodynamic simulations of the LMC's interactions with the Milky Way circumgalactic medium that accurately reproduce observations of the LMC's magnetic field structure. We recover a simulated $\gamma$-ray luminosity that is close to the observed luminosity for a wide range of cosmic ray transport models, while our control simulations of an isolated LMC not interacting with the Milky Way show the same over-prediction problem as previous investigations. Simulations including the CGM interaction yield lower $\gamma$-ray luminosities because in them cosmic rays escape the galaxy primarily via advection, significantly decreasing the emission from collisional processes and changing the dominant $\gamma$-ray emission mechanism from pion decay to inverse Compton emission. Comparisons of the detailed spectral shape show that our interacting LMC models match the observed $\gamma$-ray spectrum at photon energies below $\sim 10$ GeV, but still slightly overestimate the flux at higher energies, suggesting either more strongly energy-dependent transport than the models we have explored, or a cosmic ray injection spectrum steeper than the $p^{-2.2}$ that we adopt.
Read the full abstract

All authors Taaseen Islam, Hilay Shah, Mark R. Krumholz, Roland Crocker

Solar and Stellar Astrophysics

Collisionless stationary states of a stratified plasma in an expanding magnetic tube with stochastic heating

We investigate the collisionless kinetic structure of the upper solar atmosphere in an expanding magnetic field. Building on established velocity-filtration models, we extend the stochastic multi-temperature framework developed in previous works by including magnetic-field expansion and magnetic-moment conservation. Stochastic heating generates a non-Maxwellian boundary distribution through the superposition of particle populations associated with different temperatures. We consider a stationary two-component plasma confined within an expanding magnetic flux tube and subject to gravity, self-consistent electrostatic interactions, the Pannekoek–Rosseland electric field, and magnetic-moment conservation. Starting from the Vlasov equation, we derive fully analytical expressions for the distribution functions, density, and parallel, perpendicular, and total temperature profiles. The combined conservation of energy and magnetic moment generates a loss-cone distribution, reducing the density relative to the unmagnetized case and producing temperature anisotropy. For a single-temperature boundary, the competition between gravity and magnetic-moment conservation produces a maximum in the parallel temperature, for which we derive and numerically validate analytical scaling laws. With stochastic heating, gravitational filtering enhances the contribution of hotter populations at coronal heights, while magnetic-moment conservation amplifies the velocity-space anisotropy. Our analytical solution provides a collisionless benchmark for future kinetic models incorporating more realistic magnetic-field geometries, Coulomb collisions, and turbulent particle scattering.
Read the full abstract

Instrumentation and Methods for Astrophysics

PROVIDENCE: a ground station for space observation

PROVIDENCE is an ONERA-led project to install a new-generation optical ground station at the Haute-Provence Observatory (OHP) in South of France. Built around a 2.5-meter aperture, the station addresses five core science cases: space domain awareness, astronomy, laser activities and communications, atmospheric characterization, and instrumental prototyping. The project is organized into three parallel subprojects covering the telescope (Providence T), the building (Providence B), and the instruments (Providence I), the latter including the facility's first-light adaptive optics instrument, INTERSTELLAR. This paper presents the scientific motivation, the project architecture, the telescope technical specifications, the civil-engineering challenges associated with replacing the existing T152 equatorial telescope, and the station layout designed to host multiple co-active instrument teams. The project roadmap targets First Light in 2029, timed to track the close-approach flyby of asteroid Apophis.
Read the full abstract

All authors P.-L. Mayeur, J.-F. Sauvage, T. Fusco, N. Védrenne, C. Petit, D. Barbon-Dubosc, G. Leplat, G. Guerrin, A. Vincent-Randonnier, V. Michau, F. Dolon, M. Ferrari, F. Huppert, P.-E. Blanc, S. Basa, J. Floriot, B. Neichel, T. Guerini, J. Mingot, T. Jolivet, R. Buchard, N. Robbe, J. Dezeuze, S. Idri, F. Rosso, C. Boyer, S. Thomas, S. Egron

Instrumentation and Methods for Astrophysics

HARMONI at ELT: optical design of the spectrograph sub-system

HARMONI is an adaptive optics assisted, near infrared integral field spectrograph for the European Extremely Large Telescope. It covers a spectral range from 750 nm to 2450 nm with resolving powers from 3000 to 7000 and spatial sampling of 25 mas and 6 mas. It can operate in two adaptive optics modes, SCAO and MCAO. The spectrograph represents the last optical sub-system of the HARMONI instrument optical train that disperses the light and forms a two-dimensional spectrogram on its scientific detectors. Changes in the top-level goals and specifications have led to a significant revision of the spectrograph optical design. The current baseline design consists of a two-mirror unobscured collimator, a set of six volume-phase holographic grisms, and a two-mirror unobscured camera. The design is more compact than the previous version, it uses anamorphic freeform mirrors to correct aberrations, and is expected to reach a minimum transmission 57% with improved accessibility of the detector unit and reduced scattered light.
Read the full abstract

All authors Eduard Muslimov, Matthias Teczaa Edgar Castillo-Dominguez, James Kariuki, Liam Bolanda, Miriam Cisneros-Gonzalez, Kieran McCall, Sophie Paszynska

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

Laser Remelting for Reduced Porosity on Additively Manufactured Aluminium Mirrors

Additively manufactured (AM) AlSi10Mg mirrors are fabricated through laser powder bed fusion (LPBF), allowing the use of complex geometries such as lattices and organic structures that enable high mass reduction while maintaining mechanical stiffness. Micron-sized pores that cause optical scatter may form during LPBF as a consequence of deviations from the optimal processing window, particularly from laser energy input and scan strategy. This work proposes a laser remelting strategy aimed at reducing porosity; standard LPBF build steps automatically alternate with laser remelting passes, where previously deposited material is remelted during fabrication. Laser remelting is evaluated through fabricating 10 mm proof-of-concept cubes. Following single point diamond turning (SPDT), optical measurements characterised surface roughness and identified surface artefacts. The best-performing AlSi10Mg remelted cube exhibited no pores within sampled regions and achieved 6.4 nm average surface roughness, comparable to a conventionally manufactured RSA 6061 control cube (5.8 nm). Driven by these results, AM 52 mm diameter secondary sandwich mirrors were manufactured using LPBF and laser remelting. These incorporate an optimised diamond TPMS lattice to achieve a 50% mass reduction while accommodating design for AM considerations. Unlike the cube study, the optical surface of the remelted mirror after SPDT exhibited residual porosity and 11.8 nm average surface roughness. These results show that while the proof-of-concept confirmed the viability of laser remelting in reducing porosity within simple geometries, optimisation of the LPBF and SPDT processes are required to translate the benefits of laser remelting to lightweight AlSi10Mg AM mirrors.
Read the full abstract

All authors Joshua West, Valentina Oyarzun, Marcell Westsik, Younes Chahid, Magdalena Kraus, Scott McPhee, William Brzozowski, Samuel Tammas-Williams, Nicola Cayzer, Fraser Laidlaw, Sameer Dayanand Meshram, Daniel Ordnung, Berk Baris Celik, Michel Smet, Mirko Sinico, Wenjuan Sun, Brecht Van Hooreweder, Michael Harris, Stephen James, Carolyn Atkins

Cosmology and Nongalactic Astrophysics

Contribution to the 2026 Cosmology session of the 60th Rencontres de Moriond: Galactic moments: Understanding polarized foregrounds complexity in the quest for CMB primordial $B$ modes with LiteBIRD

Accurate modeling of polarized Galactic emission has become a major challenge for current and next-generation cosmic microwave background (CMB) $B$-mode experiments. Ignoring the spectral complexity of thermal dust and Galactic synchrotron emission when integrating along the line of sight and over large sky fractions inevitably leads to biases in CMB polarization analyses. In this work, we review how the future LiteBIRD satellite, which will benefit from an increased number of bands and sensitivity with respect to past CMB experiments such as Planck, will exploit its broader frequency coverage to characterize the spectral properties of interstellar medium emission across the three-dimensional structure of the Milky Way. We show that the canonical description of foreground spectral behavior will reach its limits for LiteBIRD, and that this challenge can be addressed using the moment expansion formalism.
Read the full abstract

All authors S. Vinzl, J. Aumont, L. Vacher, R. T. Génova-Santos, D. Adak, A. Rizzieri, for the LiteBIRD Collaboration

Cosmology and Nongalactic Astrophysics

Hubble tension: the shape wall

The standard "no-go theorem" against late-time solutions to the Hubble tension is essentially a normalization wall, since Baryon Acoustic Oscillation (BAO) measurements constrain the product $H_0r_d$, with $r_d$ the sound horizon at baryon drag. However, late-time solutions (which keep $r_d$ fixed) are tightly constrained not only by the BAO normalization $H_0r_d$, but also by the shape of the expansion history, i.e. the dimensionless expansion rate $E(z) \equiv H(z)/H_0$. We show that, if $r_d$ and the acoustic angular scale $\theta_s$ are fixed, an increase in $H_0$ needs to be matched by an equal fractional increase in the dimensionless distance integral $I \equiv \int dz/E(z)$: $\delta H_0/H_0 \simeq \delta I/I$. We use this to quantify the "shape wall" set by relative distance constraints on $E(z)$, which we reconstruct nonparametrically, using Gaussian Processes and the latest unanchored Type Ia Supernovae (SNeIa) and BAO data. In our most conservative analysis using PantheonPlus SNeIa and DESI DR2 BAO data, we find a maximum fractional increase in $H_0$ of $\lesssim 2\%$, falling well short of the $\gtrsim 8\%$ required to solve the tension. This shape wall holds even in the presence of early-time new physics, and limits the maximum increase in $H_0$ which can be contributed by late-time modifications to $E(z)$ at fixed $\theta_s$. Therefore, late-time modifications, whether invoked alone or alongside early-time new physics, face not only the well-known normalization wall, but also a stringent percent-level shape wall.
Read the full abstract

All authors Miguel A. Sabogal, Luis A. Escamilla, Davide Pedrotti, Edvig Selimaj, Sunny Vagnozzi

Solar and Stellar Astrophysics · Cross-listed from physics.atom-ph

Plasma screening and configuration interaction effects induced large enhancement on L-shell photoionization cross sections and opacity

An opacity model that incorporates improved treatments of both plasma screening and configuration interaction (CI) effects is proposed, and a 25-30% enhancement on the iron L-shell opacity is predicted at solar interior temperatures. It is originated from the plasma screening induced 14-17% enhancement on the photoionization cross sections and the CI induced 10-20% enhancement on photoexcitation and photoionization cross sections for open L-shell ions. These explain the long-standing discrepancy between theoretical and experimental iron opacity [Nature 517, 56], and the relatively weaker enhancements on chromium and nickel opacity [Phys. Rev. Lett. 122, 235001] due to the sensitivity of these effects to the different L-shell electron population and plasma temperature/density. This letter provides the systematic interpretation of L-shell opacity measurements at solar interior temperatures, and advances the accurate simulation of opacity and radiative transport in high-energy-density plasma.
Read the full abstract

All authors Fuyang Zhou, Shengbo Niu, Simei Lu, Chuangying Li, Xiang Gao, Yong Wu, Yizhi Qu, Jianguo Wang

Astrophysics of Galaxies

High-Resolution [CI] $({^3\mathrm{P}_2}\rightarrow{^3\mathrm{P}_1})$ and CO $(J=7\rightarrow6)$ Observations of Circumnuclear Disks in Nearby Seyfert Galaxies

We present [CI] $({^3\mathrm{P}_2}\mathrm{-}{^3\mathrm{P}_1})$, CO $(J=7\mathrm{-}6)$, and 800-GHz continuum observations of the circumnuclear disks (CNDs) in the Seyfert galaxies NGC 613 and NGC 1808 with the Atacama Large Millimeter/submillimeter Array. The images reveal the distributions of neutral gas and dust in nearby galaxy nuclei at a resolution of $0.2''$ (10 pc) which is unprecedented at these frequencies. Both lines and continuum exhibit peaks at the position of molecular tori. The maximum [CI]$({^3\mathrm{P}_2}\mathrm{-}{^3\mathrm{P}_1})/({^3\mathrm{P}_1}\mathrm{-}{^3\mathrm{P}_0})$ integrated-intensity ratio is $\approx1.1$ $(\mathrm{K~km~s}^{-1}$ scale) in the region of the torus, yielding an excitation temperature of $T_\mathrm{ex}\approx60~\mathrm{K}$ under the approximation of local thermodynamic equilibrium (LTE) and optically thin emission. A non-LTE radiative transfer analysis of the [CI] line ratio in the CND constrains the density and kinetic temperature of molecular gas to $n_\mathrm{H_2}>10^3~\mathrm{cm}^{-3}$ and $T_\mathrm{k}>50~\mathrm{K}$, respectively. To estimate the C/CO abundance ratio and physical conditions, we simultaneously modeled the [CI] and multiple CO lines, and obtained solutions with densities and temperatures of $n_\mathrm{H_2}\sim10^{3\mathrm{-}4}~\mathrm{cm}^{-3}$ and $T_\mathrm{k}\sim80\mathrm{-}250~\mathrm{K}$ in the central 63 pc. The C/CO ratios are $\approx2$ in NGC 613 and $\approx0.6$ in NGC 1808. C/CO tends to increase with $T_\mathrm{k}$ and decrease with $n_\mathrm{H_2}$ in the CND, which is consistent with models of photodissociation regions. Based on the kinematics of gas traced by [CI] $({^3\mathrm{P}_2}\mathrm{-}{^3\mathrm{P}_1})$ and CO $(J=7\mathrm{-}6)$, the mass of the supermassive black hole at the center of NGC 613 is estimated to be $\sim2\times10^7~M_\odot$.
Read the full abstract

All authors Dragan Salak, Suphakorn Suphapolthaworn, Yusuke Miyamoto, Naomasa Nakai, Masumichi Seta, Kazuo Sorai

Solar and Stellar Astrophysics

Comparative 3D Asymmetric Expansion and Angular Widths Evolution of Fast and Slow Coronal Mass Ejections

The radial and lateral dimensions of coronal mass ejections (CMEs) influence their duration and probability of encounter at Earth. These properties are linked to the expansion speed of CMEs in different radial and lateral directions; however, most earlier studies modeled CME evolution using a projected full ice-cream cone geometry, which does not distinguish between radial and lateral expansion. Our study investigates the asymmetric expansion (relative radial and lateral components) and kinematics of seven fast and seven slow CMEs within coronagraphic heights, using the Graduated Cylindrical Shell model. Our study confirms that CMEs expand asymmetrically, with lateral expansion exceeding radial expansion in both CME populations. This asymmetry limits the accuracy of the full ice-cream cone model. For both fast and slow CMEs, higher leading edge speeds are associated with higher expansion speeds. At a height of 10 Rs, slow CMEs with larger expansion speeds (lateral and radial) have larger angular widths (face-on and edge-on), whereas fast CMEs exhibit a negative correlation between lateral expansion speed and face-on angular width. We find that the expansion and propagation speeds of slow CMEs exhibit a two-phase evolution, whereas those of fast CMEs display more diverse trends. Overall, this study suggests that fast and slow CMEs evolve differently and should not be treated as a single population in statistical estimates of their physical parameters. Our study highlights the importance of estimating CME angular widths and expansion speeds along different directions, and beyond standard coronagraphic heights, to capture their complete physical evolution.
Read the full abstract

All authors Anjali Agarwal, Wageesh Mishra, Jie Zhang, Soumyaranjan Khuntia, Manuela Temmer

Astrophysics of Galaxies

The Stellar Mass Function of Gas-Rich Galaxies and the Underlying $M_{\rm HI}-M_{\rm star}$ Scaling Relation in the Local Universe

We estimate the Galaxy Stellar Mass Function (GSMF) of HI gas-rich galaxies using the 100% ALFALFA ($\alpha$) catalog, $\sim 98\%$ of which have optical counterparts in the Sloan Digital Sky Survey (SDSS) and a subset of them have counterparts in GALEX SDSS WISE Legacy Catalogue-2(GSWLC-2). We use the mass estimates from this subset which combines UV, optical and IR bands with individual dust corrections to recalibrate optical stellar mass estimates. We use a non-parametric method to estimate the GSMF of these gas-rich galaxies. The resulting, HI-selected GSMF is consistent with a single Schechter function with best-fit parameters $\left\{\phi_* (10^{-3}\, h_{70}^{3}\,\mathrm{Mpc}^{-3}\,\mathrm{dex}^{-1}), \log_{10} (M_*/M_{\odot}) + 2\log_{10} h_{70}, \alpha \right\} = \left\{2.30^{+0.12}_{-0.12}, \,10.83^{+0.01}_{-0.01},\, -1.14^{+0.02}_{-0.02}\right\} $. Additionally, the red and blue populations are each well described by a single Schechter function. After correcting for selection effects, we find that the red population accounts for only $\sim18\%$ of gas-rich galaxies by number, yet contributes $\sim54\%$ of the total stellar mass, with the blue population accounting for the rest. Using an optically selected sample and a joint optical-HI sample, we find gas-rich galaxies represent $\sim 33\%$ of the total stellar mass density and $\sim 39\%$ of the total galaxy number counts in the local Universe. We use the GSMF and the HI mass function (HIMF) of the HI-selected sample to obtain the $M_{\rm HI}-M_{\rm star}$ relation, which is free from selection bias.
Read the full abstract

All authors Tanya Tripty, Nishikanta Khandai, Saili Dutta, Khushi Sharma

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

MeV Electrophilic Axion-like Particles from Sun

This work explores the production of an MeV-scale electrophilic axion-like particles (ALPs) by utilizing the monochromatic 5.5MeV photon resulting from the nuclear fusion processes in the Sun. These 5.5MeV photons can undergo the Compton-like scattering with the ambient electrons in the solar matter to produce a substantial flux of MeV ALPs. Upon reaching the Earth, such ALPs can be detected via the same electron coupling, offering a new opportunity for the dark matter (DM) direct detection experiments to probe the previously unexplored parameter regions. We show that the existing data of LZ, PandaX-4T, and Borexino can attain the sensitivities $g_{ae} \lesssim 3.7 \times 10^{-6}$, $g_{ae} \lesssim 3.7 \times 10^{-6}$ and $g_{ae} \lesssim 1.7 \times 10^{-6}$, respectively, for $m_a \lesssim 1$MeV. An optimistic 200 tonne$\times$year exposure by PandaX-xT can reach $g_{ae}\lesssim 1.6 \times 10^{-6}$ for most of the mass window $m_a < 1$MeV and even $g_{ae} \lesssim 1.5 \times 10^{-7}$ with $m_a$ approaching 1MeV. Despite the stringent constraints from different laboratory experiments and astrophysical observations, our obtained limits from LZ, PandaX-4T, and Borexino can probe new parameter regions, specifically in the mass window $0.4\,{\rm MeV} \lesssim m_a \lesssim 1$MeV.
Read the full abstract

Instrumentation and Methods for Astrophysics

Infrasound Newtonian Noise Estimation at the Einstein Telescope Candidate Site Sos Enattos

We investigate the seasonal variability of atmospheric infrasound and its contribution to Newtonian noise (NN) at the Sos Enattos site, a leading candidate location for the Einstein Telescope (ET). Infrasound data recorded at three stations, SOE0 (surface), SOE1 ($-84$ m) and SOE3 ($-160$ m), are analysed over multiple seasons to characterise both the temporal variability and the depth dependence of the acoustic field. At the surface, the median amplitude spectral density (ASD) at 1 Hz shows no consistent seasonal modulation, and its variability is driven primarily by wind; the apparent dependence on sea state largely reflects its correlation with wind speed. Underground, where local wind noise is suppressed, a winter maximum of 4–7 dB above the quietest season is observed in both years of data. Using the measured sound-pressure spectra, we estimate the corresponding NN contribution for an underground cavern with the dimensions foreseen for ET, assuming an isotropic field of acoustic plane waves. The surface pressure field, acting through the atmosphere above a test mass at 200 m depth, exceeds the ET-D design sensitivity in median conditions between 1.6 and 7.1 Hz, by up to a factor of 4.7. Inside the cavern, the median NN remains a factor of 1.4 below ET-D with the pressure measured at SOE1, but exceeds it between 2.5 and 8.7 Hz, by up to a factor of 3.5, with the pressure measured at SOE3. The attenuation is not monotonic with depth: the deeper station SOE3 shows higher levels than SOE1 above a few hertz, pointing to local, anthropogenic sources in the underground infrastructure. These results show that underground placement is essential but not sufficient on its own, and that the acoustic emission of the underground infrastructure and future laboratory equipment will require careful evaluation and mitigation to keep atmospheric NN below the ET design sensitivity.
Read the full abstract

All authors Wathela Alhassan, D. Rosinska, M. Suchenek, E. Fenyvesi, T. Bulik

Astrophysics of Galaxies

What's Missing in AGN Feedback? Lessons learnt from Magneticum, IllustrisTNG and Simba

Accurately balancing gas reservoirs, star formation, and feedback across cosmic time remains a central challenge for galaxy formation models in modern hydrodynamical simulations. While different feedback prescriptions reproduce selected local galaxy properties with varying success, the most pronounced discrepancies emerge in predictions for the hot gas content of dark matter halos. We examine three state-of-the-art cosmological simulations: Magneticum, IllustrisTNG, and SIMBA, which struggle to simultaneously reproduce observed galaxy and halo gas properties in the local Universe. We confront their predictions with spatially resolved galaxy data from MaNGA and recent constraints on the hot gas mass fraction-halo mass (fgas-Mh) relation from eROSITA and Sunyaev-Zel'dovich (SZ) measurements. Reproducing the observed fgas-Mh relation requires strong active galactic nucleus (AGN) feedback. However, such feedback often leads to excessive quenching in simulated galaxy populations. Magneticum and SIMBA match the observed gas fraction relation but predict an overabundance of quenched galaxies. In contrast, IllustrisTNG implements weaker AGN feedback, yielding more realistic star-forming fractions but systematically overpredicting hot gas masses in massive groups and poor clusters. Overall, these tensions indicate current feedback models remain incomplete, not only in the total energy injected but also in the timing, location, and coupling of this energy to the surrounding gas. Our results therefore highlight the need to revisit subgrid feedback prescriptions and develop more self-consistent models capable of simultaneously regulating galaxy growth and the thermodynamic properties of halo gas. Motivated by this discrepancy, a companion study will explore whether the feedback strengths required to match halo gas constraints inevitably lead to overquenching and distorted galaxy demographics.
Read the full abstract

All authors D. T. Mazengo, P. Popesso, I. Marini, L. M. Valenzuela, N. de Isídio, V. Toptun, Mirjana Pović, Petri Väisänen, J. M. Sunzu, R.-S. Remus, V. Biffi, K. Dolag, R. Davé, A. Fraser-McKelvie, A. Dev, S. Vladutescu-Zopp, A. Merloni, G. Ponti, C. Aydar, S. Shreeram, L. A. Kahinga, J. O. Chibueze, P. Privatus

Instrumentation and Methods for Astrophysics

Fast Inference on Astronomical Time Series with Trans-Dimensional Flow Matching Posterior Estimation

The analysis of time series plays an important part in the study of (fast) transient events, including gamma-ray bursts, magnetar bursts, fast radio bursts, and solar flares. A common approach is to decompose the time series into pulses and study the pulse characteristics, such as location and amplitude, in order to constrain physical models of the source and its environment. However, estimating both the number and characteristics of these pulses presents a trans-dimensional inference problem that traditional sampling methods such as Markov Chain Monte Carlo (MCMC) and Nested Sampling struggle to solve efficiently. Simulation-based inference methods, often incorporating machine learning techniques, provide an alternative approach when traditional approaches are insufficient. Here, we introduce trans-dimensional Flow Matching Posterior Estimation (t-FMPE) implemented on a transformer architecture capable of efficient, amortized trans-dimensional inference on uniformly sampled univariate time series data. In this initial study, we apply the method to three test cases: simulated time series with known ground-truth parameters, observational data of Fast Radio Bursts and observations of X-ray bursts from magnetars. We show that t-FMPE achieves qualitative agreement with MCMC reference posteriors, successfully reproducing parameter correlations, as quantified through classifier two-sample tests. The trained network performs inference several orders of magnitude faster than MCMC and nested sampling, reaching sampling rates of 100 posterior samples per second for an 80-dimensional parameter space. The results demonstrate potential of t-FMPE for large-scale analysis of time series datasets when traditional sampling methods become infeasible, and also enable inferring unbiased posteriors in the presence of observational biases such as dead time.
Read the full abstract

All authors Nina van der Meulen, Tin Hadži Veljković, Daniela Huppenkothen, Benjamin Kurt Miller, Christoph Weniger

Cosmology and Nongalactic Astrophysics

Inferences for f(R) Models from Late-Time Megamaser Observational Data

In this work, we study three widely used models of f(R) gravity, namely HuSawicki, Starobinsky and ArcTanh along with the standard cosmology model ($\Lambda$CDM). For this, we employ the megamaser angular diameter distance and velocity measurements from the Megamaser Cosmology Project, which provide a purely geometric determination of the Hubble constant. We constrain the parameters using the Markov Chain Monte Carlo method. Our results show that values of the Hubble Constant, $H_{0}$, obtained for all four models are in concordance with its value obtained from other late-time observational data such as SNe Ia. The constraints on $H_{0}$ in all the models under study are restrictive and the marginalized estimates lie close to 73 $\mathrm {km s^{-1} Mpc^{-1}}$. The marginalized estimates of the deviation parameter, b, for the three f(R) gravity models lie close to zero. This late-time dataset, thus predicts that f(R) models mimic $\Lambda$CDM. However, the matter density, $\Omega_m$, remains weakly constrained for all the models with its marginalized estimate close to 0.5. Further, comparison of the four models (f(R) models and $\Lambda$CDM) using information criteria such as Akaike Information Criterion and Bayesian Information Criterion shows that within current uncertainties, the dataset finds f(R) models statistically indistinguishable from $\Lambda$CDM. This is consistent with the fact that the favoured value of b for each of the f(R) models lies close to zero.
Read the full abstract

Astrophysics of Galaxies

MIGHTEE: The Host-Galaxy Associated Catalogue of the Radio Sources in MIGHTEE Continuum Data Release 1

Radio continuum surveys provide samples of active galactic nuclei (AGN) and star forming galaxies (SFGs) to high redshifts, free of biases due to dust obscuration. However, radio detected sources require multi-wavelength counterparts to understand their intrinsic properties (e.g. redshift, stellar mass) and to study the evolution of star formation and AGN activity. In this work we present host galaxy counterparts for the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey continuum Data Release 1 in regions with the best ancillary data (totalling 7.5 sq. deg). We combine statistical cross-matching and visual inspection to identify Ks-band selected host galaxies, and additionally combine multiple radio components into single physical objects, where needed. This results in a combined radio catalogue of ~66 000 sources, with host counterparts and redshifts identified for ~95 per cent of sources in the COSMOS field, ~91 per cent in XMM-LSS and ~90 per cent in CDFS-DEEP. This includes a significant fraction of sources with spectroscopic redshifts within the COSMOS field (~50 per cent), with ~30 and ~20 per cent in the XMM-LSS and CDFS-DEEP fields respectively. Using the cross-matched catalogue, we make an initial identification of radio-excess and star forming galaxies based on comparisons of the radio luminosities to host star formation rates. Using this split as a proxy for radio loud AGN or SFGs, we present expectations for the redshift distributions of these sources, finding broad agreement with those from deep radio luminosity functions and simulated catalogues.
Read the full abstract

All authors C. L. Hale, J. Hamlett, I. H. Whittam, M. J. Jarvis, M. J. Hardcastle, S. L. Jung, N. Stylianou, R. G. Varadaraj, M. Vaccari, L. Barchiesi, P. N. Best, K. K. L. Charlton, E. Malefahlo, E. Moravec, J. P. Moss, M. Rakototafika, T. F. Rarivoarinoro, F. Sinigaglia, M. N. Tudorache, S. I. Loubser, L. Heino, A. Saintonge, A. Mazumder, N. Netshiavha, D. J. B. Smith, L. Verdes-Montenegro, J. Delhaize, H. Pan, F. Pozzi, Z. Randriamanakoto, L. Stockenstroom, A. A. Vărăşteanu, E. Vardoulaki, W. L. Williams, R. A. A. Bowler, L. Marchetti, A. Matthews

Solar and Stellar Astrophysics

Transport of angular momentum and chemical elements by the MRI dynamo in stellar radiative zones

The question of angular momentum transport by magnetic fields has recently been revived by the detection of magnetic fields in the deep interior of red giant stars. We aim at characterizing the efficiency of the transport of angular momentum and chemical elements in a stellar radiative zone subject to the magneto-rotational instability, in situations where hydrodynamical instabilities are not triggered. We use a large set of zero-net flux shearing-box simulations modelling a portion of a stellar radiative zone located close to the equatorial region. The shear is imposed in the direction of gravity, mimicking a radial differential rotation and rotation is perpendicular to gravity. We aim at establishing scaling laws between the transport efficiency and key parameters such as stratification and rotation. We obtain simulations where the MRI is triggered on the self-consistently built longitudinal field and then produces a state of self-sustained turbulence that we associate with dynamo action. We find that in the parameter range explored, both rotation and stable stratification strongly affect the vertical transport of angular momentum which is always dominated by the Maxwell stress component. In a similar way, we find that the transport of chemical composition is even more strongly affected by the stratification but less severely by rotation. Scaling laws are established and tentative extrapolations to stellar values are discussed. Transport by the stratified MRI dynamo is very promising to reconcile stellar evolution models and asteroseismic inversions of internal rotation rates of stars.
Read the full abstract

All authors Laurène Jouve, François Lignières, Jérôme Guilet, Alexis Vanbesien

Astrophysics of Galaxies

Measuring Outflow Distances in NGC 5548 Using Absorption-Line Variability Diagnostics

AGN-driven outflows serve as a key channel through which the energetic central engine influences host galaxy evolution. Among the physical properties of outflows, their radial distance from the galactic nucleus is particularly important for assessing AGN feedback. In this study, we investigate the UV outflow components in NGC 5548 by analyzing the variability of C IV absorption troughs in optical spectra obtained through multiple HST observations during 2013 and 2014. We construct a set of variability-based diagnostic events, labeled G1 and G2, which are sensitive to the recombination timescale ($t_r$) of ionized gas. By combining these with mock light curves generated using a damped random walk (DRW) model, we numerically establish a mapping between the G1 event probability and $t_r$. This approach allows us to constrain the radial distances of outflow components 1 and 6, whose absorption variability is primarily driven by changes in the incident ionizing continuum, to be $0.77^{+0.10}_{-0.10}$ and $1.72^{+1.74}_{-1.72}$ pc, respectively. These results are consistent with those obtained using a different method in our previous study, as well as with values reported in the literature.
Read the full abstract

Astrophysics of Galaxies

A RELHIC twin candidate near the galaxy M51

We report the discovery of a pair of H I clouds near M51 (NGC 5194) using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These clouds have no optical counterparts and are potential candidates for Reionization-Limited H I Clouds (RELHICs). We search for compact H I sources in deep FEASTS observations using SoFiA and remove objects with optical counterparts through cross-matching with the DESI Legacy Imaging Surveys. The remaining candidates are modelled as hydrostatic H I structures embedded in Navarro-Frenk-White dark matter haloes and compared with RELHIC predictions and TNG50 simulations. We identify two H I clouds, Cloud S and Cloud N, at projected distances of 70–90 kpc from M51. Each cloud has an H I mass of approximately 10^6.5 solar masses, a velocity dispersion of about 20 km/s, and no detectable optical counterpart down to a g-band surface brightness limit of approximately 27.5 mag arcsec^-2. Their stellar luminosities are constrained to be below 10^5 solar luminosities. Their H I properties are consistent with RELHIC predictions, corresponding to host halo masses of 3.7 +- 0.4 * 10^9 solar masses. Cloud S and Cloud N are promising but not definitive RELHIC candidates. A tidal origin remains possible in the interacting M51 system, especially because the clouds are unresolved by FAST and Cloud N may show a velocity gradient. Future high-resolution interferometric observations will be crucial for distinguishing between starless dark matter haloes and tidal debris.
Read the full abstract

All authors Qingze Chen, Jie Wang, Yingjie Jing, Chen Xu, Tiantian Liang, Zerui Liu, Zhipeng Hou, Yiwei Xu, Cheng Cheng, Wei Du

Astrophysics of Galaxies

Dissecting the multiple-component outflow in NGC 5548 with absorption-line Variability

AGN-driven outflows are routinely invoked as a key agent of supermassive black holes to regulate the evolution of galaxies. The radial distance from the central engine is a crucial parameter for evaluating the impact of these outflows on the host galaxy. In this work, we estimate the radial distances of ultraviolet (UV) outflow components in NGC 5548 using the most up-to-date absorption-line variability method, combined with multi-epoch HST/COS spectroscopy from the 2014 AGN STORM campaign and archival data observed in 2013. The recombination timescale (tr) of the absorbers are measured by analyzing the detection rate curves of absorption-line variability. In particular, the detection rate curves of the absorption troughs showing blended multiple velocity components are featured by distinct “multi-step' profiles, allowing for measuring tr for individual components. Among the 6 identified outflow components, four are found to be a few pc from the center and two are 30-40 pc away. Our results agree well with the more reliable results in the literature on components 1 and 4, and show overall consistency with previous works, demonstrating the power of our new methodology especially when it is aided by densely sampled HST spectra.
Read the full abstract

High Energy Astrophysical Phenomena

The Extended Ultrahigh-energy Gamma-Ray Emission in the Vicinity of PSR J2238+5903

We present a comprehensive analysis of the recently discovered TeV gamma-ray source, LHAASO J2238+5900. Based on data collected from the LHAASO, our fitting results suggest that the source is significantly extended with an angular extension of 0.54{\deg} \pm 0.01{\deg} and is spatially coincident with the pulsar PSR J2238+5903. Its spectrum is characterized by a power-law with a cutoff at 41.0\pm 3.5 TeV. Additionally, the source exhibits a significant signal of 7.9\sigma above 100 TeV, implying that it is a PeVatron candidate. While the gamma-ray emission is consistent with a pulsar wind nebula (PWN) scenario, the relatively large extension size also allows for a halo interpretation, potentially caused by electron-positron pairs escaping from the PWN.
Read the full abstract

All authors Zhen Cao, F. Aharonian, Y.X. Bai, Y.W. Bao, D. Bastieri, X.J. Bi, Y.J. Bi, W. Bian, J. Blunier, A.V. Bukevich, C.M. Cai, W.Y. Cao, Zhe Cao, J. Chang, J.F. Chang, E.S. Chen, G.H. Chen, H.K. Chen, L.F. Chen, Liang Chen, Long Chen, M.J. Chen, M.L. Chen, Q.H. Chen, S. Chen, S.H. Chen, S.Z. Chen, T.L. Chen, X.B. Chen, X.J. Chen, X.P. Chen, Y. Chen, N. Cheng, Q.Y. Cheng, Y.D. Cheng, M.Y. Cui, S.W. Cui, X.H. Cui, Y.D. Cui, B.Z. Dai, H.L. Dai, L.X. Dai, Z.G. Dai, Danzengluobu, Y.X. Diao, A.J. Dong, X.D. Duan, J.H. Fan, Y.Z. Fan, J. Fang, J.H. Fang, K. Fang, C.F. Feng, H. Feng, L. Feng, S.H. Feng, X.T. Feng, Y. Feng, Y.L. Feng, S. Gabici, B. Gao, Q. Gao, W. Gao, C. Ge, M.M. Ge, T.T. Ge, L.S. Geng, G. Giacinti, G.H. Gong, Q.B. Gou, M.H. Gu, W.M. Gu, F.L. Guo, J. Guo, K.J. Guo, X.L. Guo, Y.Q. Guo, R.P. Han, O.A. Hannuksela, M. Hasan, H.H. He, H.N. He, J.Y. He, X.Y. He, Y. He, S. Hernández-Cadena, C. Hou, X. Hou, H.B. Hu, S.C. Hu, D.H. Huang, F. Huang, J.J. Huang, X.L. Huang, X.T. Huang, X.Y. Huang, Y. Huang, Z.J. Huang, A. Inventar, X.L. Ji, H.Y. Jia, K. Jia, H.B. Jiang, K. Jiang, X.W. Jiang, Z.J. Jiang, M. Jin, S. Kaci, M.M. Kang, I. Karpikov, D. Khangulyan, D. Kuleshov, W.H. Lei, Cheng Li, Cong Li, D. Li, F. Li, H.B. Li, H.C. Li, Jian Li, Jie Li, K. Li, L. Li, R.L. Li, T.Y. Li, W.L. Li, X.R. Li, X.Y. Li, Y. Li, Zhe Li, Zhuo Li, E.W. Liang, Y.F. Liang, S.J. Lin, B. Liu, C. Liu, D. Liu, D.B. Liu, H. Liu, J. Liu, J.L. Liu, J.R. Liu, M.Y. Liu, R.Y. Liu, S.M. Liu, T. Liu, W. Liu, Y. Liu, Y. Liu, Y.N. Liu, Y.Q. Lou, Q. Luo, Y. Luo, H.K. Lv, B.Q. Ma, L.L. Ma, X.H. Ma, I.O. Maliy, J.R. Mao, Z. Min, W. Mitthumsiri, Y. Mizuno, G.B. Mou, A. Neronov, K.C.Y. Ng, S.C.-Y. Ng, M.Y. Ni, L. Nie, L.J. Ou, Z.W. Ou, P. Pattarakijwanich, Z.Y. Pei, D.Y. Peng, J.C. Qi, M.Y. Qi, J.J. Qin, H.L. Qu, A. Raza, C.Y. Ren, M.Q. Ruan, D. Ruffolo, A. Sáiz, D. Savchenko, D. Semikoz, L. Shao, O. Shchegolev, Y.Z. Shen, X.D. Sheng, F.W. Shu, H.C. Song, Yu.V. Stenkin, Y. Su, C.Y. Sun, D.X. Sun, H. Sun, J.X. Sun, M. Sun, Q.N. Sun, X.N. Sun, Z.B. Sun, N.H. Tabasam, J. Takata, P.H.T. Tam, H.B. Tan, Q.W. Tang, R. Tang, Z.B. Tang, W.W. Tian, C.N. Tong, L.H. Wan, C. Wang, D.H. Wang, G.W. Wang, H.G. Wang, J.C. Wang, J.F. Wang, J.S. Wang, K. Wang, Kai Wang, Kai Wang, L.P. Wang, L.Y. Wang, W. Wang, X.G. Wang, X.J. Wang, X.Y. Wang, Y. Wang, Y.D. Wang, Z.H. Wang, Z.X. Wang, Zheng Wang, D.M. Wei, J.J. Wei, Y.J. Wei, T. Wen, S.S. Weng, C.Y. Wu, H.R. Wu, Q.W. Wu, S. Wu, X.F. Wu, Y.S. Wu, S.Q. Xi, J. Xia, G.M. Xiang, D.X. Xiao, G. Xiao, Y.F. Xiao, B.H. Xie, Y.L. Xin, H.D. Xing, Y. Xing, D.R. Xiong, B.N. Xu, C.Y. Xu, D.L. Xu, R.X. Xu, S.S. Xu, L. Xue, D.H. Yan, T. Yan, C.Y. Yang, F.F. Yang, L.L. Yang, M.J. Yang, R.Z. Yang, W.X. Yang, Z.H. Yang, Z.G. Yao, X.A. Ye, L.Q. Yin, N. Yin, X.H. You, Z.Y. You, Y.H. Yu, Q. Yuan, H. Yue, H.D. Zeng, T.X. Zeng, W. Zeng, X.T. Zeng, M. Zha, B. Zhang, B.B. Zhang, B.T. Zhang, C. Zhang, H. Zhang, H.M. Zhang, H.Y. Zhang, J.L. Zhang, J.Y. Zhang, L.Y. Zhang, Li Zhang, P.F. Zhang, R. Zhang, R.Y. Zhang, S.R. Zhang, S.S. Zhang, S.Y. Zhang, W. Zhang, X. Zhang, X.L. Zhang, X.P. Zhang, Yi Zhang, Yong Zhang, Z.P. Zhang, J. Zhao, L. Zhao, L.Z. Zhao, X.H. Zhao, F. Zheng, T.C. Zheng, B. Zhou, H. Zhou, J.N. Zhou, L. Zhou, M. Zhou, P. Zhou, R. Zhou, X.X. Zhou, X.X. Zhou, B.Y. Zhu, C.G. Zhu, F.R. Zhu, H. Zhu, K.J. Zhu, Y.F. Zhu, Z.F. Zhu, Y.C. Zou, X. Zuo

Instrumentation and Methods for Astrophysics

GCR Spectra Reconstructed with Neutron Monitor Yield Function and Artificial Neural Networks: Comparison of Two Methods

We present a framework that reconstructs time-resolved galactic cosmic-ray (GCR) proton and helium energy spectra from the global neutron monitor network, providing data about GCR flux without direct satellite observations. Two methods are utilized and compared: a calibrated yield function plus force-field scheme and artificial neural networks trained on multi-station neutron monitor count rates coupled with heliophysical indices. The reconstructed spectral time series reproduce both large-scale solar-cycle modulation and short-term disturbances and extend to periods lacking daily spacecraft data, including 2006-2011 (consistent with PAMELA) and 2019-2022 (consistent with AMS-02 Bartels rotation averages). Artificial neural networks deliver excellent performance across energies, with markedly lower mean absolute percentage error and $\chi^2/\mathrm{dof}$ near unity. A thorough validation confirms robustness and establishes neutron monitors as an effective real-time GCR spectrometer that can be utilized for various purposes.
Read the full abstract

All authors Stepan Siruk, Vladislav Alekseev, Victor Kuzminov, Rustam Yulbarisov, Andrey Mayorov

Cosmology and Nongalactic Astrophysics

Reviving the Affleck-Dine Curvaton Scenario with Compensated Isocurvature Perturbations

We revisit the Affleck-Dine curvaton scenario in the framework of supersymmetry, taking Q-ball formation into account. In this scenario, a scalar field carrying baryon number, namely the Affleck-Dine field, generates both the curvature perturbations and the baryon asymmetry of the Universe. However, correlated baryon isocurvature perturbations are inevitably generated, resulting in excessively large correlated matter isocurvature perturbations. We propose a mechanism to suppress them, in which the correlated baryon isocurvature perturbations are canceled by anti-correlated dark matter isocurvature perturbations, thereby realizing compensated isocurvature perturbations. We show that the scenario can simultaneously account for the baryon asymmetry and the observed curvature perturbations of the Universe without generating excessively large matter isocurvature perturbations.
Read the full abstract

Solar and Stellar Astrophysics

HARPS Abundances with Korg I: 22 Element Abundances for 426 Red Giant Stars

Large stellar surveys have revealed the global abundance structure of the Milky Way, but small high-fidelity spectral samples offer a critical complement of nucleosynthetic depth. We aim to access the encoded information in an ensemble of abundances by leveraging highest-quality spectra. We used HARPS spectra (R=115,000) to determine (Teff, log(g), [M/H], vmic and vsini) and 22 element abundances (Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Fe, Ni, Zn, Sr, Y, Zr, Mo, Ba, La, Ce, Nd and Eu) for 426 red giant stars at a median internal precision of $\sim$0.02 dex evaluated from analysing repeat observations of a subset of stars. Stellar parameters and line-by-line abundances were obtained using the modern spectral synthesis code Korg – the first time it has been used for HARPS. Comparisons with the literature reveal good overall agreement. A minor 0.1 dex offset in metallicity and specific discrepancies in individual element abundances are attributed to local thermal equilibrium assumptions and inaccuracies in atomic data. We show that 22 individual elements can be collapsed into a generative 6-parameter latent-variable model of shared enrichment patterns expressed in different per-star fractions; this model accurately generates the abundances with a median $\chi_{reduced}^2 = 6$. We report element gradients with respect to selected elements from different nucleosynthetic families. These gradients are a measure of inter-element production efficiencies and indicate multiple r-process production sites. Our analysis shows that abundances occupy a low-dimensional subspace, but joint (gradient-based) information encodes nucleosynthetic signatures. We have developed a Korg-pipeline to apply across evolutionary states on high-resolution spectra to provide our precision catalogue to serve as empirical constraints on chemical evolution and as a set of benchmark red giant abundance measurements.
Read the full abstract

All authors Sarah E. Aquilina, Melissa Ness, Adam J. Wheeler, Sven Buder

Astrophysics of Galaxies

Accurate Extragalactic Magnetic Fields from Faraday Rotation with Optimal Dispersion Measure Estimators

Faraday rotation measures (RMs) are one of our few observational tools for measuring magnetic field strengths in extragalactic systems, but converting an RM to a magnetic field estimate requires knowledge of the electron column density – the dispersion measure (DM) – to the RM source. Because DMs are difficult to measure for most extragalactic radio galaxies, observers have adopted a range of strategies to estimate them from more easily measured quantities, but the accuracy of these approaches is poorly known. To address this, we carry out simulated observations of high-resolution magnetohydrodynamic simulations of a range of galactic environments to explore the performance of various possible DM estimators. We obtain the best results using an estimator $\mathrm{DM} \propto \mathrm{EM}^{\alpha} \ N_\mathrm{Hi}^{\beta}$, where EM is the emission measure and $N_\mathrm{Hi}$ is the atomic hydrogen column density, with exponents $\alpha \approx 0.2-0.4$ and $\beta \approx 0-0.1$ depending on galactic environment (e.g., galaxy centres versus outskirts, and dwarfs versus spirals). We show that the variation of these exponents with environment can be understood in terms of simple physical arguments. Based on our tests, we provide recommended best practices for extracting galactic magnetic fields from RM data as a function of galactic environment and of proxy data availability, and show that using these methods one can obtain field measurements that are accurate to a few tenths of a dex. This work therefore represents an important step toward making use of RM data from next-generation surveys with the SKA.
Read the full abstract

All authors Hilay Shah, Mark Krumholz, Naomi McClure-Griffiths, Zipeng Hu

Solar and Stellar Astrophysics

On the numerical reliability of nonadiabatic pulsation solutions in AGB envelopes

We re-examine integration methods of linear nonadiabatic radial pulsations in the envelopes of AGB and post-AGB stars using the subspace shooting formulation. In the article the accuracy of the direct integration, the Riccati based method and a continuous renormalization nonlinear integration method are examined. The last method is not widely used in stellar pulsation studies and is therefore of interest in the present comparison. A tracking transformation of the two independent solution vectors is introduced. It helps to improve the accuracy and stability of the solution and can be readily implemented for any of these integration methods. Minimum singular value maps are used to examine the structure of the pulsation spectrum over the complex-frequency plane and also to assess the fidelity with which the integration methods transport the solution subspace. These maps provide a practical alternative to the customarily used maps based on the determinant. We use these maps together with the eigenfrequency spectra to establish mesh conditions required for the reliable transport of the solution subspace, especially in the nonlinear methods. To measure the quality of the subspace transport a spillover measure is introduced. Using this measure it is also possible to verify the placement of the inner boundary, which in the case of radial modes, cannot be guided by propagation criteria based on the Lamb frequency. Visual inspection of the eigenfrequency spectra, the singular value maps, or the eigenmode amplitudes and behavior in the deep envelope does not necessarily imply preservation of the transported-solution subspace. The spillover measure provides a much clearer indicator. This approach provides explicit quantitative diagnostics for assessing the reliability of strongly nonadiabatic pulsation calculations in AGB/post-AGB stellar envelopes.
Read the full abstract

High Energy Astrophysical Phenomena

Relativistic Oblique Shocks at Finite Temperature: Detachment Angle, Shock Polars, and the Turning Parameter

Oblique shocks are ubiquitous in high-energy astrophysical environments, yet a systematic analytical treatment of how finite upstream temperature influences the maximum deflection angle has been lacking. We address this problem by developing a unified thermodynamic framework based on a novel dimensionless quantity, the turning parameter, which encapsulates the equation of state, upstream Mach number, and thermal state of the flow into a single variable. Starting from the relativistic Rankine-Hugoniot conditions and the Taub adiabat, we derive a compact turning relation and a first-order perturbative expansion in the upstream thermal parameter. We show that any finite upstream temperature monotonically suppresses the maximum deflection angle relative to the cold-fluid limit, implying that cold models systematically overestimate shock attachment. In the combined ultra-thermal and ultra-relativistic limit, the turning parameter saturates to a universal value, yielding an asymptotic detachment angle that depends only on the equation of state. Numerical shock-polar calculations validate the analytical results and reveal a non-monotonic dependence of the detachment angle on the Mach number at intermediate temperatures, arising from the competition between thermal pressure and bulk kinetic energy-a distinctly relativistic thermal effect absent in both the cold and ultra-hot limits. As an illustrative astrophysical application, we apply the framework to the Crab pulsar wind nebula, demonstrating how finite-temperature effects modify the termination-shock morphology and the observed torus geometry.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

Charged Higgs Search at Future Neutrino Telescope and Higgs Factory

We investigate the discovery potential of the charged Higgs boson at future lepton colliders and neutrino telescopes within two simplified benchmark scenarios with universal Yukawa couplings. We demonstrate that both the muon-track and cascade events at neutrino telescopes can be exploited to search for or constrain the charged Higgs boson through resonant neutrino scattering process, and we compare the sensitivities obtained under different astrophysical neutrino flux models. We study the effect of the detector volume of neutrino telescopes on the discovery potential of the charged Higgs boson. We study the signal of charged Higgs boson at future lepton colliders over different kinematic regimes using optimized reconstruction strategies. We compare the discovery potential of the two experimental approaches over the relevant parameter space and find that future lepton colliders generally provide better sensitivity in most of the parameter space, while neutrino telescopes with very large detector volume can offer competitive and complementary sensitivity in the heavy-mass region of the charged Higgs boson.
Read the full abstract

All authors YiCheng Dai, Wei Liao, Yi-Song Lu, Qi-Shu Yan