High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: Observations of elliptical galaxies suggest that black holes (BHs) might serve as dark energy candidates, coupled to the expansion of the Universe. According to this hypothesis, the mass of a BH could increase as the Universe expands. BH low-mass X-ray binaries (LMXBs) in the Galactic disk were born several gigayears ago, making the coupling effect potentially significant. In this work, we calculate the evolution of BH binaries with a binary population synthesis method to examine the possible influence of cosmologically-coupled growth of BHs, if it really exists. The measured masses of the compact objects in LMXBs show a gap around $\sim 2.5-5~{\rm M_\odot}$, separating the most massive neutron stars from the least massive BHs. Our calculated results indicate that, considering the mass growth seem to (partially) account for the mass gap and the formation of compact BH LMXBs, alleviating the challenges in modeling the formation and evolution of BH LMXBs with traditional theory. However, critical observational evidence like the detection of intermediate-mass black hole binaries is required to test this hypothesis.

Abstract: Previous theoretical works using the pre-merger orbital evolution of coalescing neutron stars to constrain properties of dense nuclear matter assume a gravitational wave phase uncertainty of a few radians, or about a half cycle. However, recent studies of the signal from GW170817 and next generation detector sensitivities indicate actual phase uncertainties at least twenty times better. Using these refined estimates, we show that future observations of nearby sources like GW170817 may be able to reveal neutron star properties beyond just radius and tidal deformability, such as the matter composition and/or presence of a superfluid inside neutron stars, via tidal excitation of g-mode oscillations. Data from GW170817 already limits the amount of orbital energy that is transferred to the neutron star to <2x10^47 erg and the g-mode tidal coupling to Qmode<10^-3 at 50 Hz (5x10^48 erg and 4x10^-3 at 200 Hz), and future observations and detectors will greatly improve upon these constraints. In addition, analysis using general parameterization models that have been applied to the so-called p-g instability show that the instability already appears to be restricted to regimes where the mechanism is likely to be inconsequential; in particular, we show that the number of unstable modes is <<100 at <~100 Hz, and next generation detectors will essentially rule out this mechanism (assuming that the instability remains undetected). Finally, we illustrate that measurements of tidal excitation of r-mode oscillations in nearby rapidly rotating neutron stars are within reach of current detectors and note that even non-detections will limit the inferred inspiralling neutron star spin rate to <20 Hz, which will be useful when determining other parameters such as neutron star mass and tidal deformability.

Abstract: The Extragalactic Background Light (EBL) that spans the UV-IR band originates from direct and dust-reprocessed starlight integrated over the history of the Universe. EBL measurements are very challenging due to foreground emission like the zodiacal light and interplanetary dust emission. Indeed, some optical/NIR direct measurements overpredict EBL models based on galaxy counts. On the other hand, there is some debate on possible additional components of the Optical-NIR photon density: e.g., population-III stars, axion-photon decay, direct collapse of black holes, intra-halo light etc. Owing to the absorption of Very High Energy (VHE) $\gamma$ rays by interaction with EBL photons, we study the prospects of accommodating an additional population of EBL sources in the Optical-NIR band on top of the standard galaxy-count--based component. To this aim we use 105 VHE spectra of 37 blazars with known redshifts, $0.03<z<0.94$. We correct the observed spectra for absorption by our model EBL. By requiring the intrinsic spectra to be non-concave and with a VHE spectral index $>$1.5, we estimate, at different wavelengths, upper limits to the additional low energy photon fields which would contribute to the absorption of $\gamma$-rays. Considering these limits we suggest that there is room for photons from Pop III stars and axion-like particle (ALP) annihilation. However, these additional hypothetical photon fields are bound to fall significantly below direct published EBL measurements by several instruments, and therefore our limits are either in tension or even inconsistent with such measurements.

Abstract: 4U 1543-47 underwent its brightest outburst in 2021 after two decades of inactivity. During its decay phase, AstroSat conducted nine observations of the source spanning from July $1^{st}$ to September $26^{th}$, 2021. The first three observations were performed with an offset of 40 arcmin with AstroSat/LAXPC, while the remaining six were on-axis observations. In this report, we present a comprehensive spectral analysis of the source as it was in the High/Soft state during the entire observation period. The source exhibited a disk-dominated spectra with a weak high-energy tail (power-law index $\geq2.5$) and a high inner disk temperature ($\sim$0.84 keV). Modelling the disk continuum with non-relativistic and relativistic models, we find inner radius to be significantly truncated at $>$$10~R_g$. Alternatively, to model the spectral evolution with the assumption that the inner disk is at the ISCO, it is necessary to introduce variation in the spectral hardening in the range $\sim$1.5-1.9.

Abstract: A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at $\sim$7 MeV cm$^{-2}$ (time period 1000~s) and $\sim$35 MeV cm$^{-2}$ (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total gravitational wave energy. These are the first limits on UHE photons from GW sources.

Abstract: We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics and gravity, and non-relativistic (NR) hydrodynamics and gravity, in an idealized model setting, we vary the initial radius of the shock and, by varying its mass and radius in concert, the proto-neutron star (PNS) compactness. We investigate two regimes expected in a post-bounce core-collapse supernova (CCSN): one meant to resemble a relatively low-compactness configuration and one meant to resemble a relatively high-compactness configuration. We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the high-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the high-compactness suite. We discuss implications of our results for CCSN simulations.

Abstract: The prediction of spin magnitudes in binary black hole and neutron star mergers is crucial for understanding the astrophysical processes and gravitational wave (GW) signals emitted during these cataclysmic events. In this paper, we present a novel Neuro-Symbolic Architecture (NSA) that combines the power of neural networks and symbolic regression to accurately predict spin magnitudes of black hole and neutron star mergers. Our approach utilizes GW waveform data obtained from numerical relativity simulations in the SXS Waveform catalog. By combining these two approaches, we leverage the strengths of both paradigms, enabling a comprehensive and accurate prediction of spin magnitudes. Our experiments demonstrate that the proposed architecture achieves an impressive root-mean-squared-error (RMSE) of 0.05 and mean-squared-error (MSE) of 0.03 for the NSA model and an RMSE of 0.12 for the symbolic regression model alone. We train this model to handle higher-order multipole waveforms, with a specific focus on eccentric candidates, which are known to exhibit unique characteristics. Our results provide a robust and interpretable framework for predicting spin magnitudes in mergers. This has implications for understanding the astrophysical properties of black holes and deciphering the physics underlying the GW signals.

Abstract: The discovery of ultra-high-energy (UHE) neutrinos has the potential to offer unique insight into fundamental questions. To capitalize on the upcoming opportunity provided by new UHE neutrino telescopes, we provide state-of-the-art forecasts of the discovery of a diffuse flux of UHE neutrinos over the next 10-20 years, focusing on neutrino radio-detection in the planned IceCube-Gen2 detector. We use state-of-the-art flux predictions and detector modeling. We find that, even under conservative analysis choices, most benchmark UHE neutrino flux models from the literature may be discovered within 10 years of detector exposure, with many discoverable sooner, and may be distinguished from each other. Our results demonstrate the transformative potential of next-generation UHE neutrino telescopes.

Abstract: GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta.