2026-09-10 | | Total: 33
Binary black holes (BBHs) have been theorised to form through various different channels, but current LIGO-Virgo-KAGRA (LVK) sky locations make the direct study of formation channels a challenge. Association of strongly lensed BBHs to their hosts has been shown to be possible with current detector networks at upgraded sensitivity and survey telescopes like Euclid, with localisations that can reach sub-galactic scales in the source plane. With Hubble Space Telescope-like high-resolution imaging, we report that strong-lensing localisation of BBHs allows us to confidently associate the BBH with specific structures within their host galaxies, specifically, galactic centres (GCs). We demonstrate that a BBH localised to the GC of the galaxy can reject light-tracing BBH populations ($p < 5\%$) in a single observation, and in most cases demonstrate high confidence rejection ($p < 1\%$), providing a conservative basis to exclude light-tracing spatial distributions associated with field-like or non-central formation channels. Such association therefore directly supports formation channels confined to GCs, such as active galactic nuclei and nuclear stellar clusters.
Over the coming years, the number of tidal disruption events (TDEs) is expected to substantially increase with observations from the Vera Rubin Observatory (g and r band) and {\it ULTRASAT} (near UV) wide-field surveys. These future samples have great promise to characterize the bottom end of the massive black hole mass function, but existing detections of intermediate mass black hole TDEs are primarily in X-rays, leaving their optical/UV emission largely unexplored. We present a time- and angle-dependent analysis of the outflow produced by dissipation near pericentre in a three-dimensional end-to-end radiation-hydrodynamics simulation of a TDE by a $10^4 M_\odot$ black hole with realistic parameters, run with the code RICH. We find that outflow anisotropy produces viewing-angle-dependent observables. Towards the poles and the pericentre region, mass-loss rates are low and bolometric luminosities reach $\sim2$--$3$ times the Eddington luminosity. Towards the stream, the properties show a stronger dependence on latitude: the mass-loss rate increases and the bolometric luminosity decreases as the line of sight approaches the orbital plane. These denser regions favour H$α$ and H$β$ emission. Despite these variations, all viewing directions show a common spectral evolution, with an initial soft X-ray flare followed (around $1.25t_{\rm fb}\approx3$~days) by the reprocessing of shock-powered emission into the UV and optical bands. Although the optical/UV luminosities we predict for this TDE are likely too dim for past surveys (e.g. ASAS-SN, ZTF), they are within the detection capabilities of LSST and ULTRASAT to horizons of $\sim 790$ and $\sim 340$ Mpc, respectively, for the brightest viewing directions.
Globular clusters are efficient factories of recycled pulsars, but searches toward high-dispersion-measure (DM) clusters can be strongly limited near 1 GHz by dispersive smearing and interstellar scattering. We present the first MeerKAT S-band (nu ~ 2.4 GHz) pulsar survey of 14 globular clusters. High time- and frequency-resolution observations were searched using segmented acceleration and jerk techniques, followed by candidate folding and targeted folding with available timing ephemerides. We re-detected 39 known pulsars and discovered four new millisecond pulsars in Glimpse-C01: J1848-0129C, D, E, and F. Multi-epoch follow-up enabled preliminary Keplerian orbital fits for J1848-0129C and J1848-0129D. J1848-0129C is an eclipsing MSP in a ~5 d orbit, placing it among long-period eclipsing systems known as huntsman binaries, while J1848-0129D is in a ~3.4 d nearly circular orbit with a massive white-dwarf companion of about 1 solar mass. Two Glimpse-C01 pulsars show large DM offsets from the cluster average. Comparison with other Galactic globular clusters indicates that intracluster DM spreads tend to increase with foreground DM, implying that narrow DM search windows may be sub-optimal for high-DM clusters. From detections and non-detections, we infer a practical single-epoch MeerKAT S-band tied-array detectability scale of about 10-20 microJy. These results demonstrate the value of high-frequency searches for pulsars in strongly dispersed and scattered cluster environments.
There is growing evidence for cosmic-ray acceleration in massive stellar clusters. Furthermore, extended gamma-ray emission suggests that particle transport in the vicinity of their sources is influenced by physical processes markedly different from large-scale diffusion in the Milky Way. We characterize extended gamma-ray emission in the direction of the Westerlund 1 stellar cluster and Kes 41 supernova remnant using > 16 years of data from the Fermi Large Area Telescope (LAT) at energies > 0.8 GeV. We test whether clusters of gamma-ray sources not associated to multiwavelength counterparts are better described by extended emission components. We report the detection of three new extended emission components with soft spectra towards regions of high gas column density in the Galactic plane. One extended component associated with the natal cloud of Kes 41 is statistically preferred over the point source previously reported towards the supernova remnant shell. The other two extended components overlap with neutral gas within ~100 pc from the edge of the Westerlund 1 superbubble. The extended emission may be explained either by mismodeled gas in the interstellar background model or by the local injection of particles. Under the latter hypothesis, explaining the component associated with Kes 41 requires converting $\lesssim 5\%$ of the supernova remnant energy into accelerated particles, while accounting for the two components near Westerlund 1 requires converting $10^{-4}$ of the cluster wind mechanical power into gamma rays. Nevertheless, our results strengthen the evidence for gamma-ray emission structures arising at intermediate spatial scales between isolated objects and the large-scale diffuse emission from the interstellar medium. This could explain a part of the soft unassociated Galactic sources detected by the Fermi LAT. (Abridged)
We investigate whether cosmogenic gamma rays and neutrinos from ultra-high-energy cosmic rays (UHECRs) produced in hadronically active blazars can significantly contribute to the observed spectral energy distribution (SED), and examine the potential of the Cherenkov Telescope Array Observatory (CTAO) to detect and discriminate this component. For four nearest neutrino bright blazars modeled by Rodrigues (2024), specifically AP Librae, TXS 1700+685, PKS 2326$-$502, and PKS 2345$-$16, we constrain the UHECR luminosity by requiring that the total emission does not overshoot available multi-wavelength data, assuming acceleration in the optically thin large scale jet, propagate UHECRs through realistic Galactic and extragalactic magnetic fields, and forecast CTAO performance. We find that cosmogenic gamma rays contribute up to about 7.5% of the total gamma-ray emission at the highest energies ($10^{20}$ eV), and thus remain subdominant in the SED. In contrast, the associated neutrino emission, at approximately $10^{44}$ to $10^{46}$ erg s$^{-1}$, offers a robust propagation signature of hadronic acceleration. Our simulations of 50 hour CTAO exposures show that the observatory can accurately reconstruct the intrinsic gamma ray spectrum, enabling strong constraints on the hadronic origin and propagation of UHECRs. These results demonstrate the power of combining CTAO gamma-ray measurements with neutrino observations to probe blazar hadronic processes.
We present the results of the first systematic search for refractive scattering substructure in active galactic nuclei (AGN) using ground-based very long baseline interferometry (VLBI) data. We selected 1074 AGNs showing a steep frequency dependence of their angular sizes and compared the observed long-baseline visibility amplitudes with the predictions of the refractive substructure model for scatter-broadened images. Nine sources showed consistency with the expected level of refractive noise, including the quasar TXS 2005+403 (J2007+4029), which confirms previously reported results for this source. Remarkably, the majority of these AGNs are located in the Cygnus region -- one of the key ``hotspots'' of strong Galactic scattering. For five selected AGNs, we identified at least two epochs with sufficiently high data quality showing long-baseline detections, suggesting the presence of scattering substructure. We showed that refractive scattering substructure in AGN is a rare effect, requiring high sensitivity, sufficient source brightness, and lines of sight with particularly strong turbulence in the interstellar medium. Nevertheless, AGNs can serve as powerful probes of the turbulent interstellar medium. We discuss strategies for improving the sensitivity of VLBI observations and identify the most promising sources for reliable detection and investigation of refractive substructure in future observations.
Motion of plasma in pulsar magnetospheres essentially follows bent magnetic field lines which may result in radiation beams sweeping past the observer's sightline. We show that under certain conditions such a passing beamlet can produce two orthogonal polarization modes (OPMs) solely from averaging of radiative contributions, without involving birefringence. This requires that radiation from different stages of the passage is contributed incoherently, hence that the timescale for the coherence is much shorter than the typical sampling time. The model produces two OPMs of comparable amount, with a wide range of (single-pulse) polarized fractions, which is consistent with general observed pulsar properties. Under the specific conditions required (essentially incoherent summation of signal at appropriate time scale) the averaging effects alone can thus lead to the appearance of two orthogonal polarization tracks with OPM jumps.
Bow-shock pulsar wind nebulae (PWNe) are synchrotron sources formed when the outflow of supersonic pulsars is confined by the surrounding interstellar medium. These sources are dominated by freshly injected particles, thus providing a unique laboratory for studying particle acceleration in relativistic outflows. The Mouse is a prototypical bow-shock PWN and is bright in the radio and X-ray bands, enabling detailed multi-wavelength spectral modeling. Using data from 17 telescopes, the Mouse is detected from 118 MHz to 353 GHz, covering over 3 decades in frequency, and in X-rays from 0.5 to 30 keV. A clear spectral break is identified at 3.7 GHz. The radio spectrum exhibits a rising trend at low frequency, peaks at the break, and then declines. We find a change in spectral index $Δα=0.65\pm0.05$ across the break. We attribute this to an intrinsic cutoff or steepening in the injected particle distribution. This could indicate the characteristic energy of leptons leaving the pulsar magnetosphere, or energy dissipation after the particles cross the termination shock.
Long-term pulse-profile variability has emerged as a potentially important source of uncertainty for high-precision pulsar timing, as changes in the integrated pulse shape can bias pulse times of arrival and affect the sensitivity of pulsar timing arrays to nanohertz gravitational waves. We present a systematic search for long-term pulse-profile variability in the six-year MeerKAT Pulsar Timing Array data, analysing observations of 84 millisecond pulsars using a two-dimensional Gaussian-process framework to identify coherent temporal and phase-dependent profile evolution. We detect pulse-profile variability in 18 pulsars. Using the frequency dependence of the observed profile evolution together with scattering simulations, we classify six pulsars as exhibiting chromatic variability consistent with interstellar scattering, while the remaining twelve display predominantly frequency-independent behaviour suggestive of intrinsic changes in the pulsar emission process. Correcting for the S/N-dependent detection bias, we infer that intrinsic profile variability is substantially more common than implied by the detected sample, with at least 40 per cent of millisecond pulsars exhibiting intrinsic profile variability at the 95 per cent credible level. Motivated by the association between pulse-profile variability and spin-down-rate ($\dotν$) switching in canonical pulsars, we search for corresponding $\dotν$ variations but find no statistically significant detections. Simulations indicate that the expected variations remain below the sensitivity of current datasets, although unresolved changes may contribute to weak achromatic timing noise. These results suggest that long-term pulse-profile variability may be common among millisecond pulsars and should be explicitly modelled in future high-precision pulsar-timing experiments.
Star-forming regions (SFRs), including HII regions, have recently attracted significant attention as potential sites of particle acceleration capable of producing PeV cosmic rays, i.e. as potential Galactic PeVatrons. Yet, SFRs known as candidate PeVatrons are still very rare. Here we report a highly likely physical association between the UHE gamma-ray source 1LHAASO J1852+0050u and coincident star-forming regions. Using $\sim$4 years of LHAASO data, the extended source 1LHAASO J1852+0050u (called J1852 in this study) was reanalyzed with a significance of $\sim$13$σ$ and 39%-containment radius of 0.35$°$. This extended source projectively covers dozens of HII regions and an energetic pulsar PSR J1853+0056. We reanalyzed $\sim$16 years of Fermi-LAT data and found a GeV source (named SrcAG) coincident with the LHAASO source J1852. The CO-line observation shows a molecular clump which coincides with the gamma-ray sources and is related with the ultra-compact HII region G34.26+0.15 at the distance of $\sim$3.3 kpc. Fitting to the spectral energy distribution indicates that PSR J1853+0056 cannot afford the observed gamma-ray flux either for J1852 or SrcAG. Although the contribution from hidden pulsars to the gamma-ray emission could not be excluded, it is most likely that the GeV-100TeV gamma-ray emission is dominantly powered by SFR through p-p hadronic interaction between the protons accelerated in the SFR and MCs: the GeV emission is ascribed to the protons accelerated by protostars' activities, while the TeV emission ascribed to the protons accelerated by massive stars in the SFR.
The origin of the optical emission in tidal disruption events (TDEs) remains one of the major outstanding questions in the field, in part due to the limited number of nearby events with high-cadence monitoring to track their evolving photometric and spectroscopic properties. We present multi-wavelength observations of the nearby ($z=0.01368$) TDE\,2025aarm, including near-daily spectroscopic coverage prior to the optical peak. Its proximity makes it one of the brightest TDEs discovered, reaching a peak magnitude of $m_r\sim15.5$ ($M_r\sim-18$). The light curve deviates from a smooth evolution, exhibiting multiple rebrightening episodes visible in both the individual filter light curves and the bolometric luminosity. Blackbody modelling reveals that these rebrightenings are associated with an increase in temperature of $> 5,000-10,000$\,K, while the inferred photospheric radius remains approximately constant. Simultaneously, the H$α$ line not only increases in blueshift but also broadens, suggesting a link between the continuum rebrightenings to changes in the kinematics of the line-forming gas. We identify a persistent absorption component at $\sim-3900$\,km\,s$^{-1}$ in multiple Balmer lines, providing further evidence for outflowing material. The H$α$ profile also exhibits excess flux compared to a Gaussian on both sides of the line, inconsistent with simple scattering-dominated outflow models. Disk-profile modelling provides evidence for the emergence of a disk-like component least $\sim20$ days after peak, with substantial changes in the disk properties between $\sim50$ and 60 days. These observations highlight the complexity of TDE emission processes and demonstrate how dense multi-wavelength monitoring can disentangle the roles of accretion, reprocessing, and outflows in shaping TDE emission.
We present a temporal and time-resolved spectral analysis of GRB 250920C using \textit{Fermi}/GBM and \textit{Swift}/BAT data. The prompt emission consists of two distinct episodes, EI and EII, separated by a significant quiescent interval. We perform Bayesian spectral fitting with empirical, thermal, composite, and physical synchrotron models. In EI, the spectra show clear evidence for an additional thermal component. The BB+PL model is preferred in most bright time bins, and joint \textit{Fermi}/GBM+\textit{Swift}/BAT fits further support the presence of this component. The blackbody temperature generally decreases with time, the blackbody flux follows the pulse profile, and the thermal flux fraction remains high. Fireball-parameter estimates give a photospheric Lorentz factor of a few hundred and an initial radius of $r_{0}\sim10^{9}$--$10^{10}$ cm, with $1+σ_{0}$ of order unity and $η\gg1$, supporting a thermally dominated baryonic outflow. In contrast, EII is dominated by non-thermal emission. Its low-energy photon indices do not significantly exceed the synchrotron line of death, and the spectra are well described by fast-cooling synchrotron radiation in a decaying magnetic field. The magnetization constraint gives $σ_{\min}\sim1.1$--$4.3$. Since these values exceed unity, they suggest that EII may be Poynting-flux dominated. These results therefore suggest a possible transition in GRB 250920C from a fireball-dominated EI phase to a Poynting-flux-dominated EII phase.
With the development of large-scale photometric surveys, an increasing number of supernova candidates are being discovered, leading to a rapidly growing demand for supernova spectra. In addition to equipping photometric surveys with follow-up spectroscopic facilities, archival spectra from large multi-object spectroscopic surveys can be mined to provide spectroscopic classifications for candidates and to find supernovae missed by previous surveys. In this work, we combine Principal Component Analysis (PCA), the Local Outlier Factor (LOF) algorithm, and the supernova classification tool SNID to search for Type Ia supernovae among 1,757,303 galaxy spectra from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1). We finally obtain 247 Type Ia supernovae and 17 supernovae of other types. Among these, 202 supernovae lack classification records in the Transient Name Server (TNS) and represent newly identified SNe. These results demonstrate the potential of multi-object spectroscopic surveys to supplement supernova samples, particularly for transients missed by traditional photometric surveys.
The present thesis is a work in progress on improving the simulation of integrated pulse profiles of pulsars. Pulsars, highly magnetized rotating neutron stars, serve as precise cosmic clocks useful for studying gravity and the interstellar medium. Although periodic pulses from many pulsars are observed and modeled at differ- ent radio frequencies, a robust realistic simulation of their profiles utilizing their parameters remain underdeveloped. The present study focuses on developing a physically consistent model to reproduce observed pulse shapes across multiple frequencies. This study will be useful in coherently understanding the various as- pects of the pulsar emission mechanism. These integrated pulse profiles are the time-averaged properties of pulsar emission and are obtained after averaging indi- vidual pulses over a few thousand rotations. Observed to be highly stable, these are indicative of the global properties of the pulsar magnetosphere. By refining simulation techniques and incorporating the propagation effects of the interstellar medium, this work will prove to be helpful in generating synthetic, yet realistic models. The results are expected to contribute to the interpretation of observa- tional data and to the development of improved timing and emission models for pulsars
We present a multiwavelength radio study of the changing-look AGN 1ES 1927+654 following its renewed X-ray brightening since mid-2022, combining VLBA imaging with multi-band and time-resolved VLA observations. Our main results are: (i) VLBA observations continue to reveal the bipolar jet structure first reported in Meyer et al. (2025), with new X-band detections confirming the morphology at larger separation and K-band imaging revealing a bridge of emission consistent with a continuous outflow; (ii) the VLBA core exhibits a GHz-peaked spectrum, and VLA/VLITE observations from 340 MHz to 45 GHz show emission from an inhomogeneous synchrotron source, self-absorbed below ~3 GHz and optically thin above; (iii) VLA and VLBA flux densities are consistent above 8 GHz, indicating emission dominated by a compact core; (iv) linear polarization is detected above ~3 GHz, with fractional polarization increasing with frequency from ~0.2% to ~6% and EVPAs exhibiting non-linear lambda^2 behavior; (v) long-term VLA X-band monitoring reveals a large, smooth rotation of the EVPA of ~137 degrees over ~81 days in 2026, followed by a partial return, at stable total intensity and suppressed fractional polarization <~1.4%, consistent with a propagating disturbance crossing the compact emitting region; and (vi) circular polarization is not detected at any epoch. These observations demonstrate that 1ES 1927+654 has transitioned from an X-ray-dominated, non-jetted state to a radio-loud AGN (Laha et al. 2025) hosting a newly launched relativistic jet, a unique laboratory for studying jet birth and early magnetic field evolution.
Terminal mass-loss in massive stars is commonly reconstructed from observations of circumstellar interaction in core-collapse supernovae (CCSNe); however, such signatures can be easily missed depending on the phase when they appear. In this paper, we present a multi-wavelength study of SN\,2022yyz, a nearby CCSN selected by its luminous late-time mid-infrared (MIR) excess in NEOWISE at $> 500$ days after the discovery. Classified as a Type II SN at discovery, it reached a peak bolometric luminosity of $1.3^{+0.7}_{-0.4}\times10^{43}$ erg s$^{-1}$ before fading away with a light-curve drop at $\sim 300$ days. We present late-time HST/WFC3 ultraviolet (UV) imaging and Keck/LRIS spectroscopy at $\sim 1200$ days, which reveal clear signs of circumstellar-medium (CSM) interaction via a UV excess consistent with a strong MgII contribution and a flat-topped H$α$ profile. We show that the bolometric light curve and the late-time UV luminosity can be explained using a semi-analytic CSM interaction model involving a broken power-law density profile: with $\sim 2.7~M_\odot$ of nearby, dense CSM within a radius of $\sim 10^{16}$ cm surrounded by wind-like outer CSM with a mass-loss rate of $\sim 10^{-4} {\rm ~M_\odot~yr}^{-1}$. The broken density profile indicates $\gtrsim 100\times$ mass-loss rate enhancement starting a few hundred years before explosion, pointing to a possible common-envelope ejection in a binary system. Despite being overlooked for follow-up at such proximity, our results demonstrate that IR photometric selection methods provide a powerful way of mapping terminal mass-loss in massive stars - laying the groundwork for the Roman Space Telescope surveys.
We consider main-sequence stars orbiting close to a supermassive black hole (SMBH), which are potential gravitational wave (GW) sources for LISA if their orbital periods are of the order of an hour. At such a short orbital period, mass transfer from the star to the SMBH occurs. The evolution of the semi-major axis and GW frequency depends on the evolution of the stellar mass and radius. We use MESA to study stars that transfer mass to SMBHs similar to Sagittarius A* starting on the zero-age main sequence (MS). We identify 4 evolutionary phases. (I) Stars initially with mass $>2 M_\odot$ remain on the MS as their mass and radius decrease. (II) Below $2 M_\odot$, the separation is sufficiently small so that the GW timescale is too short for the stars to maintain thermal equilibrium. They evolve adiabatically off the MS and shrink rapidly as they lose their high-entropy envelope. (III) Below $0.5-1 M_\odot$, depending on the initial mass, the uniform low-entropy core is exposed and stars expand adiabatically. (IV) Below $0.15-0.4 M_\odot$, the thermal and GW timescales become comparable, and stars cool and shrink while maintaining this balance. Overall, mass transfer causes the orbit to expand, slowing down the orbital evolution and leading to GW emission at lower frequencies making the GW signal harder to detect. If located around Sagittarius A*, mass transferring stars spend most of the time relatively close to the LISA sensitivity curve, with maximal SNR reaching around 600 during the transition between stages II and III.
The generation of synthetic observations is required to compare simulations of astrophysical phenomena with real sources. We present a new open-source toolkit cuDART which allows meshed simulation data to be rapidly rendered into synthetic observations of optically thin emission, automatically accounting for a range of relativistic and geometric effects such as Doppler boosting. Unlike previous visualisation schemes which generally consider a single snapshot in time, cuDART is capable of accounting for the finite time delay across astrophysical distances by rendering data from multiple simulation epochs simultaneously. We demonstrate how failing to account for a finite light time delay can result in synthetic observations with incorrect observed motion, morphologies and fluxes. The code does not require a static pattern frame for integration, relaxing standard assumptions required for analytical observational estimates. The rendering processes is accelerated through hardware (using GPUs) and software (using the 3D differential digital analyser algorithm), allowing for rapid generation of observations from a large number of viewing orientations. In this paper we describe the physical motivation, document the code operation and test the output against theoretical relativistic and geometric expectations. cuDART is publicly available via GitHub and documented using ReadTheDocs.
Circumbinary disks around supermassive black hole binaries (SMBHBs) are expected to form across a broad range of inclinations, with retrograde configurations potentially a common occurrence. Here we present the first grid-based hydrodynamical simulations of retrograde circumbinary disks around eccentric SMBHBs, solving an energy equation that balances viscous and shock heating against blackbody radiative cooling. We investigate different initial disk Mach numbers $\mathcal{M}_a \in \{10, 20, 40\}$ and consider binary eccentricities $e_\mathrm{b}\in[0.0,\,0.8]$, finding that multiple stable states exist for the same binary eccentricity and Mach number. These states differ by the sense of rotation of their minidisks; both retrograde ($\downarrow\downarrow$), both prograde ($\uparrow\uparrow$), or one of each ($\uparrow\downarrow$), a property set by the initial conditions. Our findings indicate that each state produces qualitatively distinct orbital evolution: $\downarrow\downarrow$ and $\uparrow\downarrow$ minidisks drive circular inspirals whereas $\uparrow\uparrow$ minidisks drive eccentric inspirals with potentially observable eccentricity in the LISA band. We measure the electromagnetic emission produced by the disk, finding that a binary of mass $M_\mathrm{b}=8\times10^6\mathrm{M}_\odot$ at redshift $z=1$ would be detectable by both current and upcoming optical and UV instruments. We demonstrate that the time- and azimuthally-averaged disk profiles are well described by 1D models, which naturally set a cavity radius within which angular momentum transport is dominated by Reynolds stresses rather than viscosity.
Neutron star-bearing compact-object mergers can create a kilonova, a transient powered by the radioactive decay of material synthesized via rapid neutron capture ($r$-process). During the merger, a fraction of the ejecta is launched at high speeds ($\gtrsim$ 0.4c) that can lead to neutrons evading capture onto seed nuclei, resulting in free neutrons. The free neutrons then decay and inject additional energy that power early-time ($\lesssim$1 day) emission. Here, using \texttt{Sedona}, we present a grid of the first non-local thermodynamic equilibrium and frequency-dependent opacity radiative transfer simulations to predict early X-ray/UV/optical emission from free neutron decay that span the expected theoretical range of free neutron mass, mixing with $r$-process material, and extent of the high-velocity ejecta tail. The emission can be characterized by an SED that rapidly shifts from X-rays of $\sim\rm{few} \times10^{41}-10^{42}$ erg s$^{-1}$ in the first $\sim$ minutes to a far-UV and near-UV peak on the timescale of $\sim$ minutes to hours, followed by enhanced optical and IR emission for sufficiently large free neutron masses. The properties of the free neutron ejecta are most distinguishable, in principle, at extreme-UV and far-UV wavelengths, with SED peak flux and wavelength determined by the maximum velocity of the ejecta and the mixing. In the far-UV and near-UV, free neutron ejecta masses as small as 10$^{-7}$ M$_\odot$ exhibit a unique bump compared to neutron-free models, demonstrating the importance of upcoming UV missions like UVEX and ULTRASAT to constraining the nucleosynthetic environment of neutron star mergers.
Magnetars became favored engines of fast radio bursts (FRBs) following the discovery of a luminous radio burst coincident with a hard X-ray flare from a Galactic magnetar. Several repeating FRB sources exhibit time-variable rotation measures and compact, spatially coincident persistent synchrotron radio emission, consistent with energetic-particle nebulae and confined by young supernova ejecta. Secular changes in the dispersion measure (DM) of repeating FRBs have also been observed, offering a complementary probe of their local environments; for example, the DM of FRB 121102 rose until 2019 before declining in recent years. Although rising DM evolution has been attributed to shock ionization, shocked ejecta can cool efficiently through metal-line emission and may recombine, especially if it mixes turbulently with cooler gas. Here we argue that DM variations of the observed magnitude and timescale instead arise from changes in the ionization state of supernova ejecta irradiated by X-rays from a time-variable central engine. The dominant rapidly variable contribution comes from the dense, weakly ionized shell swept up by the expanding nebula, whose ionization and recombination times are much shorter than those of the more extended ejecta. The model predicts that enhanced FRB activity should be accompanied by rising DM, with a response smoothed over the ionization/recombination time, superposed on a slower secular decline from the expanding ejecta; in FRB 121102, the DM maximum occurred close to the burst-rich 2018--2019 activity episodes, including the 2019 FAST burst storm. If the recently reported renewed activity of FRB~121102 is sustained, the model predicts that its DM decline should flatten and may reverse into a fresh rise over the coming years.
Neutron stars (NSs) provide a unique laboratory for testing gravity in the strong-field regime and for searching for deviations from General Relativity (GR). In this work, we investigate the effects of Energy-Momentum Squared Gravity (EMSG) on NS structure and examine whether its signatures can be identified from observable stellar properties using supervised machine learning (ML). We solve the modified Tolman-Oppenheimer-Volkoff equations for approximately $10^{4}$ nuclear equations of state (EOSs) for EMSG coupling parameters $α=\{-5.01,-2.50,0,+2.50,+5.01\}\times10^{-38}\,\mathrm{erg}^{-1}\mathrm{cm}^{3}$, and calculate the gravitational mass $M$, radius $R$, dimensionless tidal deformability $Λ$, and fundamental $f$-mode oscillation frequency for each stellar configuration. Imposing observational constraints on $M$, $R$, and $Λ$, we split our datasets into train and test sets, and we employ Random Forest (RF), K-Nearest Neighbors (KNN), Support Vector Machine (SVM), Logistic Regression (LR), and Gaussian Naive Bayes (GNB) to classify the representative sectors $α=\{-5.01,0,+5.01\}\times10^{-38}\,\mathrm{erg}^{-1}\mathrm{cm}^{3}$ using $(M, R,Λ,f)$. The RF classifier performs best, achieving an accuracy of approximately $99.85\%$ with precision, recall, and F1-scores exceeding $99.8\%$, while KNN also achieves accuracy above $99\%$. The nearly diagonal confusion matrices demonstrate that the observationally viable NS configurations associated with different EMSG sectors remain highly separable in the multidimensional observable space. Our results show that NS observables retain robust signatures of EMSG even after observational filtering, establishing ML-assisted NS observations as a promising complementary approach for probing modified gravity with current and future multi-messenger observations.
Rapid, precise localization of gamma-ray bursts (GRBs) requires continuous wide-field imaging at hard X-ray and soft gamma-ray energies. We develop a fixed (non-rotating) modulation-collimator architecture to address this challenge. A bigrid collimator, a pair of aligned, periodic slit grids in front of a simple, non-imaging photon-counting detector, transmits an off-axis source through a triangular fringe pattern whose local angular position, the source phase, is recovered from simultaneous count rates in several grid pairs at fixed relative phase offsets, rather than by mechanically rotating a single grid pair as in the classical rotating modulation collimator. Four-phase demodulation of these count rates robustly cancels detector background and grid leakage without requiring rotation. Because the fringe pattern repeats across the sky, a phase measurement corresponds to many candidate source positions; we resolve this ambiguity over an arbitrarily wide field of view using a cascade of progressively coarser-pitched grid stages, generalizing the two-grid vernier-scale technique originally used to localize Sco X-1. This vernier cascade technique reduces the required photon counts and grid-calibration accuracy required to achieve fine angular resolution over a broad field of view. Generalized to two dimensions with orthogonal cascades, the resulting architecture achieves wide-field localization with no moving parts, no time modulation, and only simple photon-counting detectors, a promising candidate for simple, inexpensive, continuous-monitoring instruments supporting time-domain and multi-messenger astrophysics.
Phenomenological deviation parameters used in gravitational-wave tests of general relativity need not represent universal fundamental properties of the sources. Instead, a deviation parameter can be a derived quantity that depends on source properties. Apparent violations of general relativity arising from waveform systematics or data-quality issues can also depend on source properties. In addition to diagnosing individual events, it is helpful to identify systematic trends at the population level. We introduce a nonparametric framework that uses mutual information to test whether the deviations depend on source properties. We validate the framework using simulations with massive graviton waveforms. Assuming we do not have prior knowledge of the underlying theory, massive graviton effects can be captured by a parameterized post-Einstein coefficient that is expected to depend on the luminosity distance and source mass. The results show that the mass and distance dependencies become clear for graviton masses of $\geq 1\times10^{-22}\,\mathrm{eV}$ at detector sensitivities representative of the current fourth observing run, while the controlled general relativity simulations yield null results. We then analyze public real-data products for tests of general relativity from GWTC-4.0 and 4-OGC, including phenomenological parameterized tests, modified dispersion tests, and parity violation tests. The phenomenological deviation parameters have tentative correlations with effective spin, and correlations of modified-dispersion or parity-violation parameters with source mass are visible. These results demonstrate that catalog-level dependence can identify patterns and thus be used to diagnose the origin of deviations in testing general relativity. The framework therefore complements single-event analyses and parameterized hierarchical population tests.
Neutron star-black hole (NSBH) mergers are promising multimessenger sources, but only a subset is expected to produce detectable electromagnetic (EM) counterparts, depending on the binary mass ratio, black hole (BH) spin, and neutron star (NS) equation of state (EoS). We investigate the prospects for detecting kilonova counterparts to NSBH mergers in the LIGO-India era using the Vera C. Rubin Observatory (Rubin). We simulate GW230529-motivated NSBH populations with two effective tidal deformability ranges, estimate gravitational-wave detections for the LHV (LIGO-Livingston, LIGO-Hanford, Virgo) and LHVA (LHV with LIGO-India at Aundha) detector networks, and model the detectability of the associated kilonova emission with Rubin. We find that adding LIGO-India approximately doubles the number of EM follow-ups, primarily by increasing the duty cycle and reducing the median sky-localization area, and notably increases the number of such detections at large distances. For fixed-exposure strategies, the EM detection rate saturates at exposure times of a few hundred seconds, reflecting the trade-off between depth and sky coverage. We therefore introduce an event-specific exposure-time optimization strategy based on the expected counterpart brightness and sky-localization area. With this strategy, the expected NSBH GW+EM detection rate per year is 4.19(-3.58,+9.66) for the optimistic and 0.79(-0.68,+1.83) for the conservative population with the LHVA network, compared to 1.46(-1.25,+3.37) and 0.27(-0.23,+0.62), respectively, for LHV. The required Rubin follow-up time remains within the expected Target-of-Opportunity allocation. These results show that in the next decade, a few such NSBH follow-ups may take place, with a fraction of them in the redshift range of about 0.1-0.2, which will not only be useful for measuring the Hubble constant but can also help us in probing the Hubble parameter.