2026-05-12 | | Total: 32
On November 12th, 2025 the LIGO--Virgo--KAGRA (LVK) collaboration reported gravitational waves (GWs) from a compact object merger candidate (S251112cm) with at least one sub-solar mass component. Using the Dark Energy Camera (DECam), the Fraunhofer Telescope at Wendelstein Observatory (FTW), and the Zwicky Transient Facility (ZTF), we surveyed $56\%$ of the GW localization region beginning $2.4$~hours after the GW alert. We find no kilonova (KN) counterpart, and use radiative-transfer models to rule out $42\%$ (ZTF), $68\%$ (DECam), and $92\%$ (FTW) of the KN models as possible emission from this GW candidate. Within the recently proposed disk-fragmentation (``superkilonova'') model for generating sub-solar mass neutron star mergers from stellar core-collapse, the delay between the supernova explosion time and the GW merger time is estimated to be less than a few days. Searching this time window prior to the GW event, we identify and spectroscopically classify a IIb supernova (SN~2025adtq), with a spatial association odds ratio of $\log_{10}\mathcal{I} \approx 4.8$, a chance coincidence probability of ${\sim}2$--$9\%$, and an estimated explosion time ${\sim}2$ days prior to S251112cm. SN~2025adtq is the second Type~IIb supernova found in spatial and temporal coincidence with a sub-solar mass GW candidate, following the previously reported S250818k/SN~2025ulz association; jointly, we measure an odds ratio that favors the association hypothesis over the null, however, when conditioned on finding a coincident supernova by chance, the odds ratio disfavors association. Together, these results provide suggestive but inconclusive evidence for the superkilonova formation channel.
Distance measurements are fundamental to radio pulsars' use as astrophysical probes of General Relativity and the interstellar medium. One of the primary methods for determining pulsar distances is HI kinematics, which leverages the radial velocities of HI absorption and emission features detected along pulsar lines-of-sight. This method necessarily assumes a model for Galactic rotation, our knowledge of which continues to evolve in both accuracy and precision. In this research note, we derive kinematic distances for 66 pulsars with archival HI radial velocity measurements using a state-of-the-art Galactic rotation curve. The results and software are provided in an online repository. Our kinematic distances differ by $<1σ$ from published parallaxes for nearly all pulsars in the sample that have both types of distance measurement available. Comparison to the NE2025 Galactic electron density model shows general consistency between measured and predicted distances.
At the center of our Galaxy lies the bright γ-ray point-like source HESS J1745-290, which is compatible in position with Sgr A star, although an association between the two remains uncertain. Using data obtained between 2004 and 2019 with the High Energy Stereoscopic System (H.E.S.S.) on the Galactic center region, we studied the variability of HESS J1745-290 over 353 hours of observations collected over 16 years, representing the largest dataset gathered yet on this region at TeV energies. We performed a 3D maximum-likelihood analysis of the central source and the diffuse γ-ray emission in the Galactic center region. This analysis allowed us to extract the spectral and morphological intrinsic behavior of the two components. By performing this analysis on an annual basis, we derived the light curve of HESS J1745-290 and the diffuse emission over the past 16 years. The 3D maximum-likelihood analysis method allowed us to separate the central source from the overlapping diffuse emission, enabling a recalibration of the former by the latter and alleviating some of the systematic effects. We find no long-term or yearly variability. We also provide an estimate of the sensitivity of H.E.S.S. to variation of this specific source over 16 years. We rule out any yearly gamma-ray flux variation of this source larger than 30 percent, as well as any linear flux variation exceeding 30% over this time period.
We summarize a set of constraints on a proposed population of primordial, ultra-massive, electrically charged compact objects, which we call astrons. The analysis combines charge generation, charge saturation, persistence of the charge in an ionized medium, screening by the intergalactic plasma, the Reissner--Nordström geometry of highly charged compact objects, lensing, and the cosmological implications of a sparse charged population. We also discuss the possible relation to the early structures revealed by the James Webb Space Telescope: if astrons are relevant there, they would be primordial dark seeds rather than luminous objects directly observed at high redshift. The resulting scenario is sharply constrained. Ordinary accretion saturation gives charges far below the large-charge phenomenological benchmark, screening is a serious plasma-physics issue, and a large charge can place the exterior geometry deep in the super-extremal regime. As expected at the level of a homogeneous Friedmann--Lemaître--Robertson--Walker (FLRW) description, the interaction energy of a population of charged objects scales as \(a^{-4}\), so the simplest perfect-fluid reduction does not generate asymptotically late-time acceleration; any acceleration era tied to that homogeneous component can only be transitory. The astron scenario should be regarded as a constrained framework whose viability depends on plasma physics and on a cosmological treatment beyond the homogeneous approximation.
We investigate whether reconnection of small-scale current sheets in transrelativistic supermassive black hole (SMBH) coronae can supply the nonthermal protons needed for high-energy neutrino emission, using NGC 1068 as a test case. We model the corona as a strongly turbulent, low-$β$, collisionless hydrogen plasma with characteristic size $r_{\rm co}$, magnetic field strength $B$, proton density $n_p$, and radiation energy density $u_{\rm rad}$. Combining the observed IceCube-band neutrino luminosity with the X-ray luminosity and Thomson optical depth reduces these coronal quantities to a one-parameter family. Across this family, the proton magnetization $σ_p \equiv B^2/(4πn_p m_p c^2)$ is transrelativistic with $σ_p \sim 0.3$. In this regime, we show that repeated encounters with intermittent reconnecting current sheets can energize suprathermal protons up to tens of PeV before photomeson cooling limits further acceleration. These injected particles may then be further processed by stochastic interactions with the turbulent cascade. Motivated by PIC simulations of strong turbulence at comparable magnetization, we adopt a nonthermal proton spectrum with an independently specified index and find that the predicted TeV spectral shape is broadly consistent with NGC~1068 without fitting the proton spectral slope.
Relativistic jets launched from stellar-mass compact objects embedded in the accretion disk of an active galactic nucleus (AGN) can produce nonthermal emission upon successfully breaking out of the disk. In this paper, we present a comprehensive study of the long-term propagation dynamics and broadband nonthermal radiation signatures of such jets in a realistic AGN environment, explicitly modeled as wind outflows. Our modeling reveals two distinct features imprinted by the high-density AGN medium: rapid deceleration of the jet ejecta, accompanied by a prompt downshift of the emission spectral energy distribution, and persistently strong synchrotron self-absorption, giving rise to a prominent quasi-thermal hump in the emission spectrum. Crucially, both gamma-ray burst jets and jets powered by accreting binary black hole merger remnants can produce detectable multi-wavelength emissions that substantially outshine the AGN background. Moreover, the short time delays between gravitational wave triggers and these electromagnetic counterparts--typically less than 106 s--greatly facilitate secure multi-messenger associations. Besides, our findings highlight that interaction-induced radiation from AGN-embedded jet systems offers a powerful diagnostic probe of the spatial distribution,density structure, and physical properties of the AGN medium.
Neutrino flavor evolution in dense astrophysical environments such as core-collapse supernova (CCSN) is influenced by collective effects. While the Fast Flavor Instability (FFI) and the Collisional Flavor Instability (CFI) are recognized as key drivers of rapid flavor conversion, their non-linear competition with collisional damping in spatially varying environments remains poorly understood. Motivated by recent findings that FFI and resonance-like CFI co-occur in the post-bounce phase in CCSN, we scrutinize their dynamic competitions and asymptotic states. To this end, we perform numerical simulations of the quantum kinetic neutrino transport, incorporating both spatial advection and the collision terms. We demonstrate that the interplay between these coexisting neutrino flavor instabilities and collisions leads to rich dynamics. Rather than merely inducing simple decoherence, collisional damping can substantially alter the overall dynamics of collective flavor oscillations, driving the system through complex evolutionary pathways. In all cases where flavor instability develops, we find that the system converges to the same flavor-equilibrated asymptotic state, despite the diversity of intermediate dynamics. Our results suggest that this dynamic competition could alter the widely accepted picture of collisionless FFI, highlighting the need to incorporate realistic collisional effects into studies of flavor conversions in CCSN models.
Pulse profiles probe the emission geometry of accreting X-ray pulsars, but their observed shapes may depend on instrumental response and observational setup. The pulsed fraction spectrum provides a compact spectro-timing observable that can both trace localized spectral features and serve as a quantitative cross-calibration diagnostic. We assess the consistency of energy-resolved pulse profiles obtained with simultaneous XMM-Newton/EPIC-pn and NuSTAR/FPM observations of Vela X-1, and investigate the broadband pulsed fraction spectrum as a diagnostic of spectral features from 1 to 70 keV. We construct energy-phase matrices for both instruments and derive pulsed fraction spectra after carefully accounting for instrumental and observational effects. We quantify the residual systematics in the overlapping 3-10 keV band. We then model the broadband pulsed fraction spectra phenomenologically and search for timing signatures of spectral features. After correcting for instrumental effects, the pulsed fraction spectra derived strictly over the common exposure intervals of the two instruments agree within 5% in their overlapping 3-10 keV range. Remaining discrepancies larger than 5% are confined to the iron-line region and can be attributed to the different energy resolutions of the two instruments. The broadband pulsed fraction spectrum reveals significant localized features corresponding to known emission lines in the soft band and to cyclotron resonant scattering features. Orbital-phase-resolved modeling of the EPIC-pn pulsed fraction spectrum shows that the soft-band features depend strongly on the equivalent absorption column, with emission-line signatures becoming progressively suppressed during highly absorbed intervals. The pulsed fraction spectrum serves both as a quantitative cross-calibration diagnostic and as a powerful spectro-timing diagnostic.
Context: Recent gamma-ray observations of young star clusters revealed that stellar wind termination shocks accelerate particles, with the energy reservoir provided by the mechanical power of massive-star winds. Aims: Our goal is to identify promising targets for future gamma-ray studies of stellar clusters and associations powered by massive stars. As the wind power of a single Wolf-Rayet (WR) star can rival the cumulative wind power of the most massive clusters, we also investigate isolated WR stars, many of which are indeed isolated. Methods: We ranked a large sample of stellar clusters and associations according to the number of member WR stars divided by the distance squared, a quantity proportional to the expected gamma-ray signal, and searched for spatial correlations with known gamma-ray sources. We repeated the same procedure for individual WR stars with known wind mechanical powers and distances. Results: We found a hint ($\lesssim 3 σ$ confidence) for a correlation between WR-hosting clusters and unidentified GeV gamma-ray sources, and identified new spatial associations for 11 clusters. We also found spatial coincidences between 4 isolated WR stars (WR110, WR114, WR111, and WR14) and unidentified gamma-ray sources. Although no significant correlation is found for isolated WR stars as a population, these 4 objects exhibit particularly large wind-power-to-distance-squared ratios, a necessary condition for detectability with current instruments. Assuming the gamma-ray emission is powered by WR winds, it can be interpreted as arising from interactions between particles accelerated at the wind termination shock and ambient matter or radiation fields. Conclusions: Since the wind power of an individual WR star can rival that of an entire stellar cluster, we provide a ranking of stellar clusters and isolated WR stars that may constitute potential gamma-ray emitters.
We study test-particle dynamics in the equatorial plane of a Kerr-Bertotti-Robinson black hole (BH) immersed in a cloud of strings (CS), with mass M , rotation a, magnetic parameter B, and string parameter α. Using the Hamilton formalism we recover the effective potential Ueff and the conditions for circular motion, and we compute the specific energy E and specific angular momentum L together with the radial, vertical, and azimuthal epicyclic frequencies νr , νθ , νφ. Going beyond the analytic setup, we provide the first numerical mapping of the innermost stable circular orbit (ISCO) for this background and tabulate rISCO, EISCO, LISCO, and the accretion efficiency η = 1 - EISCO for both co- and counter-rotating motion across a wide (a, B, α) grid. The CS parameter pushes the ISCO outward and raises η from 0.057 in Schwarzschild to above 0.25 for α = 0.30 at a = 0.9. We then connect the model with observed twin-peak high-frequency quasi-periodic oscillations (QPOs) in three microquasars (GRO J1655-40, XTE J1550-564, GRS 1915+105) using the relativistic-precession (RP) model and find \{chi}^2-minimum fits with α < 0.13. A general-relativistic hydrodynamical (GRH) study of Bondi-Hoyle-Lyttleton (BHL) accretion completes the picture: the CS contribution sustains shock-cone instabilities, redistributes power-spectral-density (PSD) peaks, and produces low-frequency QPO-like components that distinguish KBR+CS from pure Kerr or KBR.
We present a comprehensive analysis of Vela X-1 using two new NuSTAR observations, placed in the context of four earlier datasets obtained between 2012 and 2020. The energy-resolved pulse profiles demonstrate a significant transformation from an asymmetric structure at low energies to distinct double peaks above 12 keV, whereas the pulse fraction escalates with photon energy but decreases with flux. Broadband spectra validate the Fe K alpha emission line and disclose both fundamental and harmonic cyclotron resonant scattering characteristics (CRSF). We observe no substantial link between CRSF energies and luminosity, contrary to previous findings; rather, the photon index and folding energy demonstrate distinct anti-correlations with flux, aligning with sub-critical accretion and increased Comptonization in the accretion column. Our results provide the first clear evidence that the harmonic CRSF in Vela X-1 does not follow the long-term decay previously claimed. The fundamental line energy also displays an irregular evolution, without a clear monotonic trend. Notably, the harmonic-to-fundamental energy ratio departs from the canonical value of two, suggesting that the line-forming regions are located at different heights within the accretion column. These results provide new constraints on the accretion geometry and magnetic field topology of Vela X-1, highlighting the importance of continued monitoring with current and future X-ray observatories.
Wind is a key mechanism for supermassive black hole (SMBH) feedback to their host galaxies. In tidal disruption events (TDEs), black holes spend most of their time accreting at highly sub-Eddington rates, implying that feedback from persistent sub-Eddington winds could be significant. We investigate the effects of black hole mass, viscosity parameter and stellar debris temperature on the properties of winds from hot accretion flows in TDEs. We find that more massive black holes yield a higher accreted fraction and launch faster winds, while the debris temperature has a negligible influence on the accretion flow. For $α=0.1$, the mildly-relativistic unbound winds ($\sim 0.1c$) are launched predominantly from the outside of the accretion flows along the equatorial plane, with a kinetic energy of $\sim10^{-4}L_\mathrm{Edd}$. In contrast, convective bound outflows dominate for $α=0.01$, which differs from the true winds typically seen in active galactic nuclei and X-ray binaries. Potential applications for explaining delayed radio brightening in TDEs at $\sim10^3$ days and for searching for intermediate-mass black holes through radio and X-ray surveys are also discussed.
We investigate the thermal evolution of isolated neutron stars containing hyperon--mixed kaon--condensed matter, focusing on the role of proton superconductivity. The equation of state utilized for cooling calculation is based upon the minimal relativistic mean--field framework supplemented by chiral SU(3) dynamics for kaon condensation with an additional component on the three-baryon force, which ensures stiffness at high densities enough to meet astrophysical constraints on neutron-star masses and radii. We show that the nucleonic direct Urca processes operate at relatively low stellar masses ($M \gtrsim 1.3\,M_\odot$), erasing any observable signature of strangeness in the absence of superfluidity. However, if the proton $^1{\rm S}_0$ superconductivity works, because of suppression of fast neutrino cooling processes, the cooling scenario could become relevant with the strangeness, depending on the density regions of the pairing gap. In particular, if the proton superconductivity is so strong in high-density regions ($T_{c,p}\sim10^{10}~{\rm K}$), the nucleon and hyperon direct Urca processes shut down, which makes the kaon-induced Urca processes dominant in massive neutron stars. This scenario is in good agreement with several cold isolated neutron stars identified recently. Hence, we suggest that strong proton superconductivity can render kaon condensation observationally visible through cold neutron-star observations, providing a potential signature of strangeness in dense matter.
Neutron star (NS) low-mass X-ray binary (LMXB) systems depict a variety of X-ray spectral and timing features, which can be useful to probe the accretion-ejection mechanism in the strong gravity regime. Here, we study the relatively unexplored and faint NS LMXB 1A 1246-588, which is also an ultra-compact X-ray binary (UCXB) with a white dwarf donor. We investigate its temporal and spectral behavior using pointed NICER and AstroSat observations, supported by long-term MAXI/GSC monitoring. The MAXI light curve shows modest, recurrent outburst-like enhancements, providing the long-term flux context for interpreting the pointed observations. During the AstroSat observations in 2017, the source exhibits an absorbed 0.4-20 keV flux of $(1.18 \pm 0.02)$ x $10^{-10}$ $erg$ $cm^{-2}$ $s^{-1}$, while during the NICER observations in 2019, it spans an absorbed 0.5-10 keV flux range of $(0.7-3.7)$ x $10^{-10}$ $erg$ $cm^{-2}$ $s^{-1}$ and traces an atoll-like pattern in the hardness-intensity diagram. Broadband spectral modeling shows that the emission is well described by a soft blackbody and a hard Comptonized component, with no statistically required multicolor disk contribution. The blackbody temperature increases from 0.28 to 0.39 keV, with an emitting radius consistent within 6.9-13.8 km, while the Comptonization photon index varies from 1.8 to 2.3. We find that the observed spectral-state evolution is driven by a redistribution of accretion power between thermal emission from the NS boundary layer and Comptonized emission, consistent with atoll-type behavior. These results provide the first quantitative, multi-epoch view of accretion-state evolution in 1A 1246-588, revealing systematic changes in the thermal boundary-layer emission and the Comptonizing region in this UCXB system.
We investigate the optical appearance of a Kerr-Bertotti-Robinson (Kerr-BR) black hole illuminated by a geometrically and optically thin accretion disk. Instead of using a phenomenological power-law emissivity, we adopt a magnetically driven synchrotron emissivity proxy coupled to the local electromagnetic environment. With a backward ray-tracing framework, we examine the effects of the spin $a$, magnetic parameter $B$, and observer inclination $θ_O$ on the ray-classification maps, redshift distributions, and specific-intensity images. We show that the ISCO position is modified by both $a$ and $B$, and that rapidly rotating prograde configurations can develop an additional model-dependent inner cutoff when the magnetically dominated approximation underlying the emissivity prescription ceases to be applicable. High-resolution one-dimensional intensity profiles further separate the direct image, the $n=1$ lensing-ring contribution, and the higher-order $n\geq 2$ photon-ring subimages, while quantifying the Doppler-induced brightness asymmetry. Retrograde disks exhibit a wider emission-depleted central region because of the outwardly shifted ISCO, making the higher-order lensed components more clearly distinguishable from the direct emission. These results show that the disk inner boundary and the magnetic-field-dependent emissivity can substantially influence the observable appearance of Kerr-BR black holes.
The origin of the highly ionized iron emission (Fe XXV at $6.7\,\mathrm{keV}$) characterizing the Galactic ridge X-ray emission (GRXE) remains a fundamental puzzle in high-energy astrophysics. Although the GRXE continuum is largely resolved into discrete populations of cataclysmic variables and coronally active stars, these sources exhibit Fe XXV equivalent widths significantly lower than that of the total GRXE, leaving the intense iron line emission unexplained. In this work, we cross-correlated the XMM-Newton survey of the inner Galactic disk with Gaia DR3 astrometry to identify and characterize hard X-ray sources ($>2\,\mathrm{keV}$) with reliable stellar counterparts. We selected 107 X-ray sources located within the red giant branch of the color-magnitude diagram, many of which are verified long-period variables. These sources exhibit high X-ray luminosities ($L_{\mathrm{X}} \approx 10^{31}$--$10^{33}\,\mathrm{erg~s^{-1}}$), significantly exceeding the typical coronal saturation levels of single giants. Their X-ray spectra are notably harder than those of quiescent stellar coronae, with plasma temperatures reaching up to $kT \approx 6\,\mathrm{keV}$ and a prominent emission feature at $\sim 6.7\,\mathrm{keV}$. The combination of high $L_{\mathrm{X}}$, hard spectra, and intense Fe XXV emission identifies this population as accretion-powered binaries associated with late-type giants. Our analysis demonstrates that this population contributes $\sim 20\%$ of the total GRXE continuum and $\sim 40\%$ of its iron line emission, providing a key component to resolving the Galactic X-ray background puzzle.
Assuming that the compact source at the Galactic center, Sagittarius A*, is a massive boson star, we fit the near-infrared flare astrometry data. We consider 12 discrete boson star configurations and model the flare as a hotspot on a circular equatorial orbit. The analysis is performed in a Bayesian framework using nested sampling, yielding the marginal posterior distributions of all parameters as well as the Bayesian evidence for each model. For comparison, the same procedure is applied to a Schwarzschild black hole. The resulting Bayesian evidence values differ only marginally between the boson star and black hole cases, and the well-determined mass of Sgr~A* (${\sim}4.296\times 10^6\,M_\odot$) falls within the 68\% highest density interval in every configuration. We conclude that, under current near-infrared astrometric constraints and within the considered parameter ranges, a massive boson star and a Schwarzschild black hole remain statistically indistinguishable as the compact object at the Galactic center.
In this Paper we consider twisted accretion disks in supermassive binary black hole by analytical and numerical means. It is assumed that the disk orbiting around the more massive rotating component and that the disk rings are inclined with respect to the orbital plane. We use orbital parameters of the binary often employed in the precessing massive (PM) model of the well-known blazar OJ 287. Unlike our previous investigation of a similar problem, here we consider disks with both small and relatively large relative thicknesses $δ=h/r$, where $h$ is the disk's height at a typical radius $r$, as well as a range of values of the viscosity parameter, $α$, including the cases when $α\lesssim δ$. Similar to our previous results, we find that the twisted disk relaxes to a quasi-stationary state in the frame precessing with the Lense-Thirring frequency of the orbit. However, its shape is qualitatively different from that corresponding to the case of $δ=10^{-3}$ and $α=0.1$ considered in our previous work. In a disk with $δ=10^{-3}$ but $α\le 2\cdot 10^{-2}$, we find the new effect of generation of a twisting spiral wave near the resonance between a forcing frequency associated with the presence of the secondary and the Lense-Thirring frequency of a particular disk ring defined in the precessing frame. We propose an analytic theory of it, which is in a good agreement with our numerical results. This effect leads to multiple crossings of the orbit with the disk per one orbital period, which contradicts the PM model. When $δ\gtrsim 0.1$, a typical disk's inclination within the orbit of the binary is smaller than that of the orbit which results in only two crossings of the orbit with the disk per one orbital period. We suggest that the additional heating of the disk gas by the secondary-disk collisions may result in $δ\sim 0.1$.
The formation environments of merging binary black holes remain uncertain. While hierarchical assembly in dense stellar clusters has been widely explored as an explanation for black holes exceeding the stellar-mass limit, growth through gas accretion in active galactic nucleus (AGN) disks is an alternative that has received less observational scrutiny. Here we search for an accretion-origin subpopulation using only spin magnitudes, fitting a three-component mixture model to 166 binary black hole mergers from LIGO--Virgo--KAGRA with component shapes fixed from theoretical predictions and only the mixing fractions inferred from the data. We find strong evidence ($ln B = 5.7$) that $\sim 10\%$ (90% credible interval $[1\%, 14\%]$) of detected mergers belong to a subpopulation with primary spins clustered near $a_1 \approx 0.9$, consistent with the theoretical prediction for accretion spin-up. The hierarchical-merger prediction of $a_1 \approx 0.7$ is decisively disfavored as the location of the high-spin subpopulation ($ln B = 5.7$). Post hoc validation reveals that the accretion candidates have systematically higher masses (median $m_1 = 58\,M_\odot$) and aligned spins (median $χ_{\rm eff} = 0.33$, vs. $0.04$ for standard-dominated events). The accretion subpopulation is not limited to systems above the pair-instability mass gap: GW190517 ($m_1 \approx 39 M_\odot$) is among the top candidates, demonstrating that accretion spin-up operates across a range of masses. GW190521, previously interpreted as a hierarchical merger, shows comparable support for an accretion origin. These results provide the first population-level observational evidence for an accretion-origin subpopulation in black hole mergers.
Radiative cooling can strongly influence the structure and dynamics of black hole accretion disks. Here, we perform general relativistic magnetohydrodynamic (GR-MHD) simulations of magnetically arrested disks (MADs) around a non-spinning black hole. Radiative cooling is consistently included in the simulations and its intensity is scaled by the mass accretion rate ranging from $10^{-7}$ to $10^{-4} \dot{M}_{\mathrm{Edd}}$. Considering synchrotron and bremsstrahlung emission, we quantify how radiative losses modify the disk structure and the accretion dynamics. In the inner MAD disk regions, accumulation of magnetic field regulates gas accretion, enforcing the gas into a discrete interchange-driven filamentary structure. We identify, both analytically and numerically, a transition mass accretion rate above which radiative cooling becomes faster than the heating, which is assumed to occur via local coupling to the magnetic field. Above this mass accretion rate, cooling substantially reduces the gas thermal pressure, leading to considerably thinner and denser accretion filaments, and a substantial increase in radiative efficiency, relative to lower accretion rates. We show that under these conditions, conventional measures of the disk scale height become misleading in MAD flows. We therefore introduce an alternative definition based on the polar position of the density maximum, which more robustly characterizes the filamentary structure of the disks in the presence of strong magnetic fields and cooling.
Blazars are promising candidates for high energy neutrino sources, yet the physical origin of their neutrino emission remains uncertain. In this work, we extend our previous study by modeling the broadband spectral energy distributions (SEDs) of 103 neutrino blazar candidates (NBCs) within a hadronic framework. To estimate the maximum possible neutrino output, we adopt an assumption in which the high energy emission is dominated by p gamma interactions and the contribution from leptonic inverse Compton scattering is strongly suppressed. From the SED modeling, we constrain nine key parameters describing the emission region and particle energy distributions. We perform a partial correlation analysis to investigate the relationship between neutrino luminosity and electromagnetic emission, and we found a weak or moderate correlation between optical R band and neutrino emission. Our model predicts prominent proton synchrotron emission peaking in the MeV band for most sources, with 99 out of 103 NBCs exhibiting proton synchrotron peaks within 0.1 to 100 MeV, highlighting the MeV band as a key window for distinguishing between leptonic and hadronic scenarios. Based on the model-predicted maximum neutrino fluxes, we find that three NBCs are potentially detectable by IceCube, while up to 22, 45, and 62 sources may be detectable by KM3NeT, NEON, and TRIDENT, respectively. These results provide testable predictions for future multi-messenger observations and offer new insights into the composition and radiation mechanisms of blazar jets.
We propose a non-parametric approach to construct the statistical equation of state (EOS) continuously from the nuclear crust to the asymptotic-freedom regime. Driven by the observationally required stiffening to support two-solar-mass neutron stars (NSs) with relatively small radii for low-mass NSs, this global thermodynamic constraint suggests a clear peak of squared sound speed ($c_s^2$) in massive NSs. To prevent overshooting perturbative QCD (pQCD) energy-density bounds, this early stiffening must be actively compensated by an extended density range of softening, with $c_s^2$ not approaching $1/3$ until $\sim\!30\,n_{\rm sat}$. Consistently, the trace anomaly $Δ\equiv 1/3 - p/ε$ becomes positive beyond NS densities and approaches the pQCD limit from above. This natural emergence of $Δ> 0$ organically aligns with some anticipated microphysics, likely arising from the pressure dilution in a quark-hadron mixed phase, non-conformal pQCD corrections to quark-gluon interactions, or the symmetry-breaking effects of finite strange quark mass. By measuring the degree of this non-monotonic behavior in the posterior, we find evidence for a hadron-quark phase transition in the cores of the most massive neutron stars. This indicates that the non-perturbative quark matter is intrinsically soft, fundamentally distinguishing it from the stiff scenarios associated with the quark-star picture.
We present 30 newly discovered repeating fast radio burst (FRB) sources from the second catalog of bursts detected by the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB). These repeaters have extragalactic dispersion measures (DMs) spanning $99.4-1446.0\ \text{pc cm}^{-3}$ and burst rates between $10^{-5.7}$ and $10^{-0.5}$ hr$^{-1}$ scaled to a fluence threshold of 5 Jy ms. We report evidence of monotonic, linear DM variations in four repeaters on years-long timescales. The newly discovered sources bring CHIME/FRB's total number of observed repeating FRBs to 80, 79 of which were discovered by CHIME/FRB, between 2018 July 25 and 2023 September 15. In the full CHIME/FRB sample, only 2.4$\pm 0.4\%$ of sources have been observed to repeat, and we do not find evidence for significant evolution of this value over the duration of the experiment. We find no substantial evidence for bimodal populations of one-off and repeating FRBs in their burst rate distributions; the distribution of upper limits on repeat rates implied from observations of as-yet one-offs is entirely contained within the observed range of repeater burst rates and the distributions do not appear inconsistent. Similarly, using the population analysis framework of C. W. James (2023), we find that our observations of repeating and yet-one-off FRBs are equally well fit assuming a power-law distribution of repeat rates with 50$-$100% of the population repeating.
Perseus is the brightest X-ray strong cool-core (SCC) cluster, with a bright central radio and $γ$-ray source plus low-frequency radio mini and giant halos. It is the archetype of the cooling flow (CF) problem, with X-rays implying mass cooling rates orders-of-magnitude larger than observed in other channels. Recent work suggested that ancient ($\gtrsim$\,Gyr-old) cosmic ray (CR) halos (ACRHs), injected by the central source, would produce thermal-like soft X-ray inverse-Compton (CR-IC) emission 'boosting' the CC and alleviating the CF problem. We examine Perseus and show that a simple model of CRs injected by NGC 1275 (+satellites) simultaneously accounts for the excess CF luminosity and minihalo. The models reproduce Perseus's soft X-ray surface brightness and X-ray inferred density/temperature/pressure/metallicity/cooling time/mass deposition rates; $γ$-ray spectra; extended hard X-rays; and radio surface brightness and spectral index data, from kpc-Mpc. These also reproduce independent constraints on magnetic field strengths and mass/potential models. The evolution of the minihalo spectral index and surface brightness are predicted by an aging population of CRs boosting the apparent SCC luminosity via CR-IC, and match well the observed hard X-ray slopes. The 'giant' low-frequency halo can be predicted by the sum of ACRHs around satellites distributed throughout the cluster, dominating diffuse synchrotron at $\gtrsim 100\,$kpc. Re-acceleration is neither needed nor important in these models, and implied CR transport speeds are consistent with buoyant advection. Previous claims of upper limits to non-thermal X-rays and CR pressure relied on strong assumptions which are not valid at the CR energies of interest, e.g. a power-law spectrum of CRs. This could resolve many historical puzzles about Perseus, and makes new predictions for future observations.
We present the Swift Analysis Pipeline for Lightcurve Extraction (SAPLE), a semi-automated pipeline to extract the Swift-UVOT and Swift-XRT data products and spectral information (magnitudes, photon indices, and fluxes) for a set of observations of any point source of interest. This pipeline is not meant to substitute, but to complement the tools the Swift team has already set up. Specifically, SAPLE provides a Swift-UVOT semi-automated pipeline that also returns the absorption corrected specific fluxes for any observation and filter of interest, a tool which to our knowledge is not publicly available to the community yet. Moreover, for Swift-XRT, SAPLE enables the user to extract a lightcurve of both flux and photon index (with associated uncertainties), assuming a redshifted powerlaw spectrum. The main codes are available through a GitHub repository (L. Marcotulli & N. Torres-Albà 2026), and the following paper summarizes the main steps of the analysis.