High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: To investigate the impact of matter mixing on the formation of molecules in the ejecta of SN 1987A, time-dependent rate equations for chemical reactions are solved for one-zone and one-dimensional ejecta models of SN 1987A. The latter models are based on the one-dimensional profiles obtained by angle-averaging of the three-dimensional hydrodynamical models, which effectively reflect the 3D matter mixing; the impact is demonstrated, for the first time, based on three-dimensional hydrodynamical models. The distributions of initial seed atoms and radioactive $^{56}$Ni influenced by the mixing could affect the formation of molecules. By comparing the calculations for spherical cases and for several specified directions in the bipolar-like explosions in the three-dimensional hydrodynamical models, the impact is discussed. The decay of $^{56}$Ni, practically $^{56}$Co at later phases, could heat the gas and delay the molecule formation. Additionally, Compton electrons produced by the decay could ionize atoms and molecules and could destruct molecules. Several chemical reactions involved with ions such as H$^+$ and He$^+$ could also destruct molecules. The mixing of $^{56}$Ni plays a non-negligible role in both the formation and destruction of molecules through the processes above. The destructive processes of carbon monoxide and silicon monoxide due to the decay of $^{56}$Ni generally reduce the amounts. However, if the molecule formation is sufficiently delayed under a certain condition, the decay of $^{56}$Ni could increase the amounts through a sequence of passes instead compared with the case with lower efficiencies for the destructive processes above.

Abstract: A comparative analysis of the individual bursts between FRB 20190520B and FRB 20121102A is presented by compiling a sample of bursts in multiple wavelengths. It is found that the peak frequency ($\nu_p$) distribution of the bursts of FRB 20190520B illustrates four discrete peaks in $\sim1-6$ GHz and their spectral width distribution can be fitted with a log-normal function peaking at 0.35 GHz. The discrete $\nu_p$ distribution and the narrow-banded spectral feature are analogous to FRB 20121102A. The burst duration of FRB 20190520B in the rest frame averages 10.72 ms, longer than that of FRB 20121102A by a factor 3. The specific energy ($E_{\rm \mu_{\rm c}}$) at 1.25 GHz of FRB 20190520B observed with the FAST telescope narrowly ranges in $[0.4, 1]\times 10^{38}$ erg, different from the bimodal $E_{\rm \mu_{\rm c}}$ distribution of FRB 20121102A. Assuming a Gaussian spectral profile of the bursts, our Monte Carlo simulation analysis suggests that a power-law (PL) or a cutoff power-law (CPL) energy function can comparably reproduce the $E_{\rm \mu_{\rm c}}$ distribution of FRB 20190520B. The derived energy function index of the PL model is $4.46\pm 0.17$, much steeper than that of FRB 20121102A ($1.82^{+0.10}_{-0.30}$). For the CPL model, we obtain an index of $0.47$ and a cutoff energy of $7.4\times 10^{37}$ erg. Regarding the predicted $\nu_p$ distribution in 1-2 GHz, the CPL model is more preferred than the PL model. These results indicate that FRB 20190520B and FRB 20121102A shares similar spectral properties, but their energy functions are intrinsically different.

Abstract: The observed spectral shapes variation and tentative bimodal burst energy distribution (E-distribution) of fast radio burst (FRB) 20121102A with the FAST telescope are great puzzles. Adopting the published multifrequency data observed with the FAST and Arecibo telescopes at $L$ band and the GBT telescope at $C$ band, we investigate these puzzles through Monte Carlo simulations. The intrinsic energy function (E-function) is modeled as $dp/dE\propto E^{-\alpha_{\rm E}}$, and the spectral profile is described as a Gaussian function. A fringe pattern of its spectral peak frequency ($\nu_{\rm p}$) in 0.5-8 GHz is inferred from the $\nu_{\rm p}$ distribution of the GBT sample. We estimate the likelihood of $\alpha_{\rm E}$ and the standard deviation of the spectral profile ($\sigma_{\rm s}$) by utilizing the Kolmogorov--Smirnov (K-S) test probability for the observed and simulated specific E-distributions. Our simulations yields $\alpha_{\rm E}=1.82^{+0.10}_{-0.30}$ and $\sigma_{\rm s}=0.18^{+0.28}_{-0.06}$ ($3\sigma$ confidence level) with the FAST sample. These results suggest that a single power-law function is adequate to model the E-function of FRB 20121102A. The variations of its observed spectral indices and E-distributions with telescopes in different frequency ranges are due to both physical and observational reasons, i.e. narrow spectral width for a single burst and discrete $\nu_{p}$ fringe pattern in a broad frequency range among bursts, and the selection effects of the telescope bandpass and sensitivity. The putative $\nu_{p}$ fringe pattern cannot be explained with the current radiation physics models of FRBs. Some caveats of possible artificial effects that may introduce such a feature are discussed.

Abstract: We study the evolution of a non-relativistically expanding thin shell in radio-emitting tidal disruption events (TDEs) based on a one-dimensional spherically symmetric model considering the effect of both a time-dependent mass loss rate of the disk wind and the ambient mass distribution. The analytical solutions are derived in two extreme limits: one is the approximate solution near the origin in the form of the Taylor series, and the other is the asymptotic solution in which the ambient matter is dominant far away from the origin. Our numerical solutions are confirmed to agree with the respective analytical solutions. We find that no simple power-law of time solution exists in early to middle times because the mass loss rate varies over time, affecting the shell dynamics. We also discuss the application of our model to the observed radio-emitting TDE, AT2019dsg.

Abstract: We present a relativistic disk-corona model for a steady state advective accretion disk to explain the UV to X-ray spectral index $\alpha_{\text{OX}}$ evolution of \textbf{four} tidal disruption event (TDE) sources XMMSL2J1446, XMMSL1J1404, XMMSL1J0740, \textbf{and AT2018fyk}. The viscous stress in our model depends on gas ($P_g$) and total ($P_t$) pressures as $\tau_{r\phi} \propto P_g^{\mu} P_t^{1-\mu}$, where $\mu$ is a constant. We compare various steady and time-dependent sub-Eddington TDE accretion models along with our disk-corona model to the observed $\alpha_{\text{OX}}$ of TDE sources and find that the disk-corona model agrees with the observations better than the other models. We find that $\mu$ is much smaller than unity for TDE sources XMMSL2J1446, XMMSL1J1404, and XMMSL1J0740. We also compare the relativistic model with a non-relativistic disk-corona model. The relativistic accretion dynamics reduce the spectral index relative to the non-relativistic accretion by increasing the energy transport to the corona. We estimate the mass accretion rate for all the sources and find that the observed luminosity follows a nearly linear relation with the mass accretion rate. The ratio of X-ray luminosity from the disk to the corona increases with the mass accretion rate. The observed $\alpha_{\text{OX}}$ shows positive and negative correlations with luminosity. The disk-corona model explains the negative correlation seen in the TDE sources XMMSL1J0740, XMMSL2J1446, and XMMSL1J1404. However, TDE AT2018fyk shows a positive correlation at higher luminosity and shows a better fit when a simple spherical adiabatic outflow model is added to the relativistic disk-corona model. Even though the disk luminosity dominates at a higher mass accretion rate, we show that the accretion models without a corona are unable to explain the observations, and the presence of a corona is essential.

Abstract: The paper presents the verification of previously published fast radio bursts (FRB) from the work of V.A. Fedorova and A.E. Rodin, detected in the monitoring data of the Large Phased Array (LPA) radio telescope using a search algorithm based on the convolution of data with a scattered pulse pattern. The same 6-channel data (channel width 415 kHz) were used for verification, in which FRBs were detected with dispersion measures of 247, 570 and 1767 pc/cm3. Additional verification of the published FRB was also carried out in 32-channel data (channel width 78 kHz). We have not been able to confirm any published FRB on the signal-to-noise ratios stated in the original paper. The main errors are caused by incorrect identification of the baseline and incorrect estimation of the standard deviations of noise.

Abstract: We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th resulted in a strong 13$\sigma$ detection with a differential photon spectral index, $\Gamma$ = 3.8 $\pm$ 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target of opportunity observations with Swift, optical photometry, polarimetry and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of $\tau_{obs}$ = 6.2 $\pm$ 0.9 hours, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cut-off. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model.

Abstract: Galactic cosmic ray transport relies on the existence of turbulence on scales comparable with the gyration radius of the particles and with wavenumber vector oriented along the local magnetic field. In the standard picture, in which turbulence is injected at large scales and cascades down to smaller scales, it is all but guaranteed that turbulence on the relevant scales may be present, either because of anisotropic cascading or because of the onset of damping processes. This raises questions on the nature of cosmic-ray scattering, especially at energies $\gtrsim 1$ TeV, where self-generation is hardly relevant. Here, by means of numerical simulations of charged test-particles in a prescribed magnetic field, we investigate particle diffusion in a situation in which turbulence is mainly present at large scales, with the possible presence of a smaller power on small scales, and discuss possible implications of this setup for cosmic-ray transport phenomenology.

Abstract: We examine the effect of a jet transversal structure from magnetohydrodynamic semi-analytical modelling on the total intensity profiles of relativistic jets from active galactic nuclei. In order to determine the conditions for forming double- and triple-peaked transverse intensity profiles, we calculate the radiative transfer for synchrotron emission with self-absorption from the jets described by the models with a constant angular velocity and with a total electric current closed inside a jet. We show that double-peaked profiles appear either in the models with high maximal Lorentz factors or in optically thick conditions. We show that triple-peaked profiles in radio galaxies constrain the fraction of the emitting particles in a jet. We introduce the possible conditions for triple-peaked profiles under the assumptions that nonthermal electrons are preferably located at the jet edges or are distributed according to Ohmic heating.

Abstract: Cosmic-ray (CR) antiparticles have the potential to reveal signatures of unexpected astrophysical processes and even new physics beyond the Standard Model. Recent CR detectors have provided accurate measurements of the positron flux, revealing the so-called positron excess at high energies. However, the uncertainties related to the modelling of the local positron flux are still very high, significantly affecting our models of positron emission from pulsars and current dark matter searches. In this work, we report a new set of cross sections for positron and electron production derived from the {\tt FLUKA} code. We compare them with the most extended cross-section data-sets and show the impact of neglecting the positron production from heavy CRs. Then, we review the most significant sources of uncertainties in our current estimations of the secondary positron flux at Earth and examine for the first time the impact of considering the spiral arm structure of the Galaxy in these estimations. Finally, we provide state-of-the-art predictions of the local positron flux and discuss the limitations of our dark matter searches with positrons and difficulties to determine the contribution from pulsars to the positron flux at low energies.