Peer Revenue Play credits. Real papers.

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics

Monday 1 June

Earth and Planetary Astrophysics

Nereid as a Regular Satellite of Neptune

Matthew Belyakov, M. Ryleigh Davis, Ian Wong et al. (+2 more)

Nereid, Neptune's third largest moon, is considered to be a captured irregular satellite due to its highly eccentric orbit. However, among irregular satellites, Nereid is an outlier: it is the largest, the closest to its host planet, and the most eccentric. We present James Webb Space Telescope near-infrared spectroscopy of Nereid that demonstrates that its composition is inconsistent with its suggested captured origin. We then simulate Nereid's early orbital history subsequent to Triton's capture to demonstrate a plausible dynamical pathway for a regular satellite formed in-situ around Neptune to evolve to Nereid's present-day orbit. Based upon the available spectroscopic and dynamical evidence, we propose that Nereid is not a body captured from the Kuiper belt, but rather the sole surviving intact regular satellite of Neptune.
Read the full abstract

All authors Matthew Belyakov, M. Ryleigh Davis, Ian Wong, Konstantin Batygin, Michael E. Brown

Solar and Stellar Astrophysics

First Detection of HC5N in a Class II Disk around TW Hya

Steven C. Wampler, Ryan A. Loomis, Amina Diop et al. (+9 more)

Over the last decade of ALMA's operation the molecular inventory of protoplanetary disks has expanded rapidly, revealing a diverse set of nitrogen-bearing organics and carbon-chain molecules that trace both prebiotic chemistry and fundamental disk physics. Despite this progress, detections of larger species such as cyanopolyynes have remained limited, leaving larger carbon-chain chemistry in Class II disks largely unconstrained. Here, we report the first detection of HC5N toward the TW Hya protoplanetary disk, representing the largest cyanopolyyne identified to date in a Class II system. We derive a HC5N column density for two rotational transitions J = 41-40 and J = 37-36, N_T ~ 10^12 cm^-2 for assumed T_rot = 20-50 K and optically thin emission in LTE. We compare HC5N and HC3N formation mechanisms and analyze the HC3N/HC5N ratio. We use a chemical model to estimate the expected abundance and emitting layer of HC5N in a TW Hya-like disk. Although HC5N emission is spatially unresolved, measured column densities suggest an origin in the warm molecular layer where CN-based pathways are active. This detection extends the known carbon-chain chemistry in Class II disks and demonstrates that long cyanopolyynes can form and persist in planet-forming environments.
Read the full abstract

All authors Steven C. Wampler, Ryan A. Loomis, Amina Diop, L. Ilsedore Cleeves, Yuri Aikawa, Romane Le Gal, Viviana Guzman, Charles J. Law, Karin Öberg, Jamila Pegues, Richard Teague, Catherine Walsh

Astrophysics of Galaxies

Topology of molecular clouds in the PHANGS-JWST catalogue in relation to the physical characteristics of their galaxies

Ignacio Morales-Gil

Molecular clouds are affected by dynamical processes such as cloud merging and stellar feedback that alter their topology. For this reason, the molecular clouds in the PHANGS-JWST catalogue are studied under the light of topological data analysis. This study shows correlations between physical characteristics of galaxies and topological summaries of molecular clouds.
Read the full abstract

Instrumentation and Methods for Astrophysics

Neutrino Fingerprints: Image-Based Encodings of IceCube Events for CNN Direction Reconstruction

Floriano Tori, Brecht F. Verbeken and Vincent Ginis

Reconstructing the direction of incoming neutrinos in the IceCube Neutrino Observatory is an important problem in astrophysics. The public IceCube–Neutrinos in Deep Ice Kaggle competition provided 140 million simulated events to benchmark reconstruction techniques. To address this challenge from a novel perspective we introduce neutrino fingerprints compact $72 \times 72 \times 3$ images in which each pixel represents a single detector, with pulse timing and charge statistics encoded as color channels. This representation transforms sparse, irregular pulse data into dense images suitable for convolutional processing. Our ResNet18 model achieves a mean angular error of $1.10$ rad, indicating that convolutional networks trained on fingerprints rival more complex architectures while offering an effective, interpretable baseline for IceCube event reconstruction.
Read the full abstract

Earth and Planetary Astrophysics

One Transit Is All You Need: Detecting Exoplanets Through Learned Stellar Behaviour with EXOVEIL

Pratik Priyanshu

I present EXOVEIL, a transit detection system that learns what a star's brightness should look like and flags when reality disagrees. Unlike existing systems that require phase-folded input, EXOVEIL operates on raw flux time series and can detect planets that transit only once.A Transformer world model, trained on 16,499 Kepler light curves with transit-masked self-supervised learning, predicts expected stellar flux. A matched-filter detector with variance weighting extracts transit signals from the prediction residuals. A learned classifier (XGBoost) separates planets from false positives, achieving AUC 0.938 on Kepler DR25. Applied to single-transit injection-recovery, EXOVEIL recovers 32% of transits at 1000 ppm depth a task where all classification-based systems score 0% by construction. A blind search of 3,737 Kepler stars yields 179 new transit-like signals not present in the DR25 TCE catalogue, including 46 monotransit candidates. Applied withoutretraining to 47 confirmed TESS planets in the PLATO LOPS2 field, EXOVEIL achieves 100% recovery, demonstrating zero-shot cross-mission transfer. At PLATO's 25-second cadence, detection reaches 100 ppm – approaching the Earth-analog regime. I provide the first application of conformal prediction to transit detection (95.9% empirical coverage) and release the system as pip install exoveil with pretrained weights and a candidate catalogue.
Read the full abstract

Astrophysics of Galaxies

Reduced incidence of little red dots at z < 3 based on number density and halo mass evolution

Chenxuan Zhang, Huanian Zhang, Qingwen Wu et al. (+2 more)

An intriguing puzzle in extragalactic astronomy is the scarcity of Little Red Dots (LRDs) at $z < 3$, compared to their higher abundance at earlier cosmic epochs. To investigate this, we measure the overdensity for 98 specpically confirmed LRDs at $3<z<7$, and find that these LRDs predominantly reside in under-dense regions at $z > 4$ but shift to more typical galaxy environments at $z \sim 3.5$. Concurrently, cross-correlation analyses show that their dark matter halo masses grow rapidly, from $\lesssim 10^{10.1} \, M_{\odot}$ at $z \sim 7.5$ to $\sim 10^{11.3} \, M_{\odot}$ at $z \sim 3.5$, approaching the halo masses of normal galaxies at lower redshifts. Applying an empirical stellar-to-halo mass scaling relation, we find that LRDs still host over-massive black holes relative to their stellar masses at $z > 4$, yet converge toward the local BH-stellar mass relation at lower redshifts. The coherent evolution of LRDs' large-scale environments and halo masses toward those of normal galaxies provides a plausible explanation for their declining abundance at $z < 3$, even though the underlying small-scale physical mechanisms remain elusive.
Read the full abstract

All authors Chenxuan Zhang, Huanian Zhang, Qingwen Wu, Luis C. Ho, Jian-Min Wang

High Energy Astrophysical Phenomena

Elemental cosmic ray spectra reveal two populations of Galactic sources and an immediate transition to an extragalactic component after the knee

Timur A. Dzhatdoev and Anatoly A. Semenov

The energy spectra for individual elements and/or for groups of elements in cosmic rays (CR) in the energy range between 100 $\times$ Z GeV and 10$^{3}$ $\times$ Z PeV (where Z is the charge number of the nucleus) have a number of features, including two steepenings ("knees") with the rigidity-dependent energies $E_{k1} \approx$ 15 $\times$ Z TeV and $E_{k2} \approx$ 3 $\times$ Z PeV and three hardenings ("ankles") at $E_{a1} \approx$ 500 $\times$ Z GeV; for protons $E_{a2-p} \approx$ 150 TeV and $E_{ a3-p} \approx$ 50 PeV. While the values of $E_{a1}$ for different nuclei are rigidity-dependent, the values of $E_{a2}$ (and probably of $E_{a3}$) are not: $E_{a2-He} \approx$ 1 PeV for Helium. The recent advances in precision measurements of the elemental CR spectra in the DAMPE and LHAASO experiments, and, to some extent, in IceTop and other experiments, make it possible, for the first time, to clarify the origin of the aforementioned spectral features. We show that the elemental CR spectra are reasonably well described with a sum of three components: 1) a low-energy Galactic component with a convex spectral shape reflecting the accelerated particle spectrum in the source; this component peters out after the TeV knee, 2) a high-energy Galactic component including the PeV knee, and 3) an extragalactic component. There is no need for any third, additional component of Galactic cosmic rays in the energy range between 10 PeV and 1 EeV.
Read the full abstract

Instrumentation and Methods for Astrophysics

Bridging the UV Gap: The HST Ultraviolet Foundation for Star Formation Science in the Era of Roman, Euclid, and HWO

F. Z. Majidi, A. Bayo, K. Biazzo et al. (+4 more)

As we enter the 2030s, the astronomical landscape will be dominated by large-scale infrared (IR) and optical surveys led by JWST, Euclid, and the Nancy Grace Roman Space Telescope. While these facilities provide unprecedented views of the dusty environments of nearby star-forming regions, they are fundamentally limited in their ability to probe the high-energy physics of accretion, magnetospheric activity, and disk photoevaporation. This white paper argues for the critical continued use of the Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) and Cosmic Origins Spectrograph (COS) to bridge the "UV Gap." We demonstrate that UV spectroscopy is the only direct method for characterizing the feedback mechanisms that determine planet habitability and stellar maturation, serving as a mandatory scientific bridge toward the Habitable Worlds Observatory (HWO). The study of star formation stands at a critical intersection of multiple scientific disciplines, linking the high-energy physics of stellar birth to the chemical evolution of protoplanetary disks and the eventual habitability of exoplanets. As such, it represents one of the most compelling and essential science cases for the continued allocation of HST resources. Ensuring that HST provides high-resolution UV spectroscopic data now is a fundamental requirement for the success of future flagship missions, as these data provide the unique physical context that infrared observations alone cannot achieve.
Read the full abstract

All authors F. Z. Majidi, A. Bayo, K. Biazzo, J. M. Alcalá, K. France, E. Gaidos, M. G. Guarcello

Astrophysics of Galaxies

Signature of Bursty Star Formation in the High-Redshift Galaxies Detected with JWST

Rupam Sarkar and Saumyadip Samui

Recent JWST observations reveal an unexpectedly slow evolution in ultraviolet luminosity functions (UV LFs) at redshifts $z > 10$. To investigate this phenomenon, we develop a semi-analytical model of the UV LF, calibrated against well-constrained measurements at $z \sim 2-10$. Our analysis identifies a transition in star formation modes across cosmic epochs: at $z \lesssim 5$, a longer characteristic star formation timescale with nearly constant star formation efficiency ($f_\star$) dominates, whereas at $6 \lesssim z \lesssim 10$, shorter timescales prevail without requiring an increase in $f_\star$. For $z > 10$, the slow UV LF evolution is best explained by a shift toward even shorter star formation timescales without changing the star formation efficiency. Dust-free conditions or a top-heavy initial mass function (IMF) alone cannot reproduce the observations at $z\sim 14$. By combining UV LF with stellar mass estimates from Prospector-based SED fitting, we try to break degeneracies between IMF variations and star formation histories. Our results indicate that evolving star formation timescales rather than IMF or dust changes are the primary drivers of the observed high-redshift UV LF evolution, reflecting changing physical conditions during the earliest phases of galaxy assembly. Additionally, we show that moderate AGN activity could further boost UV luminosities at $z \sim 14$, potentially explaining the observed UV LF without changes in stellar parameters.
Read the full abstract

Earth and Planetary Astrophysics

Formation of multiple dust rings and gaps in protoplanetary discs by a single migrating planet II: radiative discs and observational signatures

Kim M. Weiskopf, Sören C. Meiners, Alexandros Ziampras et al. (+1 more)

Dust structures in protoplanetary discs have been widely observed and their creation remains an active field of research. Several possible origins have already been explored, including magneto-hydrodynamics, shadows and planets-disc interactions. The goal of this paper is to investigate whether a single migrating planet in a low-viscosity disc, including radiative processes, is capable of generating observable dust structures. We aim to examine both the lifetime of such structures and potential asymmetries within them. We perform a set of high-resolution, two-dimensional hydrodynamic simulations of migrating planets using three different equations of state: isothermal, constant $\beta$-cooling and an adaptive $\beta$ model. Dust is included in all simulations and the resulting dust density profiles are then post-processed to create radiative transfer images. For all equations of state considered, the planet undergoes one or several migration jumps, each producing dust rings and gaps. The lifetime of these structures depends on the phase of slow migration preceding and occurring between jumps, but in all cases they remain visible for at least 400 kyr. We find that cooling has a deciding effect on the migration behaviour and the number of jumps, but no measurable influence on the lifetime of the dust structures. The structures exhibit relatively few asymmetries, and large-scale vortices persist for an average of only 90 kyr. Our models highlight the capacity of planets to open multiple gaps while migrating, and stress the importance of a realistic cooling model. Care should be taken when interpreting and comparing such models directly to observations.
Read the full abstract

All authors Kim M. Weiskopf, Sören C. Meiners, Alexandros Ziampras, Cornelis P. Dullemond

Solar and Stellar Astrophysics

Precision constraints on stellar physics from main sequence detached eclipsing binaries

Mitchell T. Dennis, Harry Desmond and Jeremy Sakstein

We present a Bayesian framework to constrain {mass ($M$), metallicity ($Z$), stellar age ($\tau$), and the convective mixing length parameter ($\alpha_{\rm MLT}$)} in main-sequence (MS) detached eclipsing binaries (DEBs). These systems provide precise values of stellar mass and radius, offering stringent tests of stellar evolution models. We combine these with broadband magnitudes in the $B$ and $V$ bands and Gaussian priors on spectroscopic mass and metallicity, and perform Markov Chain Monte Carlo inference using a fast machine-learning surrogate for one-dimensional stellar evolution models computed with Modules for Experiments in Stellar Astrophysics. To make this approach computationally feasible, we implement an active learning strategy that adaptively selects new stellar models to evaluate, concentrating training data in regions of parameter space where the surrogate is most uncertain. Applying this framework to 38 stars in DEB systems, we recover ages more precise than previous isochrone-based determinations and obtain bounds on $\alpha_{\rm MLT}$ for a subset of lower-mass stars ($M \lesssim 1.5 M_\odot$), where convective envelopes provide sensitivity to the mixing length parameter. For several stars, the inferred $\alpha_{\rm MLT}$ values lie below the Solar-calibrated value, supporting previous indications that a universal mixing length parameter may not adequately describe convection across the main sequence. The active learning methods developed here provide a scalable route to Bayesian inference with stellar evolution models, with clear applications to additional stellar physics parameters and other precisely characterized stellar systems.
Read the full abstract

Astrophysics of Galaxies

The effects of yields from binary massive stars as functions of metallicity

E. Pepe, E. Spitoni, F. Matteucci et al. (+1 more)

Massive stars in binary systems that undergo mass transfer during their lifetime have a different evolution from that of single stars, possibly affecting their chemical yields. While massive stars produce most of the metals in the Universe, only few studies have investigated the effects of massive binary stars on the chemical evolution of the Milky Way. Following the most recent studies on massive binary-stripped star yields as functions of metallicity, we aim at improving previous results based on single-metallicity model grids. Here, by adopting a detailed model of chemical evolution for our Galaxy, we compute the evolution of 22 chemical species including C, N, O, $\alpha$-elements and Fe-peak elements, adopting novel prescriptions for single and binary massive star yields. Our main results can be summarised as follows: (i) consistently with previous predictions, we observe very small differences in both the predicted solar abundances and [X/Fe] vs [Fe/H] relations even when including massive binary yields depending on metallicity; (ii) when adopting the new set of stellar yields for massive single stars, as computed by Farmer et al. (2023), we are able to reproduce both the K solar abundance as well as the [K/Fe] vs [Fe/H] relation, without invoking ad hoc assumptions on nucleosynthesis prescriptions; (iii) our model adopting Farmer's yields both for single and binary massive stars is able to better reproduce the [X/Fe] versus [Fe/H] relation for both Mg and Ca, as compared with standard nucleosynthetic yields adopted in chemical evolution models; iv) we find that no models can well reproduce the [C/Fe] and [Ti/Fe] vs [Fe/H] when adopting the new yields as functions of metallicity.
Read the full abstract

All authors E. Pepe, E. Spitoni, F. Matteucci, M. Palla

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

Probabilistic Causality from Graviton Fluctuations

Giordano Cintia, Federico Piazza and Samuel Ramos

We compute the commutator of a scalar field minimally coupled to gravity at leading order in $G_N$. The commutator is operator-valued, with terms involving derivatives of Dirac deltas supported on the Minkowski light cone. When evaluated on classical/coherent graviton states, these terms “bend" the support of the commutator in precisely the way required to recover standard causality on a classical curved spacetime. However, these terms are also associated with a variance and are thus a source of uncertainty in the causal relations between events. We quantify this effect for a thermal state of gravitons at temperature $T$. In the large-time limit, the uncertainty of the lightcone position is given by $$ \text{Var}(\delta x) =\frac{4G_NTt}{3}\, , $$ where $\delta x$ is the displacement from the classical lightcone at $|\vec x| = t$. We obtain this result by subtracting a universal vacuum contribution that is logarithmically UV divergent and subleading at late times.
Read the full abstract

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

Deconstructing the Extra-Dimensional Axion

Shihwen Hor, Yuichiro Nakai, Motoo Suzuki et al. (+1 more)

We present a four-dimensional deconstruction of the extra-dimensional axion arising from a $U(1)$ gauge theory in a five-dimensional orbifold, where the axion is identified with the Wilson line of the $U(1)$ gauge field and its coupling to QCD is generated by a 5D Chern-Simons (CS) term. We construct the corresponding 4D moose (quiver) gauge theory with link scalar fields, in which the axion emerges as a collective pseudo-Nambu-Goldstone boson. The axion-gluon coupling is described by a gauged Wess-Zumino-Witten term, providing the 4D counterpart of the 5D CS term. We further analyze non-perturbative effects from zero-mode and “fractional” instanton configurations. While the latter is exponentially suppressed in the regime corresponding to the 5D description, ensuring consistency with the higher-dimensional picture, we point out that this suppression can break down for smaller instantons whose inverse size exceeds the 5D cutoff scale, leading to a potentially significant effect. We also study axion potentials induced by bulk matter fields and boundary-localized symmetry-breaking operators, reproducing the characteristic nonlocal suppression associated with propagation in the extra dimension. Our construction provides a renormalizable 4D framework with a transparent understanding of the axion shift symmetry and its quality.
Read the full abstract

All authors Shihwen Hor, Yuichiro Nakai, Motoo Suzuki, Junxuan Xu

High Energy Astrophysical Phenomena

Accretion of Primordial Black Holes in Stellar Interiors

Matteo Cantiello, Ore Gottlieb, Cameron Norton et al. (+2 more)

We study spherical accretion onto primordial black holes (PBHs) embedded in the core of a solar-type star. We compute the radiative efficiency self-consistently for the first time across the optically thin range ($10^{-16.5}$-$10^{-10}M_\odot$) with time-dependent simulations, and follow the growth up to $10^{-2}M_\odot$ using an analytical photon-trapping prescription above $5\times 10^{-13}M_\odot$. Near the Schwarzschild radius ($r_{\rm S}\sim 10^{-11}$cm for a $10^{-16}M_\odot$ PBH), gas compressed to $T\sim 10^{11}$K radiates through microphysical processes that fundamentally alter the classical adiabatic Bondi solution. We solve the time-dependent spherical Euler equations with an implicit cooling source term, determining $\dot M$, $\eta = L/\dot M c^2$, and the flow structure self-consistently. We identify three regimes for spherical accretion: a Hot Bondi regime ($M_{\rm BH}\lesssim 10^{-14}M_\odot$) in which bremsstrahlung cooling is dynamically negligible; a bremsstrahlung-cooling regime ($10^{-14}$-$5\times 10^{-13}M_\odot$) driving the flow toward isothermal with $\eta\approx 10^{-2}$; and a photon-trapping regime above $5\times 10^{-13}M_\odot$, in which the Bondi sphere is optically thick and the accretion rate remains close to the Bondi value. Cooling enhances $\dot M$ by a factor of $\sim$2-7, keeping growth super-exponential throughout the spherical regime. The radiative efficiency is an order of magnitude lower than previously assumed, and the critical initial PBH mass required to consume a solar-mass star within a Hubble time is $M_{\rm 0,crit}\sim 10^{-16}M_\odot$.
Read the full abstract

All authors Matteo Cantiello, Ore Gottlieb, Cameron Norton, Matthew Kleban, Ken Van Tilburg

High Energy Astrophysical Phenomena

Can current models predict the local black hole merger rate?

Lumen Boco, Michele Bosi, Cecilia Sgalletta et al. (+2 more)

After four observational runs, the Ligo-Virgo-Kagra collaboration estimated a local binary black hole (BBH) merger rate density of $R_{0,\textrm{LVK}}\simeq 14-26\,\textrm{Gpc}^{-3}\,\textrm{yr}^{-1}$ within the 90% credible interval. Some previous studies already pointed out that, when a realistic evolution of the metallicity-dependent cosmic star formation rate density (SFRD) is adopted, theoretical models predict a local BBH merger rate density that exceeds the observed value by at least a factor of $\sim 10$. In this paper, we confirm and strengthen this claim by constructing an empirical model for the SFRD and metallicity evolution that includes a correction accounting for iron abundance. The adopted metallicity relation is flexible, enabling us to bracket the wide range of observational uncertainties. We show that, even under the most conservative assumptions regarding both the SFRD and the metallicity relation, the local BBH merger rate density is overestimated by a factor $> 10$. Attempts to reconcile the predicted and observed merger rates by modifying only the metallicity-dependent SFRD would require unrealistically high metallicities ($Z>Z_\odot$) even in low-mass galaxies at high redshift. This finding indicates that revisions to the treatment of stellar and binary evolution are necessary to achieve consistency between theoretical predictions and observations. We suggest that even a modest steepening of the delay-time distribution could help alleviate this tension.
Read the full abstract

All authors Lumen Boco, Michele Bosi, Cecilia Sgalletta, Amedeo Romagnolo, Michela Mapelli

Astrophysics of Galaxies

A Cosmic Archipelago of lensed metal-poor galaxies at $z\sim6$

A. Bolamperti, M. Messa, A. Zanella et al. (+21 more)

The Cosmic Archipelago is an ensemble of galaxies, strongly lensed by the cluster MACSJ0416, showing extreme physical properties at $z\sim6.14$. We combine JWST/NIRCam with deep VLT/X-Shooter and JWST/NIRSpec IFU to perform a joint spectrophotometric analysis from the far ultraviolet to red optical rest-frame. We focus on CA4, a UV-faint ($M_{UV}=-17.7$), compact ($r_e=81\pm11$ pc) galaxy at $z=6.1446$, magnified by a factor $\mu=3.73$. CA4 is a young, low-mass ($M_\star =4.3\times10^6$ M$\odot$), star-forming (${\rm SFR}=0.46$ M$\odot$/yr), and metal-poor ($Z\sim0.02$ Z$\odot$) galaxy, and an efficient producer of ionizing photons ($\log(\xi_{ion}/{\rm erg^{-1} Hz})\sim25.5$). Its properties place CA4 at the poorly explored interface between massive stellar clusters and dwarf galaxies during the epoch of reionization. Moreover, CA4 shows large Ly$\alpha$ ($f_{esc}^{\rm Ly\alpha}\sim43\%$) and Lyman-continuum ($f_{esc}\sim47\%$) escape fractions, consistent with its small Ly$\alpha$ velocity offset ($\Delta v\sim100$ km/s) and extremely blue UV-continuum slope ($\beta=-3.10$). These characteristics suggest that such UV-faint, metal-poor galaxies may contribute significantly to cosmic reionization. We also confirm five additional systems at the redshift of the Cosmic Archipelago, magnified by factors up to 12.5. They are all young (mass-weighted ages $<11$ Myr) and metal-poor ($Z<0.05$ Z$_\odot$), spanning a wide range of stellar masses and SFRs. Given the large number of these bursty star-forming galaxies in a small cosmic volume, we estimate that the currently known members of the Cosmic Archipelago result in a significant overdensity at $z\sim6$ ($\Delta z\sim0.08$), with $\delta_{gal}=12.3^{+6.6}_{-4.6}$. These results highlight the Cosmic Archipelago as an unprecedented laboratory for studying the earliest groups of low-mass, low-metallicity galaxies during the epoch of reionization.
Read the full abstract

All authors A. Bolamperti, M. Messa, A. Zanella, E. Vanzella, P. Bergamini, F. Loiacono, A. M. Koekemoer, J. Vernet, R. A. Windhorst, A. Adamo, F. Annibali, F. Calura, M. Castellano, J. M. Diego, C. Grillo, M. Gronke, E. Iani, M. Meneghetti, A. Mercurio, K. Nakajima, S. Ravindranath, M. Ricotti, P. Rosati, H. Yan

Cosmology and Nongalactic Astrophysics

Beyond Self-Similarity: Reconciling X-Ray Scaling Relations in Galaxy Clusters and Groups

S. Ettori

Scaling relations hold among observed quantities that describe the thermodynamic properties of the gas in galaxy clusters and groups. However, observed data show systematic departures from the self-similar model's baseline predictions, particularly in lower-mass systems. I show that the observed departures from self-similar predictions can be efficiently described by two physical quantities modeled with power laws: the gas mass fraction ($f_g \sim T^{f_1} E_z^{f_z}$) and the temperature variation ($f_T \sim T^{t_1} E_z^{t_z}$). Using a large variety of published X-ray scaling relations, this study proceeds with an MCMC-based meta-analysis to constrain the temperature- and redshift-dependence of the meta-parameters $f_g$ and $f_T$ to calibrate the model. These calibrations indicate that, while the gas mass fraction ($f_g$) does not show significant evolution with cosmic time ($f_z = -0.11 \pm 0.03$), it decreases significantly with decreasing halo mass ($f_1 = 0.50 \pm 0.01$). On the other hand, the temperature variation shows a mild positive increase with both mass and redshift. Overall, modeling the departures from the self-similar model with $\{f_g, f_T\}$ drastically improves predictive accuracy, reducing the number of scaling relations in $>3\sigma$ ($>5 \sigma$) tension from 49 (36) percent under the self-similar scenario to just 11 (3) percent (four and one out of 39, respectively) that might be identified for their peculiarity. Moreover, the modelization through the generalized form allows me to present an extended discussion of the expected slopes and redshift evolution for several X-ray scaling laws, including the new quantity $Y_{LGT0} = L^{-1} M_g^2 T^{1/2}$, a proxy for the cluster's volume which does not depend on $f_g$ and $f_T$ by construction, and is predicted to relate directly to the mass without any redshift evolution: $M \sim Y_{LGT0} f_g^0 f_T^0 E_z^0$.
Read the full abstract

Astrophysics of Galaxies

The Radio-FIR Correlation in the Context of Deep Radio Source Counts

Tirth D. Surti, Vikram Ravi, Allison Matthews et al. (+1 more)

Increasingly deep, confusion-limited radio surveys have pushed direct radio source-count measurements down to tens of $\mu$Jy at 1.4 GHz. Confusion-noise P(D) analyses extend the statistical counts to below $1\,\mathrm{\mu Jy}$. Radio source counts have allowed for constraints on the radio-derived star formation rate density (SFRD) history through models of the backwards evolution of the local radio luminosity function, using the radio-FIR correlation, $q \propto \log(L_{\mathrm{FIR}}/L_{1.4})$, to convert radio luminosities to FIR luminosities and hence star-formation rates. Recent deep radio source counts from MeerKAT suggest a potential tension in the SFRD history between radio and UV/IR measurements at $1\lesssim z\lesssim 2$. This corresponds to a $\gtrsim 3\sigma$ discrepancy between the predicted and measured source counts near the star forming galaxy source count $S^2n(S)$ peak of ${\sim}30\,\mathrm{\mu Jy}$ under both a pure luminosity (PLE) and combined luminosity and density evolution (LADE). We consider what the requirement of agreement between radio source counts and the observed UV/IR SFRD indicates about the redshift evolution of the radio-FIR correlation and its intrinsic scatter out to $z=3$. We introduce a radio-luminosity based parameterization to $q_{\mathrm{FIR}}(z)$ based on changing thermal radio fractions alone that agrees with the observed stellar-mass dependent $q_{\mathrm{FIR}}(z)$ better than a non-evolving or decreasing $q_{\mathrm{FIR}}(z)$. Despite this, we find that a decreasing $q_{\mathrm{FIR}}(z)$ at fixed radio luminosity provides better agreement between source counts and the observed SFRD, while a $q_{\mathrm{FIR}}(z)$ that breaks down due to cosmic ray losses requires an intrinsic scatter up to $\sigma_q\approx 0.3\,\mathrm{dex}$.
Read the full abstract

All authors Tirth D. Surti, Vikram Ravi, Allison Matthews, Viviana A. Rosero

Cosmology and Nongalactic Astrophysics

Constraints on the gravitational potential from DESI DR2 BAO and its implications for the local void scenario

Indranil Banik, José Antonio Nájera and Harry Desmond

We constrain the difference in gravitational potential between our location and sources at $z \gtrsim 0.3$ using datasets at those redshifts. Our motivation is that the Hubble tension might be caused by a local void, as suggested by galaxy number counts. This would increase the redshift through outflow and gravitational redshift (GR). Only the latter is important at high redshift, where a void contributes a fixed additional GR contribution of $z_0$ due to our location on a potential hill. This $z_0$ model has various subtle effects that were not previously considered, including a hotter CMB and reduced BAO scale $r_{\rm d}$. We test whether $z_0$ can have the previously expected value of 0.84%, which was based on fitting void parameters to galaxy number counts and local $H_0$ measurements. Combining BBN, CMB, BAO, and CC datasets at $z > 0.5$, we find that $z_0 = -0.4^{+0.8}_{-0.9}\%$, which rises to $-0.1 \pm 0.7\%$ when extending our analysis down to $z > 0.29$. Although the results prefer the standard value of $z_0 = 0$, the best-fitting model with $z_0 = 0.84\%$ fits the data almost as well as $\Lambda$CDM, with $\Delta \chi^2 < 2$. We find that $\Lambda$CDM faces a $3.07\sigma$ BAO anomaly in the standard $(H_0 r_{\rm d}, \Omega_{\rm m})$ parameter space, where different regions are preferred by BAO and non-BAO datasets from $z > 0.29$. Fixing $z_0 = 0.84\%$ reduces this to $2.79\sigma$. This suggests that a local void large enough to solve the Hubble tension cannot be ruled out by higher-redshift datasets despite its novel impacts on them.
Read the full abstract

Astrophysics of Galaxies

Lyman-alpha Pressure Strongly Enhances Pre-Supernova Feedback at Cosmic Dawn: The First Multi-Dimensional Lyman-alpha Radiation Hydrodynamics Simulations

Olof Nebrin, Aaron Smith, Garrelt Mellema et al. (+2 more)

The dynamical role of Lyman-$\alpha$ (Ly$\alpha$) radiation pressure feedback has been debated for nearly a century, with recent analytical and 1D numerical studies highlighting its potential dominance over other stellar feedback processes at Cosmic Dawn. Despite this, no multi-dimensional Ly$\alpha$ radiation hydrodynamics (RHD) simulations have been performed to date. In this paper, we present the first 2D Ly$\alpha$ RHD simulations using Lydion, an RHD code with a novel M1 moment method for Ly$\alpha$ transfer, and self-consistent dust dynamics. Lydion yields a $\sim \mathcal{O}(100) \,\times$ speed-up compared to Monte Carlo radiative transfer in simple benchmarks, making 2D Ly$\alpha$ RHD feasible. We perform simulations of star clusters and isolated stars embedded in dense, metal-poor ($Z/Z_\odot \leq 0.01$) clouds, and find that Ly$\alpha$ feedback dramatically boosts outflows and dominates over feedback from direct and infrared radiation pressure. Ly$\alpha$ leakage through lower-column density channels, Doppler shifts, and Ly$\alpha$ photon destruction, while important, cannot prevent the build-up of strong Ly$\alpha$ radiation pressure in H II regions, leading to radiative forces $\sim (2 - 16) \times L_{\rm bol}/c$, and Ly$\alpha$ force multipliers $M_{\rm F} \sim 10-60$. Ly$\alpha$ feedback may not preclude efficient star formation, but raises the threshold gas surface density for this to occur. We conclude that nearly all galaxy and star formation simulations are currently missing the strongest source of radiation pressure feedback in dense and metal-poor environments.
Read the full abstract

All authors Olof Nebrin, Aaron Smith, Garrelt Mellema, Kevin Lorinc, Daniele Manzoni

Earth and Planetary Astrophysics

A Massive Hot-Jupiter Companion that Disfavors Giant Planet Formation Beyond the Water-Ice Line

Eritas Yang, Tiger Lu, Daniel A. Yahalomi et al. (+1 more)

We report evidence for a brown-dwarf companion with mass $34^{+30}_{-11}~M_{\rm J}$ in the KELT-20 system, in which an ultra-hot Jupiter transits an A2-type star. The companion's properties are inferred from a joint analysis of astrometric accelerations and transit timing variations, and its present-day orbit imposes dynamical limits on where the hot Jupiter could have formed. Given the star's current luminosity, the water-ice line is expected at $\sim$8-15 au, but the companion's inferred pericenter distance of a few au would lead to orbit crossing or long-term instability for any planet formed at such distances. If the companion formed early and remained near its current orbit over the system's lifetime, the proto-hot Jupiter must have formed within $\sim$3.7 au to avoid orbit crossing, and within $\sim$1.5 au to remain dynamically stable over the system's lifetime. These results disfavor formation beyond the ice line and point instead to formation at smaller orbital distances followed by inward migration.
Read the full abstract

All authors Eritas Yang, Tiger Lu, Daniel A. Yahalomi, Joshua N. Winn

Astrophysics of Galaxies

Euclid: Early Release Observations – Internal kinematics and the convective-transition gap of NGC 6397

M. Griggio, M. Libralato, R. Gerasimov et al. (+130 more)

We present a 'multiple-pass' data-reduction tool designed for Euclid, based on software developed for the Hubble Space Telescope (HST), which improves the astrometric and photometric precision for faint sources and in crowded fields. In this work, we apply it to Euclid Early Release Observations of the Galactic globular cluster NGC 6397. By combining our new catalogue with archival HST data, separated by a time span of approximately 20 years, we were able to measure high-precision proper motions and investigate the radial variations in the energy equipartition and velocity anisotropy of the cluster. The combination of deep and wide-field observations also allowed us to derive the present-day local mass function of NGC 6397 and to study the radial dependence of mass segregation and binary fraction. Finally, we report the discovery of a subtle under-density of stars in the colour-magnitude diagram of NGC 6397 around a stellar mass of 0.35 M$_\odot$ with a more than 5$\sigma$ confidence level. This feature is consistent with the Gaia M-dwarf gap discovered in Galactic field stars, but it has never previously been observed in a globular cluster. The gap is caused by the onset of full convection in stellar interiors. We demonstrate that the properties of the gap provide tight constraints on the distance to NGC 6397 and its intrinsic metallicity dispersion, offering a new benchmark for stellar evolution models.
Read the full abstract

All authors M. Griggio, M. Libralato, R. Gerasimov, L. R. Bedin, A. Bellini, J. Anderson, E. Dalessandro, E. Vesperini, H. Baumgardt, D. Massari, M. Cadelano, R. E. Ryan Jr., I. McDonald, F. Annibali, E. Balbinot, J.-C. Cuillandre, D. Erkal, A. M. N. Ferguson, M. Kluge, P. B. Kuzma, T. Saifollahi, M. Schirmer, K. Voggel, B. Altieri, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, A. Balestra, S. Bardelli, P. Battaglia, A. Biviano, E. Branchini, M. Brescia, S. Camera, V. Capobianco, C. Carbone, J. Carretero, M. Castellano, G. Castignani, S. Cavuoti, A. Cimatti, C. Colodro-Conde, G. Congedo, C. J. Conselice, L. Conversi, Y. Copin, F. Courbin, H. M. Courtois, M. Cropper, G. De Lucia, F. Dubath, X. Dupac, M. Farina, R. Farinelli, S. Ferriol, E. Franceschi, S. Galeotta, K. George, C. Giocoli, J. Gracia-Carpio, A. Grazian, F. Grupp, S. V. H. Haugan, H. Hoekstra, W. Holmes, F. Hormuth, A. Hornstrup, K. Jahnke, M. Jhabvala, S. Kermiche, A. Kiessling, B. Kubik, M. Kümmel, H. Kurki-Suonio, A. M. C. Le Brun, S. Ligori, P. B. Lilje, V. Lindholm, I. Lloro, G. Mainetti, O. Marggraf, M. Martinelli, N. Martinet, F. Marulli, R. J. Massey, E. Medinaceli, S. Mei, M. Meneghetti, E. Merlin, G. Meylan, A. Mora, M. Moresco, L. Moscardini, R. Nakajima, S.-M. Niemi, C. Padilla, S. Paltani, F. Pasian, K. Pedersen, W. J. Percival, V. Pettorino, G. Polenta, M. Poncet, L. A. Popa, F. Raison, A. Renzi, J. Rhodes, G. Riccio, E. Romelli, M. Roncarelli, R. Saglia, Z. Sakr, D. Sapone, B. Sartoris, P. Schneider, A. Secroun, P. Simon, C. Sirignano, G. Sirri, L. Stanco, P. Tallada-Crespí, S. Toft, R. Toledo-Moreo, F. Torradeflot, I. Tutusaus, J. Valiviita, T. Vassallo, Y. Wang, J. Weller, F. M. Zerbi, E. Zucca, V. Scottez

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

Quintessential $\alpha$-attractors fit DESI

Alessandro Borys, Joaquin Masias and Marco Scalisi

We study quintessence in $\alpha$-attractor models in light of recent DESI indications for dynamical dark energy. We show that the knee of the attractor potential provides an excellent approximation to the axion-like quintessence model used as a DESI benchmark. This leads to a simple relation between the axion decay constant $f_a$ and the attractor parameter $\alpha$, allowing the experimental constraints to be translated into a preference for $\alpha=\mathcal{O}(1)$, in agreement with string-motivated expectations. We solve the background dynamics numerically and find good agreement with the DESI-preferred evolution of $w(z)$ up to $z\sim\mathcal{O}(1)$. More generally, we point out that the agreement between axion-like and attractor potentials reflects a common requirement imposed by the data: today's potential energy and slope are both of order the Hubble scale in Planck units. We finally comment on the origin of the required initial conditions, which can naturally arise in multifield attractor scenarios.
Read the full abstract

High Energy Astrophysical Phenomena

The Life and Death of Stars That Capture Primordial Black Holes

Ore Gottlieb, Matteo Cantiello, Cameron Norton et al. (+2 more)

Primordial black holes (PBHs) in the asteroid mass window ($10^{17}-10^{23}\,{\rm g}$) remain viable dark matter candidates and can be captured by stars. We develop the first global framework for the evolution of stars that capture PBHs, combining analytic calculations, stellar evolution models, 3D general-relativistic magnetohydrodynamic simulations, and Monte Carlo population synthesis. We find that the fate of these systems bifurcates: PBHs that form an accretion disk before consuming the host drive explosive disruption, whereas PBHs captured too late or growing too slowly consume the star quietly. Capture is dominated by three-body interactions with planetary or stellar companions. For a solar-type host with a Jupiter analog, inspiral within a main-sequence lifetime requires $M_{\rm BH}^{\rm crit}\gtrsim 10^{22}\,{\rm g}$, while lighter PBHs generally require tighter companions. Once deposited at the center, the PBH grows through inefficient quasi-spherical Bondi accretion; if it reaches the angular-momentum threshold before consuming the host, the inflow circularizes into a disk. Our Monte Carlo calculations yield sizable quiet-consumption and explosive-disruption populations, with final PBH masses $M_{\rm BH}\sim0.01-1\,M_\odot$ and disk-forming PBH spins $a_\ast\approx0.8$. Disk formation is the point of no return: disk winds and relativistic jets of $\sim10^{45}-10^{50}\,{\rm erg\,s^{-1}}$ disrupt the star within minutes. The resulting transients may include a $\sim$day-long UV/blue signal, radio afterglow, and, if the jet escapes, an X-ray-flash/low-luminosity gamma-ray-burst (XRF/llGRB) signal. For an $O(1)$ PBH dark matter fraction and optimistic capture assumptions, the event rate can reach that of llGRBs. The low-mass, high-spin remnants offer a complementary PBH probe and possible source for subsolar BH mergers.
Read the full abstract

All authors Ore Gottlieb, Matteo Cantiello, Cameron Norton, Ken Van Tilburg, Matthew Kleban

Astrophysics of Galaxies

A New Record Census of Dwarf AGN and a Bimodal $M_{\rm BH}$-$M_{\star}$ Scaling Relation with DESI DR1

Ragadeepika Pucha, S. Juneau, M. Mezcua et al. (+48 more)

Using the first spectroscopic data release from the Dark Energy Spectroscopic Instrument (DESI DR1), we conduct the largest systematic search for active galactic nuclei (AGN) signatures in $\approx$1.7 million low-redshift ($0.001 \le z \le 0.45$) line-emitting galaxies. Based on the [NII]/H$\alpha$ versus [OIII]/H$\beta$ emission-line ratio diagnostic, we identify AGN in 314,245/1,211,573 (25.9%) high-mass ($\log (M_{\star}/M_{\odot}) > 9.5$) and 9648/467,214 (2.1%) dwarf ($\log (M_{\star}/M_{\odot}) \le 9.5$) galaxies. The AGN fraction in line-emitting galaxies increases monotonically with stellar mass, from $\sim$1.4% at the low-mass end to $\sim$93.3% at the high-mass end. Among these AGN, 17,949 exhibit broad H$\alpha$ emission (BL-AGN), which we use to estimate their black hole (BH) masses via single-epoch virial methods. Leveraging this large BL-AGN sample, we extend the $M_{\rm BH} - M_{\star}$ scaling relation down to $\log (M_{\star}/M_{\odot}) \approx 7.8$ and $\log (M_{\rm BH}/M_{\odot}) \approx 4.4$. We find a bimodal distribution of galaxies in this space, characterized by two distinct empirical relations. The relation for galaxies with high stellar and BH masses, extrapolates to overlap with the high-redshift overmassive BH candidates. Our results suggest at least two broad evolutionary pathways for these overmassive BHs and their host galaxies across cosmic time. With this paper, we release the EmFit value-added catalog, containing emission-line flux and width measurements for $\sim$7.4 million galaxies, the largest catalog with emission-line decomposition into narrow, broad, and outflow components to date. This work significantly expands upon the early DESI results and provides a statistical sample for probing the galaxy$-$BH connection in the low-mass galaxy regime.
Read the full abstract

All authors Ragadeepika Pucha, S. Juneau, M. Mezcua, Arjun Dey, Y.-Y. Mao, D. M. Alexander, C. Circosta, V. A. Fawcett, Wei-Jian Guo, J. Moustakas, S. Panda, M. Siudek, Z. Yu, J. Aguilar, S. Ahlen, D. Bianchi, D. Brooks, T. Claybaugh, K. S. Dawson, A. de la Macorra, P. Doel, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, E. Gaztañaga, Satya Gontcho A Gontcho, G. Gutierrez, C. Hahn, K. Honscheid, R. Joyce, R. Kehoe, T. Kisner, A. Kremin, M. Landriau, L. Le Guillou, M. Manera, A. Meisner, R. Miquel, S. Nadathur, W. J. Percival, F. Prada, I. Pérez-Ràfols, G. Rossi, E. Sanchez, D. Schlegel, M. Schubnell, J. Silber, D. Sprayberry, G. Tarlé, B. A. Weaver, H. Zou

Astrophysics of Galaxies

Toward Unbreaking the Universe: MINERVA Measurements of Color Gradients in Massive Quiescent Galaxies Can Help Ease Too-Early Star Formation Tensions

Sam E. Cutler, Luke Robbins, Danilo Marchesini et al. (+41 more)

The discovery of a population of massive, ancient quiescent galaxies within the first 2 Gyr of the Universe's history has led to significant tensions with models of galaxy formation. However, these analyses are often based on slit spectroscopy, which typically captures only the center-most region of these galaxies and, crucially, assumes these cores are representative of the entire galaxy. To illustrate the varying stellar populations present throughout these galaxies, we present an analysis of color gradients in four $z>3$, $\log(M_\star/M_\odot)>11$ quiescent galaxies which previous works have argued are in tension with models. Using medium-band photometry from MINERVA JWST observations, we measure resolved photometry in a series of elliptical annuli out to $0.7^{\prime\prime}$ ($\sim4~R_e$). We find negative color gradients in three galaxies, and for the most extreme color gradient ($\Delta(U-V)/\Delta R=-0.126\pm0.030~{\rm mag~kpc^{-1}}$), we find the stellar mass is 0.1 dex lower when compared to photometry measured within NIRSpec slits. In the limiting case where these color gradients are entirely driven by age, we find lessened tensions with extreme value statistics models out to $z\sim9.5$, though different stellar population modeling choices also contribute significantly. Ultimately, these findings highlight the need for integral field unit spectroscopy. Spatially-resolved spectra can provide the evidence needed to break the age–dust–metallicity degeneracy, and reliably separate the effects of the observed color gradients from the effects of different physical modeling assumptions on the formation histories of these galaxies.
Read the full abstract

All authors Sam E. Cutler, Luke Robbins, Danilo Marchesini, Katherine A. Suess, Adam Muzzin, Gabriel Brammer, Yoshihisa Asada, Nicholas S. Martis, Stacey Alberts, Jacqueline Antwi-Danso, Aidan P. Cloonan, Ivo Labbé, Tim B. Miller, Ikki Mitsuhashi, Alexandra Pope, Anna Sajina, Ghassan T. E. Sarrouh, Monu Sharma, Mauro Stefanon, Edgar P. Vidal, Chris J. Willot, Rachel Bezanson, Maruša Bradač, Olivia R. Cooper, Robert Feldmann, Ben Forrest, Karl Glazebrook, Jenny E. Greene, Valentina La Torre, Jamie Lin, Michael V. Maseda, Ian McConachie, Themiya Nananyakkara, Gaël Noirot, Richard Pan, Kesha A. Patel, Veronica Pratt, Marcin Sawicki, David J. Setton, John R. Weaver, Arjen van der Wel, Katherine E. Whitaker, Yunchong Zhang, Kumail Zaidi

Earth and Planetary Astrophysics

HERMES: HiERarchical Modelling for Exoplanet Science

Wasi M. F. Naqvi and Nicolas B. Cowan

ESA's Ariel Space Mission will characterise the atmospheres of approximately 1000 exoplanets to quantify population-level trends. We present HERMES (HiERarchical Modelling for Exoplanet Science), a multidimensional Bayesian framework that probes population-level correlations across multiple axes of diversity. The specific use case we present is the multidimensional relation between stellar metallicity, planetary mass, and atmospheric metallicity. Starting from the Ariel Mission Candidate Sample (Edwards & Tinetti 2022), we select confirmed planets with available masses and stellar metallicity, inject plausible multidimensional trends and demonstrate successful parameter recovery. Simulated surveys are generated with a variety of leverage and sample size, in the presence of intrinsic astrophysical scatter and measurement noise. By fitting independent Bayesian models to each survey, we confirm that survey leverage remains a reliable predictor of trend precision even in multiple dimensions and in the presence of intrinsic astrophysical scatter. For an Ariel Tier 2 transit survey of at least 400 planets, HERMES robustly recovers the correlation between stellar and planetary metallicity despite intrinsic scatter in planetary abundances as large as 1.2 dex. These results establish HERMES as a practical tool for survey design and science yield forecasting in preparation for Ariel and other surveys probing population-level trends.
Read the full abstract

High Energy Astrophysical Phenomena

The Role of Stellar Spin in Repeating Partial Tidal Disruption Events

Ananya Bandopadhyay, Benjamin Amend, Eric R. Coughlin et al. (+3 more)

The repeated tidal stripping of a star by a supermassive black hole, known as a repeating partial tidal disruption event (rpTDE), can give rise to a transient that rebrightens months to years after the first outburst. Among the rpTDE candidates so far observed, some exhibit dimmer peak luminosities during each successive outburst, which is a trend that has not been reproduced from theoretical models when the star survives more than one encounter with the black hole. Here we suggest that this trend can be recovered if the partially disrupted star is initially (i.e., prior to its first mass-stripping event) rapidly rotating, which is expected if the star was placed on its orbit through the Hills breakup of a tidally locked and tight binary. We test this hypothesis with hydrodynamical simulations of high-mass ($\geqslant 1 M_{\odot}$) main sequence stars repeatedly partially disrupted by a $10^6 M_{\odot}$ black hole, and demonstrate that successively dimmer outbursts are indeed recovered for high (tens of percent breakup) and prograde (i.e., aligned with the orbital angular momentum) stellar spins. Our results provide strong indirect evidence for the operation of the Hills mechanism in seeding the stars in rpTDEs.
Read the full abstract

All authors Ananya Bandopadhyay, Benjamin Amend, Eric R. Coughlin, C. J. Nixon, Dheeraj R. Pasham, T. Wevers

High Energy Astrophysical Phenomena

Systematic Error in Approximate Models of the GRB Early Afterglow

Benjamin Amend, Eric R. Coughlin and Jonathan Zrake

Gamma-ray burst (GRB) afterglows are thought to arise when relativistic ejecta launched by a compact central engine drive a blast wave into the surrounding circumburst medium, producing broadband synchrotron emission. We present a rigorous assessment, based on high-resolution special relativistic hydrodynamics simulations, of a widely adopted `two-zone model' for approximating the dynamics of the early afterglow phase. Before the onset of the Blandford-McKee (BMK) self-similar solution, the outflow generally produces two emission components, associated with the forward-shocked circumburst medium and the reverse-shocked ejecta. The subsequent evolution depends on whether the reverse shock significantly decelerates the ejecta as it crosses the shell, separating the so-called relativistic and Newtonian reverse shock regimes. We show that when the reverse shock is Newtonian, it crosses the ejecta shell long before BMK self-similarity is established, leaving a prolonged interval that can span $\sim$ hours in observer time in which the true hydrodynamic evolution is not captured by standard semi-analytic prescriptions. We demonstrate that this mismatch can, for representative afterglow parameters, substantially overpredict the reverse-shock emission from radio through ultraviolet frequencies, or overpredict the forward-shock emission at X-ray frequencies, depending on how the transition away from the two-zone model is prescribed.
Read the full abstract