Papers for Monday, May 18 2026

Modern particle-in-cell (PIC) codes have become an integral tool in plasma astrophysics. As most plasma phenomena grow from initially small instabilities, it is important to ensure PIC codes can suppress noise and ensure that any growing instability is indeed physical. Therefore, we introduce our efforts to implement higher-order methods for the current deposit and field interpolation as well as generalized field stencils for the field solver in the PIC code \texttt{Entity}. Our updated current deposit scheme allows for up to $11^\mathrm{th}$-order accurate interpolation, while the generalized stencils for the field solver can be tuned to suppress numerical dispersion. We perform extensive tests to ensure high accuracy of the implemented schemes for charge conservation, stabilization against numerical heating, improved energy conservation, and suppression of numerical Cherenkov effects. To supply a benchmark on performance impact, we demonstrate the scaling of the higher-order current deposit and discuss the possible performance balance between higher-order interpolation, numerical resolution, and the inclusion of additional current filtering.

Little Red Dots (LRDs), candidate high-redshift supermassive black holes accreting in dense gas, remain undetected in X-rays. In previous work, we provided the first quantitative models that reproduce the optical and near-infrared spectra of LRDs with the Sirocco radiative transfer code, thereby constraining the properties of the surrounding gas. Here, we use those constraints to predict the X-ray attenuation produced by dense gas cocoons, and explore its dependence on Balmer-break strength, metallicity, intrinsic X-ray SED, and observed bandpass as a function of redshift. We find that the X-ray constraints are very tight, requiring both extinction by a Compton-thick gas column $N_{\rm H}\sim10^{25}\,{\rm cm}^{-2}$ with moderate metallicity ($0.05$-$0.1\,Z_\odot$) and intrinsically weak X-ray emission (bolometric to X-ray luminosity ratio, $k_{\rm bol,X}\gtrsim 30$) as observed in high accretion rate, narrow-line AGN, to make LRDs sufficiently faint to evade detection. Intrinsically bright X-ray emitters as seen in typical broad-line AGN would be detected even behind the typical modest metallicity, Compton-thick gas columns inferred from the optical spectra. Very low metallicity objects could be detected in X-rays even with low intrinsic X-ray luminosities, suggesting that LRDs are not (currently) chemically pristine.

The Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields offer an exceptional combination of depth, spatial resolution, and area for identifying a shear-selected sample of dark matter overdensities. We present the first near-infrared (NIR) weak-lensing (WL) analysis of the 0.23 square degrees covered by the HST CANDELS fields: COSMOS, UDS, EGS, GOODS-N, and GOODS-S. Leveraging the high sensitivity of HST NIR imaging to distant galaxies, we achieve a WL source galaxy density of $\sim170$ galaxies arcmin$^{-2}$. Our analysis identifies 12 shear-selected overdensities spanning masses from $M_{200}=(0.2$--$2.2)\times10^{14}\ M_\odot$, with a median mass of $M_{200}=5.5\times10^{13}\ M_\odot$, demonstrating the strong capability of NIR WL for measuring low-mass systems. The systems lie in the redshift range $0.22<z<0.9$, with a mean redshift of $z=0.68$. We utilize multiwavelength data to confirm the nature of the overdensities. Seven of the overdensities have diffuse X-ray emission reported in the literature, with X-ray centroids that are spatially consistent with our WL peaks, confirming their nature as collapsed structures. We find that our WL detections broadly follow the expected X-ray luminosity--WL mass scaling relations. By stacking the tangential shear of all detections, we determine the average radial mass density profile and find that it is well fit by an NFW model with fitted concentration and mass of $4.9\pm2.1$ and $M_{200}=1.3\pm0.3\times10^{14}\ M_\odot$, respectively. These results serve as a precursor to NIR WL science with the Roman High Latitude Wide Area Survey.

Eccentric binary black holes (BBHs) formed through dynamical interactions can significantly contribute to gravitational wave (GW) detections. In this work, we present a simulated catalog of dynamically-formed, stellar-mass BBHs in the local universe, incorporating contributions from the Galactic field (flyby interactions), Galactic nucleus (eccentric Kozai-Lidov evolution), and globular clusters (N-body interactions). Our results predict a wide, highly eccentric BBH population in the Milky Way (MW), with source counts of $\sim 36, 13, 4.7, 2.3, 1.0$ (for $\mathrm{SNR} > 1, 3, 8, 20, 50$, respectively) during a 10-yr LISA observation. Extending this model to cosmological populations, we show that different dynamical channels can produce distinct eccentricity distributions in the LVK band and can contribute hundreds of additional low-SNR mHz sources. Specifically, our model yields a merger rate of $\Gamma \sim 9 \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$ and $\sim 490$ extragalactic mHz BBHs with $\mathrm{SNR} > 1$. However, due to the lower mass and weaker GW signals of stellar-mass BBHs, this number declines sharply at higher detection thresholds (e.g., $\sim 1$ for $\mathrm{SNR} > 8$). We further highlight the impact of eccentric BBH signals on the LISA global fit, showing that their individual harmonics can be independently detected in the Milky Way, and may mimic circular binaries with systematically biased chirp masses. Lastly, we show that post-Newtonian waveforms converge reliably for eccentric BBHs with masses of $\lesssim 10^3 M_\odot$ in the mHz band. Overall, eccentric BBHs represent a prevalent and promising target for future space-based GW observatories. The simulated catalog and the LISA Eccentricity Astrophysics Package (LEAP) developed in this work are publicly available at this https URL.

The classical picture for the formation of stellar bars -- key dynamical drivers of the evolution of galaxies -- is through secular evolution of instability in gas poor, stellar-dominated disks. The detection with the James Webb Space Telescope (JWST) of stellar bars and spiral arms in galaxies at early cosmic times has thus challenged LambdaCDM-based expectations, which recent studies reconcile by suggesting that these galaxies are baryon-dominated and have already consumed most of their gas. Yet, a paradox arises, as early galaxies are expected to be increasingly rich in gas, which is generally considered to prevent or slow down stellar bar formation. Here, we show the detection of a stellar bar in GN20, a gas-rich star-forming disk galaxy at a redshift of z=4.055, only 1.5 billion years after the Big Bang. Simultaneous observations of the stars, gas, and dust reveal that GN20 is indeed baryon-dominated (over dark matter; 70+/-30%), but the baryonic mass is largely in the form of gas (75+/-25%). This discovery demonstrates that gas-rich disks do support rapid stellar bar formation in the early Universe, motivating a new theoretical perspective on bar formation in gas-rich systems, and providing a potential new mechanism for very early galaxy assembly and quenching.

Observations reveal extended \OVI\, reservoirs in the circumgalactic medium (CGM) of dwarf galaxies, yet current simulations systematically underpredict \OVI\, column densities. Utilizing two suites run with different simulation codes, the \MM\, simulations (Marvelous Massive Dwarfs and Marvel-ous Dwarfs) and the publicly available FIRE-2 simulations, we explore the role of subgrid models and the resulting CGM phase in shaping \OVI\, production. By comparing observationally derived \OVI\, masses to the mass of oxygen produced over the galaxies' star formation history, we find evidence for an underproduction of oxygen for low-mass simulated galaxies. Despite clear differences in feedback implementation, CGM structure, and metal mixing, we find that \OVI\, in both suites generally self-selects cool/warm ($\rm log\, T\,/K \sim 4.5$), diffuse ($\rm log\,n_{gas}\,/cm^{-3} \sim -5.0 $), and moderately metal-enriched ($\rm log\, Z/Z_{\odot} \sim -1 $) material at large radii from the galaxy. We show that neither the choice of ultraviolet background nor plausible variations in CGM thermal structure can close the gap with observations. Taken together, our results point to a possible underproduction of oxygen in low-mass galaxies. Feedback prescriptions contribute via insufficient metal transport to large radii. Hence, the \OVI\, deficit may motivate an investigation of current modeling choices for supernova yields, star formation, and feedback in low-metallicity environments.

Predictions for observables involving the cosmological 21-cm background require calculations of spatial correlations of star formation rate densities (SFRDs) which have a nonlinear dependence on the baryon-dark matter relative velocity. Prior work derived these SFRD correlations with a simplifying assumption that neglected the difference between the correlations of the components of the velocity parallel and perpendicular to the separation between the two points being correlated. Here we calculate the full joint PDF of the squares of the peculiar velocity at two different points. The error that arises in predictions for 21-cm fluctuations if this subtlety is overlooked is generally less than a few percent, but it can be larger for some values of wavenumber $k$ and redshift $z$ if there are cancellations between different contributions to the total signal. The correct expression is easily implemented and increases the run time of the code by only a few percent.

When a rapidly-rotating, highly magnetized white dwarf (WD) approaches the Chandrashekhar limit through mass accretion, it can undergo an accretion-induced collapse (AIC) to form a proto-neutron star or protomagnetar. The protomagnetar can drive a magnetically-dominated relativistic outflow, whose low entropy can lead to efficient formation of heavy nuclei. In this work, we propose that such relativistic outflows from AIC of WDs can contribute as sources of ultra high-energy cosmic rays (UHECRs). We model the acceleration of heavy nuclei in these relativistic outflows, and show that AICs can dominantly power the observed UHECRs, if a majority of them host relativistic outflows. Accounting for uncertainties in the acceleration mechanisms and AIC rates, AICs can contribute $\sim$ a few $10^{43} - 10^{45}\ {\rm erg \ Mpc}^{-3} {\rm yr}^{-1}$ in UHECR energy generation rate density, assuming iron-like nuclei.

Understanding how galaxies and active galactic nuclei (AGN) jointly drive the reionization of the intergalactic medium (IGM) across cosmic time remains a major challenge in cosmology. We present Lumina, a large-volume radiation-hydrodynamic simulation that self-consistently follows the coupled evolution of the intergalactic medium, galaxies, and AGN through HI, HeI, and HeII reionization down to redshift $z=3$. Lumina evolves a cosmological volume of comoving side length $L_{\mathrm{box}}=500\,\mathrm{cMpc}$ with $2\times 6000^{3}$ resolution elements, corresponding to baryonic and dark-matter mass resolutions of $3.6\times 10^{6}\,\text{M}_{\odot}$ and $1.9\times 10^{7}\,\text{M}_{\odot}$, respectively. The simulation uses the moving-mesh code AREPO, combining the IllustrisTNG galaxy-formation model with a GPU-accelerated M1 radiation-transport solver in six frequency bins. The initial conditions employ separate transfer functions for baryons and dark matter and include their relative streaming velocity. Lumina predicts a late, predominantly stellar-driven hydrogen reionization, with the median sub-volume fully ionized by $z\approx 5.2$ and residual neutral HI patches persisting until $z\approx 4.75$. HeII reionization is driven self-consistently by AGN and is nearly complete by $z=3$. The simulation yields a Thomson-scattering optical depth in excellent agreement with Planck, an IGM thermal history and photoionization background broadly consistent with observational constraints, and a clear late-time thermal boost associated with HeII reionization. Its galaxy population remains consistent with the original IllustrisTNG project, while the larger volume improves statistics for rare objects, large-scale environments, and cosmic variance, enabling forward modelling of observables linking HI and HeII topologies to the evolving galaxy and AGN populations.

A tight correlation between nuclear millimeter (mm) and X-ray emission has recently been found in nearby ($z < 0.01$) and low-Eddington ratio ($\rm \lambda_{Edd} < 0.1$) radio-quiet (RQ) Active Galactic Nuclei (AGN), suggesting a common origin in the hot X-ray corona. We test this relation in nine more distant RQ AGN ($z \sim 0.06-0.16$) with higher bolometric luminosities ($\log(L_{\rm bol}/\mathrm{erg\,s^{-1}})=45.3-46.3$), Eddington ratios ($\rm \lambda_{Edd} = 0.19-0.85$), and X-ray bolometric corrections ($\kappa_{2-10}=29-194$), selected from the Burst Alert Telescope (BAT) survey. We obtained quasi-simultaneous observations with Swift at 2-10 keV and the Atacama Large Millimeter/submillimeter Array (ALMA) at 100 GHz and with high angular resolution ($<0.14$"). We find that these high-luminosity AGN lie above the mm/X-ray correlation defined by lower-luminosity sources. A joint fit to both samples yields a second-degree polynomial with an intrinsic scatter of 0.32 dex. Furthermore, the mm emission correlates linearly with both the UV disk luminosity and $L_{\rm bol}$, with intrinsic scatters of 0.45 and 0.35 dex, respectively. We propose that the deviation from the linear mm/X-ray relation arises from a two-component coronal electron population: thermal electrons that produce X-rays, but become less efficient at higher luminosities, and non-thermal electrons that produce mm emission and remain tied to $L_{\rm bol}$. Additional mm emission from outflow-driven shocks may also contribute, though SED modeling and spectral index studies favor a coronal origin.

The present understanding of C-rich AGB chemistry largely relies on observations of the archetypal carbon star IRC+10216. Current molecular abundance estimates for carbon stars are based either on single-dish spectra sampling a range of excitation conditions, or on interferometric mapping of a few lines. We aim to estimate the circumstellar abundances of SiO, SiS, and their most abundant isotopologues ($^{29}$SiO, $^{30}$SiO, $^{29}$SiS, $^{30}$SiS, and Si$^{34}$S) for a sample of five carbon stars. We derived molecular abundances using detailed 1D non-local thermodynamic equilibrium (non-LTE) radiative transfer (RT) modelling, constrained by both morphological and excitation information obtained from spatially resolved ALMA maps and single-dish observations. We further compared the derived abundances to chemical modelling results. We obtain good fits to the SiO and SiS line profiles, and derived well-constrained abundance profiles and reliable isotopic ratios for all sources except AFGL 3068. While the SiS peak abundances are very similar across the sample (2.0$\times$10$^{-6}-4.7\times$10$^{-6}$), we find that the SiO peak abundances of the rest of the stars are a factor of $\sim$5 larger than that of IRC+10216. The $e$-folding radii ($R_\mathrm{e}$) are in the range 1.3$\times$10$^{16}-7.0\times$10$^{16}$ cm for SiO and 6.0$\times$10$^{15}-1.0\times$10$^{17}$ cm for SiS. The $R_\mathrm{e}$ increases with gas density for both SiO and SiS. Chemical models reproduce the derived SiO abundance profiles well, while over-predicting the SiS $R_\mathrm{e}$ values. Our models highlight the necessity of having spatially resolved observations across a broad range of excitation conditions, while also making evident the limitations inherent in 1D RT modelling using simplified (circum)stellar models. We find that the currently assumed SiS photodissociation rate in chemical models is underestimated.

Measuring galaxy sizes is fundamental to understanding how galaxies grow and evolve. Traditional methods to measure sizes either trace the concentration of light (i.e., effective radius) or are limited by the depth of the survey (isophotal methods). With the advent of deep, wide surveys, a new physically motivated definition of size has emerged: the edge of the galaxy, defined as the most distant location where star formation has occurred or is still occurring. In this work, we take advantage of the extraordinary depth and spatial resolution of the Hubble and James Webb Space Telescopes to perform an accurate study of galaxy edges at $z=1$. Using 22 photometric bands, we derive radial age and metallicity profiles for two disk galaxies in the GOODS-South field with stellar masses of around $4\times10^{10}\ M_\odot$. The age profiles display a characteristic U-shape, while the metallicity profiles steadily decrease with galactocentric distance. The turnover in the age profile occurs near the galaxy edge, suggesting that stellar migration is responsible for the stars beyond the edge of these galaxies. Comparison with $z=0$ disk galaxies suggests that galaxies at $z=1$ grow inside-out, with little or no increase in mass within the inner 8 kpc, but a significant increase (approximately 300\%) in the outer regions.

Ocean worlds are considered as targets for life detection missions because they meet several key requirements for habitability. However, identifying potential life on other worlds requires observing clear and unambiguous biosignature signals above the existing abiotic baseline. Consequently, this necessitates evaluating uncertainty and variability in the abiotic baseline, including processes that can overlap, attenuate, or obfuscate biosignatures before they are observed. This article develops a quantitative framework for holistically evaluating abiotic baselines on ocean worlds to guide life detection strategies. Using Enceladus as an example, we assess the potential of using: i) CH$_{4}$ isotopes and their relationship with CO$_{2}$, and ii) amino acid chirality as biosignatures, demonstrating that uncertainties in abiotic processes currently prevent hypothetical future ${\delta}^{13}$C$_{\mathrm{CO2}}$ and ${\delta}^{13}$C$_{\mathrm{CH4}}$ measurements from definitively inferring a biosphere on Enceladus. Additionally, our results quantitatively show that neglecting the abiotic baseline risks false negative life detection claims for both isotopic and chiral biosignatures. Interpreting these and other alternative biosignatures on Enceladus, Europa, Titan, and similar planetary bodies therefore requires complimentary geophysical observations such as constraining internal temperatures to within $\sim$10-100$^{\circ}$C, and improving characterisation of the target's rheology, lithology, initial abiotic organic inventory and ocean transport timescales.

Pulsar wind nebulae (PWNe) constitute the largest population of Galactic very-high-energy (VHE; $E > 100$ GeV) $\gamma$-ray sources and are key laboratories for studying particle acceleration and pulsar--supernova remnant (SNR) interactions. However, realistic population-level predictions have so far lacked any detailed treatment of the reverberation phase, when the nebula is compressed by the SNR reverse shock, significantly altering its dynamics and radiative spectrum. We employ the hybrid \texttt{TIDE+L} framework, which combines a thin-shell dynamical model with a Lagrangian treatment of the SNR structure during reverberation, allowing self-consistent evolution of thousands of PWNe across all stages up to $10^5$ yr. Each source is evolved under distributions of pulsar spin-down, SNR, and environmental properties, and the resulting $\gamma$-ray fluxes are used to estimate the detectability by current and next-generation $\gamma$-ray observatories while accounting for their sensitivity and sky coverage. The model predicts that the upcoming Cherenkov Telescope Array Observatory (CTAO) will detect an order of magnitude more PWNe than those firmly detected in the TeV range, confirming its dominant contribution to the forthcoming TeV population census. Our results demonstrate that realistic modeling of reverberation is important for predicting the Galactic TeV PWNe population.

The origin of high-energy astrophysical neutrinos remains a key open question in multimessenger astrophysics. A correlation between unabsorbed hard X-ray emission and high-energy neutrino luminosity has been reported in a sample of six active galactic nuclei with the highest individual IceCube significances, linking neutrino production to compact, photon-rich environments near supermassive black holes. In this work we study whether the threshold-near IceCube excesses associated with seven NuSTAR-observed blazars are statistically consistent with that established relation. Calibrating the relation between the neutrino and hard X-ray luminosities as $\log L_\nu = \alpha + \beta \log L_\mathrm{hX} + \mathcal{N}(0, \sigma_{\rm int}^2)$ on the six published sources via a Bayesian regression with errors on both axes, the recovered slope is consistent with $\beta = 1$, and the intrinsic scatter is $\sim 0.6$\,dex. All seven new blazars are posterior-predictively consistent with this calibration ($\chi^2_7 = 1.58$, $p = 0.98$) under the working hypothesis that the published IceCube best-fit neutrino numbers $\hat{n}_s$ values reflect the signal. A null-injection test confirms that, given the present calibration sample size, the consistency test does not by itself adjudicate between signal and selected-background origins. A distance-free $L_\mathrm{hX}/L_\nu$ ratio diagnostic places both populations within the photohadronic prediction band, statistically indistinguishable. A flux-space permutation test on the 13-source joint sample, with construction-controlled $d_L^{\,2}$ distance bias, rejects random pairing $L_\mathrm{hX}$--$L_\nu$ with $p = 6.3 \times 10^{-4}$ ($3.23\,\sigma$). We interpret these results as a conditional consistency check; a detection-level statement requires either an enlarged calibration set or an X-ray-weighted IceCube stacking likelihood with internal data.

At high redshifts, star formation in galaxies is more often concentrated in clumps than in spiral arms. Although clumps are actively studied, it is rarely done considering spiral arms as objects for study as well, and the connection between clumps and spirals remains understudied. We used a sample of 159 spiral galaxies at $0.1 \leq z \leq 3.3$ observed by HST and JWST. Using the residual images from photometric models with spiral arms constructed before, which was not done previously, we have done identification of clumps and measured their properties with photometric decomposition, finding 3003 clumps in overall, and performing SED fitting for a fraction of them. We examined the overall properties of clumps, focusing on the properties of spiral structure. We found that clumps luminosities, masses and sizes are smaller than commonly accepted in literature, either for the reason spiral arms were modelled separately, or because clumps in spiral galaxies are different compared to clumpy ones. We demonstrate the connection between clumps and spirals, in particular clumps are spatially concentrated towards spirals and various parameters of clumps and spirals correlate. Also, clumps in spiral arms tend to be smaller but brighter compared to clumps in inter-arm area, but their colours are similar. There are also some differences between clumps and spirals such as colour, emphasizing the importance of their separate analysis. Our results probably indicate that spiral arms stimulate the formation of clumps, although star formation properties of clumps in spiral arms is not changed compared to inter-arm region.

The intriguing circumpolar cyclone pattern at Jupiter's poles raises fundamental questions about how these systems are organized vertically and, further, how the planet's internal heat shapes and sustains them in the absence of solar insolation. We report recent close-in observations of Jupiter's north pole acquired by NASA's Juno Microwave Radiometer (MWR), which achieved comprehensive microwave mapping of the region at an unprecedentedly high resolution. Using six-channel measurements from eleven perijove passes (PJ51-PJ61) poleward of 75N, we derive polar-mean nadir brightness temperatures and limb-darkening spectra that together point to two equally plausible atmospheric scenarios: (1) a dry-adiabatic profile with slightly depleted ammonia gas at a few bars, or (2) a moist-adiabatic profile with uniform ammonia. Markov chain Monte Carlo retrievals yield a deep ammonia abundance of 354.8+12.0/-11.0 ppmv (3+/-0.1 x solar) and a water abundance of 1.8+1.5/-1.1 x 1000 ppmv (2.1+1.8/-1.3 x solar), resembling previous estimates at lower latitudes. Remarkably, the north pole is found to be 6-7 K warmer than the equator at the 1-bar level, although the inferred difference is close to the 1-sigma uncertainty level. If confirmed, this result would suggest an enhanced internal heat flux toward the poles, which is consistent with the more intense lightning activity observed at high latitudes.

We generalize lower limits on the dark matter (DM) particle mass $m$ derived from Milky Way (MW) satellite galaxy abundances to scenarios in which DM is an ultralight scalar field produced with a field power spectrum peaked at a subhorizon wavenumber $k_*$. In these models, the DM field free-streams similar to warm dark matter while also exhibiting significant small-scale wave interference effects. The resulting dimensionless density power spectrum shows two effects: (i) free-streaming suppression at $k_{\rm fs}\sim k_{\rm eq}/[(k_*/a_{\rm eq}m)\ln(a_{\rm eq}m/k_*)]$; (ii) Poisson-like enhancement related to wave interference, at $k\gtrsim10^{-2}k_*$, which saturates near the Jeans scale $k_{\rm J}\sim k_{\rm eq}/(k_*/a_{\rm eq}m)$. Comparing these predictions with established constraints on a free-streaming cutoff in the linear matter power spectrum from the MW satellite population, we obtain $m>6\times10^{-18}\,{\rm eV}\,(k_*/10^4\,{\rm Mpc}^{-1})$ for $k_*>10^4\,{\rm Mpc}^{-1}$ at 95\% confidence. For smaller $k_*$, Poisson-noise enhancement on MW satellite scales weakens the constraint, yielding $m>6\times10^{-18}\,{\rm eV}\,(k_*/10^4\,{\rm Mpc}^{-1})^2$ for $k_*<10^4\,{\rm Mpc}^{-1}$ at 95\% confidence.

The Nova Synthetic Data Base (NSDB) is presented as the first publicly available database of synthetic spectra for classical nova shells, spanning an unprecedented range of physical parameters (e.g., ejecta mass, chemical composition, temperature, and luminosity of the white dwarf) at several post-eruption ages. Generated using detailed 3D photoionization models, this homogeneous database enables a systematic exploration of spectral features in novae. In this work, we introduce a principal component analysis/AI-based framework to derive time-dependent proxies for retrieving the physical properties of novae from limited spectral data. By analyzing the correlations between the eigenspectra and the grid's variables, a reduced set of diagnostic spectral lines is derived, paving the way for robust multiregressor machine-learning algorithms with a minimal effort observational set. The prediction capability of the method is high and robust to data noise. The results establish a proof of concept for the use of model grids combined with physically controlled AI as a tool to interpret novae observations in the context of the large number of events expected from future wide-area surveys.

Galaxy cluster abundance provides a powerful probe of the $\Lambda$CDM model and enables precise constraints on cosmological parameters. Millimeter-wavelength surveys detect clusters through the Sunyaev-Zeldovich (SZ) effect, and are particularly effective at high redshifts. However, the SZ signal can be significantly contaminated by emission from Active Galactic Nuclei (AGN), particularly AGN within the Central Galaxies (CGs). This contamination reduces the SZ signal strength at the frequencies most accessible from the ground, which reduces detection significances or converts cluster detections to non-detections, thereby diminishing survey completeness and introducing biases in cosmological analyses. In this work, we analyze three clusters that host bright AGN in their CGs using 30 and 90 GHz observations from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). In each case the AGN emission overwhelms the cluster SZ signal, resulting in non-detections in the Atacama Cosmology Telescope (ACT) survey. We present signal to noise ratio (SNR) estimates for the clusters after subtracting the AGN signal from 90 GHz ACT maps using the CARMA measurements, demonstrating high SNR cluster detections once this contaminating emission is removed. Using cluster pressure profiles derived from Chandra X-ray data, we subtract the expected SZ signal from the 150 GHz ACT maps to estimate the flux density of the AGN in that band. Leveraging the time-asynchronous CARMA observations, we also assess temporal variability in the AGN emission, and find low fractional variability for our sample. Finally, we discuss the importance of modeling and mitigating AGN contamination in SZ cluster surveys.

The explodability of massive stars, namely whether they undergo core-collapse supernovae (CCSNe) or form black holes (BHs), strongly influences galactic chemical evolution (GCE). Details of the explodability are still controversial, but realistic predictions including metallicity-dependence are becoming available through stellar-evolution and explosion calculations. In the present work, we implement recently-proposed metallicity-dependent explodability prescriptions into a GCE framework. We show that the physics-motivated explodability prescriptions reproduce the key observed abundance trends. Further, within uncertainties of the explodability models, the GCE model provides strong constraints on the region of the BH formation in the mass-metallicity space. Guided by these findings, we further construct a simplified form of the metallicity-dependent explodability designed to alleviate the red supergiant (RSG) problem and explore its compatibility with GCE constraints. We find that such a solution exists, if (1) the net outflows from the system are negligible/absent, and (2) the transition of the explodability takes place at sub-solar metallicity. These results demonstrate that GCE can provide meaningful constraints on massive-star explodability and that explodability prescriptions capable of addressing the RSG problem can be constructed without violating chemical-evolution observables. We also show that a metallicity-dependent initial mass function can improve agreement with observations; this effect becomes importance once coupled with the metallicity-dependent explodability.

We report the discovery of an extended GeV {\gamma}-ray source, 4FGL J0330.7+5845e, associated with the star-forming region AFGL 490 using 17 years of Fermi-LAT data. The emission is spatially coincident with a dense molecular cloud and centered near the massive protostar AFGL 490. Its spectral energy distribution shows a distinct high-energy cutoff. Both leptonic and hadronic models can fit the {\gamma}-ray spectrum, but energetic arguments rule out stellar winds as the primary accelerator. Instead, the protostellar jet driven by AFGL 490 is identified as a plausible site for particle acceleration, and the derived timescales and maximum particle energies are consistent with theoretical predictions for such jets.

The recent 230 GHz observations from the Event Horizon Telescope collaboration have successfully imaged the supermassive black hole shadow of the M87 galaxy. However, the relatively high radiative efficiency observed in the hot accretion flow suggests that radiative cooling is non-negligible and should be considered when calculating the electron temperature. In this study, we compare accretion models without and with radiative cooling across a range of mass accretion rates, $\dot{M}_{\mathrm{BH}} = (1.0 - 10) \times 10^{-6}\,\dot{M}_{\mathrm{Edd}}$, aiming to assess the impact of cooling on the disk structure, electron temperature distribution (eDF), black hole shadow morphology, broadband spectral energy distributions (SEDs), and flux variability. We performed general relativistic radiative transfer (GRRT) calculations on two-temperature, radiative, general relativistic magnetohydrodynamic (GRMHD) simulations, employing different electron heating prescriptions and nonthermal eDFs, analyzing the radiation transfer due to synchrotron emission at 230 GHz with inclination angle of $163^\circ$. These simulations are targeted toward M87$^{*}$. By comparing density profiles, eDFs, GRRT images, SEDs, and time variability between models, we find that the radiative cooling sharply decreases the electron temperature in the dense inner disk around the equatorial plane ($r\lesssim 10\,r_\mathrm{g}$), while slightly reducing jet sheath temperature. Cooling leads to a dimmer disk, more extended and brighter jets, and reduced total flux. For a given accretion rate, cooling reduces the high-frequency flux. Time variability originates primarily from the midplane in both non-cooling and cooling cases and decreases as accretion rates rise. Although currently below the dynamic range of EHT observations, the features identified in this study could be resolved by next-generation arrays such as the ngEHT.

We report the detailed jet properties of eROSITA Final Equatorial Depth Survey (eFEDS) J084222.9+001000 (hereafter ID830), a radio-loud super-Eddington quasar at $z=3.4351$, revealed by Very Long Baseline Array (VLBA) observations at 1.6 GHz, 4.9 GHz, and 8.2 GHz. Thanks to the high spatial resolution of the VLBA, we successfully resolve a parsec-scale core-jet structure of ID830, and find a well-collimated jet extending over $\approx 745$ pc, making it the most distant and one of the very few currently known radio-loud quasars with a resolved jet associated with super-Eddington accretion. The physical scale and evolutionary track of ID830 differs markedly from the low-$z$ analogues, such as nearby radio-luminous high-Eddington narrow-line Seyfert 1 galaxies, suggesting that this source represents a distinct high-$z$ population compared to previously known samples, with important implications for AGN feedback in early galaxy evolution. We also find that the jet has a relativistic speed of $v \gtrsim 0.19c$ and a modest viewing angle of $\phi \lesssim 79^\circ$ to the line of sight, although its emission is not significantly Doppler-boosted ($\delta \sim 1$). This provides the first evidence that such a relativistic and collimated jet can be produced over several hundred parsecs in the super-Eddington phase, lasting for at least $10^{3\text{-}4}$ yr. Our results call for further theoretical and numerical studies to understand the physical processes required to sustain such large-scale collimation in super-Eddington accretion, which remains a missing piece.

We report the detection of strong, highly circularly polarised, and rotationally modulated radio pulses from the early-mid M dwarf ASKAP J181335-604720, based on strictly simultaneous radio and optical observations with the Australian Square Kilometre Array Pathfinder (ASKAP) and the Transiting Exoplanet Survey Satellite (TESS). The ASKAP data reveal recurrent broadband radio pulses across 800-1088 MHz, with peak circular polarisation fractions approaching 100%. A dominant period of P = 5.607 +- 0.003 h is derived from the TESS light curve using a Lomb-Scargle analysis, which we interpret as the stellar rotation period. When phase-folded on this period, the radio emission is confined to narrow phase intervals and recurs at fixed rotational phases, consisting of a dominant pulse and a weaker secondary component. No contemporaneous optical flares are detected at the epochs of the radio pulses in the simultaneous TESS data. Even under conservative assumptions, the inferred brightness temperature exceeds Tb > 1.8 * 10^12 K, ruling out incoherent emission mechanisms. Combining the observed characteristics, we interpret the emission as coherent electron cyclotron maser (ECM) radiation arising from the stellar magnetosphere, with the observed frequencies implying local magnetic field strengths of at least several hundred Gauss in the radio-emitting region. This work provides a clean, well-constrained, and strictly simultaneous radio-optical case, demonstrating that coordinated radio and optical observations offer a powerful means of distinguishing magnetospheric coherent radio emission from flare-associated coronal activity in M dwarfs.

The study of the origin of life on Earth has been broadened due to panspermia models that suggest that early life may have been transferred between planets. Mars likely once had conditions that could support life, and it is interesting therefore to consider the question of early interplanetary transfer of life from Mars to the Earth. Endospore forming bacteria are ideal candidates for these studies as they can withstand harsh environmental conditions. For this reason, the idea that early life could have been delivered to Earth on Martian ejecta in the late Hadean period has gained considerable interest. To assess this, we have performed a series of both biological and astrophysical experiments. We exposed endospores shielded by a lysed colony of bacteria to extended UVC irradiation under a variety of rotation regimes, to simulate interplanetary exposure on ejecta with a variety of rotation periods. We also performed detailed n-body simulations of particles ejected from Mars at both perihelion and aphelion, finding that Martian ejecta can reach the Earth on timescales of just a few years - suggesting that, with ejection at the ideal time, transfer could occur within one year. Taken together, this study suggests this interplanetary transfer of biologically viable material from Mars to Earth is plausible under favourable conditions.

We present a search for new spider pulsar candidates through multi-wavelength cross-matching, including $\gamma$-ray, X-ray, and optical data. A search for sinusoidal-like optical modulations in TESS data of 183 eROSITA X-ray sources coincident with unassociated Fermi-LAT gamma-ray sources led to the identification of four promising spider pulsar candidates. We found optical variability periods ranging from 5 to 13 hours. All candidates display smooth sinusoidal-like phase light curves, similar to what can be expected from ellipsoidal variation; one shows double-peaked profiles indicative of harmonics. The absence of sharp minima, which are often found in black widow systems due to irradiation, together with their optical magnitudes of about G~14, suggests these sources are more likely redback-type binaries. One of the Fermi-LAT counterparts is included in a machine-learning catalog of unassociated gamma-ray sources, with relatively high pulsar probabilities. We also identify potential Gaia counterparts for several sources and estimate their distances and luminosities where parallax measurements are available. Future observations, including further spectroscopic and multi-wavelength studies, are needed to fully characterize these systems.

We investigate the sub-kiloparsec (sub-kpc) molecular ISM structure and its relation to the galactic environment and star formation in the barred spiral galaxy M83 (NGC 5236). We employ the gas density histogram (GDH), which quantifies molecular gas surface density within $550~\mathrm{pc}\times550~\mathrm{pc}\times100~\mathrm{km~s^{-1}}$ cells. The GDHs are well described by one or two log-normal components, corresponding to the lower and higher-surface-density molecular components, referred to as L-LN and H-LN, respectively. The L-LN mass ($M_{\rm L}$) is relatively uniform across the disk, whereas the H-LN mass ($M_{\rm H}$) is highly structured and traces spiral arms. The fractional contribution of the H-LN component ($f^{\prime}_{\rm H}$) shows coherent structures across the disk and is enhanced along spiral arms, consistent with our previous Milky Way results. Moreover, while the L-LN correlates only weakly with star formation rate surface density ($\Sigma_{\rm SFR}$) and shows a steep Kennicutt-Schmidt (KS) relation with surface-density saturation reminiscent of atomic gas, the H-LN exhibits a tighter, nearly linear correlation similar to the conventional molecular KS relation. These results provide direct evidence that the molecular gas in M83 consists of multiple components. Star formation is more closely linked to the H-LN component, whereas the L-LN component appears to represent a more spatially extended molecular gas. Overall, our results suggest that galactic environments control the relative contribution of the two LN components, and that enhanced H-LN contribution is associated with elevated star formation activity.

We reveal the effect of the spacetime parameters on the accretion morphology formed through the BHL mechanism around a slowly rotating disformal Kerr black hole. Thus, we investigate the measurable signatures of these parameters on the hydrodynamical morphology and the timing behavior of the accreting flow. It is shown that even weak disformal deviations from the Kerr solution modify the shock-cone structure, enhance the density in the post-shock region, and produce coherent oscillations in the accretion rate. The Kerr model produces coherent peaks at 42.99 Hz and 68.13 Hz, and these frequencies are consistent with the high-frequency QPOs observed from the source GRS 1915+105. In the models where the deviations from the Kerr solution are weak, low-frequency QPOs are produced and found to be coherent. These frequencies also fall within the frequency range observed in Galactic black-hole binaries. On the other hand, the models with large deviations from Kerr can be used to explain observational results that are more irregular, broad-band, and contain multiple peaks. In addition, by using inverse-mass scaling in this work, the numerically calculated frequencies are also compared with observations of intermediate-mass and supermassive black holes. In particular, the disformal black-hole models are found to be consistent with the observational results obtained from the sources M82 X-1, NGC 5408 X-1, and RE J1034+396. This comparison also allows the possible black-hole mass range of observed sources to be inferred from the relation between simulated and observed frequencies. This makes BHL accretion in disformal Kerr geometry a powerful framework for connecting modified-gravity black-hole spacetimes with observable QPO phenomenology.

The 4-m International Liquid Mirror Telescope (ILMT) offers a unique opportunity to detect transients in a narrow strip of sky. We explore ILMT's potential to detect astrometric and photometric transients at various ecliptic and galactic latitudes. We inspected CCD frames observed at both low and high ecliptic and galactic latitudes during the commissioning phase and the November 2023 - May 2024 observation cycle, respectively. We analysed these images using both visual inspection and the ILMT's transient detection and candidate classification pipeline. In the low ecliptic and galactic latitude field, we detected more than 500 transient candidates. We cross-matched these with the Minor Planet Checker (MPC) database, identifying 504 catalogued asteroids, all with predicted V-magnitudes brighter than 24 mag, representing a total of 152 distinct asteroids. We performed the same steps on the high ecliptic and galactic latitude field, detecting 30 MPC-catalogued asteroids, and one newly discovered photometric transient, named AT 2024fxn. We present the positions, trajectories, and magnitudes of the detected asteroids observed in the SDSS g', r', and i' spectral bands and compare results from both fields. We explore the lightcurve of AT 2024fxn, which shows partial compatibility with a supernova (SN) hypothesis, while the data invites further insights.

Neutron monitors (NMs) are ground-based devices designed to measure cosmic-ray count rates by monitoring atmospheric neutrons from cosmic-ray showers. We present results from new electronics that have recorded cross-counter time delay histograms for the {\it Princess Sirindhorn Neutron Monitor} (PSNM) at the summit of Doi Inthanon, Thailand. From these histograms, we have extracted the cross-counter leader fraction ($L$) and corrected it for atmospheric effects. For large counter separation, we measure nearly constant $L\approx0.997$, implying that 0.3\% of counts in one counter are temporally associated with later counts on a given distant counter. Monte Carlo simulations confirm that individual secondary particles cannot account for the associated counts at large counter separation, which instead requires a contribution from multiple secondary particles in the same cosmic ray shower that is apparently independent of distance over 3 to 7.5 m. We infer that $\approx$4.5\% of PSNM counts are associated with a later count in at least one of its 18 counters from a different secondary particle in the same shower. Monte Carlo simulations of atmospheric showers and NM yield functions can be validated using our measurements of neutron multiplicity across counters and the contributions of single and multiple secondary particles. These measurements also improve understanding of the single-counter $L$, which has been used for precise tracking of cosmic-ray spectral variations and extending the range of NM observations to higher energies.

The redshift range z=2-4 marks a critical transition in large scale structure formation, where the dynamically unrelaxed progenitors of local massive clusters undergo rapid stellar mass assembly. We report the discovery and physical characterization of a highly significant protocluster, PC J1001+0214, at z=2.96 within the COSMOS field. Leveraging the multi-wavelength COSMOS2025 catalog with exceptional photometric precision in conjunction with JWST/NIRCam wide-field slitless spectroscopy (WFSS) from the COSMOS-3D program, we robustly identify a cosmic overdensity with $\delta=2.66$. The structure comprises 131 member galaxies, including 21 spectroscopically confirmed members (identified primarily via He I $\lambda10830$ emission) and 110 high-fidelity photometric members. A comparative analysis against a mass-complete coeval field sample reveals a statistically significant +0.2 dex shift in the stellar mass distribution of protocluster members, signaling accelerated mass assembly in the dense environment by $z\sim3$. While the protocluster population broadly follows the star-forming main sequence, low-to-intermediate mass galaxies ($\log_{10}(M_{*}/M_{\odot})\le9.7$) exhibit a measurable star formation rate (SFR) enhancement of +0.11 to +0.15 dex. Crucially, the quiescent fraction remains extremely low and indistinguishable from the field sample, implying that environmental quenching mechanisms have not yet become dominant. Furthermore, a preliminary size-mass analysis hints at elevated morphological compactness among protocluster members at fixed stellar mass. These results suggest that PC J1001+0214 represents a growth-dominated protocluster phase in which environmental effects are already detectable, primarily through accelerated stellar mass assembly and ongoing growth rather than through strong quenching.

Mildly recycled pulsars are neutron stars partially spun up through relatively short mass-transfer phases, typically with massive carbon-oxygen (CO) or oxygen-neon-magnesium (ONeMg) white dwarf companions. PSR J2338+4818, a mildly recycled pulsar, was discovered with the Five-hundred-meter Aperture Spherical Telescope (FAST). As a pilot study on the formation and evolutionary pathways of mildly recycled pulsars, we present the updated timing solution for PSR J2338+4818 and examine its single pulses and scintillation properties. Aided by the sensitivity of FAST, the single pulses of PSR J2338+4818 were systematically studied. 27,228 single pulses with S/N > 7 have been detected in our observations. For the FAST ultra-wideband observation on MJD 61045, the receiver was still in the technical commissioning phase, and then only a preliminary single-pulse search was performed. Pulse nulling was examined using a Markov Chain Monte Carlo (MCMC) method, but no evidence for nulling was found. The possible long-term nulling reported by previous studies did not occur in any of our observations in either the 1.0 to 1.5 GHz band or the 300 to 600 MHz band. Interstellar scintillation is evident in our observations. The measured scintillation timescales and bandwidths range from 2.93 to 25.26 minutes and 1.68 to 27.41 MHz, respectively. In all observations, no clear scintillation arc was found in the secondary spectra of PSR J2338+4818.

We analyze PBH populations using state-of-the-art modeling of PBH binaries, deriving the strongest bounds on PBH abundance in the $0.6-100 M_\odot$ range from LIGO-Virgo-KAGRA O4a data and demonstrating sensitivity in the $10^{-4}-10^4 M_\odot$ range. The constraints are dominated by resolvable PBH mergers, while the associated gravitational wave background provides complementary but weaker limits. Allowing PBHs to account for a subset of the cataloged events slightly relaxes these bounds. However, a joint fit with astrophysical black holes shows no compelling evidence for a PBH contribution.

Accurate forecasting of megaelectron-volt (MeV) electrons in the outer Earth's radiation belt, which can pose significant risks to satellites, is essential for risk mitigation and spacecraft operations. We develop a machine-learning-based pipeline for forecasting 1-MeV electron flux variations, focusing first on a 6-hour forecast horizon. Using precipitating electrons measured by POES NOAA-15, near 1-MeV electron flux measured by GOES, solar wind measurements near L1, and geomagnetic activity indices as inputs in 2013-2023, we train algorithms including linear regression, 1-D convolutional and long short-term memory neural networks, and Transformer-Encoder to forecast 1-MeV electron flux in McIlwain's L-shells between 2.8 and 6.0 with 0.1 bin resolution. Particularly, we exploit the timeseries foundation model TimesFM for (1) a zero-shot prediction and (2) a hybrid application involving the ridge regression on the past dynamic covariates combined with the TimesFM inference on the residuals. Using data from January-June 2024 as an out-of-sample test, we find that the hybrid application of TimesFM, named TimesFM+Cov, yields the best results with an average R2 of 0.9 across L-shells, compared to an average R2 under 0.78 for all other models. The R2 of TimesFM+Cov remains above 0.9 for L-shells between 2.8 and 4.7 and drops to 0.77 at L=6.0, indicating improvements of 12% at the lowest L-shell and 48% at the highest L-shell compared to our second-best models. Our work offers an alternative perspective on how a pretrained foundation model could be adapted for space weather forecasting.

Massive stars, the progenitors of neutron stars and black holes, play a crucial role in shaping the chemical and radiative properties of entire galaxies through their winds and explosive deaths. Stellar pulsations are a common phenomenon in massive stars and asteroseismology -- the study of such pulsations -- provides crucial constraints on the physics of massive star interiors. The excitation of heat-driven pulsations in massive stars is expected to depend on a star's metallicity, but this remains largely uncalibrated in evolution models due to a lack of a sufficient observations. While TESS has dramatically improved the statistics for Galactic massive stars, obtaining TESS light curves for low-metallicity massive stars beyond the Milky Way is challenging, due to their faintness and heavy crowding. In this paper, we present a novel point-spread function (PSF) based light curve extraction method called {\sc Lemons}, which overcomes these challenges. We also demonstrate the limitations of the often-used simple aperture photometry (SAP) method that can provide heavily contaminated light curves. With this new technique, accurate light curves of 91 SMC massive stars in the BLOeM sample are extracted. They reveal a variety of variability types including indications of binarity (e.g. eclipses and ellipsoidal modulation) and stellar pulsations. They also enable us to investigate stochastic low-frequency (SLF) variability for massive stars in the SMC. Furthermore we demonstrate how the morphology of SLF variability probes a star's location in the Hertzsprung--Russell diagram, which appears similar to Galactic massive stars thus indicating that the underlying physical mechanism could be insensitive to metallicity.

Radio detection of extensive air showers enables near-continuous observation and precise measurements of the shower geometry and the depth of shower maximum, $X_{\rm max}$. Beyond $X_{\rm max}$, the longitudinal shower development follows a Universal Shower Profile (USP), whose shape parameters contain information on primary mass composition and hadronic interaction models. While previous studies have focused on event-by-event reconstruction of profile parameters, in this work we investigate the reconstruction of the average USP using radio interferometric techniques. Using Monte Carlo simulations, we reconstruct the average longitudinal profile directly from radio data and extract the Gaisser-Hillas shape parameters $(R,L)$. We find that the radio-derived average profiles provide enhanced separation between primary masses and hadronic interaction models compared to that obtained from fluorescence-equivalent longitudinal profiles. These results demonstrate that radio interferometry can access higher-order information on the longitudinal shower development and that the use of the average USP significantly improves the sensitivity to composition and hadronic interaction studies at ultra-high energies.

In this work, we study the effects of accretion on the primordial black hole (PBH) reheating scenario. PBHs could form from primordial fluctuations. If they have the right mass and abundance, they could dominate the Universe and complete the reheating entirely through Hawking radiation. We find accretion effects on the BH can not only increase the BH mass, but also prolong such early matter domination. The consequence of the accretion is further investigated using isocurvature induced gravitational waves (GWs), which are generated right after the sudden evaporation of the BHs from the oscillation of the gravitational potential. Big Bang nucleosynthesis limits on the energy density of the GWs put important constraints on the PBH domination scenario. Inclusion of accretion shifts such constraints significantly towards smaller formation mass and smaller initial abundance. Furthermore, the PBH could undergo mergers leading to extended mass functions. We find similar shifts in the allowed parameters with the inclusion of accretion for the merger constraint. We find the constraints from isocurvature GWs typically stronger than the constraints from mergers.

XSPECT instrument on-board XPoSat mission is a soft X-ray spectrometer sensitive in the energy band 0.8$-$15 keV. XSPECT has observed several bright neutron star low mass X-ray binaries since launch. Two well known sources, Aql X-1 and 4U 1608-52 which are soft X-ray transients, were observed by XPoSat during the decay phase of their recent outbursts in September 2024 and February 2025 respectively. During XSPECT observations, 4U 1608-52 exhibited a superburst which is a long duration thermonuclear burst, believed to be triggered by carbon burning. We carry out a detailed spectro-temporal analysis of the superburst, tracing its onset, rise, and decay over the next several hours. Using time-resolved spectroscopy, we probe the spectral evolution of the source and find that the persistent emission is suppressed during the superburst and the emission can be described by a gradually cooling blackbody component. The superburst was preceded by a precursor burst which is a normal type-I X-ray burst. We also observe a type-I burst $\sim 5$ days after the superburst, indicating resumption of burst activities which is typically quenched after a superburst. Aql X-1 also exhibited a type-I burst during XSPECT observations. The persistent emission of both the sources can be fitted using a combination of blackbody and disk blackbody emission or, alternatively, using a disk Comptonized by an optically thick plasma. Using the latter model, we find a clear flux dependence of the Comptonization parameters, with both the sources exhibiting harder spectra at higher accretion rates.

Blazars are conventionally classified into BL Lac objects and flat-spectrum radio quasars (FSRQs), commonly associated with radiatively inefficient and efficient accretion onto supermassive black holes, respectively. Increasing evidence, however, suggests that this division is not strictly discrete, with several sources occupying intermediate or transition-like spectral states. Since the photon index traces the evolution of the non-thermal emission spectrum, it provides a useful probe of jet energetics, particle acceleration, and radiative processes across blazar populations. Using multi-mission observations from Swift-XRT, Swift-BAT, NuSTAR, ROSAT, Chandra, XMM-Newton, NICER, AstroSat, TeVcat, and archival VizieR compilations, we investigate the distribution and long-term evolution of X-ray and TeV photon indices across multiple blazar subclasses. We identify broad overlap regions in photon-index space linking EHBL, HBL, IBL, LBL, and FSRQ-like populations, suggesting a continuous rather than sharply separated spectral distribution. Across X-ray instruments, the intermediate regime is concentrated near $\Gamma_{\mathrm{X}} \approx 2$, typically spanning $\Gamma_{\mathrm{X}} \sim 1.5$--$2.2$. Multi-epoch observations reveal substantial intra-source spectral evolution, including stochastic variability in Mrk~421 and state-dependent transitions in OJ~287. Several sources exhibit spectral-index changes of $\Delta \Gamma_{\mathrm{X}} > 0.5$, while occupancy of the intermediate regime reaches $\sim22$--$43%$ depending on instrument. Multiple objects repeatedly traverse spectral regions connecting traditionally distinct subclasses, supporting their interpretation as candidate transition-like blazars. Overall, the results favor a framework in which blazar subclasses represent overlapping and evolving spectral populations driven by long-term changes in jet emission and radiative processes.

Simulations and observations now probe comparable redshift regimes with unprecedented accuracy, enabling direct consistency tests through forward modeling. In a previous work, we identified a faint-end discrepancy between observed and simulated near-infrared galaxy counts in CANDELS GOODS-South. Here we investigate whether this tension originates from the forward-modeling procedure or from limitations of the underlying simulations, and we characterize the galaxy populations responsible for the tension. Using the FORECAST forward modeling code, we generated ten independent lightcone realizations and mock CANDELS images from the TNG100 and EAGLE simulations. We compared both the intrinsic lightcone catalogs and the mock-image detections with observations, testing dependencies on field and redshift, and validating the pipeline through stellar mass and multi-band analyses. The faint-end deficit is present in all CANDELS fields and appears at z>3 in both simulations. GOODS-South counts corrected for completeness exceed intrinsic simulation counts already at the 50% completeness limit, indicating that the missing population is not simply hidden below the detection threshold. Increasing the depth of mock images recovers the counts near the peak but overpredicts the faintest sources, showing that depth alone cannot resolve the discrepancy. Structural analyses reveal that compact galaxies with bright central cores observed in GOODS-South are underproduced in simulations, which instead favor diffuse low-surface-brightness systems. We conclude that the discrepancy arises both from detection losses of diffuse galaxies and, more fundamentally, from the inability of current hydrodynamical simulations to produce enough faint compact galaxies at z>3. This tension points to the need for improved modeling of early star formation, feedback, and dust treatment.

The saturated, boosted, and disrupted magnetic braking (SBD MB) model is an empirical prescription that has recently gained support from close-binary observations. Different boosting ($K$) and disruption ($\eta$) parameters appear necessary for different systems, but their physical origins remain uncertain. We aim to identify the mechanisms that boost magnetic braking (MB) and cause its disruption at the fully convective boundary in cataclysmic variables (CVs). We modelled CV evolution with MESA and compared the results with observed CV properties. We computed the convective turnover time ($\tau_c$) directly from the donor's structure rather than adopting empirical relations. We also included irradiation from the accreting white dwarf, which heats the donor's outer layers and can drive additional winds that enhance MB. The structure-based $\tau_c$ calculation reveals a pronounced spike as the donor approaches full convection, which drives the disruption parameter $\eta$ and initiates the period gap in CVs. The outcome of irradiation is sensitive to the accretion, irradiation, and wind efficiencies, all of which are poorly constrained from observations. Despite these uncertainties, plausible parameter choices allow irradiation-driven winds to provide the required boost $K$ during accreting phases. We refer to the combined prescription as the i$\tau$SBD MB model and find that it yields evolutionary tracks broadly consistent with the main CV properties. Our i$\tau$SBD MB framework offers a physically motivated interpretation of the empirical boost and disruption factors in SBD MB for CV evolution. We suggest that the convective turnover time spike at the fully convective boundary may drive MB disruption for fast-rotating stars in the saturated regime, while irradiation-driven winds may be the dominant mechanism boosting MB in accreting binaries and other strongly irradiated close systems.

We propose a new way to observe cosmic-ray-induced air showers by imaging the radio emission. With simulations we demonstrate key features for imaging the radio-wave emission from air showers, which show similarities to the well-established atmospheric imaging Cherenkov technique in gamma-ray astronomy. In addition, we find that imaging the emission with a camera, consisting of multiple antennas, resolves emission that is not accessible to a single antenna. Pursuing this technique, with a camera operating in the GHz frequency domain, might be beneficial ultra-high-energy gamma-ray astronomy and other studies that include detailed observations of air showers.

We report the discovery of a correlated energy-dependent time lag and degree of polarization for Swift J1727.8-1613 during its 2023 outburst. The energy-dependent time lag is measured around the type-C quasi-periodic oscillations (QPO) observed by IXPE on 2023-09-07, while the degree of polarization is obtained from energy-resolved polarimetric measurements. The Spearman correlation coefficient was found to be 0.8, with a null hypothesis probability of 4.2\%. Furthermore, the correlation value drops as the quality factor, or Q value, of the observed QPO frequencies decreases. The spectral properties of Swift J1727.8-1613 are analyzed using simultaneous Insight/HXMT data. Thereafter, we present model-independent theoretical arguments to show that processes other than inverse Comptonization also contributes to both the observed polarization and time lags. This correlation may therefore point to additional mechanisms contributing to the connection between the spectral, temporal, and polarimetric properties of black hole binaries in their hard state.

Stacking (i.e. averaging) the value of a given astrophysical field around sources allows us to detect new cosmological signatures, such as the kinematic Sunyaev-Zel'dovich, and gain insight on the astrophysical properties of galaxies and their environment. Further information may be gained by orienting these stacks along preferred axes defined by a local directed field, such as the transverse galaxy velocities, galaxy shapes, or the local tidal forces. Examples of this are searches for the moving lens effect, the detection of dipole signatures, or the study of cosmic filaments. Here we show that all directional stacking signals may be reconstructed, without loss of information, in terms of the cross-power spectrum between the quantity of interest and the $E$ and $B$ modes of the spin field used to define the preferred axes weighted by the local galaxy density. The power spectrum approach has several practical advantages, in terms of speed, finite-resolution effects, data visualisation, and combination with other cosmological probes. We also argue that, in some cases, such as stacking using velocities or tidal forces reconstructed from the density field, the recovered signal may be dominated by information that is already present in the cross-spectrum between the target field and the galaxy overdensity itself.

The gamma-ray emitting binary stellar system LS I +61 303 was observed with the Imaging X-ray Polarimetry Explorer (IXPE) on two successive orbits over orbital phases of 0.74 to 1.05. Polarization is detected at a significance of 4.2$\sigma$ with an average polarization degree of $13.1\% \pm 3.0\%$ in the 2-8~keV band after background subtraction. This is the second detection of polarization of the X-ray synchrotron emission from a gamma-ray binary and, again, suggests that the magnetic field in the particle acceleration region has a significant ordered component. The orbital motion on the sky of LS I +61 303 is not well determined, which leads to ambiguity in interpretation of the X-ray electric vector polarization angle (EVPA) measurement. Use of orbital elements determined via radial velocity measurements combined with radio imaging of variable nebular emission, suggests an offset between the X-ray EVPA and the compact object-massive star axis on the order of ~30$^{\circ}$. Such an offset could be produced by Coriolis forces due to binary motion. Use of two different sets orbital elements determined via optical polarimetry suggest either no offset or a perpendicular orientation, but require an unexpectedly high inclination. Use of orbital elements derived from modeling of the keV/TeV light curves suggest good alignment between the X-ray EVPA and the compact object-massive star axis. Such alignment was found for the gamma-ray binary PSR B1259-63. If the same physical situation holds for LS I +61 303, that would favor the orbital elements derived from the keV/TeV light curves.

The Milky Way in-situ halo, also known as the Splash, consists of old (Age $>$ 10 Gyr), metal-rich ([Fe/H] $> -0.7$), high-$\alpha$ stars, i.e., thick disk-like chemistry, on halo-like orbits (eccentricity > 0.6). Its origin is linked to stars formed in the disk and dynamically heated by either internal or external agents. In this work, we investigate its low-$\alpha$ counterpart, the low-$\alpha$ Splash, motivated by recent findings of an old thin disk population. We conjecture that any mechanism capable of heating disk stars should affect both of present-day high- and low-$\alpha$ old populations. Using data from the APOGEE DR17 spectroscopic catalog, we identify metal-rich low-$\alpha$ stars with halo-like kinematics similar to those of the classical high-$\alpha$ Splash. We investigate their possible heating mechanisms using the GASTRO suite of simulations, which allows us to explore the effects of star-forming clumps as well as a major merger in the proto-disk of a Milky Way analog galaxy. Our main results show that only clumpy Milky Way models are able to produce Splash populations through scattering by clumps in the early Galaxy, including the low-$\alpha$ counterpart, whereas the model including only the merger and without an early clumpy phase fails to produce these populations. In the models, the low-$\alpha$ Splash corresponds to a subset of the old thin disk that was dynamically heated by the same mechanism responsible for the formation of the high-$\alpha$ Splash.

The cosmic distance-duality relation (CDDR), expressed as $ D_L/D_A(1+z)^{-2}=1 $, is a fundamental relation in cosmology connecting luminosity distance ($ D_L $) and angular diameter distance ($ D_A $). Any departure from this relation would indicate new physics such as photon non-conservation, cosmic opacity, or non-metric gravity. We perform a stringent, model-independent test of the CDDR using a matched sample of 38 galaxy clusters from the Bonamente et al. compilation and Type Ia supernovae from the Pantheon+ sample. Employing the matched-pair technique, we simultaneously constrain the CDDR-violation parameter $ \eta $ and a possible redshift evolution of the SNe Ia absolute magnitude, parameterized as $ M_B(z)=M_0 + \varepsilon z $. We assess the robustness against matching tolerance and further supplement the analysis with DESI 2024 BAO measurements. Our results yield $ \eta = 0.050^{+0.348}_{-0.307} $ and $ \varepsilon = -0.184^{+0.724}_{-0.574} $ (68% CL), showing no statistically significant evidence for CDDR violation or SNe Ia evolution. The conclusions remain unchanged with stricter matching criteria and the inclusion of DESI BAO data. We also derive a cosmology-independent calibration of $ M_0 = -19.460^{+0.126}_{-0.124} $ mag. The standard cosmological model remains robust under this model-independent scrutiny.

The Milky Way hosts a flattened nuclear stellar disc (NSD) that dominates the gravitational potential in the inner few hundred parsecs. Whether the NSD is purely axisymmetric or contains a nuclear bar remains an open question. We test for the presence of a nuclear bar using kinematic diagnostics by combining line-of-sight velocities from the KMOS NSD survey with proper motions from VIRAC2 to construct the $ (v_\ell, v_\mathrm{los}) $ velocity ellipse. After applying strict quality cuts to minimise contamination from large-scale bar stars, we measure the vertex deviation $ l_v $ and anisotropy $ \beta $ for several subsamples. For our primary sample ($ |\ell| < 0.9^\circ $, $ -0.4^\circ < b < 0.25^\circ $, $ \mathrm{[Fe/H]} > -0.3 $), we find a significant negative vertex deviation $ l_v = -54.8^{+13.1}_{-14.8}\,^\circ $ with moderate anisotropy $ \beta = 0.16^{+0.08}_{-0.05} $. A subsample restricted to the innermost four fields yields an even stronger signal with $ l_v = -64.3^{+12.1}_{-12.2}\,^\circ $ and $ \beta = 0.38^{+0.12}_{-0.07} $. The direction of maximum velocity dispersion is oriented along Galactic longitude, opposite to that observed in large-scale bar-dominated samples. These signatures are robust against extinction-driven incompleteness, primary-bar contamination, and the choice of metallicity threshold. They are inconsistent with an axisymmetric NSD or one oriented orthogonally to the primary bar, but match expectations for a nuclear bar oriented at $ \alpha \approx 60^\circ $-$75^\circ$ to the Sun-Galactic-Centre line with its near side pointing toward positive Galactic longitude. While definitive confirmation awaits larger and more precise samples from upcoming surveys, our results provide the first kinematic indication of a possible nuclear bar in the Milky Way.

We perform a detailed study of the energetics for a well-observed solar eruption and flare that occurred on 28 October 2021. This event included a GOES class X1.0 flare, a global EUV wave, and a coronal mass ejection that reached speeds of >2000 km/s. The event was observed from a variety of spacecraft in NASA's Heliophysics System Observatory, including multiple missions near Earth, STEREO-A off the Sun-Earth line, and Solar Orbiter, near the Sun-Earth line at about 0.8 au. Using remote sensing, in situ observations, and in some cases scaling laws based on previous observations, we characterize the following quantities: free magnetic energy, energy in non-thermal electrons, energy in non-thermal ions, bolometric energy, energy deposited in the chromosphere, thermal energy radiated in the flare loops, energy dissipated by the EUV wave, CME kinetic and gravitational potential energy, CME energy flux in the heliosphere, and the energy partition in the CME shock. We find that the total energy released during the event is consistent with estimates of the pre-event stored magnetic energy, and the CME kinetic + potential energy dominates the energy partition.

Neutron star mergers are today considered a major production site for rapid neutron capture elements. While the bulk of the matter escapes at fast, but non-relativistic velocities (${\sim} 0.2\,c$), a small amount of the dynamically ejected mass reaches mildly relativistic velocities (${\gtrsim}0.6\,c$). It has been suggested earlier, that in such ejecta parts neutrons may avoid being captured and that their decay could power an early blue precursor to the main kilonova event. Here we study in detail the nucleosynthesis in such fast ejecta with nuclear network calculations along both parametrized and numerical relativity trajectories. We find that the nucleosynthesis can be divided into three channels, in one of which a substantial amount of free neutrons survives when the main r-process has frozen out. We provide a (semi-)analytical model for surviving free neutrons which agrees very well with the network calculations. If the mass fraction of the free neutrons exceeds ${\sim} 0.05$, their $\beta^-$-decay dominates the nuclear heating rate between ${\sim} 100$ and ${\sim} 10^4$ seconds. This dominance leads to a pronounced kilonova precursor that should for plausible ejecta parameters be visible for ULTRASAT out to ${\sim}200\,\rm Mpc$. Since at low electron fractions free neutrons can survive even for moderate velocities, mergers with large tidal ejecta, such as asymmetric neutron star mergers or favorable neutron star black hole mergers, may produce particularly bright blue precursors to their subsequent kilonovae.

The possible interaction between the dark components of the Universe (dark matter and dark energy) stands as an attractive alternative to the standard $\Lambda$CDM cosmological model. In this work, we present a novel analysis of three sign-changeable interaction models whose coupling term $Q$ depends explicitly on the deceleration parameter $q$ and is proportional to different energy densities: dark matter, dark energy, and total energy density. To constrain these models, we combine strong gravitational lensing data on two complementary scales: a sample of early-type galaxies acting as lenses and the galaxy cluster Abell 1689. Our results show that the interaction strength $\beta$ depends on the choice of the coupling term $Q$, with all models yielding negative values of $\beta$, indicative of a dark interaction scenario. The $\beta$ values obtained in this work are significantly larger in magnitude than those previously reported using Type Ia supernovae, CMB, and BAO. The strong-lensing constraints indicate a transition to cosmic acceleration at earlier redshifts ($z_t \sim 1.8-2.1$) than that predicted by the $\Lambda$CDM model, while remaining consistent with cosmic chronometer measurements within the reconstructed confidence regions. Therefore, our study shows that strong gravitational lensing data provide an independent and competitive cosmological probe capable of testing interacting dark energy scenarios. The sensitivity of lensing observables to the expansion history enables access to complementary information about dark-sector dynamics beyond standard cosmological probes.

We present the Data Quality Report Builder toolkit, DQRbuild, a suite of data quality tools that have been developed to vet gravitational-wave events in preparation for the fourth LIGO-Virgo-KAGRA observing run. We explain the main functionality and the many scientific tests that we support. To validate the performance of the tools included in the toolkit, we run a series of tests on all significant candidates shared as public alerts in the third observing run to compare against what was manually reported using human intervention. We find that these automated tools can now identify 96% of the problems identified by humans during this previous observing run, with a 24% false alarm rate. We conclude with a commentary on the prospects and potential challenges for fully automating the process of vetting the data quality for gravitational-wave events identified in future observing runs.

Accurate estimation of the covariance matrix of cosmic shear statistics is essential for cosmological analyses using current and upcoming wide-area weak lensing surveys. In this work, we investigate analytical methods for computing the Gaussian covariance matrix of the cosmic shear two-point correlation function (2PCF), taking into account the effects of finite survey geometry. We compute the covariance of 2PCF based on the improved Narrow Kernel Approximation (iNKA), with a projection using the Legendre transformation. We also consider other analytical covariance estimators, the $f_{\mathrm{sky}}$ approximation and the weighted quartic-counts method. We evaluate the accuracy of those analytical methods using the convergence fields with the HSC Year 3 survey mask as a test case. We find that the covariance of the 2PCF obtained by using the iNKA does not reproduce the covariance measured directly from Gaussian simulations. Although the iNKA accurately models the diagonal structure of the harmonic-space covariance, residual inaccuracies in the off-diagonal components propagate through the Legendre transformation and significantly affect the real-space covariance. In contrast, the weighted quartic-counts method shows better agreement with the simulations. Our results demonstrate that accurate modeling of the off-diagonal structure of the harmonic-space covariance is crucial for obtaining reliable covariance estimates of real-space weak lensing statistics in the presence of survey window effects.

We perform hydrodynamical simulations with radially varying resolution to study the effects of stellar feedback on the radial inflow of gas from the Central Molecular Zone (CMZ, $R\sim200$ pc) to the Circumnuclear Disk (CND, $R\sim5$ pc) of the Milky Way. The simulations include a realistic Milky Way barred gravitational potential, a cooling function coupled to a non-equilibrium chemical network, gas self-gravity, star formation, supernova feedback, and radiation feedback from massive stars computed via on-the-fly radiative transfer. Our main findings are as follows: 1) Stellar feedback drives a radial inflow that decreases monotonically with decreasing Galactocentric radius. The time-averaged inflow rate in our fiducial SNRad simulation, which includes both supernova and radiation feedback, declines from $\langle \dot{M} \rangle\sim5\times10^{-3}$ Msun/yr at $R\sim100$ pc, to $\langle\dot{M}\rangle\sim10^{-4}$ Msun/yr at $R\sim10$ pc, to $\langle\dot{M}\rangle\sim10^{-6}$ Msun/yr at $R\sim1$ pc. 2) The total inflow rate can be broken down into two components driven by two distinct mechanisms. First, feedback-driven turbulence redistributes the angular momentum of gas clouds, producing a smooth (secular) transport of mass inward, similar to a Shakura-Sunyaev viscous accretion disk. This component contributes inflow rates that vary from $\dot{M}\sim5\times10^{-4}$ Msun/yr at $R\sim100$ pc to $\dot{M}\sim10^{-7}$ Msun/yr at $R\sim1$ pc. Second, episodic inflow events can transiently increase the inflow rate by several orders of magnitude, reaching $\dot{M}\sim10^{-3}$ Msun/yr over timescales of $\Delta t\sim3$-$5$ Myr at $R=10$ pc. 3) The stellar feedback model significantly affects the episodic inflow but has little impact on the smooth component. Simulations including radiation feedback produce substantially more episodic events than those with supernova feedback alone.

Many searches for continuous gravitational waves face significant computational challenges due to the need to explore large parameter spaces characterized by unknown parameters. Parameter-space metrics are used to predict the relative loss of signal power when the searched parameters differ from the true signal parameters. In this paper we present generalized parameter-space metrics for the $\mathcal{F}$-statistic (a detection statistic used in many searches) that improve upon previous idealized metrics by incorporating realistic effects such as data gaps and varying noise floors. We derive a new marginalized $\mathcal{F}$-statistic metric that is more accurate than the previous averaged $\mathcal{F}$-statistic metric, especially for short coherent segments. We also derive a more accurate semi-coherent metric that properly accounts for the signal-power variability over segments. We provide numerical tests illustrating that the new generalized metrics provide more accurate mismatch predictions than previous expressions. More accurate metrics can result in a reduced number of templates needed for a given search, a feature that could improve the sensitivity of future searches.

Gravitational wave backgrounds from strong first-order cosmological phase transitions are key observational targets predicted by many SM extensions and might be observed by current and future observatories like LISA, the Einstein Telescope or pulsar timing arrays (PTAs). Still, their precise forecast given a specific model remains a challenge. In this article, we present TransitionListener v2.0, a Python framework for precision studies of cosmological phase transitions and their associated gravitational wave (GW) signals. The code provides an end-to-end pipeline from a user-defined scalar potential to GW spectra and signal-to-noise ratios, enabling both benchmark studies and large-scale parameter scans. Version 2 introduces a self-consistent treatment of the transition dynamics, including the evolution of the true-vacuum fraction and its backreaction on the Hubble expansion, as well as a consistent description of reheating during percolation. A direct computation of the mean bubble separation allows to faithfully map to the GW spectral templates from bubble collisions, sound waves, and turbulence stemming from state-of-the-art simulations. TransitionListener includes built-in sensitivity curves for space- and ground-based detectors and PTAs, interfaces to PTA likelihoods, and wrappers for Bayesian model inference and high-dimensional parameter scans. Compared to existing public tools, TransitionListener v2.0 improves the physical consistency and numerical stability of GW predictions across a wide range of models, with particular emphasis on the strongly supercooled and ultraslow transition regime where conventional approximations break down and the most promising GW signals are expected.

Magnetic fields between a supernova (SN) and Earth convert axions into gamma rays. The absence of such a signal in coincidence with SN 1987A neutrinos, using the coherent Milky Way field, provides well-studied constraints on $g_{ap}\times g_{a\gamma}$ (axion-proton times axion-photon couplings) and on $g_{a\gamma}$ alone. We show that the small-scale power of the turbulent magnetic field component boosts axion-photon conversion and, crucially, extends sensitivity to larger masses. The turbulent field components of the Milky Way and of the Large Magellanic Cloud (hosting SN 1987A) yield improvements of up to two orders of magnitude in $g_{ap}\times g_{a\gamma}$. Turbulence should likely impact the sensitivity of other searches based on other axion-photon conversion sites, such as starburst galaxies.

Merging binary black holes exhibit a ringdown phase in which they primarily emit gravitational waves in the shape of damped sinusoids corresponding to quasi-normal modes of the Kerr remnant. In general, each mode carries four degrees of freedom encoding amplitude and phase information. When the modes are excited with equatorial reflection symmetry, as is the case for black hole mergers with spins (anti)aligned to the orbital angular momentum, the symmetry constrains two degrees of freedom. As a result, the relationship between polarization amplitudes and phases in each mode is fixed by the viewing (inclination) angle to the equatorial plane. We use such a constrained model to fit the ringdown signals of both non-precessing and precessing systems such as GW150914 and GW190521, respectively. We show that we can measure the degree of circular polarization and handedness of ringdown signals like those of GW150914, even if only the two LIGO detectors are available; such a polarization measurement can be translated into an inferred source inclination assuming the reflection symmetry above, again using the ringdown signal alone. On the other hand, the constrained polarization model is insufficient to capture the polarization structure of signals from precessing systems, leading to biases in the inferred mode frequencies and amplitudes. We explore the magnitude of this effect by fitting GW190521-like injections with the restricted model, finding weaker predictive accuracy relative to the arbitrary-polarization model and potentially significant systematic biases. As our detectors continue to improve, using the correct polarization model is increasingly important to avoid biased ringdown measurements.

The growing demand for extensive data processing is now a standard in many scientific fields. Efficiently distributing data to processing sites and enabling seamless sharing has become crucial. The Open Science Data Federation (OSDF) builds on the success of the StashCache project to establish a global data distribution network. By expanding StashCache, OSDF integrates additional data origins and caches, enhancing accessibility and performance (20 origins and 30 caches), new access methods, and monitoring and accounting mechanisms. Additionally, the OSDF has become essential to the US national cyber-infrastructure landscape due to the sharing requirements of recent NSF solicitations. One use case for the OSDF is the data access to the Big Bear Solar Observatory (BBSO). Integrating the BBSO data into the OSDF provided standard and reliable data access. Moreover, the OSDF caches provide local data worldwide. Using the OSDF and the BBSO data, creating a pipeline to apply image processing techniques to all images from BBSO anywhere on the planet was possible.

We introduce a new series decomposition of the waveform constructed in the spirit of Debye expansions in scattering theory, and we use this to analyse the time-domain response of compact, horizonless bodies to scalar-field perturbations on curved spacetimes. The Debye decomposition separates out direct exterior propagation, surface reflection, and successive transmissions through the interior of a compact body, and it provides an intuitive interpretation of the waveform in terms of geodesic trajectories. By analysing the quasinormal-mode (QNM) content of individual Debye terms, we set out a Debye-QNM description that is complementary to the standard QNM description. With this framework, we examine a scalar field propagating on two illustrative `Schwarzschild star' compact-body spacetimes: a neutron-star-like model \(R>3M\) and an ultracompact object \(R<3M\). We show that the Debye reconstruction matches well with the exact waveform, and that (unlike the standard QNM reconstruction) it converges even at early times, giving an accurate description of all waveform features including the prompt response. In the neutron-star case, the low-order Debye terms mainly describe the ringdown and a non-modal component associated with the sub-threshold branch cut. In the ultracompact case, the Debye series organizes the waveform into a prompt/ringdown contribution followed by a succession of individually resolved echo-like wavepackets. The new Debye-QNM expansion and the standard QNM expansion have complementary spectral interpretations: the former identifies modes in individual propagation channels, whereas the latter describes collective resonances that are resummations of the former. This distinction clarifies how echo-like structures emerge from repeated interior propagation, and how pole and branch-cut contributions enter the time-domain signal.

We investigate systems containing objects with negative mass (NMOs). In a system consisting of one object with positive mass and one NMO, a bound state exists even though the force exerted by the NMO on the object with positive mass is repulsive. Unlike a standard system consisting of two objects with positive mass, the gravitational waves emitted from this system exhibit a decrease in frequency and amplitude over time. We propose a model of the time evolution of the Ellis-Bronnikov wormhole, along with a formulation that eliminates the ghost that appears when constructing the Ellis-Bronnikov wormhole, a candidate for an NMO. Furthermore, numerical simulations are performed to obtain the optical appearance of such NMOs. The observed luminosity is also compared with the Schwarzschild black hole and with the Simpson-Visser wormhole, leading to clear differences in the photon ring substructure around the central object.

Inflaton decay inevitably emits gravitons through bremsstrahlung during reheating. We show that the soft part of this emission amplitude, fixed by Weinberg's soft-graviton theorem, becomes an irreducible stochastic gravitational-wave (GW) background after accounting for cosmological evolution. The theorem fixes the infrared branch of the spectrum, $\Omega_{\rm GW}\propto f$, independently of the microscopic operator responsible for inflaton decay, while the normalization is controlled by the hard inflaton decay rate and by a phase-space factor. We carry this out for inflaton $n$-body decays, including the phase-space integrals, finding that the maximum of the spectrum scales as $2/n$ relative to the $n=2$ case. The signal can reach $\Omega_{\rm GW}h^2\sim \mathcal O(10^{-17})$ at frequencies above the GHz scale. This predicts a stochastic graviton floor from perturbative reheating: a larger signal would require either other processes beyond perturbative bremsstrahlung or inflationary scenarios beyond conventional single-field slow roll.

We introduce a Stern--Gerlach (SG) spin-kinetic control parameter for magnetic reconnection. The fully projected branch parameter, $\Xi_0=<Z>/r_L$ compares the SG cross-sheet displacement accumulated during a diffusion-region transit with the relativistic Larmor radius. For an ensemble or partially participating population the relevant effective parameter is $\Xi_{\rm Eff}=P_{eff}\Xi _0$, where $P_{eff}$ represents the surviving branch weight or effective spin/moment projection. Evaluating $\Xi_{\rm Eff}$ across representative space and astrophysical environments reveals a robust hierarchy: SG transport is negligible in the magnetotail, solar corona, active galactic nuclei (AGN)/blazar jets, and pulsar-wind nebulae, but becomes transitional to strong in magnetar current sheets and extreme near magnetar surfaces. We further show, using electron--positron particle-in-cell simulations, that the SG force sorts particles by magnetic-moment projection into opposite sides of a Harris current sheet without measurably changing the global reconnection rate in the tested regime. This identifies magnetars as the clearest natural target for strong-field spin-kinetic reconnection ($\Xi_{\rm eff}\gg 1$) near the surface; transitional in the outer magnetosphere), while SG transport is safely negligible ($\Xi_{\rm eff}\ll 1$) in all heliophysical and jet environments considered, and provides a falsifiable framework for assessing where SG physics is relevant.