High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: We present optical photometric and spectroscopic observations of the 02es-like type Ia supernova (SN) 2022ywc. The transient occurred in the outskirts of an elliptical host galaxy and showed a striking double-peaked light curve with an early excess feature detected in the ATLAS orange and cyan bands. The early excess is remarkably luminous with an absolute magnitude $\sim -19$, comparable in luminosity to the subsequent radioactively-driven second peak. The spectra resemble the hybrid 02es-like SN 2016jhr, that is considered to be a helium shell detonation candidate. We investigate different physical mechanisms that could power such a prominent early excess and rule out massive helium shell detonation, surface $^{56}$Ni distribution and ejecta-companion interaction. We conclude that SN ejecta interacting with circumstellar material (CSM) is the most viable scenario. Semi-analytical modelling with MOSFiT indicates that SN ejecta interacting with $\sim 0.05\,$M$_{\odot}$ of CSM at a distance of $\sim 10^{14}$ cm can explain the extraordinary light curve. A double-degenerate scenario may explain the origin of the CSM, either by tidally-stripped material from the secondary white dwarf, or disk-originated matter launched along polar axes following the disruption and accretion of the secondary white dwarf. A non-spherical CSM configuration could suggest that a small fraction of 02es-like events viewed along a favourable line of sight may be expected to display a very conspicuous early excess like SN 2022ywc.

Abstract: Neutron stars accreting from massive binary companions come in a wide range of types. Systems with an OB supergiant donor are often divided between persistently and transiently accreting systems, respectively called Supergiant X-ray Binaries (SgXBs) and Supergiant Fast X-ray Transients (SFXTs). The origin of this dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, the binary orbit, or magnetic field configuration, but direct observational evidence for these hypotheses remains sparse. Here, we present the results of a pilot exploration of a novel approach to this long-standing question, turning to the mm band to probe the outer regions of the stellar wind beyond the binary orbit. Specifically, we present 100-GHz NOEMA observations of a SgXB, X1908+075, and a SFXT, IGR J18410-0535. We detect the SFXT as a point source at $63.4 \pm 9.6$ $\mu$Jy, while the SgXB is not detected. The spectrum of IGR J18410-0535 is constrained to be flat or inverted by comparing with quasi-simultaneous $5.5$+$9$ GHz radio observations, ruling out non-thermal flaring and consistent with thermal wind emission. Additional X-ray measurements further constrain the wind mass loss rate and velocity of the SgXB. We compare our targets with each other and earlier wind estimates, and reflect on future opportunities using this novel observational approach to characterize stellar winds in X-ray binaries.

Abstract: The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetised fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field amplification, or the transport of angular momentum. Sub-grid models offer a means to overcome those difficulties. In a recent paper we presented MInIT, an MHD-instability-induced-turbulence mean-field, sub-grid model based on the modelling of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors. While in our previous work MInIT was assessed within the framework of the magnetorotational instability, in this paper we further evaluate the model in the context of the Kelvin-Helmholtz instability (KHI). The main difference with other sub-grid models (as e.g. the alpha-viscosity model or the gradient model) is that in MInIT we track independently the turbulent energy density at sub-grid scales, which is used, via a simple closure relation, to compute the different turbulent stresses relevant for the dynamics. The free coefficients of the model are calibrated using well resolved box simulations of magnetic turbulence generated by the KHI. We test the model against these simulations and show that it yields order-of-magnitude accurate predictions for the evolution of the turbulent Reynolds and Maxwell stresses.

Abstract: The detection of gamma-ray emission from accreting pulsars in X-ray binaries (XRBs) has long been sought after. For some high-mass X-ray binaries (HMXBs), marginal detections have recently been reported. Regardless of whether these will be confirmed or not, future telescopes operating in the gamma-ray band could offer the sensitivity needed to achieve solid detections and possibly spectra. In view of future observational advances, we explored the expected emission above 10 GeV from XRBs, based on the Cheng & Ruderman model, where gamma-ray photons are produced by the decay of pion-0 originated by protons accelerated in the magnetosphere of an accreting pulsar fed by an accretion disc. We improved this model by considering, through Monte Carlo simulations, the development of cascades inside of and outside the accretion disc, taking into account pair and photon production processes that involve interaction with nuclei, X-ray photons from the accretion disc, and the magnetic field. We produced grids of solutions for different input parameter values of the X-ray luminosity (L_x), magnetic field strength (B), and for different properties of the region where acceleration occurs. We found that the gamma-ray luminosity spans more than five orders of magnitude, with a maximum of ~1E35 erg/s. The gamma-ray spectra show a large variety of shapes: some have most of the emission below ~100 GeV, others are harder (emission up to 10-100 TeV). We compared our results with Fermi/LAT and VERITAS detections and upper-limits of two HMXBs: A0535+26 and GROJ1008-57. More consequential comparisons will be possible when more sensitive instruments will be operational in the coming years.

Abstract: Since the release of the Gravitational Wave Transient Catalogue GWTC-2.1 by the LIGO-Virgo collaboration, sub-threshold gravitational wave (GW) candidates are publicly available. They are expected to be released in real-time as well, in the upcoming O4 run. Using these GW candidates for multi-messenger studies complement the ongoing efforts to identify neutrino counterparts to GW events. This in turn, allows us to schedule electromagnetic follow-up searches more efficiently. However, the definition and criteria for sub-threshold candidates are pretty flexible. Finding a multi-messenger counterpart via archival studies for these candidates will help to set up strong bounds on the GW parameters which are useful for defining a GW signal as sub-threshold, thereby increasing their significance for scheduling follow-up searches. Here, we present the current status of this ongoing work with the IceCube Neutrino Observatory. We perform a selection of the sub-threshold GW candidates from GWTC-2.1 and conduct an archival search for sub-TeV neutrino counterparts detected by the dense infill array of the IceCube Neutrino Observatory, known as "DeepCore". For this, an Unbinned Maximum Likelihood (UML) method is used. We report the 90% C.L. sensitivities of this sub-TeV neutrino dataset for each selected sub-threshold GW candidate, considering the spatial and temporal correlation between the GW and neutrino events within a 1000 s time window.

Abstract: The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute time scales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection layer, is a candidate solution to this problem. In this work, we employ state-of-the-art particle-in-cell simulations in a statistical comparison with observations of a flaring episode of a well-known blazar, Mrk 421, at very high energy (VHE, E > 100 GeV). We tested the predictions of our model by generating simulated VHE light curves that we compared quantitatively with methods that we have developed for a precise evaluation of theoretical and observed data. With our analysis, we can constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma, viewing angle and the reconnection layer orientation in the blazar jet. Our analysis favours parameter spaces with magnetic field strength 0.1 G, rather large viewing angles (6-8 degrees), and misaligned layer angles, offering a strong candidate explanation for the Doppler crisis often observed in the jets of high synchrotron peaking blazars.

Abstract: Short-duration gamma-ray bursts (sGRBs) are commonly attributed to the mergers of double neutron stars (NSs) or the mergers of a neutron star with a black hole (BH). While the former scenario was confirmed by the event GW170817, the latter remains elusive. Here, we consider the latter scenario in which, a NS is tidally disrupted by a fast spinning low-mass BH and the accretion onto the BH launches a relativistic jet and hence produces a sGRB. The merging binary's orbit is likely misaligned with the BH's spin. Hence, the Lense-Thirring precession around the BH may cause a hyper-accreting thick disk to precess in a solid-body manner. We propose that a jet, initially aligned with the BH spin, is deflected and collimated by the wind from the disk, therefore being forced to precess along with the disk. This would result in a quasi-periodic oscillation or modulation in the gamma-ray light curve of the sGRB, with a quasi-period of $\sim 0.01-0.1$ s. The appearance of the modulation may be delayed respective to the triggering of the light curve. This feature, unique to the BH-NS merger, may have already revealed itself in a few observed sGRBs (such as GRB 130310A), and it carries the spin-obit orientation information of the merging system. Identification of this feature would be a new approach to reveal spin-orbit-misaligned merging BH-NS systems, which are likely missed by the current gravitational-wave searching strategy principally targeting aligned systems.

Abstract: X-ray binaries (XRBs) are thought to regulate cosmic thermal and ionisation histories during the Epoch of Reionisation and Cosmic Dawn ($z\sim 5-30$). Theoretical predictions of the X-ray emission from XRBs are important for modeling such early cosmic evolution. Nevertheless, the contribution from Be-XRBs, powered by accretion of compact objects from decretion disks around rapidly rotating O/B stars, has not been investigated systematically. Be-XRBs are the largest class of high-mass XRBs (HMXBs) identified in local observations and are expected to play even more important roles in metal-poor environments at high redshifts. In light of this, we build a physically motivated model for Be-XRBs based on recent hydrodynamic simulations and observations of decretion disks. Our model is able to reproduce the observed population of Be-XRBs in the Small Magellanic Cloud with appropriate initial conditions and binary stellar evolution parameters. We derive the X-ray output from Be-XRBs as a function of metallicity in the (absolute) metallicity range $Z\in [10^{-4},0.03]$. We find that Be-XRBs can contribute a significant fraction ($\sim 60\%$) of the total X-ray budget from HMXBs observed in nearby galaxies for $Z\sim 0.0003-0.02$. A similar fraction of observed ultra-luminous ($\gtrsim 10^{39}\ \rm erg\ s^{-1}$) X-ray sources can also be explained by Be-XRBs. Moreover, the predicted metallicty dependence in our fiducial model is consistent with observations, showing a factor of $\sim 8$ increase in X-ray luminosity per unit star formation rate from $Z=0.02$ to $Z=0.0003$.

Abstract: The spin parameter of the black hole in the accreting X-ray binary LMC X-1 has been measured in a number of studies to be $a_*\gtrsim 0.9$. These measurements were claimed to take into account both statistical and systematic (model-dependent) uncertainties. We perform new measurements using a recent simultaneous observation of LMC X-1 by NICER and NuSTAR, providing a data set of very high quality. We use the disk continuum method together with improved models for coronal Comptonization. With the standard relativistic blackbody disk model and optically thin Comptonization, we obtain values of $a_*$ similar to those obtained before. We then consider modifications to the standard disk model. Using a phenomenological color correction set to 2, we find lower values of $a_*\approx 0.64$--0.84. We then consider disks thicker than the standard one, i.e., with some dissipation in surface layers, as expected if partially supported by magnetic pressure. To account for that, we assume the disk is covered by a warm and optically thick layer, Comptonizing the emission of the underlying disk. Our model with the lowest $\chi^2$ yields then a low range of the spin, $a_*\approx 0.40^{+0.41}_{-0.32}$. That last model is also in agreement with the inverse disk temperature-luminosity relation found in this source. We conclude that determinations of the spin using the continuum method is highly sensitive to the assumptions about the disk structure.

Abstract: Recent associations of high-energy neutrinos with active galactic nuclei (AGN) have revived the interest in leptohadronic models of radiation from astrophysical sources. The rapid increase in the amount of acquired multi-messenger data will require soon fast numerical models that may be applied to large source samples. We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed. LeHaMoC solves the Fokker-Planck equations of photons and relativistic particles (i.e. electrons, positrons, protons, and neutrinos) produced in a homogeneous magnetized source that may also be expanding. The code utilizes a fully implicit difference scheme that allows fast computation of steady-state and dynamically evolving physical problems. We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHE$\nu$A, a well-tested code of similar philosophy but different numerical implementation. We find a good agreement (within 10-30%) with the numerical results obtained with ATHE$\nu$A without evidence of systematic differences. We then demonstrate the capabilities of the code through illustrative examples. First, we fit the spectral energy distribution from a jetted AGN in the context of a synchrotron-self Compton model and a proton-synchrotron model using Bayesian inference. Second, we compute the high-energy neutrino signal and the electromagnetic cascade induced by hadronic interactions in the corona of NGC 1068. LeHaMoC is easily customized to model a variety of high-energy astrophysical sources and has the potential to become a widely utilized tool in multi-messenger astrophysics.