Papers

We present a multi-epoch, primarily X-ray study of the narrow-line Seyfert 1 (NLSy1) galaxy Mrk~766 by combining archival XMM-Newton observations spanning nearly two decades together with a XRISM observation obtained during its Performance Verification Phase. We analyzed the absorption and emission features in the soft X-ray spectra from RGS and EPIC-pn data, revealing a stratified multi-phase potentially ionised absorber with a column density varying at a timescale of $\sim$1 day. The covering fraction shows positive correlation with the continuum flux, however with significant intrinsic scatter, between the epochs. Mrk~766 exhibits a multi-component Fe\,K$\alpha$ emission line as revealed by fits to the XRISM data. A quasi-phenomenological model comprising a torus and a broad diskline component distinguishes the two components and constrains the inner radius of the potential broad line emitting disk substructure between $R_{in}=40-60\,r_{\rm g}$. Analysis of multiple epochs of the Fe K$\alpha$ complex across XMM-Newton observations reveals that the broad component tracks the continuum flux closely and effectively reflects the trend with respect to the continuum. The narrow component is consistent with a distant emitter, potentially the torus. Quasi-simultaneous monitoring with Swift-XRT and UVOT further indicates UV emission lagging the X-rays by $5.9^{+4.1}_{-5.7}$ hours, consistent with disk reprocessing of coronal X-ray fluctuations, and corresponds to a light travel distance of up to $10^{3}\,r_{\rm g}$. Overall, our analysis maps the extent of the dynamically evolving absorber, the iron line emitter, and the accretion disk in this highly accreting NLSy1 system.

Recent IceCube observations point to the ultra-high energy neutrino flux exhibiting a very soft spectral nature beyond tens of TeV, which is unexpected from the standard cosmic-ray neutrino connection. As a potential explanation for this behaviour, we investigate neutrino self-interactions in an extra-dimensional $ {U(1)}_{L_{\mu}-L_\tau} $ gauge theory, where symmetry breaking leads to a tower of Kaluza-Klein gauge bosons in compactified four dimensions. These multiple gauge bosons may enhance the attenuation of astrophysical neutrinos as they propagate in the C$\nu$B medium. Unlike single-mediator scenarios, for example, which arise from broken $ {U(1)}_{L_\mu - L_\tau} $ gauge symmetries in four dimensions, the presence of a Kaluza-Klein tower of mediators produces multiple closely spaced resonances whose interference gives rise to a rich energy-dependent behaviour of the scattering cross-section over a vast range of incident neutrino energies. Alongside $s$-channel resonances, off-resonant $t$- and $u$-channel contributions also become important. We explore the possibility that repeated resonant scattering between astrophysical neutrinos and those in the C$\nu$B, mediated by these new multiple gauge bosons may increasingly attenuate the former at higher energies, thereby addressing the possibility of softening its spectral nature within the consideration of a single power law.

We analyze observations of the Cassiopeia A core-collapse supernova remnant (CCSNR) and identify a pair of opposite rings, attributing their formation to unequal-energy explosion jets that also imparted the neutron star (NS) its kick velocity via the kick-BEAP (kick by early asymmetrical pair of jets) mechanism. We replicate the formation of these ring structures, which we identify as circum-jet rings, with three-dimensional hydrodynamical simulations of two opposite explosion jets in which the northern jet is much more powerful, forming the larger northern ring known as the Crown. The line connecting the centers of the two opposite rings passes through the previously-established point-symmetric center of Cassiopeia A. The momentum difference between the two jets also explains the observed southwards NS kick velocity of ~400 km/s. The kick-BEAP and point-symmetric morphology suggest powering and shaping by several pairs of jets during the explosion process, strengthening the jittering jets explosion mechanism for the Cassiopeia A CCSNR and core-collapse supernovae in general.

We present a three-dimensional visco-magnetoacoustic framework for analog gravity in magnetohydrodynamic (MHD) flows. While standard fluid dissipation typically breaks the Lorentzian signature of acoustic metrics, we demonstrate that evaluating wave perturbations in the eikonal limit alongside the Shakura-Sunyaev $\alpha$-viscosity prescription preserves a well-defined effective spacetime geometry for the fast magnetoacoustic mode. The slow-magnetoacoustic mode and Alfvén mode do not admit a non-degenerate metric. To investigate analog horizon thermodynamics, we utilize astrophysical accretion disks as background media, specifically modeling numerical magnetized advective accretion flows around rotating black holes and analytical advection-dominated inflow-outflow solutions (ADIOS). Standard self-similar ADIOS models strictly enforce a constant Mach number, precluding horizon formation. We therefore introduce a magnetic field perturbation that breaks self-similarity, generates a dynamic Mach number, and enables the formation of a visco-magnetoacoustic horizon. By evaluating the spontaneous phonon emission at these horizons, we reveal that the analog Hawking temperature is highly sensitive to the magnetic field topology; the spatial orientation of the background magnetic gradients dictates whether the Hawking radiation is amplified or suppressed. Furthermore, we find that increasing the viscosity parameter leads to a monotonic increase in the analog Hawking temperature.

Kehle and Unger have recently shown that, under suitable conditions, a massless charged classical scalar field can produce an extremal black hole in finite time. This result appears to violate the third law of thermodynamics (TLT), since extremal black holes have zero temperature. However, because the TLT is intrinsically a quantum principle, a self-consistent analysis of this problem must account for quantum effects. In particular, we point out that the electric field of the system may spontaneously polarize the vacuum, thereby creating massless charged pairs that will discharge the black hole. This mechanism guarantees the quantum validity of the TLT in the regime $|eQ|\geq\hbar/2$ (here ${Q,e}$ denote, respectively, the electric charges of the black hole and the scalar field). This leaves open the question of what happens for lower values of $|eQ|/\hbar$. Recently, Schneider derived a classical bound showing that extremal black holes cannot be formed à la Kehle-Unger in the small-charge regime $|eQ|\leq\hbar/3$. His analysis leaves open the possibility that a stronger bound may exist. We conjecture here that this bound can indeed be strengthened such that classical considerations would prevent extremal (zero-temperature) charged Reissner-Nordström black holes from forming dynamically in the regime $|eQ|<\hbar/2$, where quantum polarization effects are too weak to preserve the validity of the TLT.

A sufficiently strong hadron-to-quark first-order phase transition can give rise to a third family of stable compact stars that are commonly referred to as hybrid hadron-quark stars. Stable hybrid stars that have the same gravitational mass as neutron stars are referred to as twin stars. Although equilibrium twin stars may exist, whether they can be dynamically formed remains an open question. We investigate this problem by examining the gravitational binding energy of competing equilibrium configurations at fixed baryonic rest mass and by performing general relativistic hydrodynamical simulations of several possible transition channels. While twin stars are more gravitationally bound than neutron stars with the same rest mass, this energetic preference alone does not determine the dynamical outcome. Compression and shocks during the evolution generate thermal pressure that can prevent the system from settling on the cold twin star branch. Our simulations show that sufficiently rapid cooling can remove this thermal support and enable twin star formation, whereas slower or no cooling generally favors a neutron star remnant. Accessing the twin star branch through the formation channels considered here requires cooling on a timescale comparable to or shorter than the stellar dynamical timescale. Since realistic cooling mechanisms operate on much longer timescales, our results suggest that in channels that conserve the total rest-mass neutron stars may be dynamically favored even when a more gravitationally bound twin star configuration exists with the same rest mass. Our results demonstrate a point of principle, at least for equations of state where the quark deconfinement density does not change appreciably for temperatures up to $10-20\,\textrm{MeV}.$

The "Little red dots" (LRDs) are a population of accreting supermassive black holes (SMBHs) in the early Universe which often exhibit undermassive or even undetectable stellar hosts. Their early emergence, high space density, and extremely large black-hole-to-stellar mass ratios pose a serious challenge to conventional seeding scenarios that rely on baryon for both the formation and growth of black holes. Here we demonstrate that the above anomalies can be naturally resolved if dark matter is self-interacting. We apply a fully relativistic, non-equilibrium halo-evolution model, first developed in our earlier work, to trace the complete gravothermal evolution of self-interacting dark matter (SIDM) halos, from the initial collapse into BH seeds to the subsequent accretion of dark matter. We find that in highly concentrated halos assembled before reionization, gravothermal collapse efficiently produces stellar-mass black-hole seeds within a few hundred million years. Remarkably, and contrary to standard expectations for dark-matter accretion, heat conduction in SIDM then sustains a prolonged super-Bondi inflow that drives these seeds to supermassive scale by the LRD epoch, without baryonic assistance. The halo conditions required for completing these processes, together with the probability of avoiding major mergers that disrupt gravothermal evolution, result in an SMBH population consistent with the observed abundance and redshift distribution of LRDs. Our findings establish a pathway in which SMBHs are seeded and assembled primarily from dark matter, well before substantial galaxies form around them, thereby offering both a compelling physical explanation for LRDs and a new observational probe of dark-matter microphysics.

The ultraviolet continuum from a high-redshift accreting black hole is routinely used to infer its mass. GN-z11 offers a sharp test of whether that inference is secure. A thin-disk fit to its continuum gives an Eddington mass of $M_{\rm E}=1.12\times10^{7}\,\Msun$, an order of magnitude above the $\log(\Mbh/\Msun)=6.2\pm0.3$ inferred from broad N~iv. We show that this apparent tension is not imposed by the ultraviolet slope alone. We derive a closed-form criterion comparing the radius where observed photons are produced with the radius where a supercritical flow departs from efficient solution. For GN-z11, the criterion places the supercritical transition inside the ultraviolet-emitting region. Composite disks retaining the outer solution confirm this numerically, requiring $\Mdot_{\rm out}=3.1$–$4.4\,\Msun\,\mathrm{yr^{-1}}$ and placing $\Rsph=(1.2$–$1.8)\times10^{3}\,\rg$ inside $\Ruv=(2.6$–$2.9)\times10^{3}\,\rg$. Replacing the inner $\sim10^{3}\rg$ by a supercritical flow shifts the fitted slope by only $\Delta\betaUV=+0.08$ to $+0.11$, comparable to uncertainties from standard spectral modelling. Reversing the radial ordering requires the effective transition to move outward by a factor $1.7$–$2.1$. Within our composite models, the measured slope constrains black-hole mass only to around $10^{6}\,\Msun$, an order of magnitude below the continuum Eddington mass. Continuum Eddington masses should therefore be treated as model-dependent bounds that assume global thin-disk efficiency. The same hierarchy appears in three additional JWST sources under the disk-dominated interpretation, where separation masses are only $4$–$6$% of continuum Eddington masses. This distinction can substantially weaken the seed-mass and early-growth demands inferred from ultraviolet continua at cosmic dawn for the first massive black holes.

Thermal conduction can act as a heating mechanism for cool-core clusters, yet the same transport must erode any temperature discontinuity it crosses. Ophiuchus permits both effects to be tested in the same atmosphere. Its temperature rises from about $1$ keV in the innermost core to $9$ keV at $r\sim30$ kpc, making conductive heating especially favorable, while two sharply resolved cold fronts independently constrain the same transport. \textcolor{black}{We construct a steady radial energy balance from deprojected Chandra density and temperature profiles, with XRISM-based limits on turbulent heating and a subsonic inflow contribution.} Closing the balance with conduction requires only sub-Spitzer coefficients — \Ftwentyfive\ at 25 kpc and \Finner\ and \Fouter\ at the inner and outer cold fronts — so conduction can energetically supply the missing heat. Yet the same coefficients, applied normal to the cold-front surfaces, would broaden the fronts to their observed width limits within \Tinner\ and \Touter\ Myr, short compared to the characteristic sloshing timescales that generate and sustain such interfaces. For a representative 200 Myr residence time, the radial conductivity demanded by the energy balance exceeds the cross-front conductivity permitted by front survival by roughly a factor of three, with the required anisotropy exceeding unity in 99.9% and 98.8% of accepted models. The result leaves two broad regimes. Thermal conduction is either subdominant in Ophiuchus or contributes substantially but is strongly anisotropic, with heat transport suppressed across the cold fronts relative to the radial direction. Together they disfavor locally isotropic conduction as a dominant heating mechanism in cool cores generally.

The radio/X-ray/black-hole-mass Fundamental Plane (FP) is widely used to estimate masses of accreting intermediate-mass black hole (IMBH) candidates. Its inverse use is justified only when radio luminosity traces a compact jet core and X-rays arise from the same sub-Eddington hard-state accretion flow. We formulate an accretion-state validity criterion for inverse-FP masses using verified BH X-ray binary radio/X-ray tracks as the empirical reference. Using 424 detections from 21 BHXBs and candidates, we find that both standard and curved tracks are present and the luminosity at which radio/X-ray behavior changes is source dependent. A low X-ray Eddington ratio is therefore not, by itself, sufficient to establish inverse-FP validity. Within a scale-invariant jet/hot-flow framework, track steepening reduces the effective mass leverage of the inversion, while at much lower accretion rates uncertainties in the dominant emission mechanism and compact-core production can invalidate the underlying FP assumptions altogether. We apply the criterion to representative IMBH-relevant systems without fitting a new plane. NGC 4395 is low-Eddington but compact-radio deficient: its VLA-scale and VLBI jet-base luminosities imply formal inverse-FP mass underestimates of about 1.0 and more than 2.8 dex if the deficit is misread as a mass signal. Globular-cluster X-ray populations in NGC 1399 and NGC 4472 are source-identity limited: 71% and 67% lie below $L_{\rm Edd}(10\,M_\odot)$, and the high-luminosity tail alone does not identify an IMBH accretor. The radio non-detection of the $\omega$ Centauri IMBH candidate constrains gas supply, accretion efficiency, and compact-core production rather than excluding the black hole. The resulting test defines which SKA/ngVLA-era radio cores can be converted into FP masses, and which instead constrain gas supply, source identity, or compact-core production.

Next-generation gravitational-wave observatories like Cosmic Explorer will provide an unprecedented opportunity to probe the cosmological stochastic gravitational-wave background, providing insight into the early Universe that is inaccessible with other means. We investigate the prospects for measuring a Gaussian cosmological gravitational-wave background in Cosmic Explorer in the presence of a distinguishable, non-Gaussian foreground of binary black hole mergers, using a Bayesian search that simultaneously fits binary black hole and background parameters. We successfully recover a Gaussian cosmological background with a signal-to-noise ratio of 5 in one day of simulated data in Cosmic Explorer without contamination from unresolved signals or data exclusion. However, when introducing waveform systematic errors to our analysis, our inference on the stochastic background parameters becomes biased and can falsely infer the presence of a background when there is none. The bias is driven by signal mismatches at low frequencies (<20 Hz), especially when the recovery waveform is missing physical effects included in the injected signals. While we propose a metric for diagnosing if measurements will be strongly biased by waveform systematics, our results suggest improvements in waveform modeling or other mitigation techniques will be required to accurately recover a primordial stochastic gravitational-wave background in next-generation detectors.

We provide transition probabilities and Maxwellian-averaged collision-strengths for the electron-impact excitation of Ce ii - iv. These ions are potentially abundant in kilonovae (KNe) and feature both forbidden and allowed transitions in the infrared. These atomic data are obtained using the grasp$^0$ atomic structure code, interfaced with the non-perturbative darc R-matrix close-coupling codes. We additionally use autostructure to produce dielectronic and radiative recombination rate coefficients for Ce ii - vi. Using this data, we calculate an approximate ionization balance including non-thermal ionization rates with deposition rates from realistic models. We find Ce ii - iv to be of roughly equal abundance, subject to ejecta conditions. We find that generally we need $\gtrsim$ 10$^{-3}$ M$_{\odot}$ of Ce in total to produce emission of comparable luminosity to the features in the spectrum of AT2023vfi. In particular, it is found that [Ce iv] 4.5 $\mu$m requires both strong ionization and weak recombination, as well as a production of Ce considerably stronger than that of the solar abundance in order to have luminosity comparable to the apparent emission in this wavelength vicinity in \grb. %which is seemingly incompatible with the current understanding of the bulk composition of KNe. With such strict requirements on individual element masses, the detection of individual lanthanides on a line-by-line basis in the nebular phase seems unlikely - but lanthanides may remain important for cooling and early time analyses.

When propagating through the universe, gamma rays at very-high energy (VHE, E > 100 GeV) and ultra-high energy cosmic rays (UHECRs, E > 1 EeV) can interact with the optical, infrared and microwave photon fields that permeate the universe. These interactions result in a characteristic absorption imprint in the spectra of extragalactic gamma-ray sources at VHE, and in a change in the mass-composition of UHECRs. The study of both VHE gamma rays and UHECRs therefore requires precise knowledge of the intensity of the cosmic photon fields. In this work, we explore the impact of the current uncertainties on the optical and infrared photon fields on the propagation of astroparticles. We restrict the range of available models to those best matching the recent measurements, and compare the different reconstructions resulting from these models. We find that the knowledge on the cosmic background light is no longer the dominant source of uncertainties in understanding the phenomenology of astroparticle sources in the low redshift universe (z < 0.1), enabling robust spectral and composition inference with the next generation of both gamma-ray and UHECR measurements.

The neutron skin, defined by the difference between neutron and proton root-mean-square radii, is a characteristic manifestation of isospin asymmetry and an important probe of the isovector nuclear interaction. This focused review examines how neutron skins influence particle emission and collective dynamics in heavy-ion collisions, from the Fermi-energy regime to ultra-relativistic energies. By modifying the initial neutron and proton density profiles, the neutron skin affects the isospin composition and geometry of the participant region, pre-equilibrium emission, particle production, fragment formation, and collective flow. We review neutron-to-proton and $\rm{t}/^3\rm{He}$ yield ratios, light clusters, pion ratios, bremsstrahlung photons, isoscaling and fragment momentum distributions, and neutron-proton differential flow and momentum observables, emphasizing their interplay with the symmetry energy and transport dynamics. At high energies, neutron skins also modify the initial geometry, eccentricities, multiplicities, and anisotropic flows in isobar and heavy-nucleus collisions. We discuss the challenge of disentangling these effects from deformation, surface diffuseness, shell structure, clustering, and model dependence. Broader connections to parity-violating electron scattering, dipole responses, coherent elastic neutrino-nucleus scattering, SRC-induced proton skins in momentum space, and neutron-star observables are also explored. Finally, we highlight opportunities from radioactive beams, improved collision experiments, microscopic many-body and transport calculations, and Bayesian inference. Combining multiple reaction systems and observables with complementary nuclear-structure and astrophysical information will be essential for quantitatively constraining neutron skins and the density dependence of the symmetry energy.

The interpretation of gravitational wave sources depends on the specific choice of model used to interpret signals. Previous analyses of GWTC-3 event candidates using SEOBNRv4PHM, IMRPhenomXPHM, and NRSur7dq4 have sometimes arrived at notably different conclusions about event properties. Subtle analysis settings and code differences can also produce notable differences. We revisit 60 publicly available events from the first three observing runs, released across GWTC-1, GWTC-2, GWTC-2.1, and GWTC-3, using a consistent analysis framework and three waveform models: two state-of-the-art time-domain models with higher-order modes (SEOBNRv5PHM and IMRPhenomTPHM) and one older frequency-domain model without higher-order modes (IMRPhenomPv2). The two state-of-the-art models agree well overall, although about 20% of events have a visible difference in at least one one-dimensional marginal posterior. By contrast, the older model often arrives at qualitatively different conclusions for key events across the mass spectrum. We also identify a bug in published LVK GWTC-2.1 SEOBNRv4PHM results and in a later reprocessing of GW200105. At the adopted sampling rates, the model's Nyquist-frequency restriction excluded part of the allowed parameter space and artificially truncated the mass-ratio posterior away from equal mass for seven low-mass events. These results emphasize the need to use multiple state-of-the-art waveform models and carefully validated analysis settings to characterize uncertainty in gravitational-wave source properties.

Astrophysical neutrinos are a clear identifier of photohadronic interactions in sources of Ultra-High-Energy Cosmic Rays (UHECRs). As the Pierre Auger Observatory has shown that UHECRs exhibit a composition heavier than protons at Earth, we expect that these heavy primaries undergo nuclear cascades due to photodisintegration in the dense source environments from which they originate. In this contribution, we will derive the neutrino spectra including those originating from nuclei undergoing photohadronic interactions. This will be explored in the context of modeling the emission region of a gamma-ray burst. Photohadronic interactions will be modeled using the new framework for Cosmic Ray Stochastic Interactions for Propagation (CRISP), which employs an analytic approach to compute the underlying probabilistic description of UHECR interactions.

The paper discusses the fraction of extragalactic cosmic iron nuclei with energies of 10^19 eV that reached Earth from distant sources without fragmentation in intergalactic space. It is shown that the fraction of particles with energies E = 10^19 eV is significant (about 10%), even if the redshifts of the sources reach 0.2, and not only for nearby sources.

GRS 1915+105 is one of the most studied black-hole X-ray binaries, known for its extreme variability and rich phenomenology. Previous studies of this source with the Rossi X-ray Timing Explorer reported a transition of the phase-lags of type-C quasi-periodic oscillations (QPOs) from soft, where low-energy photons lag the high-energy ones, when the QPO frequency is higher than ${\sim} 2 ~{\rm Hz}$, to hard when the QPO frequency decreases below ${\sim} 2 ~{\rm Hz}$. The hard-lags of the QPO coincide with episodes of strong radio emission. We analyse NICER observations of GRS 1915+105 obtained between 2018 and 2020, during a period in which the source flux decreased steadily, and perform a detailed spectral-timing study of the detected type-C QPOs. We find a type-C QPO with frequencies in the range of ${\sim} 1.3$-$3.9 ~{\rm Hz}$, which displays soft lags and, contrary to the RXTE observations, shows no evidence of hard phase lags at frequencies below $2~{\rm Hz}$. Quasi-simultaneous AMI-LA radio observations show consistently low radio flux (${\lesssim} 5 ~{\rm mJy}$) during this period. These results appear to show that the hard QPO lags in GRS 1915+105 are linked to the presence of strong radio activity, suggesting that the relativistic jet is responsible for the hard phase lags, supporting a scenario in which QPO phase lags trace changes in coronal geometry and accretion-ejection coupling in GRS 1915+105.

The magnetic reconnection process near a black hole can efficiently convert magnetic energy into the kinetic and thermal energy of plasma outflows, enabling the extraction of rotational energy from a spinning black hole. However, its impact on the surrounding accretion disk remains insufficiently explored. In this work, we develop a framework to investigate the 'heating' effect of magnetic reconnection outflow plasma on thin accretion disks around Kerr black holes. We incorporate the energy and angular momentum fluxes of reconnection outflow into the disk conservation laws as additional source terms, yielding a modified radiative flux. We investigate the dependence of the modified radiative flux on the black hole spin, reconnection radius, plasma magnetization, and outflow orientation. For a fixed magnetic field prescription, the contribution of the reconnection outflow plasma to the modified radiative flux decreases as the black hole spin increases. Its dependence on the reconnection radius differs inside and outside the ergosphere: inside, a smaller radius enhances disk heating through more efficient extraction of black hole rotational energy, whereas outside, a larger radius produces stronger heating due to reduced radial dilution of the outflow. The reconnection outflow plasma contribution also increases with plasma magnetization and the radial component of the outflow velocity. Our results demonstrate that magnetic reconnection can provide an additional source of energy and angular momentum for relativistic accretion disks, establishing a connection between plasma processes near rotating black holes and the observable accretion disk properties, providing a new perspective on the strong gravitational field regime.

We derive an extension of Chandrasekhar's pulsational stability criterion for relativistic stars with a pressure boundary condition at the surface. We then apply the new criterion to hylotropes, a model describing rapidly accreting supermassive stars with a core-envelope structure. We find that confinement typically stabilises stars, delaying the onset of the general relativistic instability when the external pressure at the surface is roughly $\sim10^{-3}$ of the stars' hydrostatic pressure scale. Beyond this value, the critical mass at marginal stability increases rapidly, exceeding the value of $\sim 10^5$ M$_{\odot}$ for helium burning, isolated $n=3$ polytropes by orders of magnitude. We discuss how the required pressure boundary conditions may be produced by the combined ram and magnetic pressure supplied by the accretion flow as it assembles the supermassive star, and provide a fitting formula for the final mass before collapse. As a rough estimate assuming equipartition between magnetic and thermal pressures in a standard $\alpha$-disk, we find that the critical mass doubles at accretion rates of order $\dot M\sim100\,{\rm M_\odot}\,{\rm yr}^{-1}$, again increasing rapidly after this value. Our results show under what conditions, potentially realized in a subset of high redshift halos, non-rotating SMS may collapse into black holes beyond the typical scale of heavy seeds.

We present a deterministic machine-learning method for denoising candidate signals from all-sky searches for continuous gravitational waves from rotating neutron stars. Using the time-domain $F$-statistic search pipeline, we post-process the resulting $F(f,\dot{f})$ candidate patterns with iterative $\alpha$-(de)blending, a deterministic diffusion-type generative model implemented as a U-Net, well suited to the non-Gaussian, correlated noise of the $F$-statistic output. Two models, trained on data based on simulated 6- and 12-day time-domain segments with software-injected signals similar to the LVK hardware injections, are tested by comparing deblended images to the library of expected signal patterns via the Structural Similarity Index Measure. In general the method recovers the correct sky-position-dependent pattern for signal-to-noise ratios $\rho \gtrsim 4$, moderately below typical all-sky detection thresholds, with recovery depending strongly on pattern morphology, demonstrating that deterministic diffusion-based denoising may serve as a consistency/veto tool ahead of the semi-coherent coincidence stage. We also discuss limitations and possible improvements of the method.

Context. Swift J1727.8-1613 is an X-ray binary observed during its 2023-2024 outburst, including bright and dim hard states with similar X-ray polarimetric signatures but fluxes differing by two orders of magnitude. Aims. We aim to reproduce the spectro-polarimetric properties of these two hard states. Methods. We modelled their spectra in the JED-SAD framework and computed the X-ray polarisation properties of the best-fit configurations with the MONK general relativistic radiative transfer code. Results. Both spectra are well reproduced within the JED-SAD framework, although the dominant radiative process differs: self-Comptonised bremsstrahlung in the bright state and self-Comptonised synchrotron in the dim state. In the bright state, the simulated polarisation fraction remains below the observed value of about 4%, while the polarisation angle differs from the observations by 90 degrees because of the JED's large vertical optical depth. In contrast, the dim state's simulated polarisation fraction is about 3% and its polarisation angle is perpendicular to the disk, both broadly consistent with IXPE measurements. The dim-state results are largely insensitive to the magnetic field configuration because of Faraday depolarisation. Conclusions. The dim hard state's spectro-polarimetric properties are well reproduced, whereas those of the bright state are not. We discuss a non-thermal particle distribution and electron scattering in an accretion-disk wind as possible explanations for the observed bright-state polarisation.

The closest and brightest millisecond pulsar (MSP), PSR J0437-4715, is a high-priority pulsar for precision timing experiments, including the Parkes Pulsar Timing Array (PPTA) project, where it is used for both pulsar timing and for polarization calibration. Using approximately seven years of observations obtained with the Ultra-Wideband Low-Frequency receiver on Murriyang, CSIRO's Parkes radio telescope, we investigate the long-term spectro-polarimetric variability of PSR J0437-4715. We report the discovery of two consecutive discrete profile change events exhibiting the same phase-localised morphological variations, the first such behaviour reported in any MSP. Both events affect three phase-localised regions of the pulse profile that exhibit coupled and anti-correlated temporal evolution. We also identified, and corrected for, systematic profile variability induced by polarization calibration errors arising from long-term drift of the polarimetric properties of the receiver. The profile change exhibits frequency-dependent evolution in total intensity, together with variations in linear polarization, position angle and fractional linear polarization. These results confirm a magnetospheric origin for the events rather than propagation effects caused by the interstellar medium, and provide potential new constraints on MSP emission models.

Temporal gaps in Swift-XRT X-ray afterglows of gamma-ray bursts (GRBs) limit how precisely the standard empirical plateau function constrains three parameters: the plateau end time $T_a$, the plateau flux $F_a$, and the post-plateau decay index $\alpha$. These parameters underpin the luminosity-time and fundamental-plane correlations of GRB plateaus so that smaller measurement uncertainties can improve the statistical power of GRBs as cosmological probes. As the fifth study in a series on light-curve reconstruction, this work benchmarks four models on a parent sample of 545 Swift-XRT GRBs: a Physics-Informed Neural Network (PINN) with seven afterglow-motivated priors, ReFANN, a Siamese dual-branch network with Monte Carlo dropout, and Polynomial Quantile Regression (PQR). All four reduce the fractional uncertainties in $\log T_a$, $\log F_a$, and $\alpha$ relative to the original observations. The PINN with a broken power-law prior gives the largest reductions ($\sim$41-49%), at a higher outlier rate ($\sim$15-20%), while ReFANN, Siamese network and PQR give more modest reductions ($\sim$19-27%) with outlier fractions $\lesssim$4%, across a wider range of light-curve morphologies. A reduced-$\chi^2$ prior-selection scheme recovers a morphological classification consistent with independent labelling, though it does not match the best single-prior reduction. Together, these results provide a systematic benchmark of physics-informed and data-driven reconstruction strategies for irregularly sampled GRB afterglows.

Fast radio bursts (FRBs) are millisecond-duration extragalactic pulses whose discovery requires searches over large signal-parameter spaces and the rejection of candidate sets dominated by radio-frequency interference and noise. Integrating quantum processors into FRB searches requires a workflow connecting telescop data to hardware-compatible models. Here we develop a pipeline that converts raw search-mode PSRFITS data into dynamic-spectrum segments and identifies FRB-like signals using a hybrid Quantum Vision Transformer (QViT). On labelled FAST data, the selected QViT achieves a mean accuracy of 94.00% and recall of 98.30%, with similar performance to a compact classical Vision Transformer. Recall reaches 97.11% on an FRB source excluded from training and model selection. Real-device execution on 82 segments yields a mean simulator-hardware Hellinger fidelity of 0.9161 +/- 0.0074. The workflow connects raw telescope data to quantum-assisted identification and provides an experimental basis for assessing the measurement requirements that must be addressed before survey-scale deployment.

We report new optical spectroscopic results for the gamma-ray blazar PKS~0903$-$57. The optical properties of this source have historically been complicated to establish, due to a bright foreground star separated from the blazar by only 0.67$^{\prime\prime}$. During its strong high-energy (HE) and very high-energy (VHE) gamma-ray flaring period between 2020 and 2022, the Southern African Large Telescope and Very Large Telescope (VLT) took spectra of its optical counterpart. Photometric monitoring with the Rapid Eye Mount (REM) telescope complemented the observations. The observational results were inconclusive. A major improvement was achieved with a repeat of the VLT observations in 2024 using a 0.5$^{\prime\prime}$ slit under excellent seeing conditions. This led to the detection of five narrow emission features in the spectra, yielding a robust redshift measurement of $z$ = 0.2621 $\pm$ 0.0006. We also detected a broad symmetric H$\alpha$ emission component with a full width at half maximum of 4020 $\pm$ 30 km/s in the spectra. Non-thermal jet emission was found to dominate the continuum, supporting the classification of PKS~0903$-$57 as a flat-spectrum radio quasar. Another interesting observation was a significant velocity displacement of $\sim$1500~km/s between the broad H$\alpha$ line and its corresponding narrow line. Such a displacement may be associated with unusual dynamics in the broad-line region or accretion flow, a system containing two supermassive black holes, or a recoiling black hole produced following a supermassive black-hole merger. Building on these results, we carried out additional VLT, REM, and Swift observations towards the end of 2024 and during 2025. We present in this paper the preliminary results from the analysis of these new data.

Using nearly 17 years of Fermi-LAT data, we analyzed the diffuse gamma-ray emission associated with different spiral arms of the Milky Way in the ranges $l = 105^\circ$ to $145^\circ$ and $b = -5^\circ$ to $5^\circ$. The gamma-ray emissions from these spiral arms all exhibited clear pion-bump features, which can be well explained by a power-law spectrum of cosmic-ray protons. The spectral indices corresponding to the three distinct spiral arms – the Local Arm, the Perseus Arm, and the Outer Arm – are found to be approximately $-2.75$, $-2.55$, and $-2.80$, respectively. The energy densities of the cosmic rays above 10 GeV are estimated to be all around 0.2$\rm~eV~cm^{-3}$. We further compared our results with previous studies, which have established various propagation models predicting that cosmic ray density decreases and spectra soften with increasing Galactocentric radius. We found higher densities and harder spectral indices than predicted in the outer Galaxy, suggesting possible contributions from a larger halo, slower magnetic field decay, or a combination of these and other effects. In particular, the spectral hardening observed in the Perseus Arm may be attributed to locally enhanced cosmic-ray acceleration.

Equations of state (EoS) for dense matter are commonly provided in tabulated pressure-energy density relations, complicating numerical implementation and potentially compromising thermodynamic consistency. In this work, we propose a piecewise-continuous analytical form with a common parametrization capable of representing a broad class of dense matter EoSs. We validate this parametrization against tabulated EoSs from the CompOSE and LALSimulation repositories. Following the initial deterministic fitting, we employed two distinct Bayesian approaches: one based on synthetic EoSs generated by adding noise to the fitted EoS and another based on multi-messenger measurements of neutron-star masses and tidal deformabilities. In both approaches, the initial best-fit parameters are used as reference values. For the considered EoSs, spanning from very soft to very stiff, the resulting fits reproduce the tabulated EoSs and the associated macroscopic neutron-star observables in the repositories with sufficient accuracy. The inferred parameters exhibit strong correlations arising from the continuity conditions imposed at the segment boundaries. These correlations persist in the low- and intermediate-density segments under multi-messenger inference but become weak in the core, reflecting the limited constraining power of current observations at high densities. The ability of a single functional form to describe EoSs derived from different microscopic frameworks, together with the recurrence of similar parameter correlations, suggests that diverse dense matter EoSs may share a common underlying structure.

The hypothesis of strange quark matter (SQM) and the possible existence of strange stars have been extensively investigated within the thermodynamic bag model, typically employing the free Fermi gas approximation with or without perturbative QCD corrections at zero temperature, where quark confinement is modeled by the bag pressure. In this work, based on the perturbative inclusion of quark mass effects, we develop a thermodynamically consistent equation of state (EOS) for SQM at finite temperature, incorporating perturbative corrections up to $\mathcal{O}(\alpha_s)$ in the strong coupling constant $\alpha_s$ while ensuring full adherence to the Maxwell relations. These corrections are essential for accurately describing hot quark stars potentially formed during core-collapse events. We present results for both fixed and running couplings, using a phenomenological model to include the breakdown of the perturbative running of $\alpha_s$ at low momenta. Our results demonstrate that incorporating finite-temperature perturbative QCD corrections leads to SQM configurations that fall within the absolute stability window, with equilibrium energies per baryon lying below that of iron, the most tightly bound nucleus. Our EOS supports compact stars with masses exceeding $1.4\,M_\odot$ (and above $2\,M_\odot$ for some values of the bag constant), in agreement with current astrophysical constraints from pulsar and gravitational-wave observations.

We consider collisions between particles moving freely in the Schwarzschild background. We suggest classification of scenarios that lead to unbounded energy output in the center of mass frame. Although this is forbidden in the exterior of a black hole, the phenomenon becomes possible with account of white hole region and mirror universe in the complete space-time diagram. Also, we scrutinize collisions in the vicinity of the bifurcation point and reveal the role of resting observers in the inner region beyond the horizon. In some scenarios $E_{c.m.}$ not only diverges but becomes seemingly infinite that is, as a matter of fact, impossible according to the principle of so-called kinematic censorship. We resolve a corresponding paradox and argue, why this principle remains valid.