# High Energy Astrophysical Phenomena (astro-ph.HE)

• Gamma-ray burst (GRB) jets are narrow, and thus typically point away from us. They are initially ultra-relativistic, causing their prompt $\gamma$-ray and early afterglow emission to be beamed away from us. However, as the jet gradually decelerates its beaming cone widens and eventually reaches our line of sight and the afterglow emission may be detected. Such orphan afterglows were not clearly detected so far. Nevertheless, they should be detected in upcoming optical or radio surveys, and it would be challenging to clearly distinguish between them and other types of transients. Therefore, we perform detailed, realistic calculations of the expected afterglow emission from GRB jets viewed at different angles from the jet's symmetry axis. The dynamics are calculated using 2D relativistic hydrodynamics simulations of jets propagating into different power-law external density profiles, $\rho_{\rm ext}\propto{}R^{-k}$ for $k=0,\,1,\,1.5,\,2$, ranging from a uniform ISM-like medium ($k=0$) to a stratified steady stellar-wind like profile ($k=2$). We calculate radio, optical and X-ray lightcurves, and the evolution of the radio afterglow image size, shape and flux centroid. This may help identify misaligned relativistic jets, whether initially ultra-relativistic and producing a GRB for observers within their beam, or (possibly intrinsically more common) moderately relativistic, in either (i) nearby supernovae Ib/c (some of which are associated with long duration GRBs), or (ii) in binary neutron star mergers, which may produce short duration GRBs, and may also be detected in gravitational waves (e.g. GW$\,$170827/GRB$\,$170817A with a weak prompt $\gamma$-ray emission may harbor an off-axis jet).
• The X-ray to radio afterglow emission of GRB 170817A/GW 170817 so far scales as $F_\nu\propto\nu^{-0.6}t^{0.8}$ with observed frequency and time, consistent with a single power-law segment of the synchrotron spectrum from the external shock going into the ambient medium. This requires the effective isotropic equivalent afterglow shock energy in the visible region to increase as $\sim t^{1.7}$. The two main channels for such an energy increase are (i) \emphradial: more energy carried by slower material (in the visible region) gradually catches up with the afterglow shock and energizes it, and (ii) \emphangular: more energy in relativistic outflow moving at different angles to our line of sight, whose radiation is initially beamed away from us but its beaming cone gradually reaches our line of sight as it decelerates. One cannot distinguish between these explanations (or combinations of them) using only the X-ray to radio $F_\nu(t)$. Here we demonstrate that the most promising way to break this degeneracy is through afterglow imaging and polarization, by calculating the predicted evolution of the afterglow image (its size, shape and flux centroid) and linear polarization $\Pi(t)$ for different angular and/or radial outflow structures that fit $F_\nu(t)$. We consider two angular profiles -- a Gaussian and a narrow core with power-law wings in energy per solid angle, as well as a (cocoon motivated) (quasi-) spherical flow with radial velocity profile. For a jet viewed off-axis (and a magnetic field produced in the afterglow shock) $\Pi(t)$ peaks when the jet's core becomes visible, at $\approx2t_p$ where the lightcurve peaks at $t_p$, and the image can be elongated with aspect ratios$\;\gtrsim2$. A quasi-spherical flow has an almost circular image and a much lower $\Pi(t)$ (peaking at $\approx t_p$) and flux centroid displacement.
• We currently know about 30 magnetars: seemingly isolated neutron stars whose properties can be (in part) comprehended only acknowledging that they are endowed with magnetic fields of complex morphology and exceptional intensity-at least in some components of the field structure. Although magnetars represent only a small percentage of the known isolated neutron stars, there are almost certainly many more of them, since most magnetars were discovered in transitory phases called outbursts, during which they are particularly noticeable. In outburst, in fact, a magnetar can be brighter in X-rays by orders of magnitude and usually emit powerful bursts of hard-X/soft-gamma-ray photons that can be detected almost everywhere in the Galaxy with all-sky monitors such as those on board the Fermi satellite or the Neil Gehrels Swift Observatory. Magnetars command great attention because the large progress that has been made in their understanding is proving fundamental to fathom the whole population of isolated neutron stars, and because, due to their extreme properties, they are relevant for a vast range of different astrophysical topics, from the study of gamma-ray bursts and superluminous supernovae, to ultraluminous X-ray sources, fast radio bursts, and even to sources of gravitational waves. Several excellent reviews with different focuses were published on magnetars in the last few years: among others, Israel and Dall'Osso (2011); Rea and Esposito (2011); Turolla and Esposito (2013); Mereghetti et al. (2015); Turolla et al. (2015); Kaspi and Beloborodov (2017). Here, we quickly recall the history of these sources and travel through the main observational facts, trying to touch some recent and sometimes little-discussed ramifications of magnetars.
• We demonstrate that the shower-to-shower fluctuations of the muon content of extensive air showers correlate with the fluctuations of a variable of the first interaction of Ultra High Energy Cosmic Rays, which is computed from the fraction of energy carried by the hadrons that sustain the hadronic cascade. The influence of subsequent stages of the shower development is found to play a sub-dominant role. As a consequence, the shower-to-shower distribution of the muon content is a direct probe of the hadron energy spectrum of interactions beyond 100 TeV center of mass energies.
• We report a new Fast Radio Burst (FRB) discovered in real-time as part of the UTMOST project at the Molonglo Observatory Synthesis Radio Telescope (MOST). FRB170827 is the first detected with our low-latency ($< 24$ s), machine-learning-based FRB detection system. The FRB discovery was accompanied by the capture of voltage data at the native time and frequency resolution of the observing system, enabling coherent dedispersion and detailed off-line analysis, which have unveiled fine temporal and frequency structure. The dispersion measure (DM) of 176.80 $\pm$ 0.04 pc cm$^{-3}$, is the lowest of the FRB population. The Milky Way contribution along the line of sight is $\sim$ 40 pc cm$^{-3}$, leaving an excess DM of $\sim$ 145 pc cm$^{-3}$. The FRB has a fluence $>$ 20 $\pm$ 7 Jy ms. The burst's width at 10$\%$ maximum is 400 $\mu$s, making it one of the narrowest FRBs. The burst shows three temporal components, the narrowest of which is $\sim$ 30 $\mu$s, and a scattering of $4.1 \pm 2.7$ $\mu$s. The FRB shows spectral modulations on frequency scales of 1.5 MHz and 0.1 MHz. Both are prominent in the dynamic spectrum, which shows a bright region of emission between 841 and 843 MHz, and weaker, patchy emission across the entire band. We show the fine spectral structure could arise in the FRB host galaxy, or its immediate vicinity.
• We have explored the capabilities of dust extinction and $\gamma$ rays to probe the properties of the interstellar medium in the nearby anti-centre region. We have jointly modelled the $\gamma$-ray intensity and the stellar reddening, E(B-V) as a combination of H$_{\rm I}$-bright, CO-bright, and ionised gas components. The complementary information from dust reddening and $\gamma$ rays is used to reveal the dark gas not seen, or poorly traced, by H$_{\rm I}$, free-free, and $^{12}$CO emissions. We compare the total gas column densities, $N_{\rm{H}}$, derived from the $\gamma$ rays and stellar reddening with those inferred from a similar analysis (Remy et al. 2017) of $\gamma$ rays and of the optical depth of the thermal dust emission, $\tau_{353}$, at 353 GHz. We can therefore compare environmental variations in specific dust reddening, E(B-V)/$N_{\rm H}$, and in dust emission opacity (dust optical depth per gas nucleon), $\tau_{353}/N_{\rm{H}}$. Over the whole anti-centre region, we find an average E(B-V)/$N_{\rm H}$ ratio of $(2.02\pm0.48)\times$ $10^{-22}$~mag~cm$^2$, with maximum local variations of about $\pm30\%$ at variance with the two to six fold coincident increase seen in emission opacity as the gas column density increases. In the diffuse medium, the small variations in specific reddening, E(B-V)/$N_{\rm H}$ implies a rather uniform dust-to-gas mass ratio in the diffuse parts of the anti-centre clouds. The small amplitude of the E(B-V)/$N_{\rm H}$ variations with increasing $N_{\rm{H}}$ column density confirms that the large opacity $\tau_{353}/N_{\rm{H}}$ rise seen toward dense CO clouds is primarily due to changes in dust emissivity. The environmental changes are qualitatively compatible with model predictions based on mantle accretion on the grains and the formation of grain aggregates.
• Cyclic sedimentation has varied at several timescales and this variability has been geologically well documented at Milankovitch timescales, controlled in part by climatically (insolation) driven sea-level changes. At the longer (tens of Myr) timescales connection between astronomical parameters and sedimentation via cyclic solar-system motions within the Milky Way has also been proposed, but this hypothesis remains controversial because of the lack of long geological records. The absence of a physical mechanism that could explain the connection between climate and astronomy at these longer timescales led to the explanation of plate motions as the main driver of climate on Earth. Here we statistically show a prominent and persistent ~36 Myr sedimentary cyclicity superimposed on two megacycles (~250 Myr) in a relatively well-constrained sea-level (SL) record of the past 542 Myr (Phanerozoic eon). Given the possible link between amplitudes of the ~36 and ~250 Myr cyclicities in SL record and the potential that these periodicities fall into the frequency band of solar system motions, we suggest an astronomical origin, and model these periodicities as originating from the path of the solar system in the Milky Way as vertical and radial periods that modulate the flux of cosmic rays on Earth. Our finding of the ~36 Myr SL cyclicity lends credibility to the existing hypothesis about the imprint of solar-system vertical period on the geological record. The ~250 Myr megacycles are tentatively attributed to a radial period. However, the tectonic drivers also remain potentially plausible. The potential existence of a correlation between the modeled astronomical signal and the geological record may offer an indirect proxy to understand the structure and history of the Milky Way by providing a 542 Myr long record of the path of the Sun in our Galaxy.
• Non-thermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection, where the magnetization $\sigma$, defined as the ratio of magnetic energy density to enthalpy density, is $\sim 1$, is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. We focus on the dependence of the electron energy spectrum on $\sigma$ and the proton $\beta$ (i.e., the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is non-thermal and can be generally modeled as a power law. At $\beta \lesssim 3 \times 10^{-3}$, the slope, $p$, is independent of $\beta$ and it hardens with increasing $\sigma$ as $p\simeq 1.8 +0.7/\sqrt{\sigma}$. Electrons are primarily accelerated by the non-ideal electric field at X-points, either in the initial current layer or in current sheets generated in between merging magnetic islands. At higher values of $\beta$, the electron power law steepens for all values of $\sigma$. At values of $\beta$ near $\beta_{\rm max}\approx1/4\sigma$, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. These particles are accelerated via a Fermi-like process by bouncing in between the reconnection outflow and a stationary magnetic island. We provide an empirical prescription for the dependence of the power-law slope and the acceleration efficiency on $\beta$ and $\sigma$, which can be used in global simulations of collisionless accretion disks.
• Dwarf galaxies represent a powerful probe of annihilating dark matter particle models, with gamma-ray data setting some of the best bounds available. A major issue in improving over existing constraints consists in the limited knowledge of the astrophysical background (mostly diffuse photons, but also unresolved sources). Perhaps more worrisome, several approaches in the literature suffer of the difficulty of assessing the systematic error due to background mis-modelling. Here we propose a data-driven method to estimate the background at the dwarf position and its uncertainty, relying on an appropriate use of the whole-sky data, via an optimisation procedure of the interpolation weights. While this article is mostly methodologically oriented, we also report the bounds based on latest Fermi-LAT data and updated information for J-factors for both isolated and stacked dwarfs. Our results are very competitive with the Fermi-LAT ones, while being derived with a more general and flexible method. We discuss the impact of profiling over the J-factor as well as over the background probability distribution function, with the latter resulting for instance crucial in drawing conclusions of compatibility with DM interpretations of the so-called Galactic Centre Excess.
• We address the significance of the observed GeV emission from \textitFermi-LAT on the understanding of the structure of long GRBs. We examine 82 X-ray Flashs (XRFs), in none of them GeV radiation is observed, adding evidence to the absence of a black hole (BH) formation in their merging process. By examining $329$ Binary-driven Hypernovae (BdHNe) we find that out of $48$ BdHNe observable by \textitFermi-LAT in \textitonly $21$ of them the GeV emission is observed. The Gev emission in BdHNE follows a universal power-law relation between the luminosity and time, when measured in the rest frame of the source. The power-law index in BdHNe is of $-1.20 \pm 0.04$, very similar to the one discovered in S-GRBs, $-1.29 \pm 0.06$. The GeV emission originates from the newly-born BH and allows to determine its mass and spin. We further give the first evidence for observing a new GRB subclass originating from the merging of a hypernova (HN) and an already formed BH binary companion. We conclude that the GeV emission is a necessary and sufficient condition to confirm the presence of a BH in the hypercritical accretion process occurring in a HN. The remaining $27$ BdHNe, recently identified as sources of flaring in X-rays and soft gamma-rays, have no GeV emission. From this and previous works, we infer that the observability of the GeV emission in some BdHNe is hampered by the presence of the HN ejecta. We conclude that the GeV emission can only be detected when emitted within a half-opening angle $\approx$60$^{\circ}$ normal to the orbital plane of the BdHN.
• Enhanced angular momentum transfer through the boundary layer near the surface of weakly magnetised accreting star is required in order to explain the observed accretion timescales in low-mass X-ray binaries, cataclysmic variables or young stars with massive protoplanetary discs. Accretion disc boundary layer is locally represented by incompressible homogeneous and boundless flow of the cyclonic type, which is linearly stable. Its non-linear instability at the shear rates of order of the rotational frequency remains an issue. We argue that hydrodynamical subcritical turbulence in such a flow is sustained by the non-linear feedback from essentially three-dimensional vortices, which are generated by quasi-two-dimensional trailing shearing spirals grown to high amplitude via the swing amplification. We refer to those three-dimensional vortices as cross-rolls, since they are aligned in the shearwise direction in contrast to streamwise rolls generated by the anti-lift-up mechanism in rotating shear flow on the Rayleigh line. Transient growth of cross-rolls is studied analytically and further confronted with direct numerical simulations (DNS) of dynamics of non-linear erturbations. DNS performed in a tall box show that transition Reynolds number $R_T$ as function of shear rate accords with the line of constant maximum transient growth of cross-rolls. The transition in the tall box has been observed until the shear rate three times higher than the rotational frequency, when $R_T\sim 50000$. Assuming that the cross-rolls are also responsible for turbulence in the Keplerian flow, we estimate $R_T\lesssim 10^8$ in this case. The most favourable shear rate to test the existence of turbulence in the quasi-Keplerian regime may be sub-Keplerian and equal approximately to $1/2$.
• Graham et al. (2015a) reported a periodically varying quasar and supermassive black hole binary candidate, PG1302-102 (hereafter PG1302), which was discovered in the Catalina Real-Time Transient Survey (CRTS). Its CRTS optical light curve is well fitted to a sinusoid of an observed period of $\approx 1,884$ days and well modeled by the relativistic Doppler boosting of the secondary mini-disk (D'Orazio et al. 2015). However, the CRTS light curve from MJD $\approx 53500$ to MJD $\approx 56400$ covers only ~1.5 cycles of periodic variation, a short baseline that can be highly susceptible to normal, stochastic quasar variability (Vaughan et al. 2016). In this Letter, we present a re-analysis of PG1302, using the latest light curve from the All-Sky Automated Survey for Supernovae (ASAS-SN), which extends the observational baseline to the present day (MJD $\approx 58200$), and adopting a maximum likelihood method which searches for a periodic component in addition to stochastic quasar variability. The extended ASAS-SN light curve does not support the persistence of its periodicity, and thus the suggested binary model of PG1302 is disfavored.
• The quasinormal resonant modes of massless neutral fields in near-extremal Kerr-Newman-de Sitter black-hole spacetimes are calculated in the eikonal regime. It is explicitly proved that, in the angular momentum regime ${\bar a}> \sqrt{{{1-2{\bar\Lambda}}\over{4+{\bar\Lambda}/3}}}$, the black-hole spacetimes are characterized by slowly decaying resonant modes which are described by the compact formula $\Im\omega(n)=\kappa_+\cdot(n+{1\over2})$ [here the physical parameters $\{{\bar a},\kappa_+,{\bar\Lambda},n\}$ are respectively the dimensionless angular momentum of the black hole, its characteristic surface gravity, the dimensionless cosmological constant of the spacetime, and the integer resonance parameter]. Our results support the validity of the Penrose strong cosmic censorship conjecture in these black-hole spacetimes.
• The observed multi-GeV gamma-ray emission from the solar disk --- sourced by hadronic cosmic rays interacting with gas, and affected by complex magnetic fields --- is not understood. Utilizing an improved analysis of the Fermi-LAT data that includes the first resolved imaging of the disk, we find strong evidence that this emission is produced by two separate mechanisms. Between 2010-2017 (the rise to and fall from solar maximum), the gamma-ray emission is dominated by a polar component. Between 2008-2009 (solar minimum) this component remains present, but the total emission is instead dominated by a new equatorial component with a brighter flux and harder spectrum. Most strikingly, although 6 gamma rays above 100 GeV are observed during the 1.4 years of solar minimum, none are observed during the next 7.8 years. These features, along with a 30-50 GeV spectral dip which will be discussed in a companion paper, were not anticipated by theory. To understand the underlying physics, Fermi and HAWC observations of the imminent Cycle 25 solar minimum are crucial.
• Based on broadband SEDs, we estimate jet physical parameters of 1392 $\gamma$-ray loud active galactic nuclei (AGNs), the largest sample so far. The (SED) jet power and magnetization parameter are derived for these AGNs. Out of these sources, accretion disk luminosity of 232 sources and (extended) kinetic jet powers of 159 sources are compiled from archived papers. We find, 1) FSRQs and BL Lacs are well separated by $\Gamma=-0.127\log L_{\gamma}+8.18$, in $\gamma$-ray luminosity versus photon index plane, with a succuss rate 88.6\%. 2) Most FSRQs present the (SED) jet power larger than accretion power, which suggests that the relativistic jet launching mechanism is dominated by the Blandford-Znajek process. This result confirms previous finding. 3) We find a significant anti-correlation between jet magnetization and a ratio of the (SED) jet power to the (extended) kinetic jet power, which, for the first time, provides supporting evidence for jet energy transportation theory: high magnetization jet may be more easily to transport energy to large scale than low magnetization jet.