High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: Context. The study of Quasi-Periodic Oscillations (QPO) at low and high frequency in the variability of the high-energy emission from black-hole binaries and their physical interpretation in terms of signatures of General Relativity in the strong-field regime. Aims. To understand the nature of the 67 Hz QPOs observed in the X-ray emission of the peculiar black-hole binary GRS 1915+105 within the general classification of QPO and to determine the spin of the black hole in the system by applying the Relativistic Precession Model (RPM). Methods. Within the RPM, the only relativistic frequency that is stable in time over a large range of accretion rates and can be as low as 67 Hz (for a black-hole mass as measured dynamically) is the Lense-Thirring frequency at the Innermost Stable Circular Orbit (ISCO). In the application of the model, this corresponds to type-C QPOs. Under this assumption, it is possible to measure the spin of the black hole. We re-analysed a large number of RossiXTE observations to check whether other timing features confirm this hypothesis. Results. The identification of the 67 Hz QPO as the Lense-Thirring frequency at ISCO yields a value of 0.706 +/- 0.034 for the black hole spin. With this spin, the only two QPO detections at higher frequencies available in the literature are consistent with being orbital frequencies at a radius outside ISCO. The high-frequency bumps often observed at frequencies between 10 and 200 Hz follow the correlation expected for orbital and periastron-precession frequencies at even larger radii.

Abstract: Using the Zwicky Transient Facility (ZTF) and other observatories, we have identified three candidate Tidal Disruption Events (TDEs) in spatial and temporal coincidence with high-energy neutrinos detected by IceCube: AT2019dsg, AT2019fdr and AT2019aalc. All three of these events have been shown to be able to produce high-energy neutrinos. In these proceedings, I will give an overview of Tidal Disruption Events, outline our follow-up program with ZTF, describe the observations carried out for each of those coincident events and highlight their similarities and differences.

Abstract: The motion of the center of mass of a coalescing binary black hole (BBH) in a gravitational potential imprints a line-of-sight acceleration (LOSA) onto the emitted gravitational wave (GW) signal. The acceleration could be sufficiently large in dense stellar environments, such as globular clusters (GCs), to be detectable with next-generation space-based detectors. In this work, we use outputs of the \textsc{cluster monte carlo (cmc)} simulations of dense star clusters to forecast the distribution of detectable LOSAs in DECIGO and LISA eras. We study the effect of cluster properties -- metallicity, virial and galactocentric radii -- on the distribution of detectable accelerations, account for cosmologically-motivated distributions of cluster formation times, masses, and metallicities, and also incorporate the delay time between the formation of BBHs and their merger in our analysis. We find that larger metallicities provide a larger fraction of detectable accelerations by virtue of a greater abundance of relatively lighter BBHs, which allow a higher number of GW cycles in the detectable frequency band. Conversely, smaller metallicities result in fewer detections, most of which come from relatively more massive BBHs with fewer cycles but larger LOSAs. We similarly find correlations between the virial radii of the clusters and the fractions of detectable accelerations. Our work, therefore, provides an important science case for space-based GW detectors in the context of probing GC properties via the detection of LOSAs of merging BBHs.

Abstract: I explore a possibility to estimate an upper limit of the current iron abundance of the barion matter. The upper limit is determined by the minimal iron abundance, at which the gamma-ray background, produced by the decay of $^{56}$Ni synthesised in the Universe to date, contradicts the observational MeV gamma-ray background. I calculate the gamma-ray background from SNe~Ia and SNe~II with the gamma-ray scattering and absorption in supernova envelope. It is shown that the model background does not contradict the observed MeV background, if the present day iron abundance of the barion matter is less than 15\% of the solar abundance.

Abstract: The late-time spectra of the kilonova AT 2017gfo associated with GW170817 exhibit a strong emission line feature at $2.1\,{\rm \mu m}$. The line structure develops with time and there is no apparent blue-shifted absorption feature in the spectra, suggesting that this emission line feature is produced by electron collision excitation. We attribute the emission line to a fine structure line of Tellurium (Te) III, which is one of the most abundant elements in the second r-process peak. By using a synthetic spectral modeling including fine structure emission lines with the solar r-process abundance pattern beyond the first r-process peak, i.e., atomic mass numbers $A\gtrsim 88$, we demonstrate that [Te III] $2.10\,\rm \mu m$ is indeed expected to be the strongest emission line in the near infrared region. We estimate that the required mass of Te III is $\sim 10^{-3}M_{\odot}$, corresponding to the merger ejecta of $0.05M_{\odot}$, which is in agreement with the mass estimated from the kilonova light curve.

Abstract: Over a period of four active episodes between January 2021 and January 2022, the magnetar SGR J1935+2154 emitted a total of 343 bursts observed by the \textit{Fermi}/GBM and 82 bursts observed by GECAM-B. Temporal and spectral analyses reveal that the bursts have an average duration of 145 ms and a fluence ranging from $1.2 \times 10^{-9} \ \mathrm{erg \cdot cm^{-2}}$ to $4.1 \times 10^{-4} \ \mathrm{erg \cdot cm^{-2}}$ (8 - 200 keV). The spectral properties of these bursts are similar to those of earlier active episodes. Specifically, we find that the emission area of the Double Black Body (BB2) model shows a Log-Linear correlation to its temperature, and there is a weak relation between fluence and $E_{\mathrm{peak}}$/$\alpha$ in the CPL model. However, we note that the temperature distributions of BB2/BB models in GECAM-B are different from those in \textit{Fermi}/GBM, due to differences in the energy range used for fitting. To understand this difference, we propose a Multi-Temperature Black Body (MBB) model for analyzing thermal radiation, assuming that the BB temperatures follow a power law distribution. Our analysis shows the minimum temperature $kT_{\mathrm{min}} \sim 5$ keV of the MBB model is consistent between \textit{Fermi}/GBM and GECAM-B, and reveals the spectra of magnetar bursts tending to be soft, which may be composed of multiple BB components. The slope of the temperature distribution is steep which indicates that the majority of the BB temperatures are concentrated around the minimum temperature.

Abstract: With the breakthrough in PeV gamma-ray astronomy brought by the LHAASO experiment, high-energy sky is getting more completed than before. Lately LHAASO Collaboration reported the observation of a gamma-ray diffuse emission with energy up to the PeV level from both the inner and outer Galactic plane. In these spectra, there is one bump which is hard to explain by the conventional cosmic-ray transport scenarios. Therefore, we introduce two extra components corresponding to unresolved sources with exponential-cutoff-power-law (ECPL) spectral shape, one with index of 2.4, and 30 TeV cutoff energy, and another with index of 2.3 and 2 PeV cutoff energy. With our constructed model, we simulate the Galactic diffuse neutrino flux and find our results are in full agreement with latest IceCube Galactic plane search. We estimate the Galactic neutrino contribute of $\sim 9\%$ of astrophysical neutrinos at 20 TeV. In the high-energy regime, as expected most of neutrinos observed by IceCube should be from extra-galaxy.

Abstract: We report the discovery of Swift J221951-484240 (hereafter: J221951), a luminous slow-evolving blue transient that was detected by the Neil Gehrels Swift Observatory Ultra-violet/Optical Telescope (Swift/UVOT) during the follow-up of Gravitational Wave alert S190930t, to which it is unrelated. Swift/UVOT photometry shows the UV spectral energy distribution of the transient to be well modelled by a slowly shrinking black body with an approximately constant temperature of T~2.5x10^4 K. At a redshift z=0.5205, J221951 had a peak absolute magnitude of M_u,AB = -23 mag, peak bolometric luminosity L_max=1.1x10^45 erg s^-1 and a total radiated energy of E>2.6x10^52 erg. The archival WISE IR photometry shows a slow rise prior to a peak near the discovery date. Spectroscopic UV observations display broad absorption lines in N V and O VI, pointing toward an outflow at coronal temperatures. The lack of emission in the higher H~Lyman lines, N I and other neutral lines is consistent with a viewing angle close to the plane of the accretion or debris disc. The origin of J221951 can not be determined with certainty but has properties consistent with a tidal disruption event and the turn-on of an active galactic nucleus.

Abstract: We present a maximum likelihood (ML) algorithm that is fast enough to detect gamma-ray transients in real time on low-performance processors often used for space applications. We validate the routine with simulations and find that, relative to algorithms based on excess counts, the ML method is nearly twice as sensitive, allowing detection of 240-280% more short gamma-ray bursts. We characterize a reference implementation of the code, estimating its computational complexity and benchmarking it on a range of processors. We exercise the reference implementation on archival data from the Fermi Gamma-ray Burst Monitor (GBM), verifying the sensitivity improvements. In particular, we show that the ML algorithm would have detected GRB 170817A even if it had been nearly four times fainter. We present an ad hoc but effective scheme for discriminating transients associated with background variations. We show that the on-board localizations generated by ML are accurate, but that refined off-line localizations require a detector response matrix with about ten times finer resolution than is current practice. Increasing the resolution of the GBM response matrix could substantially reduce the few-degree systematic uncertainty observed in the localizations of bright bursts.

Abstract: We show that gas disks around the components of an orbiting binary system (so-called minidisks) may be susceptible to a resonant instability which causes the minidisks to become significantly eccentric. Eccentricity is injected by, and also induces, regular impacts between the minidisks at roughly the orbital period of the binary. Eccentric minidisks are seen in vertically integrated, two-dimensional simulations of a circular, equal-mass binary accreting from a circumbinary gas disk with a $\Gamma$-law equation of state. Minidisk eccentricity is suppressed by the use of an isothermal equation of state. However, the instability still operates, and can be revealed in a minimal disk-binary simulation by removing the circumbinary disk, and feeding the minidisks from the component positions. Minidisk eccentricity is also suppressed when the gravitational softening length is large ($\gtrsim 4\%$ of the binary semi-major axis), suggesting that its absence could be an artifact of widely adopted numerical approximations; a follow-up study in three dimensions with well-resolved, geometrically thin minidisks (aspect ratios $\lesssim 0.02$) may be needed to assess whether eccentric minidisks can occur in real astrophysical environments. If they can, the electromagnetic signature may be important for discriminating between binary and single black hole scenarios for quasi-periodic oscillations in active galactic nuclei, which may in turn aid in targeted searches with pulsar timing arrays for individual supermassive black hole binary sources of low-frequency gravitational waves.

Abstract: In this paper, we analyse sound waves arising from a cosmic phase transition where the full velocity profile is taken into account as an explanation for the gravitational wave spectrum observed by multiple pulsar timing array groups. Unlike the broken power law used in the literature, in this scenario the power law after the peak depends on the macroscopic properties of the phase transition, allowing for a better fit with pulsar timing array (PTA) data. We compare the best fit with that obtained using the usual broken power law and, unsurprisingly, find a better fit with the gravitational wave (GW) spectrum that utilizes the full velocity profile. We then discuss models that can produce the best-fit point and complementary probes using CMB experiments and searches for light particles in DUNE, IceCUBE-Gen2, neutrinoless double beta decay, and forward physics facilities at the LHC like FASER nu, etc.

Abstract: We present six epochs of optical spectropolarimetry of the Type II supernova (SN) 2023ixf ranging from $\sim$ 2 to 15 days after the explosion. Polarimetry was obtained with the Kast double spectrograph on the Shane 3 m telescope at Lick Observatory, representing the earliest such observations ever captured for an SN. We observe a high continuum polarization $p_{\text{cont}} \approx 1$ % on days +1.4 and +2.5 before dropping to 0.5 % on day +3.5, persisting at that level up to day +14.5. Remarkably, this change coincides temporally with the disappearance of highly ionized "flash" features. The decrease of the continuum polarization is accompanied by a $\sim 70^\circ$ rotation of the polarization position angle ($PA$) as seen across the continuum. The early evolution of the polarization may indicate different geometric configurations of the electron-scattering atmosphere as seen before and after the disappearance of the emission lines associated with highly-ionized species (e.g., He II, C IV, N III), which are likely produced by elevated mass loss shortly prior to the SN explosion. We interpret the rapid change of polarization and $PA$ from days +2.5 to +4.5 as the time when the SN ejecta emerge from the dense asymmetric circumstellar material (CSM). The temporal evolution of the continuum polarization and the $PA$ is consistent with an aspherical SN explosion that exhibits a distinct geometry compared to the CSM. The rapid follow-up spectropolarimetry of SN 2023ixf during the shock ionization phase reveals an exceptionally asymmetric mass-loss process leading up to the explosion.

Abstract: Magnetorotational instability (MRI)-driven turbulence and dynamo phenomena are analyzed using direct statistical simulations. Our approach begins by developing a unified mean-field model that combines the traditionally decoupled problems of the large-scale dynamo and angular-momentum transport in accretion disks. The model consists of a hierarchical set of equations, capturing up to the second-order cumulants, while a statistical closure approximation is employed to model the three-point correlators. We highlight the web of interactions that connect different components of stress tensors -- Maxwell, Reynolds, and Faraday -- through shear, rotation, correlators associated with mean fields, and nonlinear terms. We determine the dominant interactions crucial for the development and sustenance of MRI turbulence. Our general mean field model for the MRI-driven system allows for a self-consistent construction of the electromotive force, inclusive of inhomogeneities and anisotropies. Within the realm of large-scale magnetic field dynamo, we identify two key mechanisms -- the rotation-shear-current effect and the rotation-shear-vorticity effect -- that are responsible for generating the radial and vertical magnetic fields, respectively. We provide the explicit (nonperturbative) form of the transport coefficients associated with each of these dynamo effects. Notably, both of these mechanisms rely on the intrinsic presence of large-scale vorticity dynamo within MRI turbulence.