Peer Revenue Now accepting positions

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Friday 14 August

Solar and Stellar Astrophysics

JWST/NIRSpec Spectra for Three Ultracool Brown Dwarfs Detected in Extragalactic Surveys: Further Evidence for Phosphine Absorption

We present JWST NIRSpec prism spectra for three ultracool brown dwarfs discovered in extragalactic survey data, two from the JWST Advanced Deep Extragalactic Survey (JADES), and one from Public Release IMaging for Extragalactic Research (PRIMER) survey observed as part of the Mirage or Miracle (MoM) program. The spectra for these sources indicate that one is a T6 dwarf (JADES-GS-BD-11, T$_{\mathrm{eff}} = \sim 700$ K) and two are Y0-Y1 dwarfs (JADES-GS-BD-5, T$_{\mathrm{eff}} = \sim 400$ K, and MoM-239450, T$_{\mathrm{eff}} = \sim 500$ K). Model atmospheric fits with $\texttt{NIFTY}$ to the spectra are consistent with this classification, and indicate that JADES-GS-BD-5 is only $\sim 150$ pc from the Sun, MoM-239450 is $\sim 700 - 800$ pc from the Sun, and JADES-GS-BD-11 is $\sim 1$ kpc from the Sun, with these latter two more distant sources being best fit at sub-solar metallicities. JADES-GS-BD-5 has an observed spectrum with significantly weaker J and H band emission than Y dwarf atmospheric models, potentially indicating the presence of water ice clouds in the brown dwarf. The spectrum for JADES-GS-BD-11 has a feature at 4.3$\mu$m consistent with absorption from the rarely seen phosphine molecule at $2.6\sigma$ confidence. Given the low metallicity for this source ([M/H] $ = -0.7$), our finding supports the theory that detecting phosphine in brown dwarf atmospheres is tied to atmospheric metallicity. JWST/NIRSpec spectroscopy continues to be a powerful tool for understanding the properties of these distant, and very cold brown dwarfs.
Read the full abstract

All authors Kevin N. Hainline, Samuel A. Beiler, Adam J. Burgasser, Evan S. Chen, Jakob M. Helton, Jarron Leisenring, Brittany E. Miles, Mark S. Marley, Sagnick Mukherjee, Stefi Baum, Andrew J. Bunker, Stefano Carniani, Francesco D'Eugenio, Eiichi Egami, Tobias J. Looser, Pierluigi Rinaldi, Christopher N. A. Willmer

High Energy Astrophysical Phenomena

Unveiling Neutrino Nature with the Diffuse Supernova Background

The true nature of neutrinos–whether Dirac or Majorana–is a foundational, unresolved question. We demonstrate that the diffuse supernova neutrino background (DSNB) offers an untapped avenue to resolve this issue, provided neutrino magnetic moments are $\gtrsim 10^{-14}\mu_B$. The intense magnetic fields characteristic of a subset of collapsing massive stars can trigger resonant chirality flips. This flavor conversion physics alters DSNB fluxes in measurable ways that depend on the neutrino nature. A $20$~yr combined exposure at Hyper-Kamiokande loaded with gadolinium and JUNO can unravel this signature at $90\%$ ($99\%$) confidence if the fraction of magnetorotational events exceeds $12\%$ ($20\%$) of cosmic core collapses. This result holds independent of the mass ordering and establishes the DSNB as a critical gateway to unveiling the true identity of the neutrino.
Read the full abstract

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

Catastrophic Inflation in the Axiverse

In the landscape of compactifications of type IIB string theory on O3/O7 orientifolds of Calabi–Yau threefolds, we identify regions of the axion potential in specific geometries that support small-field inflation, under the assumption that the overall Calabi–Yau volume can be stabilized at certain tuned values. At the tuned volume multiple critical points of the axion potential merge in a “catastrophe”, creating a region with small gradient and small Hessian determinant — necessary conditions for slow-roll inflation. We describe an algorithm that identifies several explicit examples of these catastrophes that have small field range $f / M_{\rm Pl} \ll 1 $ yet inflate for thousands of efolds. Some examples correctly reproduce the observed spectral tilt but not the amplitude of scalar perturbations, while others reproduce the observed amplitude but with a tilt that is slightly too red. All our examples have negligible tensor perturbations and no significant non-Gaussianity generated during the slow-roll phase. Our results demonstrate that string theory axions can be the source of cosmic inflation without alignment or large field ranges, provided that the Kähler moduli can be stabilized at special values of the volume.
Read the full abstract

All authors Naomi Gendler, Oliver Janssen, Matthew Kleban, Cameron Norton

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

Blazar Boosted Dark Matter in IceCube

We study the sensitivity of IceCube to blazar-boosted dark matter in a fermionic dark matter model with a massive vector mediator coupling to quarks. To this aim, we compute the diffuse flux arising from a sample of 324 blazars with proton spectra inferred from multiwavelength observations, adopting conservative dark matter spike profiles around the central supermassive black holes and consistently accounting for attenuation effects during propagation through the Earth. The dark matter-nucleon scattering cross section is evaluated by including elastic, resonant single pion production, and deep inelastic contributions, with particular emphasis on resonant single-pion production channels in order to smoothly cover the transition between the elastic and deep inelastic regimes. Using IceCube neutrino data, we derive constraints on the parameter space of the model and show that this detection strategy can surpass the sensitivity of conventional direct-detection experiments for dark matter masses below $\sim 1$ GeV. We find that the signal is dominated by deep inelastic scattering and is therefore more sensitive to comparatively heavy mediators, while resonance processes provide a reduction of the event rate, reaching up to about $9\%$ near the experimental threshold. Our results demonstrate that IceCube constitutes a powerful probe of sub-GeV dark matter scenarios through the observation of blazar-boosted dark matter.
Read the full abstract

All authors Alberto M. Gago, Jaime Hoefken Zink, Joel Jones-Pérez, Gabriel D. Zapata

High Energy Astrophysical Phenomena

Cosmic Ray Diffusion and the Origin of Very High Energy Gamma-Ray Emission in Young Massive Stellar Clusters

The search for Galactic sources capable of accelerating cosmic rays (CRs) to PeV energies has advanced significantly in recent years. High-energy observatories such as LHAASO have detected extended gamma-ray halos around several sources, suggesting that CRs escape their acceleration sites through anomalously slow diffusion. Theoretical studies propose that magnetic mirror diffusion combined with pitch-angle scattering in turbulent flow can naturally suppress CR transport. Here, we first show how mirror diffusion combined with scattering suppresses cosmic-ray transport, leading to an energy-dependent diffusion coefficient $D(E)\propto E^{1/3}$. We then combine a 3D magnetohydrodynamic (MHD) simulation of a young massive stellar cluster (YMSC) with Monte Carlo CR propagation calculations (CRPropa). The model includes the background gas density, magnetic field, stellar blackbody and dust emission, the cosmic microwave background, and the Galactic interstellar radiation field. Using the YMSC W43 as a benchmark, we compare two CR injection geometries: a central source and a spherical shell representing the cluster's collective wind shock. We show that mirror+scattering diffusion $D(E)\propto E^{1/3}$, combined with a CR injection spectrum $E^{-2}$, reproduces the gamma-ray spectrum observed by Fermi and LHAASO. In contrast, stronger energy-dependent diffusion requires a harder CR injection spectrum, $\sim E^{-1.6}$, to match the data. The relative contributions of leptonic inverse-Compton and hadronic emission depend sensitively on the diffusion regime. Finally, the resulting spectra show little dependence on the CR injection location, aside from a lower injection luminosity in the central-source case. Overall, our results indicate that the observed gamma-ray emission is shaped primarily by the diffusive propagation regime rather than by the precise location of the CR source.
Read the full abstract

All authors Lucas Barreto-Mota, Elisabete M. de Gouveia Dal Pino, Gaetano Di Marco, Rafael Alves Batista, Stela Adduci Faria

Astrophysics of Galaxies

ATLAS. IV. A JWST+MUSE Demographic Study of Ly$\alpha$ Profiles in Little Red Dots

We present an initial demographic study of Ly$\alpha$ profiles in little red dots (LRDs) at $z=3$–9 using $R\sim1000$–4000 spectroscopy. Our sample consists of 8 LRDs observed in the VLT/MUSE Deep and Wide surveys and 23 LRDs observed with JWST/NIRSpec grating spectroscopy from JADES, CANUCS, and GO programs including SPURS. We identify Ly$\alpha$ emission in 5 MUSE LRDs and 9 JWST LRDs. Only two of them exhibit broad Ly$\alpha$ emission (FWHM $>1000\ \mathrm{km\ s^{-1}}$), both reported previously. The other Ly$\alpha$-emitting LRDs show narrow Ly$\alpha$ emission (FWHM $<1000\ \mathrm{km\ s^{-1}}$), with FWHMs mostly in the range 300–600 $\mathrm{km\ s^{-1}}$, comparable to or slightly larger than those of high-redshift star-forming galaxies (100–500 $\mathrm{km\ s^{-1}}$). We measure the fraction of broad Ly$\alpha$ emitters above a broad Ly$\alpha$ luminosity threshold of $L_{\mathrm{Ly}\alpha,\mathrm{broad}}=10^{42}\ \mathrm{erg\ s^{-1}}$, obtaining $0.10^{+0.12}_{-0.07}$ for LRDs, about five times higher than the $2\sigma$ upper limit of $<0.02$ for high-redshift star-forming galaxies. Although we reproduce the broad Ly$\alpha$ component reported in a previous stacking analysis of eight LRDs, albeit with a large uncertainty, we find no evidence for broad Ly$\alpha$ emission in either the larger JWST stack, after excluding the two individually detected broad Ly$\alpha$ emitters, or the higher-resolution MUSE stack. These results suggest that broad Ly$\alpha$ emission is not ubiquitous among LRDs. Instead, LRDs with broad Ly$\alpha$ emission appear to represent a rare population that may correspond to a particular evolutionary stage, potentially associated with unusually strong outflows or other distinctive physical conditions.
Read the full abstract

All authors Yuta Kageura, Masami Ouchi, Hiroto Yanagisawa, Makoto Ando, Yuichi Harikane, Tomokazu Kiyota, Minami Nakane, Yoshiaki Ono, Yui Takeda

Earth and Planetary Astrophysics

Hydrogen Engulfment into Sub-Neptune Cores through Magma Ocean Convective Surface Renewal

A convective surface renewal explanation for ingassing of hydrogen from H2-rich envelopes into molten cores of growing planets is presented, with particular application to sub-Neptunes. The ingassing of hydrogen occurs at a moving front comprising the upper boundary of the magma ocean that engulfs hydrogen from the envelope above as it moves. The mechanism of engulfment is surface renewal as the upper diffusive boundary layer is diluted by convection. Surface renewal results in ingassing even where the stabilizing buoyancy of the low-density hydrogen neutralizes the destabilizing thermally-driven buoyancy. It is ultimately driven by cooling. Beginning during the gaseous disk phase, when hydrogen is effectively available in infinite supply, sub-Neptunes can acquire enough hydrogen to permit global core-envelope equilibrium. The structures are therefore consistent with hydrogen-rich, supercritical magma oceans overlain by hydrogen-rich envelopes. The process of surface renewal has implications for other magma ocean - atmosphere interactions in general.
Read the full abstract

Astrophysics of Galaxies

The deepest color-magnitude diagrams for the benchmark open cluster NGC 2437 from Gaia and VVVX

Deep photometry of Galactic star clusters provides one of the most powerful tools for determining their physical properties. In particular, in low Galactic latitude regions that suffer from heavy extinction and crowding. NGC2437 is the most extended star cluster in the near-IR footprint of the VVV Extended Survey (VVVX), covering more than one degree on the sky. We aim to characterize its physical properties using Gaia DR3 in the optical and the VVVX in the near-IR. We use Gaia DR3 proper motions to select NGC2437 members in order to make optical and near-IR color-magnitude and color-color diagrams. We further exploited the newly constructed VVVX deep stack images to obtain the deepest near-IR color-magnitude diagram currently available for this cluster. We estimate the main physical parameters for NGC2437, including the mean parallax of 0.608 mas and PMs (-3.85, 0.41) mas/yr. The mean reddening E(J-Ks) = 0.059 mag and extinction of A_k = 0.034 mag for the cluster field, with no significant differential reddening spread. A distance modulus of 11.08 mag is estimated, equivalent to a distance of 1644 pc. This places NGC2437 at z=115 pc above the Galactic plane and at a galactocentric distance of 9.24 kpc. We measure the cluster structural parameters, obtaining a core radius of 10.79 arcmin. The estimated total absolute magnitudes are Mk = -4.91 mag and Mv = -3.70 mag. The cluster mean age is 350 Myr, using PARSEC-COLIBRI isochrones for solar metallicity. We measure a binary fraction of 28.6%. We also discuss the implications of the revised cluster parameters for the nearby open cluster NGC2425, the planetary nebula NGC2438, and the evolved OH/IR source OH231.8+04.2. The VVVX deep stacks increase the Ks photometric depth by 1.6 mag, nearly doubling the detected point sources and enabling significantly improved studies of stellar populations throughout the southern Galactic plane.
Read the full abstract

All authors Tali Palma, Matías Gómez, Dante Minniti, Alan Montecinos, Nicholas J.G. Cross, Javier Alonso-García, Daniel Majaess, Bruno Dias, Maren Hempel, Roberto K. Saito, Joyce B. Pullen

Astrophysics of Galaxies

OH Line Detections in Southern Galaxies of the IRAS Revised Bright Galaxy Sample

We present a systematic study of OH main-line emission and absorption in 186 southern galaxies from the IRAS Revised Bright Galaxy Sample, using archival MeerKAT snapshot data. OH features are detected in 38 galaxies, including eight with OH maser emission (three new) and 30 showing OH absorption, mostly unreported previously. Four absorption systems exhibit weak OH emission superposed on strong absorption. OH-emitting regions are generally more compact than the associated radio continuum. Most absorption profiles are well fit by two Gaussian components (1667 and 1665 MHz), with an average integrated line ratio of $\sim$1.5. LIRGs show an OH emission detection rate of ~13%, versus significantly lower rates in non-LIRGs. For sources with radio continuum flux densities >20 mJy, OH absorption detection rates reach ~36% (LIRGs) and ~27% (non-LIRGs), while no OH absorption features were detected among sources with lower radio continuum flux densities. This suggests that sufficient background continuum is likely an important factor for the detection of OH absorption. Detected OH emitters follow the empirical $L_{\rm OH}$–$L_{\rm FIR}$ relation, consistent with far-infrared pumping, while non-detections show upper limits below the relation. No significant differences are found between OH absorbers and non-detections in infrared luminosity or radio continuum compactness. Stacked spectra of non-detections reveal no significant OH features, suggesting that sensitivity and orientation alone do not fully explain the absence of absorption. In contrast, mid-infrared colors (e.g., W2–W3) and q_TIR differ between the two populations. OH absorption galaxies occupy an intermediate regime in L_HCN/L_CO between OH megamasers and non-detections, implying that OH absorption detectability is linked to dense molecular gas conditions, with extreme star formation potentially suppressing its occurrence.
Read the full abstract

All authors Han Zhao, Zhongzu Wu, Bo Zhang, Timur Mufakharov, Yongjun Chen, Zhiqiang Shen, Yulia Sotnikova

Cosmology and Nongalactic Astrophysics

How X-rays heat the IGM in different 21-cm simulation codes: a comparison between Licorice and Beorn

Any interpretation of the 21-cm signal of neutral hydrogen using Bayesian inference methods can only be as accurate as the underlying simulation code used to model the state of the intergalactic medium (IGM). 3D radiative transfer (RT) simulation codes may capture complex physics, but are computationally expensive and, therefore, faster, more approximate codes have been developed. To improve our understanding of the convergence of simulation codes in the 21-cm science community, we present a comparison of the X-ray heating of the IGM modelled in Licorice, a 3D RT simulation code, and Beorn, a 1D RT code. We use Beorn to process sources extracted from Licorice simulations, using the same physics of the sources, in order to obtain two versions of the temperature of the IGM heated by X-rays. We observe a good agreement between the luminosity fields, mean temperatures, and global 21-cm signal of the two setups, but discrepancies in the distribution of temperature and 21-cm signal, which result in a $\sim 30\%$ difference in the 21-cm power spectrum. We attempt to isolate the approximations that lead to these differences and find that common approximations used in 1D RT codes produce effects of that magnitude. Using an emulator of the Licorice power spectra in an MCMC pipeline, we translate these differences between power spectra into differences between posterior distributions over the astrophysical parameters. We observe a typical bias between 1D posteriors of $\gtrsim 1 \sigma$ (with a noise level corresponding to 100h of SKA observations).
Read the full abstract

High Energy Astrophysical Phenomena

XRISM reveals sloshing-driven gas motions in the core of Abell 2029

We investigate the velocity structure of the intracluster medium (ICM) in the core of the relaxed cool-core cluster Abell 2029 using XRISM Resolve spectroscopy. We analyze combined XRISM Resolve observations and divide the central region into several subregions. To account for photon mixing caused by the XRISM point spread function, we perform a spatial-spectral mixing analysis. We detect an ordered line-of-sight bulk-velocity gradient across the cluster core: the northern regions are blueshifted relative to the brightest cluster galaxy (BCG), while the southern regions are close to zero velocity or slightly redshifted. The maximum velocity difference is about $280~{\rm km\,s^{-1}}$. In contrast, the turbulent velocity dispersion is smaller, with measured values and upper limits of $\lesssim150~{\rm km\,s^{-1}}$, implying a non-thermal pressure fraction below $\sim2.5\%$. The velocity pattern is consistent with gas sloshing associated with the spiral structure seen in Chandra X-ray images. Averaged over all regions, the inferred turbulent heating rate is below the radiative cooling rate, indicating that turbulent dissipation alone is insufficient to offset cooling in the entire core. These results reveal that A2029 is not kinematically featureless: sloshing-induced bulk motions are present, while the observed line-of-sight velocity dispersion indicates only a limited contribution to pressure support and core heating.
Read the full abstract

All authors Yuusuke Uchida, Yuna Saito, Naomi Ota, Erwin T. Lau, Ming Sun, Eric D. Miller, Tommaso Bartalesi, Stefano Ettori, Kotaro Fukushima, Caroline Kilbourne, Lorenzo Lovisari, Kyoko Matsushita, Brian R. McNamara, Arnab Sarkar, Kazunori Suda, Irina Zhuravleva

High Energy Astrophysical Phenomena

Optical Spectroscopy of TeV-emitting BL Lac Candidates

TeV-emitting BL Lacertae (BL Lac) blazars have important implications for the study of jet phenomenology, particle acceleration, and ultra high-energy cosmic ray production. They also allow for indirect studies of the extragalactic background light, cosmic magnetic fields, and axion-like particles. The upcoming Cherenkov Telescope Array (CTA), a ground-based gamma-ray observatory, will expand the observed TeV-emitting BL Lac population and enhance the aforementioned fields of study. However, reliable redshifts are crucial for planning and interpreting observations with CTA and only about 50% of gamma-ray BL Lacs have spectroscopic redshifts. We performed medium resolution optical spectroscopy of 16 TeV-emitting BL Lac candidates with the SALT and MDM telescopes. We measured spectroscopic redshifts for the full sample, ranging from 0.059 to 0.4, including 5 new spectroscopic redshift measurements.
Read the full abstract

Astrophysics of Galaxies

The Galactic Centre G+0.633-0.0604 Molecular Cloud: A New Gold Mine for Astrochemistry

Astrochemistry is living a golden age, with more than a quarter of the ~350 molecules in the current interstellar census having been detected over the last three years. One of the sources driving this progress is the G+0.693-0.027 cloud, located in the northern part of the Galactic Centre Sgr B2 complex. In this contribution, we present the astrochemical characterisation of G+0.633-0.0604, a newly discovered chemically rich molecular cloud at the southern edge of Sgr B2. With an inventory of >120 species, G+0.633 provides robust second detections of several prebiotic molecules only reported towards G+0.693, establishing it as the first confirmed astrochemical twin of G+0.693 while demonstrating that the extraordinary chemistry of this cloud is not unique. Furthermore, G+0.633 offers an observational advantage over G+0.693 since it displays half narrower linewidths. Together, G+0.633 and G+0.693 form a unique benchmark pair for unveiling molecular complexity and prebiotic chemistry in the interstellar medium.
Read the full abstract

All authors D. San Andrés, V. M. Rivilla, L. Colzi, M. Sanz-Novo, S. Martín, I. Jiménez-Serra, S. Zeng

High Energy Astrophysical Phenomena

The Evolution of the ACIS Contamination Layer on the Chandra X-ray Observatory from 2010 to 2026

The Chandra X-ray Observatory (CXO) was launched over 27 years ago and has been delivering spectacular science over the course of its mission. The Advanced CCD Imaging Spectrometer (ACIS) is the prime instrument on the satellite, conducting over 90% of the observations. The CCDs operate at a temperature of $-$120$^\circ$C and the optical blocking filter (OBF) in front of the CCDs is at a temperature of approximately $-$60$^\circ$C. The surface of the OBF has accumulated a layer of contamination over the course of the mission, as it is the coldest surface exposed to the interior to the spacecraft. We have been characterizing the thickness, chemical composition, and spatial distribution of the contamination layer as a function of time over the mission. The contamination model has required several revisions over the course of the mission as the properties of the contamination layer have changed and our understanding of the layer has improved. In this paper, we evaluate the performance of the current contamination model (N0016 released in CalDB 4.12.3 on 16 December 2025) using the most recent calibration observations conducted from 2023 to 2026 by using the standard model spectrum for the supernova remnant 1E 0102.2-7219 (E0102) developed by the International Astronomical Consortium for High Energy (IACHEC), spectral data from the cluster of galaxies known as Abell 1795, and high resolution X-ray spectra of Mrk 421. This evaluation has been complicated by the decreasing observed counts at low energies, especially the O VII He$\alpha$ line complex and the O VIII Ly$\alpha$ line from E0102 which are no longer useful for this purpose. The analyses of the E0102, Abell 1795, and Mrk 421 data show that the current model of the contamination adequately predicts the additional absorption through mid-2026.
Read the full abstract

All authors Paul P. Plucinsky, Peter W. Ratzlaff, Akos Bogdan, Herman L. Marshall

Cosmology and Nongalactic Astrophysics

Cosmicflows-4 and the Cosmological Consistency of Large-Scale Motions

We study coherent flows in Cosmicflows-4 (CF4) and its constituent catalogs. We fit monopole and dipole moments in nonoverlapping radial bins using a log-distance radial-velocity estimator. For each sample, the Planck 2018 $\Lambda$CDM prediction accounts for the actual positions, statistical weights, velocity covariance, and measurement covariance. Across six bins extending to $300,h^{-1},\mathrm{Mpc}$, the All Individual sample is broadly consistent with $\Lambda$CDM, with dipole and joint deviations of $2.29\sigma$ and $2.36\sigma$. The strongest localized feature occurs at $120$–$160,h^{-1},\mathrm{Mpc}$, where $|V_{\rm dip}|=628\pm82,\mathrm{km,s^{-1}}$ and the dipole deviation reaches $3.40\sigma$. The excess is concentrated toward negative supergalactic X and depends strongly on the survey window. TFR and 6dFGS favor larger dipoles, while SDSS and SN favor smaller ones. Removing 6dFGS reduces the upper-bin dipole deviation to $1.56\sigma$, while removing SN raises it to $3.77\sigma$. Nevertheless, the SN and 6dFGS dipoles remain mutually consistent when their covariance is included. SDSS is consistent with $\Lambda$CDM through $160,h^{-1},\mathrm{Mpc}$ but shows a separate rise at larger radii whose significance depends on the treatment of the monopole and shared distance-scale covariance. Overall, CF4 does not show uniformly anomalous motion across its full radial range. Instead, it contains localized and window-dependent features that should be tested with selection-matched mock catalogs and correlated calibration uncertainties.
Read the full abstract

High Energy Astrophysical Phenomena

Reconnection-Driven Plasma Loading for Prompt Emission in Tidal Disruption Event Jets

A tidal disruption event (TDE) occurs when a star passes sufficiently close to a black hole (BH) to be disrupted by its tidal forces. Although several hundred TDEs have been identified to date, only a small fraction launch relativistic jets. These jets are thought to be Poynting-flux-dominated outflows powered by the Blandford-Znajek (BZ) mechanism. However, the origin of the jet plasma required to account for the observed prompt emission remains poorly understood and constitutes a key unresolved problem in the formation of relativistic TDE jets. In this Letter, we propose magnetic reconnection near the BH as a physical mechanism for injecting plasma into the BZ jet during the super-Eddington accretion phase. Once sufficient magnetic flux accumulates near the BH, an equatorial current sheet forms, in which magnetic reconnection produces high-energy gamma rays. We demonstrate that photon-photon interactions generate a sufficient particle supply to account for the observed prompt emission. We further show that our plasma loading scenario allows the jet to remain magnetically dominated out to the prompt emission region and that magnetic dissipation can consistently reproduce the observed prompt emission properties. Furthermore, our scenario predicts MeV and GeV gamma-ray flares, providing a direct means of testing this scenario with future gamma-ray observatories.
Read the full abstract

Solar and Stellar Astrophysics

Magnetoacoustic Portals in Quiet-Sun Fluxtubes Revealed by Chromospheric Spectropolarimetry with SUNRISE III/SCIP

Acoustic waves propagate into the chromosphere, contributing to energy transport and their dynamics. Their upward propagation is restricted to frequencies above the acoustic cutoff frequency. The magnetic field configuration plays a key role in determining whether acoustic waves can propagate upward because the cutoff frequency is reduced in regions where the field is inclined with respect to gravity, forming so-called magnetoacoustic portals. Previous studies linked magnetic fields and oscillations in quiet regions, but these analyses were based on photospheric magnetic field information, leaving chromospheric structure unconstrained. This study investigates the coupling between acoustic waves and magnetic topology using photospheric and, for the first time, chromospheric spectropolarimetry in a quiet region, obtained with the Sunrise Chromospheric Infrared SpectroPolarimeter (SCIP) aboard the Sunrise iii balloon-borne solar observatory launched in 2024. The SCIP sit-and-stare observations sampled magnetic features in which the line-of-sight field strength exhibits multiple sharp spatial peaks in the photosphere while becoming broader and weaker at two heights in the chromosphere, indicating expanding fluxtubes. The chromospheric velocity field in these fluxtubes exhibits strong 5-minute oscillations, while the surrounding regions show weak 3-minute oscillations. In these fluxtubes, sawtooth temporal velocity variations are associated with intensity enhancements, suggesting steepened shocks. Fluxtubes with low-frequency oscillations are identified not only in network regions but also in weak internetwork regions. These results provide observational evidence that fluxtubes expanding into the chromosphere act as magnetoacoustic portals, in both network and internetwork regions, allowing low-frequency waves to propagate upward and driving chromospheric dynamics via shocks.
Read the full abstract

All authors Takayoshi Oba, Yukio Katsukawa, Masahito Kubo, Yusuke Kawabata, Takuma Matsumoto, Ryohtaroh T. Ishikawa, Yoshihiro Naito, Toshifumi Shimizu, Hirohisa Hara, Fumihiro Uraguchi, Toshihiro Tsuzuki, Kazuya Shinoda, Tomonori Tamura, Yoshinori Suematsu, Jose Carlos del Toro Iniesta, David Orozco Suárez, Maria Balaguer Jiménez, Azaymi L. Siu-Tapia, Carlos Quintero Noda, Sami K. Solanki, Andreas Lagg, Achim Gandorfer, Pietro Bernasconi, Thomas Berkefeld, Tino L. Riethmüller, Alberto Álvarez-Herrero, H. N. Smitha, Bianca Grauf, Michael Carpenter, Alexander Bell, Valentín Martínez Pillet, Francisco Javier Bailén, Julian Blanco Rodríguez, Juan Sebastián Castellanos Durán, Edvarda Harnes, Johannes Hölken, Francisco A. Iglesias, Hanna Strecker, Dušan Vukadinović

Cosmology and Nongalactic Astrophysics

Cosmographic Reconstruction of the Quintessence Potential in Scalar-Tensor Gravity

We generalise the cosmographic reconstruction programme of Chakraborty, Dunsby & Scherrer to scalar-tensor gravity. Working in the Jordan frame with a general non-minimal coupling $F(\Phi)$, we derive exact closed-form expressions for the potential slope $\lambda_U$, the coupling slope $\lambda_F$, and their curvature parameters $\Gamma_U$ and $\Gamma_F$ in terms of the cosmographic parameters $(q,j,s)$, the Planck-mass running $\alpha_M = d\ln F/d\ln a$ and its derivatives. We express $\alpha_M$ and its derivatives in terms of a set of coefficients $g_n$ describing the time variation of the gravitational constant $G\propto F^{-1}$, casting all four quantities in fully observable form. We then evaluate the reconstruction against current data, including DESI~DR2 baryon acoustic oscillation measurements and Lunar Laser Ranging constraints on $\dot G/G$, propagating the cosmographic uncertainties through to the reconstructed potential. All expressions reduce analytically to those of Chakraborty, Dunsby & Scherrer in the minimally-coupled limit, which we verify symbolically. We find that the potential slope $\lambda_U$ is recovered at the few-tenths level while the curvature $\Gamma_U$ is essentially unconstrained. We further show that the positivity of the scalar kinetic term is equivalent to $w_{\rm eff,0}\ge-1$, so that the region of cosmographic parameter space closed to minimally-coupled quintessence is exactly the phantom region, and that the non-minimal coupling reopens it through a single coefficient $g_2$ of the gravitational tower. Because Lunar Laser Ranging forces $|g_1|\lesssim2\times10^{-4}$ while the reconstruction is sensitive to $g_2\sim O(0.1)$, the framework presupposes that $\alpha_M$ is passing through zero at the present epoch, as expected on the Damour–Nordtvedt least-coupling attractor.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Constraining Primordial Power Asymmetry from Galaxy Clustering and Peculiar Velocity Information

Primordial power asymmetry would probe departures from statistical isotropy, offering a clue to the physics of the primordial Universe. We investigate how such asymmetries can be tested with galaxy surveys and quantify how peculiar-velocity information can improve these tests. We develop a unified bipolar spherical harmonic (BipoSH) analysis of the auto- and cross-power spectra of galaxy density and line-of-sight peculiar velocity, and perform Fisher forecasts combining a Euclid-like spectroscopic galaxy survey with peculiar velocities reconstructed from Simons Observatory-like CMB maps through the kinetic Sunyaev-Zel'dovich effect. For dipolar asymmetry, we consider both scale-independent and scale-dependent modulations proportional to $k^{-0.5}$. Galaxy clustering provides most of the constraining power in both cases, while peculiar velocities add only modest information. For quadrupolar asymmetry, we account for the recently identified anisotropic galaxy-bias response and marginalize over its amplitude. We find that galaxy clustering alone suffers from a degeneracy between the primordial quadrupolar modulation and anisotropic galaxy bias. This degeneracy can be substantially broken by adding peculiar-velocity information: the density-velocity cross-spectrum provides complementary information to the galaxy auto-spectrum, while the velocity auto-spectrum constrains the primordial modulation independently of anisotropic galaxy bias. The velocity information becomes increasingly effective for more negative scale dependence, yielding constraints tighter than those from BOSS for scale dependences proportional to $k^{-1}$ and $k^{-2}$. Our results demonstrate that peculiar-velocity information provides a complementary avenue for testing primordial power asymmetry with large-scale structure, particularly for quadrupolar asymmetry in the presence of anisotropic galaxy bias.
Read the full abstract

All authors Keita Minato, Atsushi Taruya, Teppei Okumura, Maresuke Shiraishi

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

Physics of Circular Polarized Ion-Scale Waves in Hybrid Simulations of Alfvénic Fluctuations

Ion cyclotron waves (ICW) and fast magnetosonic/whistler waves (FMW) are fundamental electromagnetic modes at ion kinetic scales, yet their generation mechanisms and roles in plasma evolution remain poorly understood. We analyze a 2.5D hybrid simulation of broadband Alfvénic fluctuations, where the proton velocity distribution is modeled as a sum of two bi-Maxwellian components: a thermal core and a drifting beam. Using wavelet-based wave identification, bi-Maxwellian VDF fitting, and the PLUME linear dispersion solver, we find that ICW behave as linear modes. Growth is intermittent, occurring when core temperature anisotropy builds up, and is driven mainly by the core (the beam contributes negligibly). Poynting flux analysis shows that ICW are predominantly forward-propagating, with a net energy flux ratio of $+1$ across all frequencies, consistent with the initial condition. FMW present a stark contrast: PLUME solutions often yield very small (near-zero) linear growth/damping rates. The species decomposition breaks down when $|\gamma/\omega_r| \gtrsim 0.368$, indicating that linear theory predicts these waves to be strongly damped and not describable by linear eigenmodes. Nevertheless, FMW are clearly observed in the wavelet helicity spectrogram, indicating that they are generated by nonlinear processes (e.g., parametric decay or phase steepening) and persist despite linear damping. The net energy flux ratio for FMW is close to $+1$ at low frequencies but decreases at higher frequencies, yet never reaches zero (net energy flow remains forward). These results demonstrate that ICW are linear, core-driven waves that transfer energy to the plasma, while FMW are heavily damped, nonlinearly generated waves.
Read the full abstract

All authors Hai Yang Harry Qian, Trevor A. Bowen, Carlos A. Gonzalez, Nikos Sioulas, Alfred Mallet, Kristopher G. Klein, Daniel Verscharen, Stuart D. Bale

Astrophysics of Galaxies

A Reproducible Two-Boundary Kinematic Correction for Baryonic Rotation-Curve Reconstruction in an 84-Galaxy SPARC Benchmark

We present a reproducible computational validation and failure analysis of the empirical omega kinematic correction introduced by Flynn and Cannaliato (2025). The algorithm is deliberately minimal: one coefficient per galaxy is calculated from the innermost and outermost measured rotation-curve points and applied to the full observed radial profile before comparison with a baryonic reconstruction assembled from SPARC gas, disk, and bulge components. We preserve the predecessor's frozen 84-galaxy benchmark and publish its exact membership so that the transformation, not sample re-selection, is the object of validation. Without fitting the transformation to the baryonic residual, the primary maximum-disk reconstruction reduces the mean observed-baryonic discrepancy from 51.82 to 30.15 km/s across the frozen 84-galaxy benchmark. Bounded mass-to-light-ratio optimization further reduces the descriptive sensitivity-fit value to 25.45 km/s, while the simple Keplerian reference has a mean RMSE of 74.20 km/s in our earlier analysis [15]. Recalculation from the 84 per-galaxy records shows a resolved mass-to-light optimization benefit (delta RMSE > 0.05 km/s) in 53 galaxies and no resolved change in 31. Six galaxies do not beat the Keplerian reference; all six occur at Upsilon_max <= 0.111, whereas their omega values are not concentrated at the high end of the sample. We specify the complete deterministic workflow, native units, endpoint invariants, uncertainty propagation, and formula-level regression checks required to prevent grouping and sign errors. The complete 84-galaxy panel set, population-level error distributions, and failure diagnostics are retained as inspectable outputs. The result is a reproducible astronomical data-transformation benchmark rather than a proposed force law or replacement for dark matter or modified gravity.
Read the full abstract

Solar and Stellar Astrophysics

The B-Index as a Diagnostic of Cool Stars: Assessing Metallicity Dependence of TiO Absorption in the Visible Spectrum

Molecular absorption features constitute essential diagnostics of late-type stellar atmospheres, with titanium oxide (TiO) bands serving as sensitive tracers of effective temperature and magnetic activity. While near-infrared TiO indices are well established, visible-band diagnostics provide complementary constraints for stellar classification and parameter estimation. This study investigates the metallicity dependence of the B-index, a TiO absorption measure centered at 567 nm, through synthetic spectra generated with ATLAS9 and high-resolution observations from the HARPS spectrograph. The synthetic grid spans effective temperatures between 3500-4000 K and metallicities from [Fe/H] =-4 to +0.2, while the observational sample comprises 23 MK spectral types stars with effective temperatures between 3500-4000 K. Both theoretical and empirical analyses demonstrate that the B-index exhibits negligible sensitivity to stellar metallicity, with only marginal correlations detected. Therefore, there is no requirement to include this feature in the calculation of the B-index. The B-index's low sensitivity to metallicity reinforces its utility as an effective tool for spectral classification, stellar atmosphere modeling, and the study of magnetically active stars.
Read the full abstract

Instrumentation and Methods for Astrophysics

$\texttt{Aether.jl}$ : A High-Performance 3D MHD and Multifluid Dust Code Written in a Dynamic Language with an Interactive Human-AI Development Framework

We present $\texttt{Aether}$, a new finite-volume code for compressible hydrodynamics and magnetohydrodynamics, written in Julia and primarily designed for GPU systems. The code solves the MHD equations with constrained transport in Cartesian, cylindrical, and spherical-polar coordinates, using standard high-order Godunov methods. An arbitrary number of dust fluids can be coupled to the gas through stiff mutual drag. It was developed from scratch with interactive Human-coding agent workflow; the paper documents the framework of this workflow alongside the numerical methods. Performance-critical kernel is written through $\texttt{KernelAbstractions}$, and supports runs on CPUs and GPUs from multiple vendors. $\texttt{Aether}$ can be ran either from an interactive notebook or batch scripts, keeping prototyping, production runs, and analysis in a single language. We verify the implementation through a series of hydrodynamic, MHD, and dust tests. Although written in a dynamic language, $\texttt{Aether}$ achieves comparable or even higher single-GPU throughput than C++ code on the same hardware. In weak scaling on Frontier, parallel efficiency stays above $93\%$S on 4096 GCDs. These results show that a dynamic language now supports production astrophysical MHD simulations on exascale systems. $\texttt{Aether}$ and its Jupyter notebook example suite are publicly available.
Read the full abstract

Solar and Stellar Astrophysics

High-order Paschen emission from the quiet-Sun off-limb chromosphere

We report the detection of high-order hydrogen Paschen emission lines (Pa~15, Pa~16, and Pa~17) in the quiet-Sun chromosphere off the solar limb using the Chromospheric Infrared SpectroPolarimeter (SCIP) on board the Sunrise~iii balloon telescope. These lines reveal thread-like structures resembling spicules and exhibit systematically smaller Doppler velocities than Ca~II~854.2~nm, suggesting that they are optically thinner and more affected by line-of-sight averaging, especially near the limb. Non-LTE radiative transfer synthesis using the spherically symmetric one-dimensional code rhsphere reproduces the overall spectral properties. The observed ratios among three Paschen lines show systematic deviations from synthetic and theoretical results, suggesting that additional physical effects may influence the formation of high-order Paschen lines. The study demonstrates the potential of high-order Paschen lines as a new diagnostic of optically thin plasma in the off-limb chromosphere.
Read the full abstract

All authors Haocheng Yu, Yukio Katsukawa, Mingde Ding, Takuma Matsumoto, Sami K. Solanki, Julian Blanco Rodríguez, David Orozco Suárez, Masahito Kubo, Andreas Lagg, Achim Gandorfer, José Carlos del Toro Iniesta, Pietro Bernasconi, Thomas Berkefeld, Alex Feller, Tino L. Riethmüller, Alberto Álvarez-Herrero, H. N. Smitha, Bianca Grauf, Michael Carpenter, Alexander Bell, Valentín Martínez Pillet, Francisco Javier Bailén, Juan Sebastián Castellanos Durán, Edvarda Harnes, Johannes Hölken, Francisco A. Iglesias, Ryohtaroh T. Ishikawa, Yusuke Kawabata, Takayoshi Oba, Azaymi L. Siu-Tapia, Hanna Strecker, Dušan Vukadinović, Hirohisa Hara, Toshifumi Shimizu

Astrophysics of Galaxies

The galaxies' energy balance problem solved

We attempt to resolve the long-standing energy balance problem encountered by Radiative Transfer (RT) models, particularly in edge-on galaxies by incorporating a treatment for the clumpy structure of the Interstellar Medium (ISM). A subgrid approach is adopted, treating the quiescent dust clouds as pseudo-dust grains, with equivalent optical and thermal emission properties. Deriving key quantities such as the absorption, scattering and extinction cross-sections enables the virtualisation of the macroscopic clump into a microscopic pseudo-grain that can be included alongside the existing dust model constituents. The addition of the pseudo-grain results in a flatter extinction curve. A library of clump emission spectral energy distributions (SEDs) is constructed for radiation fields of various colours and intensities. The new clumpy model is applied to the edge-on galaxy NGC 891 and, for the first time, is able to achieve a good energy balance, simultaneously fitting both the submm and Near Infrared (NIR) data. The clumpy model is further applied to a small sample of seven galaxies of various inclinations, and the results are compared with those from the purely diffuse models. The clumpy models are characterised by a reduction in dust opacity, and therefore attenuation, compared to their purely diffuse counterparts. Thus, the maximum face-on optical depth in the $B$-band, ${\mathrm max}(\tau^{\mathrm f}_{\mathrm B})$, derived from the clumpy models is found to be lower by factors ranging from 1.3 to 2.8. Of the seven galaxies, two are found to be optically thick in their centres, two are found to be moderately optically thick, and three are found to be optically thin.
Read the full abstract

Astrophysics of Galaxies

The SAMI Galaxy Survey: Linking Tidal Features and Orbit Populations Using Schwarzschild Modelling

The evolution of angular momentum in galaxies is shaped by a combination of internal secular processes and external mechanisms such as mergers. Orbit-superposition based dynamical modelling provides a powerful means of linking the intrinsic orbital structures of galaxies to their global properties and merger histories. We construct Schwarzschild orbit-superposition models of massive ($\log(M/M_{\odot})>10$) SAMI galaxies using the DYNAMITE code, utilising deep KiDS photometry to accurately reproduce each galaxy's luminosity distribution. We find that the fractions of hot, cold, warm, and counter-rotating orbits all show significant correlations with the spin parameter proxy $\lambda_{R_e}$, with the strongest correlation arising from the combined hot plus counter-rotating fraction. When controlling for stellar mass and environment, we find that the fraction of hot and cold orbits show significant correlations with stellar age, whereas warm orbits do not. We further find that the lower values of $\lambda_{R_e}$ for young galaxies with shell merger features as compared to the full sample is driven by an excess of hot orbits and a deficit of cold orbits, with no dependence on warm orbits. We suggest that the kinematic transformation in this SAMI sample proceeds through stars transitioning directly from cold to hot orbits. As warm orbits are expected to arise from secular heating processes, these findings indicate that merger-driven heating is the dominant mechanism governing the redistribution of angular momentum and the reduction of rotational support in massive galaxies.
Read the full abstract

All authors T. H. Rutherford, J. van de Sande, S. M. Croom, A. Fraser-McKelvie, G. Santucci, G. van de Ven, S. M. Sweet, S. Brough, Y. Mai, J. Bryant, J. Bland-Hawthorn, M. Goodwin, J. Lawrence, N. P. F. Lorente

High Energy Astrophysical Phenomena

Detection of possible burst oscillation in the neutron star low-mass X-ray binary 4U 1323-62

Burst oscillations observed during thermonuclear X-ray bursts arise from asymmetric brightness patterns on the neutron star surface and provide a direct probe of the neutron star spin frequency. We present a detailed timing analysis of the neutron star low-mass X-ray binary 4U 1323-62 using 2024 observations with XMM-Newton and NuSTAR observatories. We identify nine thermonuclear X-ray bursts in the XMM-Newton/EPIC-pn data, along with eclipsing dips in the light curve. One of the XMM-Newton bursts exhibits a rare doublet structure. In addition, NuSTAR detects six bursts, four of which occur simultaneously with those observed with XMM-Newton. We identify a possible burst oscillation signal at $\sim$611.5 Hz in the XMM-Newton data. The strongest oscillation, detected during the primary burst of the doublet burst, reaches a maximum $Z_1^{2}$ power of $\sim35$. An analytical estimate accounting for the searched frequency and time intervals gives a significance of $\sim3.0\sigma$, whereas independent Monte Carlo simulations yield a more robust global significance of only $\sim2.4\sigma$. We therefore interpret the signal as a tentative detection of burst oscillation. The folded pulse profile in the 0.5-10 keV band is well described by a sinusoid, with a fractional rms amplitude of $\sim30\pm6$%. The oscillation frequency corresponds to a neutron-star spin period of $\sim$1.635 ms, suggesting that 4U 1323-62 may harbor a rapidly rotating millisecond neutron star.
Read the full abstract

High Energy Astrophysical Phenomena

X-ray Activity of the RS CVn-type Star {\sigma} Gem with the First-Year Observations of Einstein Probe

Context. Stellar flares are energetic events driven by the sudden release of magnetic energy in the stellar atmosphere. Studying these flares is crucial for understanding their impact on exoplanets, the circumstellar environment, and stellar evolution itself. The launch of the Einstein Probe (EP) offers a unique opportunity to systematically detect such events. Aims. We present a systematic analysis of the flaring activity of the active RS CVn-type binary {\sigma} Gem, utilizing the first-year monitoring data from the Wide-field X-ray Telescope (WXT) aboard EP. Our goals are to demonstrate the unique capability of EP in monitoring stellar X-ray activity and detecting flares, by identifying and characterizing extreme X-ray flares on {\sigma} Gem and estimating their occurrence rate. Methods. We developed a data-processing pipeline to select and extract EP-WXT observations, producing a background-subtracted, vignetting-corrected light curve. We employed the Bayesian Blocks method to detect significant flares in the long-term X-ray light curve. For each identified flare, we performed light curve and spectral fitting to derive the flare parameters. Results. Between October 2024 and April 2025, WXT detected 6 distinct flares from {\sigma} Gem. Their durations ranged from 21 hours to 3 days, with peak X-ray luminosities (0.5-4 keV) of 3.7 * 10^31 to 7.0 * 10^32 erg/s and total energies of 1.1 * 10^36 to 4.4 * 10^37 erg, placing them among the "superflare" class. Conclusions. Using {\sigma} Gem as a case study, we demonstrate an analysis process for flare detection and analysis with EP-WXT data, which provides new statistical constraints on its flaring behavior. Applying this methodology to the growing EP stellar archive promises to yield a vast sample of X-ray flares, which will significantly advance our understanding of stellar magnetic activity.
Read the full abstract

All authors Xuan Mao, Giuseppina Micela, Fabio Favata, Weimin Yuan, He-yang Liu, Huaqing Cheng

Astrophysics of Galaxies

How mergers shape galaxy morphology in the IllustrisTNG simulation

How galaxy mergers drive morphological evolution remains an open question. Traditional views hold that major mergers produce elliptical galaxies, while minor mergers form dispersion-dominated components such as galactic bulges. However, more recent work has challenged this simple picture, suggesting a more complex evolutionary scenario. In this study, we use the IllustrisTNG cosmological simulation to investigate how mergers shape galaxy morphology, across a broad range of galaxy masses and merger mass ratios. Our results show that the post-merger galaxy morphology is primarily determined by three factors: collision angle $\overline{\theta}$, cold gas fraction $f_\mathrm{cold\,gas}$, and pre-merger galaxy morphology $\mathrm{(B/T)_{*,pre}}$. Specifically, spiral-in mergers with large $\overline{\theta}$ increase the rotational support of the system, allowing the gas to settle into an extended disk. When the system is rich in cold gas, star formation within the newly formed gas disk can further strengthen the disk-dominated structure of the remnant. On the other hand, head-on mergers with small $\overline{\theta}$ typically disrupt ordered galactic motion, producing more dispersion-supported remnants. Overall, we interpret our results within a unified picture of merger-driven morphological transformation in galaxies.
Read the full abstract

Instrumentation and Methods for Astrophysics

On the Acceleration of Pulsar Timing computations using Normalising Flows and Parallelisation

Single-Pulsar Noise Analysis (SPNA) and Gravitational Wave (GW) searches done on Pulsar Timing Array (PTA) datasets have everlastingly suffered from the computational bottleneck arising due to high dimensionality and multi-modality of the PTA likelihood landscape, along with strong correlations amongst various single-pulsar noises and ensemble-level common noise processes. We addressed this outstanding issue by employing a Normalising Flow-based Preconditioned Monte-Carlo sampling technique implemented in the POCOMC package, for the first time on PTA-specific computations, and comparing the achieved acceleration with the widely used PTMCMCSAMPLER and DYNESTY packages. We further investigated the acceleration achieved via parallelisation over an increasing array of communicating nodes on a high-performance computing (HPC) resource, by employing the PARALLEL_BILBY architecture with DYNESTY. We tested the acceleration on realistic long baseline simulated datasets with SPNA and Common Red Noise (CRN) analysis. We found PARALLEL_BILBY to be the most efficient in parallelisation, achieving a runtime of ~10min and ~100min with 16 nodes for spatially uncorrelated and Hellings and Downs-correlated CRN searches, respectively. POCOMC outperforms in single node performance requiring only ~10h for correlated search. PTMCMCSAMPLER was found to be the least efficient. We envisage POCOMC to be of great importance for PTA analyses, without requiring any GPU or HPC support, while also performing ensemble-level GW searches within a manageable time span. These results have everlasting implications with increasing data volumes and need to incorporate more complicated models, which were otherwise beyond reach due to the associated computational costs.
Read the full abstract