Daily arXiv

Ringdown signals from perturbed black holes (BHs) offer a clean window into BH spacetime, strong-field gravity, and fundamental physics. Presently the quasi-normal modes of stellar-mass BH ringdowns have been successfully extracted in the ground-based gravitational wave (GW) observations. Looking ahead, the future space-borne observatories will listen to the ringdowns from massive BH binary coalescences more loudly and resolve multiple modes to unprecedented precision, which calls for efficient approaches to mitigate the sharply increasing computational burden. We develop a practical ringdown analysis pipeline for space-borne detectors by implementing FIREFLY, a novel acceleration algorithm validated in ground-based detectors, and for the first time demonstrate its compatibility and effectiveness with the time-delay interferometry (TDI) observables. With high fidelity, we achieve a $\sim 200$-fold speedup for a simulated ringdown signal including six modes, providing a viable and scalable route for multi-mode ringdown analysis in the space context. This new approach has sound statistical interpretation and is extensible to other GW sources in band.

We present perturbative quantum chromodynamics (pQCD) predictions for high-momentum particle yield modification in very light ion collisions - ${}^{10}\mathrm{B}+{}^{10}\mathrm{B}$, ${}^{6}\mathrm{Li}+{}^{6}\mathrm{Li}$, ${}^{4}\mathrm{He}+{}^{4}\mathrm{He}$, and ${}^{3}\mathrm{He}+{}^{3}\mathrm{He}$ - with and without medium-induced energy loss. We find non-trivial suppression in symmetric systems from ${}^{208}\mathrm{Pb}+{}^{208}\mathrm{Pb}$ to ${}^{3}\mathrm{He}+{}^{3}\mathrm{He}$ and in asymmetric $A+B$ systems, with the suppression scaling approximately as $R_{AB} \simeq (\sqrt{AB})^{1/3}$. Further, we find that ${}^{3}\mathrm{He}$ and ${}^{6}\mathrm{Li}$ offer particularly clean environments for observing final-state partonic energy loss from quark-gluon plasma (QGP) formation in extremely small systems. Finally, we show that energy loss models generically predict $v_2\{\mathrm{SP}\} \approx 0$ in small systems, indicating that the large measured $v_2 > 0$ in $p+{}^{208}\mathrm{Pb}$ is not due to energy loss.

We investigate the properties of dense quark matter and strange quark stars within a nonperturbative, Poincaré-covariant framework. Employing a symmetry-preserving vector$\,\otimes\,$vector contact interaction model, we extend the quark gap equation to the regime of zero temperature and finite quark chemical potential. From the resulting momentum-independent quark propagator, we construct the equation of state (EOS) and solve the Tolman-Oppenheimer-Volkoff (TOV) equations to evaluate the mass-radius relations and tidal deformabilities of strange quark stars. We systematically analyze the sensitivity of the EOS and the macroscopic stellar properties to the model parameters, specifically the effective interaction strength and the ultraviolet cutoff. We demonstrate that reducing the coupling constant stiffens the EOS, whereas increasing the ultraviolet cutoff softens it. By confronting our predictions with multi-messenger astrophysical constraints-including pulsar mass measurements and gravitational-wave data-we identify parameter regimes that successfully describe current observations. Specifically, we find that parameter sets with $α_{ir}=0.735π$, $Λ_{uv}=0.905\,\mathrm{GeV}$ and $α_{ir}=0.588π$, $Λ_{uv}=0.9955\,\mathrm{GeV}$, alongside a vacuum bag pressure of $B \approx (0.106\,\mathrm{GeV})^4$, yield stellar properties in excellent agreement with empirical constraints.

The IceCube Neutrino Observatory has opened a new window into the high-energy Universe, providing measurements of neutrinos over a broad energy range. This contribution presents recent results, including a follow-up on the first identification of a steady neutrino source NGC 1068, measurements of the flavor composition of the diffuse astrophysical flux, limits on prompt atmospheric neutrinos, and searches for neutrinos from dark matter annihilation in the Sun. These measurements probe neutrino production mechanisms, fundamental particle interactions, and physics beyond the Standard Model. Looking forward, the recently deployed IceCube Upgrade will enhance sensitivity to lower-energy neutrinos and reduce systematic uncertainties, while the planned IceCube-Gen2 will expand the detector volume, increase the neutrino detection rate, and extend energy reach, enabling more detailed studies of cosmic sources and high-energy particle physics.

We investigate the dynamics of transonic advective accretion flows around spinning black holes in the presence of viscosity. The spacetime of a Kerr black hole is approximated using a pseudo-potential. We study viscously driven shock oscillations over a range of black hole spin parameters. Our results show that the frequency range of quasi-periodic oscillations (QPOs) obtained from the power density spectra depends strongly on the black hole spin. Low-spin systems predominantly exhibit low-frequency QPOs, whereas rapidly rotating black holes (0.9) produce QPOs spanning a broad range from low to high frequencies, comparable to those observed in black hole X-ray binaries. We further obtain a correlation between the QPO frequency and the power-law photon index by computing the spectrum for a 10 solar mass black hole.

We investigate the thermodynamic properties of static nonlinearly scalarized black holes in Einstein-scalar-Gauss-Bonnet theory with polynomial coupling functions. Based on the scalarized solutions constructed previously, we compute thermodynamical quantities of these scalarized black holes. Moreover, we examine the first law of black hole thermodynamics and consider the phase transitions between Schwarzschild and scalarized black holes. It shows that a phase transition from Schwarzschild black hole to scalarized black hole is a first-order with non-zero latent heat.

Supernovae are one of the most promising gravitational wave sources. But, since the system of the supernovae is nearly spherically symmetric, the expected gravitational waves from them are relatively weak, compared to the case of the compact binary mergers. Thus, at least using the current gravitational wave detectors, only the gravitational waves from a supernova that occurred in our galaxy could be detected. To reliably extract information from gravitational waves originating from such a low event rate, thorough preparation is essential. However, because supernova gravitational waves strongly depend on model parameters, such as progenitor mass and the equation of state for dense matter, it may be difficult to extract physical properties even if the gravitational waves are detected. The universal relations between gravitational-wave signals and physical properties, independent of model parameters, are important for solving this difficulty. To discuss such a universal relation, in this article, we systematically examine the protoneutron-star oscillation frequencies with the linear analysis, the so-called asteroseismology, and compare them with the gravitational wave signals in the simulations.

The gravitational wave event GW231123, with component black hole masses lying within or above the pair-instability mass gap, poses a significant challenge to current stellar evolution models. In this Letter, we investigate its origin by coupling the galaxy formation model GAMESH with the cluster population synthesis code RAPSTER, and with two distinct binary population synthesis codes (SEVN and BSEEMP). This framework allows us, for the first time, to reconstruct the life cycle of GW231123-like candidates within the same cosmological simulation, enabling a self-consistent comparison between different formation channels. We find that, although both population synthesis codes can in principle produce black holes compatible with GW231123, isolated binary evolution fails to reproduce the inferred merger redshift. In SEVN, massive black hole binaries form with semi-major axes > 10^3 Rsun , preventing coalescences within a Hubble time. In BSEEMP, candidates arise only at extremely low metallicities (Z = 10^{-10}), which contribute negligibly to the star formation rate density in our overdense simulated volume. Our results instead strongly support a dynamical, hierarchical origin. The observed black hole masses are naturally reproduced through successive mergers in dense globular clusters. The high dimensionless spins reported by the LIGO-Virgo-KAGRA Collaboration are consistent with this hierarchical population. We find a local merger rate density of 0.78 Gpc^{-3} yr^{-1}, with a peak at z = 4 - 6, tracing the maximum formation rate of globular clusters in metal-poor environments (Z = 0.006). Overall, GW231123 may represent a benchmark event for a robust population of hierarchical black holes formed in the early Universe.

The ground-state properties of superheavy $Z = 122$ isotopes are investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc). Bulk properties, including binding energies, Fermi energies, nucleon separation energies, quadrupole deformations, and root-mean-square radii, are calculated. The results are compared with those obtained from the relativistic continuum Hartree-Bogoliubov (RCHB) theory. By examining the dependence on the angular-momentum cutoff and the effects of triaxial and octupole deformations, a strategy for determining the ground states is suggested. Furthermore, based on an analysis of the Fermi and nucleon separation energies, the proton and neutron drip lines for $Z = 122$ isotopes are determined within both the DRHBc and RCHB frameworks. The possible magic numbers $N=184$, 258, and 350 are also suggested. Finally, the evolution of single-particle levels, deformation, charge and neutron radii as well as average pairing gaps with increasing neutron number, is discussed.

Fast radio bursts (FRBs) are millisecond-duration extragalactic radio transients of unknown origin. Rotation measures (RMs) probe their local magneto-ionic environments and provide important clues to their nature. While RM variability has been observed in several repeating FRBs, it is typically gradual or stochastic. Recently, observations of FRB~20220529 revealed an abrupt RM excursion followed by rapid recovery on week-long timescales, termed an ``RM flare'', suggesting a potentially distinct form of RM variability associated with localized magnetized plasma. In this work, we perform a systematic search for RM flare candidates in repeating FRBs with multi-epoch RM measurements. Using a $3σ$ significance threshold, we identify two candidates with multiple observational epochs (FRB~20121102A and FRB~20201124A) and two additional single-epoch candidates (FRB~20180916B), in addition to the event in FRB~20220529A. Our results suggest that RM flares, if confirmed, may not be rare among repeating FRBs and point to highly dynamic magnetized environments local to the sources. Future high-cadence polarimetric observations, particularly following the discovery of RM excursions, will be essential for confirming these candidates and constraining their physical origin.

Neutrino telescopes have detected astrophysical neutrinos with energies up to ${O}(100)$ PeV. Several current and proposed experiments aim to observe neutrinos at even higher energies, with the goal of detecting cosmogenic neutrinos. This increase in neutrino energy makes tests of Lorentz invariance violation (LIV) particularly appealing, since the effects of higher-dimension LIV operators on neutrino propagation grow rapidly with energy. In this work, we investigate the potential of the upcoming experiments GRAND and POEMMA to probe LIV in the neutrino sector through the detection of ultra-high-energy tau neutrinos. We generate the cosmogenic neutrino flux using SimProp and interface it with a calculation of neutrino flavor transition probabilities in the presence of LIV effects. Deviations from standard flavor transition probabilities manifest as changes in the expected tau neutrino event rates at GRAND and POEMMA. We first consider the case with a single nonzero LIV operator of various dimensions, and find that the projected sensitivities exceed existing limits from lower-energy probes by orders of magnitude. We then explore scenarios with multiple nonzero LIV parameters and show that their interplay can significantly modify the sensitivities compared to the single-parameter case. Overall, we find that upcoming observations of ultra-high-energy tau neutrinos will place some of the most stringent constraints on LIV.

Gamma-ray bursts (GRBs) are generally classified as Type~I GRBs, which originate from compact binary mergers, and Type~II GRBs, which originate from massive collapsars. The traditional correspondence between short--Type~I GRBs and long--Type~II GRBs, separated by a duration of 2 seconds, has been challenged by recent observations of long GRBs associated with kilonovae (i.e., Type~I-L GRBs) and a short GRB associated with a supernova. In this paper, we focus on GRBs with precursor emission (PE) and compile 366 GRBs detected by Fermi/GBM. Applying the unsupervised machine learning methods t-SNE and UMAP, we are able to distinguish Type~I (including subclass Type~I-L) and Type~II GRBs for the first time and identify PE as a key feature for distinguishing GRBs of different origins. Inspired by results of machine learning, we propose a diagnostic parameter, the $E_{\rm p,ME}$-precursor index ($EPI$), defined as ${\rm log_{10}}(E_{\rm p,ME}^{2}/(T_{\rm 100,PE}T_{\rm 100,QE1}^{1/2}T_{\rm MVT,PE}))$, where most Type~I GRBs have $EPI > 6.2$ and most Type~II GRBs have $EPI < 6.2$. This parameter can help the community to diagnose the origin of any GRB with PE based solely on its prompt emission and rapidly plan for follow-up observations. The validation using Swift GRBs provides illustrative evidence that our method may also be applicable to GRBs observed by instruments other than Fermi.

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe-K band on second timescales. Using data from the Performance Verification phase in April 2024, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve's capability to probe plasma environments in high-mass X-ray binaries.

We investigate dense-matter equations of state (EOSs) within a Bayesian framework, with particular emphasis on whether recent small-radius compact-star candidates can be accommodated in a twin-star scenario. For the hadronic sector, we adopt a meta-modeling EOS constrained by the NICER mass--radius measurements of PSR J0030$+$0451, PSR J0437$-$4715, PSR J0614$-$3329, and the massive pulsar PSR J0740$+$6620. The hadronic inference indicates that PSR J0614$-$3329 favors a somewhat softer EOS than the other two \(\sim1.4\,M_\odot\) pulsars, while the \(\sim2\,M_\odot\) constraint prevents the EOS from becoming too soft. We then introduce a strong first-order phase transition through a constant-speed-of-sound quark-matter segment. Using HESS J1731$-$347 and XTE J1814$-$338 to constrain the phase-transition parameters, we find a preferred transition density of \(n_\mathrm{t}\sim2.7\text{--}2.8\,n_0\), a sizable energy-density jump of \(600\text{--}700\) MeV, and a relatively large post-transition sound speed of \(c_s^2/c^2\sim0.85\). Such a phase transition generates a disconnected hybrid branch with radii of about \(6\text{--}7\) km at masses around \(1.2\text{--}1.4\,M_\odot\), and strongly suppresses the dimensionless tidal deformability relative to the purely hadronic branch. This pronounced change in tidal deformability is a characteristic signature of the twin-star mechanism and may provide an important observational tool for identifying phase transitions in neutron-star matter in future multimessenger measurements. These results show that small-radius compact stars can provide direct constraints on both the strength of a first-order phase transition and the stiffness of the post-transition phase in dense matter.

There has been a recent revival in understanding the spontaneous scalarization phenomenon in scalar-tensor gravity as a phase transition. Using the tools of the Landau theory, we now know that first-order transitions where scalarization occurs in a discontinuous manner is more prominent than what had been considered in the literature, and this might lead to novel observation channels. However, the examples so far have been restricted to specific quadratic scalar coupling terms and spherically symmetric stars. Here we explore the phase transition structure of scalarization for more general couplings, considering linear as well as quadratic terms in the conformal scaling factor of the theory. Moreover, we also investigate the effect of rotation on the scalarization phase transition. Both of these considerations are natural choices since the coupling in a scalar-tensor theory can appear at all orders, and astrophysical neutron stars commonly have angular momentum. The introduction of linear coupling leads to a complex solution space which is harder to explore. However, we demonstrate that the Landau model of scalarization enables us to systematically find the branches of scalarized solutions that are commonly overlooked in numerical searches, providing a novel tool. On the other hand, the main effect of stellar rotation is shifting the stellar masses at which the phase transition occurs to higher values, but the qualitative picture remains similar to what happens under spherical symmetry.

By monitoring the times of arrival of radio pulses from millisecond pulsars, Pulsar Timing Arrays (PTAs) serve as unique gravitational wave (GW) laboratories in the nanohertz band. To date, the primary astrophysical sources of GWs targeted in this frequency range have been inspiraling supermassive black hole binaries (SMBHBs) on circular and eccentric orbits. In this work, we demonstrate that, thanks to the so-called pulsar term in the timing residual waveform of GW signals, PTAs can probe individual SMBHBs that merged before timing observations began. We refer to the latter as \emph{zombie binaries}. Using SMBHB population models consistent with current PTA constraints, we find that while the probability of detecting such systems in existing PTA datasets remains low, the Square Kilometer Array observatory is expected to achieve sufficient sensitivity to have a few zombie binaries with a signal-to-noise ratio exceeding 3 in its data. Although their confident identification might be challenging, this new class of PTA sources opens a novel window for studying the most massive SMBHBs in our local universe.

We present multi-wavelength study of the $γ$/X-ray transient EP250416a (also designated GRB 250416C), triggered by the Einstein Probe (EP) Wide-field X-ray Telescope and also by SVOM and Konus-Wind. Observations spanning the gamma-ray, X-ray, and optical bands facilitated detailed analysis of the burst's prompt emission, afterglow evolution, and physical origin. EP250416a exhibits a burst duration of 30 s in X-ray and 17.7 s in gamma-rays, with joint spectral fitting of 0.5-5000 keV data gives $E\rm_{peak}=342_{-232}^{+90}$ keV. Optical spectroscopy of the afterglow, acquired with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South, yielded a redshift of $z=0.963$. Accounting for the measured redshift, the isotropic energies are $E\rm_{X,iso}=2.7_{-0.5}^{+0.9}\times10^{50}$ erg and $E\rm_{γ,iso}=7.34_{-2.1}^{+5.1}\times10^{51}$ erg, aligning with the Amati relation for long GRBs. The fluence ratio $\rm S(25-50~keV)/S(50-100~keV)=0.78_{-0.15}^{+0.1}$ classifies EP250416a as an X-ray rich (XRR) GRB. The X-ray afterglow shows an initial shallow decay ($α\approx -0.5$) transitioning to a canonical decay phase ($α\approx -1$), with a very late jet break at $t\sim 1.5\times 10^6$ s, corresponding to a jet half-opening angle of $θ_j=10.6_{-1.8}^{+1.9}$ degrees. EP250416a is optically dark, as it shows only a faint $r$-band detection ($r=24.16$ mag) from Gemini South-GMOS and a low optical-to-X-ray spectral index $β_{\rm OX} = 0.3$. This may be attributed to significant host-galaxy extinction, with a required $A_V^{\text{host}}=5.5\ \text{mag}$ derived from the extinction curve model.

Pulsar timing arrays record gravitational waves from supermassive black hole binaries at two spacetime points: an Earth term, measured when the wave passes the Earth, and a pulsar term, measured when the wave passed each pulsar at an earlier epoch. We show that a future $μ$Hz-band detection of a nearby massive binary by a mission such as $μ$Ares would turn PTA pulsar terms into targeted probes of binary evolution. In analogy with supernova light echoes, each pulsar term acts as a gravity echo: a dated snapshot of the binary at an earlier stage of its inspiral. Together, the $μ$Hz Earth-term measurement and the nHz pulsar-term echoes provide a temporal baseline that neither detector could access alone. For a fiducial equal-mass binary with total mass $10^9\,M_\odot$ at 80~Mpc, we find a combined pulsar timing array echo signal-to-noise ratio of 33, with up to 24 pulsars individually resolving the signal among pulsars with 50-year baselines. The angular dependence of the single-pulsar echo sensitivity alone enables independent sky localization of the source to $\sim$10--100~deg$^2$, and the resolved pulsar-term frequencies directly measure the binary inspiral rate hundreds to thousands of years ago. With sufficient pulsar distance precision, a small set of anchor pulsars could additionally phase-connect the array and trace the post-Newtonian evolution coherently over kpc baselines. The source population required for gravity echoes is drawn from the same massive-end census responsible for the observed nanoHertz stochastic background.

Around 10 % of hydrogen-poor supernovae explode inside compact ($\sim 10^{15}$ cm), massive ($\sim 0.1 \ \mathrm{M_\odot}$) circumstellar material (CSM), signalling an episode of enhanced pre-explosion mass loss whose mechanism remains unclear. The extreme members of this population are considered to constitute some of the Fast Blue Optical Transients (FBOTs), which exhibit rapid rise times of $\sim$ few days and high peak luminosity $\sim 10^{44} \ \mathrm{erg}$. Recent binary evolution calculations show that the expansion of helium stars during their latest evolutionary stages can trigger a rapid but stable mass-transfer episode that can form a dense circumbinary disc (CBD) that may explain the observed dense CSM. However, a detailed, quantitative analysis of this process and the resulting CBD properties such as its mass, radius and density profile has not yet been undertaken. We present a set of models that solve the viscous evolution of such a CBD under time-dependent mass injection. We find that although the injected mass is initially sub-Keplerian, a lower ``accretion eigenvalue'' $χ$ prevents more mass from falling back onto the central binary. For our fiducial set of models, the CBD immediately prior to the explosion reaches a mass of $0.07-0.20 \ \mathrm{M}_\odot$, a half-mass radius of $640 - 4000 \ \mathrm{R}_\odot$, and an aspect ratio of $θ= H/R \sim 0.1$. We also show that the interaction between SN ejecta and the CBD can power some of the fastest-evolving interacting Type Ibc SNe that can be classified as FBOTs, such as SN 2018gep or SN 2019jc. Despite uncertainties in the model parameters, our results demonstrate that CBD formation triggered by rapid, stable mass transfer is a viable mechanism to explain the dense circumstellar environments observed around rapid, hydrogen-poor interacting SNe. (abridged)

Ultra high energy cosmic rays have been see coming from the direction of the local cosmic void. We use this fact to argue that at least some of these these cosmic rays are relatively light magnetic monopoles and that their relative fraction above 1020 eV can be found from full sky observations.

Complex astrophysical processes regulate the growth of galaxies by injecting energy and momentum into their surroundings, redistributing baryons across megaparsec scales. The clustering of matter on these scales, as measured via weak lensing and galaxy surveys, encodes critical cosmological information on the dynamical dark energy, the nature of dark matter and the sum of neutrino masses. The suppression of matter clustering due to feedback processes limits the interpretation of cosmological measurements. Multiple probes of the baryon distribution have attempted to quantify the strength of feedback via measurements of suppression in the matter power spectrum. The dispersion measures (DMs) of fast radio bursts (FRBs) have emerged as a powerful new probe of baryons, with the advantage over other probes of being unbiased with respect to density and temperature. Here, we use a sample of 109 FRBs with redshifts and DMs to directly measure the spatial fluctuations in the baryon density field, quantifying the effects of feedback on the matter power spectrum at scales of $k \sim 0.1-3$ h/Mpc, and the gas fraction in galaxy groups and clusters ($10^{13}-10^{15} M_\odot$). We use a halo-model prescription to conduct inference, and find that FRB data reduces the posterior variance at k $\sim$ 1 h/Mpc by a factor of $\sim 8$ relative to the prior. The statistical precision of inferred FRB constraints is similar to other baryon tracers, while probing a complementary redshift regime ($z \lesssim 0.3$). A comparison with several hydrodynamical simulations excludes extreme large-scale feedback scenarios at $\sim 2σ$ confidence. This work establishes FRBs as a sensitive probe of feedback-regulated structure formation. As next-generation experiments deliver orders-of-magnitude larger samples, FRBs are poised to drive the constraints on baryonic physics in the era of precision cosmology.

Astrophysical dark matter particles with masses well below GeV-scale can be difficult to detect using conventional nuclear recoil experiments due to their low velocities in our Milky Way halo. Elastic scattering with high-energy cosmic rays or thermal production inside core-collapse supernovae can accelerate sub-GeV DM to (semi-)relativistic velocities, producing nuclear recoil energies above the keV threshold that paleo-detectors can record over geological timescales. Using olivine as the target with 100$\,$g$\cdot$Gyr exposure, we compute track length distributions from such (semi-)relativistic dark matter fluxes, incorporating all major backgrounds (neutrinos, uranium-chain neutrons, thorium recoils) with a statistical analysis on an Asimov dataset. We derive 95 C.L. projected sensitivity of paleo-detectors to the DM-nucleon cross section for dark matter masses between a few MeV and hundreds of MeV. Our results show that paleo-detectors are able to probe large parameter regions that are not covered by current and near-future experiments designed to detect dark matter and neutrinos. In particular, paleo-detectors offer a unique ability to record the dark matter flux from Galactic supernova events over geological times. Such cumulative exposure enables sensitivity gains of a few orders of magnitude compared to conventional experiments.

Little Red Dots (LRDs) have emerged as one of the central puzzles of the JWST era. Their spectra increasingly require dense gas close to the source, yet the physical origin of that cocoon-like structure remains unclear. We examine whether late pulsational mass loss from supermassive stars (SMS)leads to dense gas cocoons. We analyze five accreting GENEC models at different metallicities with characteristic masses of order $10^5\,M_\odot$, following them through post-accretion evolution with radial pulsation calculations and general relativistic (GR) stability diagnostics. Mass loss during the final stages of evolution occurs not as a steady wind, but through discrete strange-mode ejection episodes. In the $Z=10^{-2}\,Z_\odot$ model, which provides the clearest LRD analogue, four late episodes last $41$--$282$ yr and eject $10$--$348\,M_\odot$ each, for a total loss of $(4.8-10)\times10^2\,M_\odot$; the final episode alone contributes $\simeq 73\%$ of that budget. Since the last episode dominates the mass-loss, it is the only event sufficiently massive enough to leave behind a compact, optically thick shell extending out to 0.4 pc that reproduces the LRD dense gas cocoon. The final ejecta are H/He dominated but chemically distinctive, with a robust nitrogen-rich composition, $\log(\mathrm{N/O})\simeq0.13$ and $\log(\mathrm{C/O})\simeq-0.23$. The SMS reaches GR instability at an age of $\sim 1$ Myr and collapses in $\sim10^4$ s, retaining $\sim 99\%$ all of its mass. Across the full metallicity range from Pop III to $10^{-2}\,Z_\odot$, this shell-ejection channel persists. Pulsational mass-loss from SMSs therefore provides a physically motivated origin for the compact cocoon-like structure implied by LRDs, while remaining the natural progenitors of the massive black hole seeds invoked in direct collapse scenario.

X-ray polarimetry is a new tool capable of probing the geometry of accretion onto weakly magnetized neutron stars. Here we present the first X-ray spectropolarimetric results from coordinated observations of the atoll source 4U 1735-44, conducted with the Imaging X-ray Polarimetry Explorer (IXPE), NICER, and NuSTAR. Over the 2-8 keV energy range, we obtained a marginal detection of polarization with the polarization degree of $1.4\%\pm0.7\%$ and polarization angle of $-29^\circ\pm14^\circ$, corresponding to a $3σ$ upper limit on the polarization degree of 3.5\%. The best-fit model to describe the spectrum comprises a thermal component associated with the accretion disk, a Comptonized blackbody component, and a relativistic reflection component. From the reflection model, we infer a disk inclination of $\sim 40^\circ$. The spectroscopic and polarimetric properties of 4U 1735-44 are consistent with those observed in other atoll sources studied by IXPE, with its low polarization likely due to its low inclination.

We propose a cosmological scenario in which, beyond matter and radiation, an additional barotropic fluid with positive equation of state $ω_s$ contributes to the cosmic energy budget, in contrast to Early Dark Energy (EDE). We investigate the theoretical implications of this framework, here dubbed the $Λ_{ω_s}$CDM model, at both the background and perturbative levels, exploring its impact on the expansion history and structure formation. We show that, while remaining subdominant at late times and therefore consistent with current observational bounds, the additional fluid modifies the early-time expansion rate, leading to a higher inferred value of the Hubble constant. Thus, we perform a full Bayesian analysis using a modified version of the \texttt{CLASS} Boltzmann code interfaced with \texttt{MontePython}, considering combinations of \textit{Planck} 2018 Cosmic Microwave Background (CMB) data, DESI DR2 Baryon Acoustic Oscillations (BAO) measurements, Pantheon Type Ia supernovae (SNe Ia), and SH0ES determinations of $H_0$. We find that the inclusion of the SH0ES prior, $H_0 = 73.04 \pm 1.04\,\mathrm{km/s/Mpc}$, leads to a preference for a nonvanishing barotropic fluid. In particular, we obtain $ω_s = 0.290^{+0.017(0.021)}_{-0.007(0.028)}$ and density $10^5Ω_s = 1.47^{+0.35(1.14)}_{-0.62(0.94)}$ for the dataset combination CMB + BAO + Pantheon + SH0ES, and $ω_s = 0.302^{+0.024(0.034)}_{-0.013(0.038)}$ and $10^5Ω_s = 1.21^{+0.31(1.10)}_{-0.65(0.86)}$ when BAO data are excluded. We further compare our scenario with the EDE framework and show that, statistically, no strong evidence is found against the $Λ_{ω_s}$CDM model. Finally, we provide a physical interpretation of our fluid in terms of matter with pressure, indicating that the standard cosmological model may be incomplete in its current minimal formulation.

We compute the abundance of cosmological gravitational waves produced during the evolution of an ultralight vector (spin-1) dark matter field. A homogeneous background vector field breaks spatial isotropy, requiring a Bianchi I geometry and inducing a mixing between the scalar, vector, and tensor perturbation sectors. We derive the perturbation equations in this background and show that, as a consequence of this mixing, scalar perturbations act as a source of tensor modes, generating a stochastic GW background. The production and cosmological evolution of these gravitational waves are implemented in a modified version of CLASS, from which we obtain the present-day spectrum.

We present formulations of effective potentials suitable for approximating general relativistic effects in Newtonian simulations of core-collapse supernovae. Assuming a spherically symmetric spacetime and a stress-energy tensor that includes both fluid and neutrino contributions, Eulerian and Lagrangian projections of the Einstein equations are made to determine general relativistic corrections to the Newtonian gravitational potential. We implement the effective potentials in both the Chimera and Flash-X codes, and perform a series of adiabatic and core collapse simulations. The results are compared to Newtonian and fully general relativistic simulations, as well as another widely used effective potential formulation. We find close agreement between our new effective potentials and the fully general relativistic results from multiple other codes.

The impact of Coulomb core polarization on the neutron density distribution of $^{48}\mathrm{Ca}$ is investigated, revealing that Coulomb-driven proton shifts toward the nuclear surface lead to neutron skin impurity, defined as the presence of protons in the neutron skin region. The Coulomb boundary radius, where protons begin to be driven outward, is located before the neutron skin region, leading to proton admixture and establishing impurity as an intrinsic feature of the density profile in $^{48}\mathrm{Ca}$. Using the $(^3\mathrm{He},t)$ isobaric analog state (IAS) reaction at 420 MeV, we quantify impurity-induced modifications of the probed neutron distribution. For a neutron skin thickness of 0.144 fm, consistent with the thin-skin scenario suggested by CREX ($0.121 \pm 0.035$ fm), Coulomb core polarization effectively enhances the neutron distribution by approximately 0.052 fm via the symmetry potential response to proton admixture. Our analysis shows that neutron skin impurity is thickness-dependent and becomes suppressed in large-skin scenarios as the symmetry potential neutralizes core polarization, thereby motivating IAS charge-exchange measurements as a complementary probe of neutron skin observables.

Einstein's Weak Equivalence Principle (WEP), the universality of free fall, is a fundamental component of general relativity and other metric theories of gravity. Its validity can be tested through the post-Newtonian parameter gamma, which quantifies the amount of spacetime curvature due to the presence of unit rest mass. In this paper, we use high-energy neutrino events detected by IceCube and associated with the gamma-ray blazars TXS 0506+056 and PKS 0735+178 to test the WEP via the Shapiro delay induced by the gravitational potential of Laniakea. We find that violation of the equivalence principle for neutrinos and photons is limited to an accuracy of 10^-6, 10^-7 and 10^-8, representing improvements of one, two, and three orders of magnitude, respectively, over previous constraints obtained from other high-energy neutrino-blazar associations and up to six orders of magnitude tighter compared to the constraints obtained with MeV neutrinos from SN1987A.

Long period radio transients (LPTs) are the slowest radio-pulsing sources ever found, with the current population spanning periods of seven minutes to over six hours. Two of the thirteen published LPTs, ILT J1101+5521 and GLEAM-X J0704--37, have been associated with an M dwarf closely orbiting a white dwarf (WD) through optical spectroscopy. Here, we present new Keck I/LRIS optical spectroscopy of ILT J1101+5521, which reveals H$α$ emission from the M dwarf and confirms an orbital period nearly matching the radio period (2.092 hr). Radio pulses in both systems arrive just after maximum M dwarf redshift, assuming the radio period matches the orbital period. Based on Gaia proper motions and systemic velocities, we find that these systems are kinematically hotter and less concentrated in the Galactic plane than other LPTs. Both systems harbor unusually massive and cool WDs, with $M_\mathrm{WD} \approx 0.84-1.0 M_\odot$ and $T_\mathrm{eff} \approx 5200-7300$ K, implying that their carbon-oxygen cores are nearly entirely crystallized. Both systems are unusually close to being face-on binaries ($i=13^\circ-28^\circ$), signaling that the production of coherent radio pulses may be a strongly inclination-dependent phenomenon. We present MESA models that show that the M dwarf in each system will fill its Roche lobe within $\sim1$ Gyr, becoming a cataclysmic variable. Finally, we place lower limits on the space density of WD + M dwarf LPTs ($ρ\gtrsim 10^{-8}\;\mathrm{pc}^{-3}$); based on the broader population of WD + M dwarf binaries, we estimate that there are 100 (2000) WD + M dwarf LPTs within 2 kpc if current radio findings are 100% (10%) complete. Current and upcoming radio surveys will be sensitive to many such systems, and M dwarf optical counterparts out to $\sim$2 kpc will be detectable with the Rubin Observatory Legacy Survey of Space and Time (LSST).

We study the gravitational collapse of a non-interacting mix of perfect fluid and a spatially homogeneous scalar field within a Chiellini-integrable framework. We choose an extended Higgs-type self-interaction potential and reduce the Klein-Gordon equation into a generalized damped Milne-Pinney class of differential equation. We derive a closed-form analytical solution for the scalar field, the scale factor and explore the collapsing branch of the same. We find that it exhibits an asymptotic collapse in which the proper volume decreases monotonically but never reaches zero at finite time. We analyze the energy conditions for the constituent elements of the collapsing sphere. While the scalar field remains canonical in nature, we find that the perfect fluid can violated the Null Energy Condition. We also study the formation of apparent horizon condition and find multiple possibilities depending on the parameter space : either no trapped surface or the formation of multiple apparent horizons. We match the interior homogeneous solution to a generalized Vaidya exterior via the Israel-Darmois junction conditions, yielding the corresponding boundary mass function, ensuring a smooth collapse scenario.

We model the multi-band lightcurves of 80 SNe Ic-BL, including 11 associated with lGRBs, using a magnetar engine model with $^{56}$Ni decay. We find that the data are all consistent with a magnetar central engine, and such a model yields high-quality fits across the sample. The medians with $1σ$ regions of the key parameters are $P_{\rm{i}}\sim2.04^{+1.84}_{-0.96}\,{\rm{ms}}$, $B_{\rm{p}}\sim3.96^{+3.28}_{-1.40}\times10^{15}\,{\rm{G}}$, $M_{\rm{ej}}\sim2.30^{+1.48}_{-1.02}\,M_\odot$, and $M_{\rm{Ni}}\sim0.18^{+0.14}_{-0.09}\,M_\odot$, with strong and statistically significant correlations observed for both $M_{\rm{ej}}-P_{\rm{i}}$ (anti-correlation) and $M_{\rm{Ni}}-M_{\rm{ej}}$ (correlation). Comparing the SN Ic-BL samples with and without lGRB association using fitting parameters, we find no significant difference between them, although the GRB-associated sample is slightly brighter, possibly due to an observational bias. Relative to ordinary SNe Ic, SNe Ic-BL have similar $^{56}$Ni and ejecta masses, suggesting comparable pre-SN progenitor properties, with differences possibly arising from the presence of a magnetar engine. In comparison with other possible magnetar-powered SESNe, including SLSNe Ic and FBOTs, we confirm a strong universal $M_{\rm{ej}}-P_{\rm{i}}$ correlation, indicating a common origin. SNe Ic-BL and SLSNe Ic have similar ejecta mass distributions, typically $M_{\rm ej}\gtrsim0.5\,M_\odot$, while FBOTs mostly lie below this value. Differences between SNe Ic-BL and SLSNe Ic may arise from magnetar properties, with SN Ic-BL magnetars rotating faster and having stronger fields. Moreover, the $P_{\rm{i}}-B_{\rm{p}}$ distribution of lGRB magnetars largely overlaps with that of SN Ic-BL magnetars. In connection with binary simulation results, we propose a unified physical classification and progenitor framework for magnetar-powered and ordinary SESNe.

X-ray flashes (XRFs) are a type of gamma-ray bursts (GRBs) with prompt emission predominantly below 30 keV poorly detected by previous missions. The advent of the SVOM mission, with its wide-field instrument ECLAIRs, provides a new way to detect soft X-ray transients such as XRFs. We present photometric and spectroscopic observations of XRF 241001A detected by SVOM, a soft, sub-luminous, and low-energetic burst located in a poorly populated region of the Amati relation. We investigate the origin of its faint, soft high-energy emission to assess its connection to the long GRB population. We analyze the SVOM/ECLAIRs prompt emission and model its afterglow emission from X-ray to-radio. We present JWST/NIRSpec and SVOM/VT observations of the associated supernova (SN 2024aiiq), which we model with an Arnett radioactive decay component and compare its properties with previously detected GRB/SNe. XRF 241001A is located at z = 0.573 and has a prompt emission dominated by photons below 20 keV with a duration of T90 = 3.14 seconds. Its spectrum can be modeled by non-thermal or thermal models, all pointing towards a low Epeak < 10 keV and Eiso ~ 8x10^49 erg. The X-ray-to-radio afterglow modeling favors an origin from a relativistic jet viewed on-axis. In the optical, XRF 241001A exhibits an early blue emission, similar to that detected in some fast X-ray transients and inconsistent with synchrotron emission. The JWST/NIRSpec observations firmly established its collapsar origin by revealing a SN Type Ic with broad lines, comparable to SN 1998bw and SN 2025kg-like events. XRF 241001A is a soft, low-luminosity collapsar event produced by a weak relativistic jet observed on-axis, supporting the view that part of the XRF population forms the low-energy tail of the long GRB population. It demonstrates the potential of SVOM/ECLAIRs to probe the soft regime of the high-energy transient population.

General relativity is a background-independent theory of a dynamical classical spacetime geometry. Quantum theory is formulated in a classical spacetime, as an intrinsically probabilistic, contextual theory of non-classical, interfering probabilities, with a fixed Born rule for computing those probabilities. We argue that the quantum nature of spacetime, which includes a non-commutative dual companion to the (observed) classical spacetime, is the reason behind an intrinsically probabilistic and contextual nature of quantum theory, with the fixed Born rule. In quantum gravity, we claim, quantum theory is gravitized into a background-independent structure with dynamical and contextual quantum probabilities. This proposal implies intrinsic triple and higher-order interference in the context of massive quantum probes, which sheds light on string theory and the observed vacuum energy as well as the masses of elementary particles.

We investigate timelike geodesics in asymptotically flat regular black holes supported by a phantom scalar field characterized by a scalar charge $A$. This parameter removes the central singularity and continuously deforms the Schwarzschild geometry while preserving asymptotic flatness. We derive the equations of motion for massive test particles and classify bounded and unbounded trajectories in terms of the conserved energy and angular momentum. We determine circular and critical orbits, including the innermost stable circular orbit (ISCO), and analyze the transition between capture and scattering. We show that the scalar charge modifies the location of the unstable and stable circular orbits, the ISCO, and the threshold angular momentum for scattering, exhibiting a nontrivial dependence on the radial coordinate. Their physical scales are naturally described in terms of the invariant areal radius $R(r)=\sqrt{r^2+A^2}$. In the weak-field regime, we compute the perihelion precession and obtain corrections proportional to the scalar charge, allowing us to constrain the scalar charge from Solar System observations. We also analyze the motion with vanishing angular momentum and show that, while the qualitative structure of the trajectories remains connected to the Schwarzschild limit $A\to 0$, the quantitative deviations encode the geometric effects of the scalar hair.

High energy neutrinos can be injected in the early Universe from the decay or annihilation of long lived primordial relics. We analyse the possibility that the ultrahigh energy neutrino event recently observed by the KM3NeT neutrino telescope could have such an origin. This possibility has the advantage of leading to a sharp spectral feature in a way that the neutrino flux can be small at all energies except at the KM3NeT event energy. Thus, along this scenario the tension with null results from other experiments is reduced with respect to the usual power law case analysed by the KM3NeT and IceCube experiments. At such energies and for an emission around the recombination time, interactions of these neutrinos with background neutrinos prove to be relevant and must be determined from the development of a dedicated code. These interactions, as well as final state radiation processes, modify the spectrum. Interestingly, it turns out that the scenario can also leave an imprint in the CMB that could be probed in the near future. Interestingly too, this scenario does not predict an associated $γ$-ray flux beyond observation. All in all we do find that the high energy neutrino could be a primordial high energy neutrino, provided it has been produced around the recombination time or later.

The interaction between a Nambu-Goto string and a Kerr black hole gives one of the methods of energy extraction from a rotating black hole. Although the properties of such processes have been well studied for rigidly rotating strings, little is known for non-rigidly rotating strings. In this paper, we study time evolution of a Nambu-Goto string on the equatorial plane of a Kerr spacetime, which sticks on the horizon and extends to spatial infinity. The time evolution is studied by the series expansion with respect to $t$ and the numerical simulations, which give reliable results for $t\lesssim 4M$ and $t\lesssim 38M$, respectively, where $M$ is the black hole mass. Since the angular velocity of the string on the horizon must coincide with the horizon angular velocity to keep the timelike property, the string is dragged into rotation and coils around the black hole. The negative energy is observed to fall into the black hole, but the positive energy follows after that, meaning that the energy extraction occurs for a short period of time. In the outside region, a wave is generated and propagates to the distant region carrying the extracted energy. After the propagation of the wave, the system approaches the time-independent configuration found by Boos and Frolov, and the total extracted energy is estimated as $E_{\rm ext}\lesssim μM$, where $μ$ is the tension of the string.

The paper develops a model to understand the effective quantum geometry of a black hole horizon and the emission of Hawking spectrum in 2+1 dimensions. We argue that one may view the black hole horizon as formed out of quantised lengths of elementary quanta of value $8π\ell_{P}\, n$, where $n\in \mathbb{N}$, and $\ell_{P}$ is the Planck length. To an observer near the black hole horizon, the entropy (or length of horizon cross-section) is related to the black hole energy. Hence, one may develop a formulation of length ensemble (similar to the area canonical ensemble of Krasnov) from which the black body spectrum may be obtained directly. To this local observer, the temperature of the Hawking spectrum is modified due to the Tolman factor.

While massive black hole binaries (MBHBs) merge at gravitational-wave frequencies above the pulsar timing array (PTA) sensitivity band, we show that they leave orphaned low-frequency contributions in the PTA pulsar term. Due to the light-propagation time between each pulsar in the array and Earth, the pulsar term acts as a time-delayed probe of a chirping merger with a specific frequency response determined by the direction of origin and intrinsic properties of the MBHB. We provide a detailed consideration of how such a multiband signal would manifest in a full PTA, demonstrate an approach to stack these orphaned pulsar terms across the array, and discuss prospects for an archival, multiband search in conjunction with MBHB mergers observed in astrometric data or spaceborne interferometers like LISA.

SDSSJ110546.07+145202.4 stands out as a unique radio changing-look Narrow-line Seyfert 1 (NLS1) galaxy that has brightened dramatically and shows an exceptionally long duration of its "on" phase. We present the first high-frequency radio observations, the first simultaneous radio spectral energy distributions (SEDs), the first optical--UV--X-ray SEDs, and the first X-ray monitoring and spectroscopy of this recently discovered event. Importantly for understanding the nature of the outburst, we show that the X-ray spectrum is soft with a photon index Gamma_X=2.5; line-of-sight absorption and extinction are low or absent; the radio SED is peaked at low frequencies ~2 GHz; and the radio outburst emission is very long-lived (t > 8 yr) and roughly constant. The softness of the X-ray spectrum, low supermassive black hole (SMBH) mass, and high Eddington ratio all corroborate the optical NLS1 classification. We discuss multiple outburst scenarios, including lensing, absorption, a binary SMBH merger, a long-duration giant-star tidal disruption, a newly ignited active galactic nucleus (AGN), and an accretion-rate change. While most of them can be either excluded or are deemed too rare and lack positive evidence so far, most or all types of these transients are expected to be detected in ongoing VLA and upcoming SKA surveys. SDSSJ110546.07+145202.4 itself is well explained by an accretion rate change that triggered the powerful radio jet emission. The low redshift and SMBH mass of this system offer a unique perspective of the physical processes of radio-jet ignition that are expected to operate in the early Universe around growing SMBHs.

We investigate the shadow, timelike geodesic structure, radiation properties of thin accretion disks, and optical appearance of a static spherically symmetric regular black hole, constructed based on the Dehnen-type density profile. Using observational data from M87* and Sgr A*, we constrain the model parameter $a$ at both $1σ$ and $2σ$ confidence levels. Based on the Page--Thorne model, we calculate the local radiative flux, redshift factor distribution, and the radiation flux received by a distant observer, systematically examining the effects of the parameter $a$ and the viewing angle on the black hole image. The results show that larger $a$ will enlarge the effective radiation area of the accretion disk and significantly enhance the asymmetry and Doppler boosting effects of the direct and secondary images at large viewing angles.

In a search for low surface brightness radio nebulae using the ASKAP-EMU survey, we discovered a faint radio shell, G289.6+5.8, and its central point radio source at the position of the soft gamma-ray source IGR J11187-5438. The central radio source is spatially coincident with a previously known low-mass X-ray binary (LMXB) with an M-type donor star. However, the newly determined Gaia DR3 distance of 267 pc and correspondingly low X-ray luminosity (3 x 10e31 erg/s) cast doubt on the LMXB classification. Neither radio nor X-ray pulsations are detected. Chance-alignments between radio shell, central radio source, optical star, gamma-ray, and X-ray sources appear unlikely. By combining all currently available evidence, we conclude that G289.6+5.8 is a remnant of a low-energy core-collapse explosion of an intermediate mass star (~8Msun) in a binary system with an M-type secondary, which remained bound after the explosion. In this scenario, G289.6+5.8 is a supernova remnant, while the central gamma- and X-ray source is associated with a young neutron star driving a pulsar wind interacting with its M-type stellar companion, making IGR J11187-5438 a nascent spider-type X-ray binary.

We use the two-flavor Linear Sigma Model with quarks as an effective description of QCD to investigate the nature of the chiral phase transition at finite baryon chemical potential and zero temperature. We work at one-loop order to set up and solve the system of self-consistent coupled equations for the particle pole masses. The chemical potential-dependent value of the chiral order parameter is obtained by minimizing the one-loop effective potential. This treatment goes beyond the conventional ring-diagram approximation and provides a description valid for arbitrary values of the chemical potential. We find that the phase transition is of first order, and occurs when the quark chemical potential reaches the value of the vacuum quark mass for the chosen set of parameters. The first order nature of the transition is signaled by the discontinuous behavior of the chiral condensate, the masses and the couplings. The thermodynamics of the system is readily implemented and in particular, we find that the square of the speed of sound exhibits a discontinuity at the phase transition and then smoothly approaches the conformal limit from below.

The acceleration of cosmic rays up to PeV energies at supernova remnant shocks requires an amplification of the ambient magnetic field. The amplification mechanism must operate upstream of the shock, to prevent the escape of particles from the system. Observational evidence of field amplification has been indeed obtained by means of X-ray observations. However, such observations constrain the magnetic field strength downstream of the shock only. Here we describe a mechanism for magnetic field amplification that operates downstream of the shock. It is based on a plasma instability triggered by the drift of charged interstellar dust grains overtaken by the shock. We compute the growth rate of the instability, we estimate the level of magnetic field amplification expected downstream of supernova remnant shocks, and we compare our results with observations. In some cases (most notably Cas~A) this mechanism might explain the presence of the X-ray filaments observed at supernova remnant shocks, without requiring any amplification of the magnetic field upstream of the shock and therefore no acceleration of CRs to ultra-high energies.

We investigate electron--positron pair creation through the Schwinger process in newborn magnetars with millisecond spin periods and surface dipole fields close to or above the QED critical field, $B_{\rm Q} = 4.414\times10^{13}\,\mathrm{G}$. In the unscreened field scenario, we derive the analytical global pair creation flux and recast it into a compact form with accurate analytic approximations. For a fiducial model with $B_{\rm p} = 10^{14}\,\mathrm{G}$ and $P_0 = 1\,\mathrm{ms}$, the Schwinger channel exceeds the classical Goldreich--Julian particle supply by many orders of magnitude and becomes the dominant source of charges at the earliest stage of the magnetar. The associated discharge removes about $90\%$ of the initial rotational energy within 30 ms, suppresses the gravitational-wave loss channel, and implies that the observable millisecond phase is extremely short in this unscreened scenario. The rapid energy release over such a short timescale may also provide a viable power source for astrophysical transients. Extending the same fiducial model to $10^4\,\mathrm{yr}$ gives spin periods of order seconds, linking newborn millisecond magnetars to the mature magnetar population.

We present a systematic study of particle production in $Ne+Ne$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.36$ TeV using the A Multi-Phase Transport (AMPT) model with string melting (SM) configuration. The analysis compares spherical and deformed configurations of ${}^{20}\mathrm{Ne}$ to investigate the influence of initial-state nuclear deformation on bulk observables. Charged-particle pseudorapidity ($\langle dN_{\mathrm{ch}}/dη\rangle$) densities, identified particle yields ($dN/dy$), transverse momentum ($p_T$) spectra, mean transverse momentum ($\langle p_{\mathrm{T}} \rangle$), and $p_{\mathrm{T}}$-differential particle ratios ($K/π$ and $p/π$) are studied as functions of multiplicity and centrality. The results show that all observables exhibit the expected dependence on event activity, including smooth multiplicity scaling, mass ordering in $\langle p_{\mathrm{T}} \rangle$, and characteristic features associated with radial flow and quark coalescence. Differences between the two configurations on bulk observables remain small across all observables, typically at the level of a 2\%--6\% percent, with slightly enhanced sensitivity observed in peripheral collisions. These findings suggest that, within the AMPT-SM framework, the collective dynamics and hadrochemical composition are primarily governed by the overall system density and interaction dynamics, while the influence of initial-state deformation is subleading. This study provides a baseline for understanding deformation effects in light-ion collision systems and highlights the limited sensitivity of bulk observables to initial nuclear geometry in transport-based approaches.

We investigate the quantum properties of gravitational waves (GWs) generated by a classical energy-momentum tensor. Treating the GW field as a quantum field coupled to a classical source, we evaluate the expectation value of the GW operator. We demonstrate that this expectation value exactly reproduces the classical retarded solution. Furthermore, we show that the mean and variance of the number of emitted gravitons are equal. This suggests that the graviton emission is a Poisson process, as expected for a coherent state. We establish a quantitative criterion for the validity of the classical wave description. By applying this criterion, we find that the classical approximation is remarkably accurate for astrophysical sources, but laboratory-scale systems may reside in a regime where the discrete nature of graviton emission becomes significant.

The recent detection of a stochastic gravitational wave background by pulsar timing arrays has opened a new window in understanding supermassive black hole binaries and in probing the universe at the early time. Recently, pulsar timing array (PTA) collaborations have been further paving the way to probe anisotropies in the stochastic gravitational wave background. This study investigates dipole-type statistical anisotropy in the primordial power spectrum within a phenomenological framework. We demonstrate that the primordial dipole induces both dipolar and quadrupolar anisotropies in the energy density spectrum of scalar-induced gravitational waves (SIGWs), without generating extra polarization modes. Based on this anisotropic spectrum, we derive the corresponding PTA overlap reduction functions (ORFs), which exhibit frequency dependence, with the anisotropies enhanced on small scales. Furthermore, owing to the non-uniform distribution of millisecond pulsars over the sky in current PTA dataset, the ORFs exhibit a morphology that explicitly depends on the preferred direction of the anisotropy. However, our bayesian analysis of the NANOGrav 15-year dataset still yields no significant evidence for a preferred direction and a weak upper limit on anisotropy amplitude $(g\lesssim0.5)$. This result arises because the observational frequency band lies below the spectral peak, where our models predict suppressed anisotropic contributions. This limitation highlights the potential of future PTA observations. Specifically, datasets with broader frequency coverage are expected to tighten constraints on dipole-type anisotropy.

We present new regular black hole solutions in general relativity (GR) within a static, spherically symmetric framework governed by a variable equation of state, following the approach of [Class. Quant. Grav. 42, 025024 (2025)]. The matter supporting these geometries is identified as a purely magnetic monopole configuration of the Maxwell-Faraday tensor in the context of nonlinear electrodynamics (NLED). We explicitly reconstruct the corresponding NLED Lagrangian and analyze the asymptotic and central behaviors of the solutions. The geometric structure is examined through the metric functions, the regularity of the Kretschmann scalar, and the profiles of energy density and pressures, including a discussion of the resulting energy conditions. Using Event Horizon Telescope observations of Sgr A$^*$, we constrain the model parameters by comparing the predicted size of the central dark region with the inferred observational images, taking into account the effective geometry experienced by photons in the presence of NLED. Finally, we investigate the dynamical stability of these configurations under scalar perturbations by computing the quasinormal mode spectrum and performing a time-domain analysis.

G310.7-5.4 is a supernova remnant (SNR) candidate identified as a faint shell in the second epoch Molonglo Galactic Plane Survey (MGPS-2), but this has not been followed up with multi-wavelength observations until now. It is an example of an SNR at high Galactic latitude showing spatially coinciding $γ$-ray emission. Here, we make the first detailed investigation of the radio emission from the G310.7-5.4 region, aiming to characterise the radio structure, polarisation measurements and the coinciding GeV emission. We used recent radio continuum observations at 943.5 MHz from the EMU and the POSSUM surveys with ASKAP, as well as 16.5 years of Fermi-LAT observations. We furthermore considered the multiwavelength context of the object by investigating observations previously conducted with other instruments, such as infrared and X-ray surveys. We confirm the SNR candidate as a new supernova remnant, dubbed Abeona. We detect the presence of a faint, extended, bilateral radio shell of the size of around 30' diameter and ASKAP radio flux density of $1.5^{+1.5}_{-0.1}$ Jy with no obvious infrared counterparts. With a radio surface brightness of about $2.4^{+2.4}_{-0.1}\times10^{-22}$ W m$^{-2}$ Hz$^{-1}$ sr$^{-1}$, this SNR is one of the faintest radio SNRs known. The northern part of the shell shows linearly polarised radio emission, characteristic of synchrotron emission in SNRs. The physical size of the SNR is estimated to be around $42^{+42}_{-21}$ pc, which would give a distance of around $4.9^{+4.9}_{-2.5}$ kpc. Furthermore, the spatially coincident $γ$-ray source 4FGL J1413.9-6705 shows an energy flux of $1.26\pm0.35\times 10^{-6}$ MeV cm$^{-2}$ s$^{-1}$ with a significance of 5.7 $σ$ between 100 MeV and 100 GeV. The SNR is also put in context with known high-latitude SNRs with $γ$-ray counterparts and compared with their observational properties.

When an accretion flow extends to the event horizon, their intersection defines the contour of the inner shadow. However, the morphological evolution of this critical feature remains largely unexplored within a torn accretion disk system, a configuration comprising distinct sub-disks formed when a tilted disk is disrupted by frame-dragging. To address this, we phenomenologically construct a torn accretion disk model and numerically simulate the inner shadow of a Kerr black hole using relativistic backward ray-tracing. We discover that the torn disk geometry profoundly alters the black hole's observational signatures, inducing severe erosion of the inner shadow and generating novel features such as bifurcated shadows, crescent-like structures, and multiple orders of shadow rings. These exotic morphologies, which are predominantly governed by the spatial discontinuity between the sub-disks and the tilt angle of the outer sub-disk, are exceedingly difficult to replicate within standard equatorial accretion paradigms. Our findings demonstrate that these distinctive shadow structures hold significant potential to serve as robust diagnostic probes for torn accretion environments, simultaneously implying that relying solely on the inner shadow to test gravity theories is fundamentally insufficient.

Massive stars commonly form binaries that can evolve into compact systems via common envelope evolution (CEE), a critical but poorly understood phase -- especially when the companion is a neutron star. Understanding the drag force exerted on a neutron star during CEE is a key to the quantitative evaluation of orbital decay, merger timescale, and compactness of the resultant binary. In this paper, we conduct general-relativistic hydrodynamical simulations under a novel strategy of multi-layer domain-decomposition to treat the vast disparity of $10^4$--$10^7$ between the neutron star radius and the accretion radius. Our 10-model survey spans diverse physical conditions that the neutron star encounters in the envelope of a massive star. We find that nested bow shocks with alternating orientations commonly form. This configuration is qualitatively different from those in the conventional picture and results in an enhancement of the drag force by one to two orders of magnitude from what the Bondi--Hoyle--Lyttleton formula predicts. Moreover, the direction of the net force can reverse depending on the envelope conditions, contrary to the standard picture in which the drag always decelerates the companion. These results will serve as a basis for improvements of the drag force prescription in CEE modeling, and have implications for binary evolution theory.

The recent O4a release from the LIGO-Virgo-KAGRA collaboration, which significantly increased the number of gravitational-wave (GW) detections, reveals features with potentially important astrophysical implications. One notable example is a hint of the so-called pair-instability mass gap. In particular, the observed decline in the number of black holes (BHs) with mass above 45 Msun, together with indications of possibly higher spins for BHs above this threshold, has been interpreted by Antonini et al. and Tong et al. as evidence for pair-instability supernovae. In this work, we investigate whether mass transfer in binary systems can produce BH components' mass distribution that mimics the pair-instability limit. We use both the population synthesis code StarTrack and a simple semi-analytical framework to highlight the impact of mass transfer efficiency on the BH masses. We find that efficient mass transfer (over 50%) during the first Roche-lobe overflow, followed by a highly non-conservative second mass transfer phase, naturally limits the mass of the first-born BH and produces a cutoff that mimics a mass gap. While the upper mass limit for the more massive BH is increased through accretion during the first mass-transfer phase and is ultimately set by the pair-instability limit, the less massive BH is limited to the stripped primary mass. As a result, the fraction of systems in which the less massive BH exceeds 45 Msun is negligible. While the pair-instability mass gap is a plausible interpretation of current GW data, similar features can naturally arise from binary evolution. Future detections will help distinguish between these scenarios. In particular, a predominance of positive effective spins or low-spin events within the gap would challenge the pair-instability interpretation and instead support a binary-interaction origin for high-mass BHs.

Current space missions primarily depend on Earth-based Guidance, Navigation, and Control (GNC) systems involving human-in-the-loop operations. X-ray pulsar-based navigation offers a promising alternative by using the very precise periodic X-ray emissions from pulsars for fully autonomous state estimation. This study presents a comprehensive analysis of pulsar selection criteria that significantly influence overall navigation performance. Observational data from the NICER mission is used to derive realistic estimates of measurement noise. Key mission-level constraints, including pulsed flux, pulsar visibility, geometric configuration, and long-term timing stability, are integrated into the pulsar selection process, addressing limitations of existing studies. An extended Kalman filter (EKF) is used for onboard spacecraft state estimation. The proposed system is evaluated in two scenarios: a Low Earth Orbit (LEO) satellite at 600 km altitude and an interplanetary transfer from Earth to Jupiter. Simulation results show that including the Crab pulsar yields position errors below 7 km in LEO and 20 km during interplanetary transfer with an instrument effective area of 200~cm$^2$; however, the Crab's limited timing stability leads to filter divergence after 20 days without timing model updates. In contrast, more stable pulsars enable long-term autonomy but with reduced accuracy. These results highlight the trade-offs involved in pulsar selection for autonomous navigation and the need to balance competing objectives. Overall, this study demonstrates the feasibility of X-ray pulsar-based navigation and marks a key step towards fully autonomous spacecraft operations.

Mrk421 displayed its highest flux state ever observed in February of 2010 with very high TeV fluxes and interesting cross-band correlations and a spectral energy distribution (SED) evolution not entirely consistent with the standard single zone leptonic synchrotron self-Compton model. The source was already in a high state in January 2010 and displayed strong variability in the days preceding the highest state. We study the temporal evolution of the spectra in January to extract information about the particle dynamics and the physical properties of the emission region. We build up on the temporal variability and correlations studied in the previous work (MAGIC collaboration - Abe et al. 2025) and attempt to improve the SED model fits with a physics oriented approach. The multi-wavelength data was processed and the SEDs were fit using JetSeT. The SED evolution and cross band correlations were modelled using leptonic log-parabola with a low energy power-law branch (LPPL) and pile-up distributions that are predicted in a stochastic acceleration scenario. A simplified temporal evolution model was developed and fit to the SEDs and the resulting trends and phenomenology were characterised in context of theoretical literature. An expanding emission region model was also tested. We find the spectral variability to be well in agreement with stochastic acceleration. Our analysis suggests that the standard LPPL distribution develops a Maxwellian pile-up component at the transition from acceleration to cooling dominated phase on 3 nights in the dataset, as also hinted by the very-high energy and X-ray light curves. The resulting phenomenology of our sequential snapshot evolution SED model agrees well with theoretical and numerical simulation studies on temporal evolution using the diffusion equation approach.

We present a practical implementation of a phase-coherent mapping technique for pulsar timing arrays that resolves the full complex polarisation state of the gravitational-wave sky as a function of direction and frequency. Unlike standard cross-correlation methods, this approach preserves the amplitude, phase, and polarisation of the signal in every sky pixel. The resulting maps constitute a compact, minimally processed summary of the data from which all subsequent analyses -- characterisation of a stochastic background, searches for anisotropy, and identification of individual sources -- can be derived within a single unified framework. Our implementation is fully compatible with established pulsar timing data analysis methods. We validate the framework through a series of realistic simulations with varying array configurations, noise properties, and signal types. We demonstrate robust recovery of source amplitudes and sky locations across different scenarios, and discuss the impact of polarisation leakage, noise, and direction-dependent array sensitivity on the recovery of astrophysical signals.

We present the cosmic expansion scenario in the framework of $f(T)$ gravity by employing a dark energy equation of state (EoS) parameter. Specifically, we proceed with the power-law form of the function $f(T) = α$$(-T)^{n}$, in conjunction with the Gong-Zhang parametrization of the dark energy EoS. We derive the expansion rate in terms of the redshift for the considered model, providing deeper insights into the underlying cosmic dynamics. The model is further utilized to explore the expansion history of the universe and the evolution of several cosmological parameters. By using the Bayesian methods based on the $χ^{2}$-minimization technique, the median values of the model parameters are determined for the cosmic chronometer (CC) and joint (CC+Pantheon) datasets. The evolution of the deceleration parameter, energy density, pressure, EoS parameter and the energy conditions for dark energy is analyzed in detail. The model captures the observed acceleration as a transient phenomenon, followed by future deceleration. Additionally, the nature of both geometrical and dynamical diagnostics robustly indicates a quintessence-like behavior at the present epoch. Finally, the thermodynamic viability of the model is confirmed through the generalized second law of thermodynamics and the estimated age of the universe further supports the model's compatibility with astronomical observations.

Superradiant instability of ultralight bosons can produce clouds around rotating black holes, whose continuous gravitational wave (GW) emission is a promising observational target. Precise predictions of the signal frequency and its evolution are essential for detecting such continuous GWs. For axions, self-interactions can populate multiple superradiant modes via nonlinear couplings, and GW emission can occur through various channels. To calculate the frequency shifts of GWs emitted through these channels, we employ relativistic perturbation theory based on a bilinear form. We apply this framework to self-interaction effects for the first time, and also revisit the treatment of the self-gravity contribution. Our results provide a simple and unified framework for calculating frequency shifts, including cases in which multiple modes are excited, and are relevant for next-generation GW observations.

Most Active Galactic Nuclei (AGN) are Radio Quiet, with radio emission that may arise from star-formation activity, AGN-driven winds, weak jets, and coronal activity. Disentangling these mechanisms is challenging and requires detailed multi-wavelength investigation, but it is crucial for quantifying AGN feedback in galaxy evolution. We present a detailed radio investigation of 21 X-ray selected AGN in the Supermassive Black Hole Winds in X-Rays (SUBWAYS) sample (log Lbol = 44.9-46.3 erg/s, z=0.1-0.5), selected to systematically search for Ultra-Fast Outflows (UFOs). UFOs are detected in 30% of the targets, making the sample particularly well-suited for investigating the role and signatures of multi-scale outflows at different frequencies. We build the radio SED of the sources complementing our proprietary data, collected with the JVLA at 1.5 and 6 GHz, with images from LoTSS and other publicly available radio surveys between 150 and 1400 MHz. We investigate the role and occurrence of the aforementioned mechanisms, with particular interest in outflows and their possible relation with UFOs. We combined information on spectral indices, luminosities, and morphologies of the radio emission with properties derived in other wavebands, such as Star Formation Rate, X-ray luminosity, Eddington ratio or the UFO kinetic luminosity. All the sources are detected and are mostly consistent with RQ AGN. For 80% of the sources the data suggest the presence of an outflow (wind or weak jet). Interestingly, our results indicate that AGN with UFOs tend to have larger radio extension and a steep radio spectrum consistent with outflows. Moreover, the radio emission of the 6 UFO hosts is consistent with predictions from wind-driven shock models, possibly indicating a direct connection between the two phases. Alternatively, this may reflect physical conditions favouring the rise of both phenomena.

Large families of confining holographic QFTs, described by Einstein-Dilaton gravity, are considered on constant-curvature manifolds in the presence of a $θ$-angle. The space of ground states of such theories is explored as a function of the UV parameters, namely the dimensionless curvature and the $θ$ angle. The free energy is computed, and the phase structure is determined. For constant negative curvature manifolds, we find solutions dual to single QFTs as well as solutions describing interfaces. The single QFTs exhibit an infinite family of saddle points, with the leading one dominating the gravitational path integral and no phase transitions present. For constant positive curvature manifolds, like de Sitter, the ($θ$-angle, curvature) phase diagram exhibits both first and second order phase transitions, as a function of the class of theories considered. We also show that when $θ=0$, a holographic Vafa-Witten-like theorem can be proven.

Building on recent work introducing the idea of Quasilocal Probability in curved spacetime, we develop its observational implications for black hole ringdown in detail. We show that horizon-induced probability flux leads to an effective non-Hermitian dynamics producing three distinctive signatures, which are correlated multi-mode deviations, weak amplitude dependence and a mismatch between waveform damping and energy accounting. These effects arise from a single boundary-flux mechanism and therefore exhibit a constrained, low-dimensional structure not expected in generic modified gravity scenarios. We demonstrate that while individual deviations may be mimicked, their combined pattern provides a robust discriminator of quasilocal probability. We further argue that upcoming gravitational wave observations can probe these signatures at meaningful precision. We also establish that black hole ringdown is a novel arena to test whether quantum mechanical Hermiticity is really a fundamental property or an emergent symmetry in quantum gravity.

We establish a quantitative relation between local spin polarization and quantum entanglement in two-qubit systems by deriving an upper bound on the concurrence at fixed local polarization, showing that increasing polarization constrains the maximum achievable entanglement. We further demonstrate that this bound is saturated by pure states in certain cases. As a concrete physical application, we consider the parity-violating process $e^+e^- \to Z^0 \to q\bar{q}$, which generates final-state spin polarization. We show that the maximal concurrence is attained in specific kinematic regions and is significantly reduced relative to the unpolarized case. These results establish a general, process-independent framework connecting local polarization, maximal entanglement, and the role of pure states.

The Vela supernova remnant (SNR), G263.9-3.3, and its pulsar wind nebula (PWN), Vela X, is one of the closest such systems, and it has been studied using observations across the electromagnetic spectrum. SNRs are known sources of gamma rays with energies from GeV to the TeV range. In the GeV band, a cluster of cataloged unidentified Fermi-LAT point sources are found across the large angular extension of the Vela SNR. We aim to search for a high-energy signature associated to the SNR. We applied two independent machine learning algorithms to classify unidentified point sources in the Vela region by comparing their properties to those of known populations of Fermi pulsars and active galactic nuclei. We analyzed LAT data and modeled the spectrum of any emission attributable to Vela using leptonic and hadronic processes typical of SNRs. We find that most of the "point sources" cataloged within the extent of Vela do not share characteristics with those of the two most common Fermi point-like source populations and that even after the emission attributed to these "point sources" is subtracted, considerable residual emission is seen throughout Vela. Morphologically, most of the GeV emission is found within the shell of the SNR. We conclude that the majority of the cataloged point sources are likely spurious, and the GeV gamma rays come from an extended source, which we argue is the counterpart of the Vela SNR. Adopting a simple morphology given by a uniform disk for the emission the resulting extension is 6.5 deg. The northeastern portion of G263.9-3.3, where the ambient density is thought to be higher, is brighter in gamma rays than the south and west. The spectrum of the emission is best fit with a hadronic model. These facts make the hadronic origin for the gamma rays more likely.

This paper presents a systematic study of the chaotic dynamics of charged test particles around purely magnetically charged black holes immersed in a uniform external magnetic field within the framework of Einstein-ModMax theory. By constructing an explicit symplectic integrator, we obtain high-precision numerical solutions of the equations of motion. Combining the observational constraints from the Event Horizon Telescope (EHT) shadow images, we further restrict the parameter ranges of the model. We apply Shannon entropy and MIPP (mutual information for particle pairs) as effective indicators to identify the chaotic behavior of charged test particles in the spacetime of this black hole. Numerical results indicate that these indicators can clearly distinguish between regular and chaotic motion of orbits in strong gravitational field systems. Further analysis reveals that, compared to the key conserved quantities that determine the global dynamical behavior of the system -- energy $E$ and angular momentum $L$, the sensitivity of the system parameters $e^{-ν}$ and $Q_{m}$ to transitions in the orbital dynamical states is significantly reduced. This study provides a new perspective for a deeper understanding of the characterization and evolution mechanisms of chaotic dynamics in strong gravitational fields.

Gravitational spectroscopy tests the coupling of gravity to matter by measuring gravitationally induced frequency shifts of quantum transitions. While modern optical clocks probe the gravitational response of electronic transitions with extraordinary precision, tests in the nuclear sector have not progressed since the Mössbauer measurements of the gravitational redshift by Pound and Rebka. Here we introduce a new approach to nuclear gravitational spectroscopy based on phase-sensitive heterodyne interferometry of time-resolved nuclear resonant scattering of synchrotron radiation. In this scheme the gravitational redshift appears as a slowly accumulating phase drift of a delayed heterodyne beat signal, converting nuclear gravitational spectroscopy from energy-domain detection to time-domain interferometry. A Fisher-information analysis supported by Monte Carlo simulations and experimentally confirmed photon count rates shows that the nuclear gravitational redshift of $^{57}$Fe can be detected within hours on a few-meter-scale vertical baseline, with percent-level precision on deviations from general relativity becoming accessible on day-scale timescales. The method thus establishes an experimentally realistic and scalable platform for precision tests of gravitational coupling to nuclear structure.

These notes give a concise introduction to General Relativity at the advanced undergraduate level, starting from the weak field limit and gravitational waves, then introducing curved manifolds and Riemannian geometry. The nonlinear gravitational action is used to derive the nonlinear field equations, with applications to black holes and cosmology. It is assumed that special relativity and electromagnetic waves have been previously studied. Some advanced topics such as Rindler and Hawking radiation are derived, and recent developments in gravitational wave detection are briefly covered. Problems are included, both those suitable for homework, and simpler ones that could be worked out by students during class sessions.

We present a new public-domain MOlecular Fitting Analysis Tool (MOFAT) designed to probe molecule-forming regions in supernovae (SNe) through analysis of molecular features in the near- and mid-infrared. MOFAT employs a novel data-driven approach to explore the physical properties of these regions using time-independent radiative transfer simulations that include multidimensional, clump-like structures, constrained by high-precision observations. Such structures are required to reproduce the flux ratio between fundamental and overtone bands, overcoming limitations of traditional one-zone forward-modeling, such as optical-depth effects and initial configurations. Our approach enables spectral fits that can reconstruct overall abundances and temperatures and determine parameterized small-scale structures associated with physical instabilities. We systematically study the relationship between physical parameters and the profiles of CO and SiO, showing that free parameters are constrained, while detection of small-scale structure requires optically thick bands. As a demonstration, MOFAT is applied to SN2024ggi at +285 and +385 days post-explosion. We find that CO formation triggers SiO formation in the inner layers of the CO-rich region previously studied. The inner edge of the SiO-emitting region recedes from velocities of v1 from 1500 to 1000 km/s, indicating continued SiO formation. The SiO mass decreases from about (2-6)E-3 Mo by roughly an order of magnitude, suggesting ongoing evaporation. SiO features indicate clumping, but most of the flux originates from optically thin regions. SiO contributes negligibly to cooling, and we find no evidence for dust formation. Finally, we discuss observational strategies to trace the evolution of molecule formation and its connection to dust formation.

NICER observations of millisecond pulsars (MSPs) suggest that non-dipolar magnetic fields are required to explain their surface X-ray hotspots. C. Kalapotharakos et al. (2021) modeled the NICER light curve of MSP J0030+0451 (J0030) using a static vacuum offset dipole-plus-quadrupole field and corresponding force-free (FF) solutions to jointly reproduce the X-ray and Fermi-LAT $γ$-ray emission. We substitute their static vacuum field model with a more realistic swept-back configuration that accounts for rotational effects. This field more closely resembles the corresponding FF solutions, making it a more physically motivated choice for future multiwavelength modeling. We adopt a centered swept-back vacuum multipolar magnetic field (SVM2F; J. Pétri 2015), expressed as a complete expansion in vector spherical harmonics, enabling flexible descriptions of arbitrary magnetic field geometries. We introduce a metric to quantify the complexity among different field prescriptions, illustrated for the static offset vacuum field. To efficiently explore parameter space, we train a neural network surrogate (G. Olmschenk et al. 2025) on SVM2F light curves including components up to the octupole, accelerating Markov chain Monte Carlo sampling by $\sim 10^3$ compared to direct physical model evaluations. Applying this framework to J0030, we constrain the field parameter space and find that a centered swept-back multipolar field including terms up to the octupole adequately reproduces the bolometric thermal X-ray light curve. Our study highlights the importance and inherent complexity of prescribing different multipolar magnetic field models for rotating stars, and can be extended to other MSPs to ultimately constrain the masses and radii of neutron stars, and hence their equation of state.

In core-collapse supernovae (CCSNe), deleptonization normally favors $ν_e$ over $\barν_e$ emission. However, lepton-number emission self-sustained asymmetry (LESA) can make the energy-integrated emitted lepton-number flux negative along some directions. We derive a simple diagnostic for this transition and test it in 33 independent 3D CCSN simulations: 25 Princeton/Fornax models ($8.1$--$100\,M_\odot$) and 8 Garching models, including non-, slow-, and fast-rotating $15\,M_\odot$ cases. Of 23 non-black-hole-forming Princeton models, 22 cross the threshold, with median onset $t_c=225\,\mathrm{ms}$, IQR $162$--$264\,\mathrm{ms}$, and cross-model scatter $\mathrm{CV}=18.6\%$. Full-sky flux-sign searches show that the threshold identifies the anti-LESA-pole transition, distinguishing the global LESA-driven crossing from early localized turbulent crossings. The fast-rotating Garching $15\,M_\odot$ model, where rapid rotation suppresses the LESA dipole, is correctly classified as a non-crosser without using any rotation parameter. Both black-hole-forming Princeton models cross near $250\,\mathrm{ms}$ post-bounce and remain above threshold for $1807$ and $2463\,\mathrm{ms}$ before collapse. Thus, in the next nearby CCSN, the emitted $\barν_e$ energy flux may exceed the $ν_e$ flux along some lines of sight. Such directions may also correlate with sustained fast flavor instability, although testing this requires local phase-space distributions or dedicated linear stability analysis. The relevant quantity here is the energy-integrated emitted flux field, i.e. a luminosity difference per steradian, not a neutrino number flux.

Flavor instabilities develop in neutrino plasmas through emission of flavomons, the quanta of flavor waves. We derive the flavomon equations of motion in slowly varying environments, notably the matter gradients of supernovae, and use them to construct a flavomon ray tracing framework. Combined with a quasi-linear description of flavomon growth, we thus develop a new approach to the global evolution of flavor instabilities. As a first application, we show that the growth of neutrino-mass-induced instabilities is slowed down, but not suppressed, by the inevitable matter gradients. Local stability analysis alone cannot gauge the impact of inhomogeneities and instead must be coupled to flavomon ray tracing.

In this paper, we investigate the quasinormal modes of two covariant effective black holes characterized by the quantum parameter $ζ$, charge $Q$, and cosmological constant $Λ$, under the scalar perturbation. By employing the pseudo-spectral method, we numerically calculate the quasinormal frequencies and analyze the influence of $ζ$ on the spectra with respect to $Q$. Our results demonstrate that while the quantum parameter $ζ$ significantly modifies the quasinormal frequency spectrum, the non-monotonic behavior and overtone outbursts persist. Notably, the impact of quantum gravity on the overtone outbursts is not merely limited to enhancement or suppression; instead, it introduces additional spectral features. Furthermore, a comprehensive analysis of the full quasinormal mode spectrum reveals rich interactions between complex and purely imaginary modes, including damping-rate crossings and merging-splitting behavior. These phenomena typically accompany overtone outbursts in near-extremal regimes, suggesting a potential connection between mode interactions and overtone outbursts. This work emphasizes the necessity of analyzing the full quasinormal frequency spectrum rather than focussing solely on fundamental modes, and provides novel insights into its underlying spectral structures.

Gravitational lensing constitutes one of the most direct observational manifestations of spacetime curvature and provides a powerful probe of compact astrophysical objects. In this work, we present a comprehensive analysis of the bending of light in curved spacetime, beginning with the fundamental aspects of gravitational lensing and the Newtonian approximation to light deflection. The relativistic formulation of lensing is then developed through the lens equation, lensing potential, and a geometrical interpretation of Fermat's principle using an effective refractive index in curved spacetime. Photon trajectories and light deflection are subsequently investigated in static, spherically symmetric geometries, followed by a detailed study of photon motion in the equatorial plane of the Kerr spacetime. Analytical expressions for the closest approach distance and critical parameters governing photon orbits are derived. Furthermore, the bending angle is examined using the Rindler-Ishak method and the Gauss-Bonnet theorem within the optical geometry. Finally, the analysis is extended to axisymmetric spacetimes using the OIA and GW-OIA formalisms, providing a unified geometrical framework for computing light deflection in both static and rotating gravitational fields.

We report the discovery of a radio-quiet to radio-loud transition in the narrow-line Seyfert 1 galaxy J1105+1452. The source has undergone a long-term evolution from a radio-quiet state in the 1990s to a persistently radio-bright state after 2017. Post-2017 flux densities in the $0.8$-$7$ GHz range cluster between $32$ and $43$ mJy, whereas the $144$ MHz flux density is only $1.94 \pm 0.23$ mJy. This indicates strong low-frequency suppression from a compact, absorbed component. Modeling the radio spectral energy distribution with a synchrotron self-absorption model yields a turnover frequency $ν_{\rm p} = 0.48 \pm 0.03$ GHz and a peak flux density $S_{\rm p} = 38.9 \pm 4.7$ mJy. These parameters classify J1105+1452 as a megahertz peaked-spectrum source, consistent with the new episode of an early-stage compact jet. Under the assumption of equipartition, we derive an intrinsic physical radius $R \sim 0.68$ pc and an average apparent expansion velocity $β_{\rm app} \approx 0.64$. The observed brightness temperature $T_b \approx 6.0 \times 10^{11}$ K necessitates a Doppler factor $δ\approx 12$, implying a relativistic jet viewed at $θ\lesssim 5^\circ$. Despite the dramatic radio evolution, the X-ray spectrum remains stable and steep ($Γ\simeq 3.0$), suggesting that the X-ray emission remains dominated by the disk-corona, while the radio band has become jet-dominated. Our results identify J1105+1452 as a rare radio changing-state NLSy1, providing a unique laboratory for studying the birth and early evolution of relativistic jets at high Eddington ratios.

We compare the correlation functions of inflationary perturbations computed either with quantum or classical dynamics. Even if they are enforced to agree at a specific time during inflation, classical and quantum correlations will differ at the end of inflation, provided that interactions are relevant. The difference between the results of the classical and quantum computations is exponentially sensitive to the number of e-folds elapsed from the time of agreement. We illustrate this finding with the tree-level bispectrum of the primordial curvature fluctuation and the one-loop power spectrum of tensor modes. We also show that classical evolution from a finite time does not imply the appearance of poles in the scalar bispectrum.

Neutron stars make a unique astrophysical test bench for our understanding of quantum physics at kilometre scales. The rotation of a neutron star features glitches, sudden spin-ups that interrupt the otherwise regular stellar spin-down, which are often attributed to the dynamics of pinned quantised vortices in one or several of the superfluid phases inside the star. Laboratory experiments probing superfluid vortices have inspired neutron star theory and simulations from the beginning. Here we argue that vortex experiments in superfluids contained in aerogels show phenomenology that offers a highly appealing but vastly unexplored analogue for neutron star physics. We build a point-vortex simulation that allows analysing experiments in a crust-like and a core-like aerogel, extracting two different regimes of pinned vortex (non-)dynamics and validating a microscopic picture of very strong vortex pinning. In the crust-like aerogel, vortices get depinned once the ambient superflow is fast enough, while in the core-like aerogel pinned vortices are never released and rotational velocity changes are accommodated by the avalanche-like production of new vortices instead. Finally, we show that these concepts should apply also in neutron stars and may thus revolutionise the analysis of neutron star observations.

The shock fronts of supernova remnants (SNRs) are believed to be significant sites of acceleration of cosmic ray particles. Previous researchers have shown that a particle distribution similar to a log-parabola can be generated when particles have an energy-dependent escape. We explore the acceleration of electrons at SNR shock fronts, and show that modeling this energy-dependent particle escape model can produce spectral energy distributions consistent with observations of two lepton-radiation-dominated SNRs: RX J1713.7-3946 and SN 1006. The model includes the evolution of both the electron distribution and photon spectra as a result of the combined effects of the SNR evolution and electron energy loss. The electron-escape energy dependence is quite weak, but the electron distribution and photon spectra turn out to be very sensitive to changes in the electron escape. We also explore how sensitive the spectra and electron distributions are to the parameters used in the log-parabola model.

In this paper, we investigate the cosmic expansion scenarios within the framework of $f(Q,T)$ gravity by using the affine equation of state (EoS) parameter. Specifically, we consider the linear form $f(Q,T)=Q+βT$, where $β$ is a free model parameter. We use Bayesian statistical methods, specifically the $χ^2$ minimization technique to constrain the model parameters using Cosmic Chronometer (CC), Pantheon+SH0ES and DESI BAO data. We further analyze the characteristics of the derived cosmological model. A comprehensive study of energy density, pressure, equation of state parameter and cosmographic parameters are carried out to understand the evolution of the Universe in this model. The determination of present age of the universe for this model is within $1σ$ with Planck results.

We study uncertainties in the equation of state of neutron stars using conformal prediction as a distribution-free and model-agnostic method that provides coverage guarantees. In particular, we apply the Conformalized Quantile Regression (CQR) method to posterior samples calculated from Bayesian inference, creating reliable uncertainty bands without assuming a specific form of the underlying distribution. We first construct CQR bands as a postprocessing step to the posterior samples of neutron star mas-radius relations provided by the NMMA collaboration and to Quantum Monte Carlo calculations of pure neutron matter. In all cases, empirical coverage studies confirm the robustness of the method.

We investigate the sensitivities of upcoming MeV gamma-ray telescopes to sterile neutrino dark matter in the mass range $(0.2-100)\,{\rm MeV}$. Sterile neutrinos in this regime can produce observable photon signals through radiative two-body decays and three-body decays with final-state radiation. We perform a Fisher forecasting analysis incorporating realistic astrophysical background modeling and detector response to derive projected constraints on the sterile neutrino decay rate. We find that future MeV instruments can improve existing limits by several orders of magnitude across a wide region of parameter space. Our results highlight the discovery potential of next-generation MeV telescopes in probing sterile neutrino dark matter.

We consider a search for the anisotropic stochastic gravitational-wave background (SGWB) that decomposes the sky map into its spherical harmonics components in order to obtain estimators of the angular power spectrum. Such a search often requires the inversion of a Fisher information matrix which contains small singular values. Rather than dealing with biases induced by regularization methods used to facilitate this matrix inversion, we opt to avoid this inversion step entirely by working in the so-called ``dirty map" space, and we introduce methodology for statistical model inference in this space. We apply our methodology to simulated model signals added to detector noise characterized by Advanced LIGO's third observing run and consider angular power spectra for both the SGWB auto-correlation search as well as a cross-correlation search between the SGWB and electromagnetic tracers of matter structure in the universe. In both cases we are able to reliably recover simulated model parameters for sufficiently strong signals up to maximum order spherical harmonic modes of $\ell_{max}=10$. We find the limitations of our methodology to arise from the computational cost of testing complex models, the assumption of Gaussianity of the angular power spectrum, and a cosmic variance-like source of uncertainty which scales with the strength of the underlying signal.

Recent LIGO-Virgo-KAGRA (LVK) analyses have revealed complex structure in the binary black hole (BBH) population, including distinct features in the primary mass spectrum and nontrivial spin-mass correlations. However, the phenomenological models used to capture these features often lack analytic transparency, making it difficult to isolate robust physical laws from modeling artifacts. To address this, symbolic regression is applied to the posterior inference products of the GWTC-4 catalog, discovering compact, closed-form analytic expressions for four key population relationships: (i) the merger-rate evolution with redshift; (ii) the mass-ratio dependence of the effective-spin distribution; (iii) the redshift evolution of the effective-spin distribution; and (iv) the conditional mass-ratio distributions associated with the 10 solar mass and 35 solar mass primary mass peaks. This framework successfully compresses both rigid and highly flexible models into differentiable phenomenological laws, dynamically recovering a consistent low-redshift merger-rate slope without assuming an a priori power-law form. The exact analytic derivatives provided by symbolic regression show that the mass ratio--effective spin and redshift--effective spin correlations are robustly driven by broadening of the posterior widths rather than shifts in the mean. Furthermore, qualitatively distinct functional forms for the mass-ratio distributions conditioned on the 10 solar mass and 35 solar mass primary mass peaks are identified. These closed-form expressions enable exact analytic gradient diagnostics and compact surrogate summaries, particularly for flexible numerical posteriors that are not otherwise available in low-dimensional analytic form. They also facilitate rapid downstream calculations for rate forecasting, formation channel comparison, and stochastic background estimation.

We investigate the leading gravitational eikonal in nonlocal $D$ dimensional theories of gravity. We analyze the simplest cases of $2\rightarrow2$ massless and massive scalar scattering at tree level, studying the effects of nonlocal form factors in the gravitational sector. We give an interpretation of our results in terms of geodesic motion in effective generalized Aichelburg-Sexl geometries for the massless case, and in smeared linearized Schwarzschild metrics for the massive case in the probe limit. Combining our results for the geometries at linearized level with general requirements about the behaviour of the solutions in the core, we propose a nonlinear completion of the geometries. The resulting spacetimes describe singularity-free, asymptotically flat deformations of the Schwarzschild solution with a de Sitter core. We also analyze the main geometric and thermodynamic features of these solutions.

Studying the content and distribution of molecular gas provides key insights into how feedback from Active Galactic Nuclei (AGN) and star formation influences galaxy evolution, since molecular gas is the primary fuel for star formation. Ultra-Luminous Infrared Galaxies (ULIRGs) are ideal candidates to study how AGN and/or starbursts affect the interstellar medium due to their intense AGN and star forming activity. We present spatially-resolved multi-phase molecular gas study of IRAS20551-4250, a nearby ($z=0.0429$) ULIRG, using JWST/MIRI-MRS and ALMA. Mid-infrared diagnostics do not rule out the presence of AGN in IRAS20551-4250. [OIII]$λ$5007 in VLT/MUSE data reveal ionised gas outflows with $w_{80}^{\rm [OIII]} \sim 790$ km s$^{-1}$ and $\dot{M}_{\rm out}^{\rm[OIII]}<0.01$ M$_{\odot}$ yr$^{-1}$. No outflows are observed in either molecular phases. JWST/MIRI-MRS data reveal several rotational transitions of warm H$_{2}$ (T$\sim500-1400$ K) within the central $\sim4\times4$ kpc$^{2}$ region. Excitation temperature maps suggest that the warm H$_{2}$ is primarily heated by UV radiation from the central source. The CO-based cold molecular component dominates the molecular gas mass, accounting for $>$95% of the total molecular gas mass. Warm H$_{2}$ maps show two tidal tails and the velocity centroid maps show disturbed, non-rotational motions and a systematic gradient across the field-of-view, similar to that of ALMA CO-based cold molecular gas and consistent with a late-stage merger. Together, our analysis indicate that the molecular gas composition in IRAS20551-4250 is consistent with ongoing star formation in the host galaxy and the outflows observed in ionised gas phase appear insufficient to expel the molecular gas or quench ongoing star formation.

The origin of intergalactic magnetic fields (IGMFs) remains one of the key open questions in cosmology. Gamma-ray pair halos produced by electromagnetic cascades from TeV-emitting blazars provide a powerful indirect probe of these fields. In this work, we present a novel search for pair halos that explicitly exploits their expected anisotropic morphology, aligning with the projected orientation of blazar jets on the sky. Using a Monte Carlo framework to model the spatial distribution of cascade emission, we identify an optimal sample of 21 high-synchrotron-peaked BL Lac objects with well-constrained jet position angles from radio interferometry. By rotating and stacking \textit{Fermi}-LAT observations of these sources along their jet directions, we enhance sensitivity to anisotropic extended emission that would be diluted in traditional orientation-agnostic analyses. Applying a likelihood analysis to the combined dataset, we find evidence for a non-zero IGMF, excluding the null hypothesis at $3.8σ$ level and obtaining a best-fit field strength of $B_0 = 2.8 \times 10^{-16}\,\mathrm{G}$, with a $99\%$ confidence interval of $0.9 \times 10^{-16}\,\mathrm{G} < B_0 < 8.9 \times 10^{-16}\,\mathrm{G}$. Our result is consistent with previous constraints from spectral, spatial, and temporal studies, while demonstrating that incorporating anisotropic information provides a significant gain in sensitivity. This approach opens a new avenue for probing intergalactic magnetism and highlights the potential of future high-angular-resolution gamma-ray observations to directly image pair halos and map magnetic fields in cosmic voids.

We report the first detection of radio emission from the Type II supernova SN 2007it, located at a distance of 12.2 Mpc in NGC 5530. The observations were obtained with the Australian Telescope Compact Array (ATCA) more than 18 yr after the explosion as part of the Rebrightening in Interacting Supernova Emission (RISE) program, which monitors nearby core-collapse supernovae for late-time interaction with dense circumstellar material. SN 2007it was detected on 2026 April 8 (08:00-12:00 UTC) at 5.5 GHz with a flux density of $3.30 \pm 0.13$ mJy and at 9.0 GHz with $3.54 \pm 0.24$ mJy. Its non-detection in publicly available 0.88 GHz ASKAP data from 2026 January 11 suggests either rapidly rising emission or significant internal absorption at lower frequencies. We assess the prospects for detection at other wavelengths and encourage coordinated follow-up observations across the radio, optical, X-ray, and $γ$-ray band

Gravitational wave (GW) detector LISA will observe near-coalescence extreme mass ratio inspirals (EMRIs), which typically form in galactic central accretion disks. Gas torques on EMRI will alter its GW-driven inspiral trajectory from the vacuum expectation, leading to potentially LISA-observable GW dephasing ($Δψ_{\rm gas}$). Most studies compute $Δψ_{\rm gas}$ for a thin, laminar disk, with negligible flow turbulence, where the disk exerts a fairly well-understood linear torque ($T_{\rm lin}$). However, these disks must be turbulent due to magneto-rotational instability in the inner regions. Hence, we present a proof-of-concept general, agnostic prescription for the turbulent torque ($T_{\rm turb}$) acting on an EMRI by modeling it as a Gaussian distribution around $T_{\rm lin}$, based on recent advances from a global hydrodynamical (HD) study. We compute $Δψ_{\rm gas}$ for the ``golden'' circular EMRI with total source mass $M=10^6~{\rm M}_\odot$ and mass ratio $q=5\times10^{-5}$ in its final four-year evolution at redshift $z=0.276$ and signal-to-noise ratio (SNR) $=50$ by varying Eddington ratio ${\rm f}_{\rm Edd}$, turbulence normalization $C$ ($=~360$ in the aforementioned HD study), disk aspect ratio $h_0$, and turbo-viscous coefficient $α$ in a reasonable parameters space. We find that for ${\rm f}_{\rm Edd}\gtrsim0.3$, $C\gtrsim300$, $h_0\gtrsim0.03$, and $α\gtrsim0.1$, gas-induced dephasings are unobservable if only considering $T_{\rm lin}$ but could become detectable ($Δψ_{\rm gas}>8/$SNR) if EMRIs exhibit chaotic migration due to turbulent gas flow. Hence, this work motivates running MHD simulations of accretion disks with embedded LISA EMRIs in the early in-spiral phase over long enough timescales to understand the evolution of their orbital elements and the imprint of the turbulent environment on their gravitational waveforms.

Optical reciprocity--the principle that light retraces the same path when source and detector are interchanged--is a foundational concept in geometric optics. In this Letter, we demonstrate that this ``symmetry-protected'' behavior can be qualitatively overturned in a rotating black hole when spontaneous Lorentz symmetry breaking introduces a nonminimally coupled background structure with a preferred direction. Through numerical ray-tracing simulations, we reveal a striking macroscopic signature: upon optical-path reversal achieved by exchanging the source and the observer, the shadow of the same black hole morphs from a quasi-symmetric rugby-ball shape into a distinct teardrop profile. This high-contrast nonreciprocity effectively turns the black hole into a cosmic-scale optical diode, offering a novel pathway to probe fundamental symmetries using current and next-generation horizon-scale imaging.

We investigate the performance of quantum algorithms for light nuclear systems by studying the deuteron (2H), triton (3H), and helium-3 (3He) nuclei within a lattice formulation of pionless effective field theory (EFT). We first compute ground-state energies using classical exact diagonalization (ED), serving as a benchmark reference for variational quantum algorithms. We then perform Variational Quantum Eigensolver (VQE) calculations using noiseless classical statevector simulations of quantum circuits, enabling a controlled assessment of algorithmic performance in the absence of hardware-induced noise. We calibrate the two-body low-energy constant using the deuteron system and fit the three-body interaction strength to the triton, then consistently apply the resulting Hamiltonian parameters to the helium-3 nucleus. Our VQE calculations employ physically motivated ansatze targeting the relevant particle-number sector, with explicit particle-number-conserving constructions implemented for the triton and helium-3 systems. The variational optimization includes an analysis of the Hamiltonian energy variance roviding additional insight into convergence behavior and the quality of the variational states. We find that the VQE results are in good agreement with the corresponding classical ED ground-state energies across all three systems, including the isospin-asymmetric helium-3 nucleus with Coulomb interactions. Overall, our study provides a transparent and reproducible benchmark for assessing the applicability of variational quantum algorithms to few-body nuclear systems. Additionally, we perform a noisy VQE simulation with a depolarizing noise model for the triton system to illustrate the impact of realistic Noisy Intermediate-Scale Quantum (NISQ)-era hardware noise on variational energy estimation.

This study investigates the visual characteristics of a rotating black hole (BH) within the fabric of $4$D Einstein-Gauss-Bonnet gravity illuminated with two illumination models, such as a celestial light sphere and a thin accretion disk. To visualize the BH shadow images, we use a recent fisheye camera model and ray-tracing method. And then, we focus on investigating the impact of the coupling parameter $α$ and the spin parameter $a$ on the shadow images. The results exhibit that the shadow radius decreases, while the shadow deviation increases with the aid of $α$. However, with respect to $a$, the shadow radius is slightly increased compared to the corresponding shadow deviation. For a celestial light sphere, the increasing values of $α$, lead to a decrease in the corresponding photon ring, while the space-dragging effect becomes more prominent with increasing $a$. For a thin accretion disk, we enhance its inner edge to the BH event horizon, and the particle motion is different in the regions inside and outside the innermost stable circular orbit. The result demonstrates that the shadow becomes progressively asymmetric with $a$, while the overall size of the inner shadow gradually decreases with the variations of $α$. Subsequently, we also investigated the distinct features of red-shift configurations on the disk for both direct and lensed images. Additionally, we used the latest observational data from M87* and Sgr A* to impose certain parameter constraints on $α$; the results depict the consistency of our considering the BH model.

Systematic errors in the parameter estimation (PE) of gravitational wave (GW) mergers can arise from various sources, including waveform systematics, noise mischaracterization, data analysis artifacts, and other unknown factors. In this study, we analyze selected events from the first three observing runs of the LIGO-Virgo-KAGRA (LVK) collaboration. We choose events that have been flagged in various studies as potentially affected by systematic errors. Here, we reanalyze these events using a couple of parametric models developed in previous work that incorporate uncertainties in both the phase and amplitude of the GW waveform. In this data-driven approach, we apply sufficiently broad priors on the uncertainty parameters to account for potential systematic errors. Our findings show that the proposed method effectively reduces systematic errors, even those arising from data artifacts, such as glitches occurring near a signal and the deglitching process in GW frame files. Similarly, inconsistent results from different waveform models become much more consistent in our framework. One noteworthy event we examine is GW191109\_010717, which is particularly interesting due to its anti-aligned spin properties. We report that, within our framework, the event still exhibits anti-aligned spin characteristics, but the inference results become consistent across raw and deglitched frame files, as well as across the waveform models used for this event (IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4). A similar trend is observed for the event GW200129\_065458, which previously yielded a high, but inconsistent precession parameter among different waveform models. In contrast, we observe a non-zero and consistent value of $χ_{p}=0.60^{+0.31}_{-0.33}, 0.58^{+0.30}_{-0.29}$ and $0.56^{+0.31}_{-0.28}$ for the IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4 waveform models, respectively.

We show that the LIGO--Virgo--KAGRA (LVK) verification of Hawking area law carries profound consequences for quantum gravity if such a law is postulated to hold exactly. The observed mergers can be produced in local Stelle gravity and in nonlocal quantum-gravity theories with entire or fractional form factors either by (i) singular Ricci-flat black holes or (ii) possibly regular classical black holes under very restrictive conditions: absence of $R^2$ and (Riemann)${}^2$ terms in the action, absence of extra real poles in the graviton propagator, and positivity of its spectral representation. To date, this is the strongest simplification of the ambiguities of this class of theories. We also prove that the classical standard black-hole entropy-area law holds as a consequence of Hawking area law, and provide a rigorous realization of Barrow's fractal black holes otherwise.

IceCube is a cubic-kilometer-scale neutrino detector located at the geographic South Pole. A precise directional reconstruction of IceCube neutrinos is vital for associations with astronomical objects. In this context, we discuss neural posterior estimation of the neutrino direction via a transformer encoder that maps to a normalizing flow on the 2-sphere. It achieves a new state-of-the-art angular resolution for the two main event morphologies in IceCube - tracks and showers - while being significantly faster than traditional B-spline-based likelihood reconstructions. All-sky scans can be performed within seconds rather than hours, and take constant computation time, regardless of whether the posterior extent is arc-minutes or spans the whole sky. We utilize a combination of $C^2$-smooth rational-quadratic splines, scale transformations and rotations to define a novel spherical normalizing-flow distribution whose parameters are predicted as a whole as the output of the transformer encoder. We test several structural choices diverting from the vanilla transformer architecture. In particular, we find dual residual streams, nonlinear QKV projection and a separate class token with its own cross-attention processing to boost test-time performance. The angular resolution for both showers and tracks improves substantially over the whole trained energy range from 100 GeV to 100 PeV. At 100 TeV deposited energy, for example, the median angular resolution improves by a factor of $1.3$ for throughgoing tracks, by a factor of $1.7$ for showers and by a factor of $2.5$ for starting tracks compared to state-of-the art likelihood reconstructions based on B-splines. While previous machine-learning (ML) efforts have managed to obtain competitive shower resolutions, this is the first time an ML-based method outperforms likelihood-based muon reconstructions above 100 GeV.

The 511 keV electron-positron annihilation feature near the galactic center has been detected for more than half a century, yet its origin remains a mystery. In this paper, we describe a concept for a balloon-borne 511 keV $γ$-ray mission called the 511-Spectrometer Mission. The mission will use Transition-Edge Sensor (TES) arrays with thick metal absorbers that are thermally coupled to the TES. The strength of the approach is a projected energy resolution of 200 eV Full Width Half Maximum (FWHM) at 511 keV, enabling detailed studies of the shape and substructure of the 511 keV emission from the galactic center region. A first mission equipped with 8,192 $γ$-ray detectors and a fully active shield and collimator could detect the galactic center with ~35 $σ$ statistical significance. We present the mission concept as well as first results obtained with a prototype detector equipped with $1.35\times1.35\times2$ mm$^{3}$ Bi absorbers. The detector has a quantum efficiency of 15% for 511 keV photons in photoelectric effect interactions. In tests with a $^{137}$Cs source, these prototype detectors show an energy resolution of 525 eV FWHM at 662 keV. We end with a discussion of follow-up missions that use coded mask imaging, or use concentrating or focusing optics to scrutinize the sources of 511 keV $γ$-rays on smaller angular scales.

We study the late-time evolution of an anisotropic Bianchi-I universe with radiation in the framework of asymptotically safe gravity. We first discuss the radiation-dominated universe for the perfect fluid with the equation of state $p=ρ/3$, and find that the classical evolution involves logarithmic terms, which lead to a slow approach toward isotropy. The quantum effects introduce subleading corrections that soften the anisotropy in the intermediate stage. Next we discuss the universe with magnetic fields. For a vanishing classical cosmological constant, we find that the universe in general evolves toward a Kasner-type regime with persistent anisotropy while the expansion rate is enhanced by quantum effects, leading to a faster decay of the magnetic field. In contrast, for a nonzero classical cosmological constant, the late-time dynamics are dominated by the cosmological constant, and the universe asymptotically approaches an isotropic de Sitter phase with exponential decay of both anisotropies and the magnetic field. Finally, we employ Hodge duality to demonstrate that these cosmological findings apply equally to environments dominated by electric fields.

We present high-sensitivity Very Long Baseline Interferometry (VLBI) observations of four ultraluminous X-ray sources (ULXs): Holmberg II X-1, IC 342 X-1, NGC 6946 X-1, and NGC 925 X-1. No compact emission was detected on milliarcsecond scales, with rms noise levels reaching approximately 5--20 $μ$Jy. The corresponding $5σ$ flux density upper limits reach $\sim 26\,μ\mathrm{Jy}$, implying radio luminosity limits $L_{\rm R} \lesssim 2 \times 10^{33}\,\mathrm{erg\,s^{-1}}$. This disfavors any persistently bright hard-state-like compact core at our sensitivity level. The previously reported VLBI core in Holmberg II X-1 exhibits significant long-term variability, broadly consistent with an overall decline over the past decades. This behavior is consistent with emission from optically-thin ejecta undergoing adiabatic expansion. The VLBI non-detections may reflect intrinsically weak/intermittent compact emission, and/or low--surface--brightness structure that is resolved out by VLBI, and/or absorption/propagation effects such as free--free absorption in dense, ionized winds.

The properties of the spectral evolution during the prompt emission phase of gamma-ray bursts (GRBs), which are closely related to the radiation mechanism (synchrotron or photosphere), are still a subject of debate. Two spectral evolution patterns (``hard-to-soft'' and ``intensity-tracking'') have been commonly observed in GRB prompt emission spectra. Here we present a well-defined sample of 20 single-pulse GRBs detected by \emph{Fermi} whose prompt emission spectra exhibit the intensity-tracking pattern. By performing a time-resolved spectral analysis, we derive $E_{\rm p}$ and the energy flux $F$ from the same time bins and introduce a matched-bin lag, $t_{\rm lag}^{\rm F} \equiv t_{\rm p}(E_{\rm p})-t_{\rm p}(F)$, where $t_{\rm p}$ denotes the time at which each quantity reaches its maximum. We find that the intensity-tracking pattern subdivides into three distinct subclasses: Type I (5/20), with aligned $E_{\rm p}$ and flux peaks; Type II (13/20), with $E_{\rm p}$ peaking before the flux; and Type III (2/20), with $E_{\rm p}$ peaking after the flux. The early-peaking Type~II subclass dominates the sample. The subclasses also exhibit systematic differences in their spectral and temporal properties. Type II bursts are systematically harder than Type I, show broader flux pulses, and more often display asymmetric rising and decaying $E_{\rm p}$-$F$ branches. Type I is consistent with tightly coupled spectral and power evolution, whereas Type II is more naturally explained by nonthermal or hybrid prompt-emission scenarios in which spectral hardening precedes the peak radiative output. Type III appears to form a rare positive-lag tail whose physical origin remains uncertain.

We analyze the effect of using the Fermi-Dirac statistics, rather than its Boltzmann approximation, in numerical simulations of perfect spin hydrodynamics of particles with spin 1/2. The system considered is boost invariant, transversely homogeneous, with corrections to the baryon current and the energy-momentum tensor that are second order in the spin polarization tensor $ω$, and the spin tensor considered is first order in $ω$. The study shows the feasibility of this approach, as the special functions defined by integrals that appear in the coefficients in the Fermi-Dirac case can be conveniently parametrized. For sets of initial conditions used in previous works, the differences in parameter evolution between the two underlying particle statistics are about one order of magnitude smaller than corrections coming from spin feedback. We also discuss when and why the numerical solutions of the equations of perfect spin hydrodynamics break down for very large values of spin polarization in one of the geometric configurations considered.

This work presents a non-equilibrium framework for thermodynamicized black holes, inspired by the entropy-functional interpretation of emergent gravity and by residue-based methods in black hole thermodynamics. The main idea is to unify three components: an entropy functional principle for selecting physical on-shell backgrounds, a Euclidean and contour-based description of the horizon temperature through simple pole singularities, and a topological residue classification of multi-horizon black hole configurations. On this basis, the paper introduces a quasi-stationary non-equilibrium partition functional in which irreversible entropy production appears as an additional contribution to the singular action. The formalism reproduces the standard equilibrium relations in the adiabatic limit, while also extending them to dynamically driven black-hole systems with matter, charge, and rotational fluxes. The framework is then applied to Kerr Newman type black holes in constant curvature f(R) gravity, where the equilibrium entropy remains weighted by the derivative of f at the background curvature, while non-equilibrium corrections arise from flux-induced deformations of the effective thermodynamic action. The analysis further shows that the outer and inner horizons carry opposite topological orientations, so the non-extremal Kerr Newman family stays in the topological class W = 0 unless a horizon bifurcation or merger changes the singularity structure. Finally, several function plots are provided to illustrate the behavior of equilibrium and non-equilibrium free energy, the Kerr Newman temperature curve, and the entropy production law.

The Gamma-Ray Monitor (GRM) is a key scientific payload onboard the Space-based Multi-band Variable Object Monitor (SVOM) satellite, designed specifically for the detection and study of gamma-ray bursts (GRBs). Launched into a 625 km low-Earth orbit on 22 June 2024, GRM serves as a large-area, wide-field-of-view instrument capable of observing the hard X-ray and soft gamma-ray emissions in the energy range of 15 keV to 5 MeV. Its primary scientific objectives include: promptly triggering and localizing GRBs (with particular sensitivity to short-hard GRBs), measuring spectral and temporal properties of bursts, monitoring charged particle fluxes in orbit. GRM successfully detected its first GRB (GRB 240627B) on 27 June 2024, and has since maintained a detection rate of more than 100 GRBs per year. Cross-instrument comparisons with detectors such as GECAM and Fermi/GBM have validated the performance and data quality of GRM. This paper provides a comprehensive overview of GRM instrument design, reliability verification through ground testing, in-orbit triggering and localization algorithms, performance calibration, and preliminary in-orbit results, demonstrating its capability as a versatile gamma-ray all-sky monitor.

\texttt{gwNRHME} is a framework that converts multi-modal (i.e., containing several spherical harmonic modes) quasi-circular waveforms into their eccentric counterparts, provided the quadrupolar eccentric mode is known, by exploiting universal eccentric modulation functions. Leveraging this framework, we combine the quasi-circular NR surrogate model \texttt{NRHybSur3dq8} with the quadrupolar, non-spinning, eccentric surrogate \texttt{NRSurE\_q4NoSpin\_22} to construct a multi-modal, non-spinning, eccentric model, denoted as \model{}, which includes nine modes: $(2,\{1,2\})$, $(3,\{1,2,3\})$, $(4,\{2,3,4\})$, and $(5,5)$. When compared against 156 eccentric SXS NR waveforms, \model{} achieves median frequency-domain mismatches (computed using the Advanced LIGO design sensitivity) of $\sim 9\times 10^{-5}$, with a standard deviation of $\sim 2 \times 10^{-4}$. To demonstrate the modularity of the framework, we further combine \texttt{NRSurE\_q4NoSpin\_22} with effective-one-body (EOB) models \texttt{SEOBNRv5HM} and \texttt{TEOBResumS-Dali} in their non-spinning limits, yielding eccentric waveforms with median mismatches of $\sim 2\times10^{-4}$ and $\sim 10^{-3}$, respectively, with standard deviation of $\sim 2 \times 10^{-3}$ and $\sim 2 \times 10^{-2}$ respectively. Finally, we provide both a surrogate model, \texttt{gwEccEvolve\_q4NoSpin\_Sur}, and an analytical model, \texttt{gwEccEvNSv2}, for the eccentricity evolution up to $2M$ before merger, based on eccentricity definitions derived from the universal modulation functions. The \texttt{gwNRHME} framework is publicly available through the \texttt{gwModels} package, and the resulting waveform models will be released via the \texttt{gwsurrogate} package.