Search

We present a complete analysis of Fermi Large Area Telescope (LAT) data of GRB 221009A, the brightest Gamma-Ray Burst (GRB) ever detected. The burst emission above 30 MeV detected by the LAT preceded by 1 s the low-energy (< 10 MeV) pulse that triggered the Fermi Gamma-Ray Burst Monitor (GBM), as has been observed in other GRBs. The prompt phase of GRB 221009A lasted a few hundred seconds. It was so bright that we identify a Bad Time Interval (BTI) of 64 seconds caused by the extremely high flux of hard X-rays and soft gamma rays, during which the event reconstruction efficiency was poor and the dead time fraction quite high. The late-time emission decayed as a power law, but the extrapolation of the late-time emission during the first 450 seconds suggests that the afterglow started during the prompt emission. We also found that high-energy events observed by the LAT are incompatible with synchrotron origin, and, during the prompt emission, are more likely related to an extra component identified as synchrotron self-Compton (SSC). A remarkable 400 GeV photon, detected by the LAT 33 ks after the GBM trigger and directionally consistent with the location of GRB 221009A, is hard to explain as a product of SSC or TeV electromagnetic cascades, and the process responsible for its origin is uncertain. Because of its proximity and energetic nature, GRB 221009A is an extremely rare event., Comment: 60 pages, 38 figures, 9 tables

While SNRs have been considered the most relevant Galactic CR accelerators for decades, CCSNe could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV gamma-rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. The multi-wavelength observation of SN 2023ixf, a Type II SN in the nearby galaxy M101, opens the possibility to constrain CR acceleration within a few days after the collapse of the RSG stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. We performed a maximum-likelihood analysis of gamma-ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. We do not find significant gamma-ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. A more sophisticated, inhomogeneous characterisation of the shock and the progenitor's environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times., Comment: Accepted in A&A. 13 pages, 12 figures, 4 tables

We analyzed the globular cluster M5 (NGC 5904) using 15 years of gamma-ray data from the Fermi Large Area Telescope (LAT). Using rotation ephemerides generated from Arecibo and FAST radio telescope observations, we searched for gamma-ray pulsations from the seven millisecond pulsars (MSPs) identified in M5. We detected no significant pulsations from any of the individual pulsars. Also, we searched for possible variations of the gamma-ray emission as a function of orbital phase for all the six MSPs in binary systems, but did not detect any significant modulations. The gamma-ray emission from the direction of M5 is well described by an exponentially cutoff power-law spectral model, although other models cannot be excluded. The phase-averaged emission is consistent with being steady on a time scale of a few months. We estimate the number of MSPs in M5 to be between 1 and 10, using the gamma-ray conversion efficiencies for well-characterized gamma-ray MSPs in the Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars, suggesting that the sample of known MSPs in M5 is (nearly) complete, even if it is not currently possible to rule out a diffuse component of the observed gamma rays from the cluster., Comment: 14 pages, 5 figures, 3 tables, published in ApJ

We recently constructed the G4Jy-3CRE, a catalog of extragalactic radio sources based on the GLEAM 4-Jy (G4Jy) sample, with the aim of increasing the number of powerful radio galaxies and quasars with similar selection criteria to those of the revised release of the Third Cambridge catalog (3CR). The G4Jy-3CRE consists of a total of 264 radio sources mainly visible from the Southern Hemisphere. Here, we present an initial X-ray analysis of 89 G4Jy-3CRE radio sources with archival X- ray observations from the Neil Gehrels Swift Observatory. We reduced a total of 615 Swift observations, for about 0.89 Msec of integrated exposure time, we found X-ray counterparts for 61 radio sources belonging to the G4Jy-3CRE, 11 of them showing extended X-ray emission. The remaining 28 sources do not show any X-ray emission associated with their radio cores. Our analysis demonstrates that X-ray snapshot observations, even if lacking uniform exposure times, as those carried out with Swift, allow us to (i) verify and/or re ne the host galaxy identi cation; (ii) discover the extended X-ray emission around radio galaxies of the intracluster medium when harbored in galaxy clusters, as the case of G4Jy 1518 and G4Jy 1664, and (iii) detect X-ray radiation arising from their radio lobes, as for G4Jy 1863., Comment: 35 pages, 17 figures, 2 tables; second paper of a series, pre-proof version

Since the early sixties, our view of radio galaxies and quasars has been drastically shaped by discoveries made thanks to observations of radio sources listed in the Third Cambridge catalog and its revised version (3CR). However, the largest fraction of data collected to date on 3CR sources was performed with relatively old instruments, rarely repeated and/or updated. Importantly, the 3CR contains only objects located in the Northern Hemisphere thus having limited access to new and innovative astronomical facilities. To mitigate these limitations we present a new catalog of powerful radio sources visible from the Southern Hemisphere, extracted from the GLEAM 4-Jy (G4Jy) catalog and based on equivalent selection criteria as the 3CR. This new catalog, named G4Jy- 3CRE, where the E stands for "equivalent", lists a total of 264 sources at declination below -5 degrees and with 9 Jy limiting sensitivity at ~178 MHz. We explored archival radio maps obtained with different surveys and compared then with optical images available in the Pan-STARRS, DES and DSS databases to search for optical counterparts of their radio cores. We compared mid-infrared counterparts, originally associated in the G4Jy, with the optical ones identified here and we present results of a vast literature search carried out to collect redshift estimates for all G4Jy-3CRE sources resulting in a total of 145 reliable z measurements., Comment: 72 pages, 35 figures, 5 Tables

We present the first X-ray observation at sub-arcsecond resolution of the high-redshift ($z=6.18$) radio-loud quasar CFHQS J142952+544717 (J1429). The ~100 net-count 0.3-7 keV spectrum obtained from $\sim 30$ ksec Chandra exposure is best fit by a single power-law model with a photon index $\Gamma=2.0\pm0.2$ and no indication of an intrinsic absorber, implying a 3.6-72 keV rest-frame luminosity $L_{\rm X}=(2.3^{+0.6}_{-0.5})\times10^{46}$ erg s$^{-1}$. We identify a second X-ray source at 30 arcsec, distance from J1429 position, with a soft ($\Gamma\simeq 2.8$) and absorbed (equivalent hydrogen column density $N_{\rm H} <13.4\times 10^{20}$ cm$^{-2}$) spectrum, which likely contaminated J1429 spectra obtained in lower angular resolution observations. Based on the analysis of the Chandra image, the bulk of the X-ray luminosity is produced within the central $\sim 3$ kpc region, either by the disk/corona system, or by a moderately aligned jet. In this context, we discuss the source properties in comparison with samples of low- and high-redshift quasars. We find indication of a possible excess of counts over the expectations for a point-like source in a 0.5 arcsec-1.5 arcsec ($\sim 3-8$ kpc) annular region. The corresponding X-ray luminosity at J1429 redshift is $4\times 10^{45}$ erg s$^{-1}$. If confirmed, this emission could be related to either a large-scale X-ray jet, or a separate X-ray source., Comment: 9 pages, 4 figures, accepted for publication in MNRAS

Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1-100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 +/- 0.4)x10^-6 photons cm^-2 s^-1. Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to kT_shock = 4 keV. The lack of a detectable 6.7 keV Fe K_alpha emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin-Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova., Comment: 20 pages, 9 figures, 3 tables. Accepted to MNRAS

We analyze radio continuum emission of early-type galaxies with dynamical measurements of central super-massive black hole (SMBH) masses, and well-characterized large-scale environments, but regardless on the exact level of the nuclear activity. The 1.4 GHz radio fluxes collected with $\sim$arcmin resolution for 62 nearby targets (distances $\lesssim$ 153 Mpc), correspond to low and very low monochromatic luminosities $L_{\rm r} \sim 10^{35} - 10^{41}$ erg s$^{-1}$. We quantify possible correlations between the radio properties with the main parameters of supermassive black holes, host galaxies, and hot gaseous halos, finding a general bimodality in the radio luminosity distribution, with the borderline between ``radio-bright'' and ``radio-dim'' populations $\log L_{\rm r} / L_{\rm Edd} \simeq -8.5$. We analyze the far-infrared data for the targets, finding that all radio-bright sources, and over a half of radio-dim ones, are over-luminous in radio with respect to the far-infrared--radio correlation. High-resolution radio maps reveal that the overwhelming majority of radio-dim sources are unresolved on arcsecond scale, while the bulk of radio-bright sources display extended jets and lobes of low- and intermediate-power radio galaxies; those jets dominate radio emission of radio-bright objects. Regarding the origin of the radio emission of radio-dim sources, we discuss the two main possibilities. One is the ADAF model, in which the radio and the nuclear X-ray radiative outputs at very low accretion rates, are both dominated by unresolved jets. The other possibility is that the radio-dim sources, unlike the radio-bright ones, are characterized by low values of SMBH spins, so that their radio emission is not related to the jets, but instead is due to a combination of starforming processes and past nuclear outbursts., Comment: revised version accepted for publication in the Astrophysical Journal

We present an analysis of multi-wavelength observations of the dark globule DC\,314.8--5.1, using data from the Gaia optical, 2MASS near-infrared, and WISE mid-infrared surveys, dedicated imaging with the Spitzer Space Telescope, and X-ray data obtained with the Swift-XRT Telescope (XRT). The main goal was to identify possible pre-main sequence stars (PMSs) and young stellar objects (YSOs) associated with the globule. For this, we studied the infrared colors of all point sources within the boundaries of the cloud. After removing sources with non-stellar spectra, we investigated the Gaia parallaxes for the YSO candidates, and found that none are physically related to DC\,314.8--5.1. In addition, we searched for X-ray emission from pre-main sequence stars with Swift-XRT, and found no 0.5--10\,keV emission down to a luminosity level $\lesssim 10^{31}$erg\,s$^{-1}$, typical of a PMS with mass\,$\ge 2 M_\odot$. Our detailed inspection therefore supports a very young, ``pre-stellar core'' evolutionary stage for the cloud. Based on archival Planck and IRAS data, we moreover identify the presence of hot dust, with temperatures $\gtrsim 100$\,K, in addition to the dominant dust component at 14\,K, originating with the associated reflection nebula., Comment: Accepted to AJ

A well-known active galaxy of the blazar type, S5 0716+714, is characterized by a particularly high variability duty cycle on short-time scales at optical frequencies. As such, the source was subjected to numerous monitoring programs, including both ground-based as well as space-borne telescopes. On closer inspection of the most recent accumulation of the data provided by the Transiting Exoplanet Survey Satellite, we have noticed several conspicuous events with `volcano-like' symmetric shape, lasting all for several hours, which closely resemble the achromatic events detected with the previous Whole Earth Blazar Telescope campaigns targeting the source. We propose that those peculiar features could be due to the gravitational micro-lensing of the innermost segments of the precessing jet in the system, by a binary lens. We study the magnification pattern of the lens with the inverse-ray shooting method, and the source trajectory parameters with the Python package MuLensModel. In this way, we were able to fit successfully all the selected events with a single lens, adjusting slightly only the source trajectory parameters for each lensing event. Based on the fit results, we postulate the presence of a massive binary lens, containing an intermediate-mass black hole, possibly even a super-massive one, and a much less massive companion (by a factor of $\lesssim 0.01$), located within the host galaxy of the blazar, most likely the central kiloparsec region. We discuss the major physical implications of the proposed scenario regarding the quest for the intermediate-mass and dual supermassive black holes in active galaxies., Comment: accepted for publication in the Astrophysical Journal

We report the Fermi-LAT gamma-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ~20-23 GeV, the highest-energy photons detected in some sub-intervals, are consistent with $\pi^{\rm 0}$-decay emission from shocks in the ejecta as proposed by Tatischeff & Hernanz (2007) for its previous 2006 outburst. The LAT light-curve displayed a fast rise to its peak >0.1 GeV flux of $\simeq$6x10^-6 ph cm^-2 s^-1 beginning on day 0.745 after its optically-constrained eruption epoch of 2021 August 8.50. The peak lasted for ~1 day, and exhibited a power-law decline up to the final LAT detection on day 45. We analyze the data on shorter timescales at early times and found evidence of an approximate doubling of emission over ~200 minutes at day 2.2, possibly indicating a localized shock-acceleration event. Comparing the data collected by the AAVSO, we measured a constant ratio of ~2.8x10^-3 between the gamma-ray and optical luminosities except for a ~5x smaller ratio within the first day of the eruption likely indicating attenuation of gamma rays by ejecta material and lower high-energy proton fluxes at the earliest stages of the shock development. The hard X-ray emission due to bremsstrahlung from shock-heated gas traced by the Swift-XRT 2-10 keV light-curve peaked at day ~6, later than at GeV and optical energies. Using X-ray derived temperatures to constrain the velocity profile, we find the hadronic model reproduces the observed >0.1 GeV light-curve., Comment: ApJ, accepted. 21 pages, 10 figures, 4 tables

Studies of the most distant AGNs allow us to test our current understanding of the physics present in radio-jetted AGNs across a range of environments. The decrease in apparent luminosity with distance is the primary difficulty to overcome in the study of these distant AGNs, which requires highly sensitive instruments. Our goal is to employ new long wavelength radio data to better parametrise the broad-band SED of GB 1508+5714, a high-redshift (z=4.30) AGN. Its high redshift, high intrinsic luminosity, and classification as a blazar allow us to test emission models that consider the efficient cooling of jet electrons via inverse Compton losses in interactions with the dense CMB photon field at high redshifts. A significant detection of this effect in GB 1508+5714 may partly explain the apparent sparsity of high-redshift radio galaxies in wide-field surveys; detections of this kind are only becoming possible with the current generation of SKA precursors. We used international LOFAR telescope to image the long wavelength radio emission around the high-redshift blazar GB 1508+5714 on arcsecond scales at frequencies between 128 MHz and 160 MHz. This allowed us to compare the spatially resolved structure with higher frequency observations, and to construct spectral index maps. The LOFAR image shows a compact unresolved core and two resolved emission regions around 2 arcsec to the east and to the west of the radio core. We find structure consistent with previous VLA observations, as well as a previously unreported emission region to the east. We interpret the arcsecond-scale radio structure of GB 1508+5714 as a FR II-like radio galaxy at a small viewing angle. Our SED modelling shows that a scenario featuring significant quenching effects caused by interaction with the CMB provides a good description of the data, and notably explains the suppressed radio emission., Comment: 11 pages, 10 figures, 2 tables

Flat spectrum radio quasars (FSRQs) are the most luminous blazars at GeV energies, but only rarely emit detectable fluxes of TeV gamma rays, typically during bright GeV flares. We explore the gamma-ray variability and spectral characteristics of three FSRQs that have been observed at GeV and TeV energies by Fermi-LAT and VERITAS, making use of almost 100 hours of VERITAS observations spread over 10 years: 3C 279, PKS 1222+216, and Ton 599. We explain the GeV flux distributions of the sources in terms of a model derived from a stochastic differential equation describing fluctuations in the magnetic field in the accretion disk, and estimate the timescales of magnetic flux accumulation and stochastic instabilities in their accretion disks. We identify distinct flares using a procedure based on Bayesian blocks and analyze their daily and sub-daily variability and gamma-ray energy spectra. Using observations from VERITAS as well as Fermi, Swift, and the Steward Observatory, we model the broadband spectral energy distributions of PKS 1222+216 and Ton 599 during VHE-detected flares in 2014 and 2017, respectively, strongly constraining the jet Doppler factors and gamma-ray emission region locations during these events. Finally, we place theoretical constraints on the potential production of PeV-scale neutrinos during these VHE flares., Comment: 34 pages, 13 figures. Accepted for publication in the Astrophysical Journal

The radio source 1146+596 is hosted by an elliptical/S0 galaxy NGC\,3894, with a low-luminosity active nucleus. The radio structure is compact, suggesting a very young age of the jets in the system. Recently, the source has been confirmed as a high-energy (HE, $>0.1$\,GeV) $\gamma$-ray emitter, in the most recent accumulation of the {\it Fermi} Large Area Telescope (LAT) data. Here we report on the analysis of the archival {\it Chandra} X-ray Observatory data for the central part of the galaxy, consisting of a single 40\,ksec-long exposure. We have found that the core spectrum is best fitted by a combination of an ionized thermal plasma with the temperature of $\simeq 0.8$\,keV, and a moderately absorbed power-law component (photon index $\Gamma = 1.4\pm 0.4$, hydrogen column density $N_{\rm H}/10^{22}$\,cm$^{-2}$\,$= 2.4\pm 0.7$). We have also detected the iron K$\alpha$ line at $6.5\pm 0.1$\,keV, with a large equivalent width of EW\,$= 1.0_{-0.5}^{+0.9}$\,keV. Based on the simulations of the {\it Chandra}'s Point Spread Function (PSF), we have concluded that, while the soft thermal component is extended on the scale of the galaxy host, the hard X-ray emission within the narrow photon energy range 6.0--7.0\,keV originates within the unresolved core (effectively the central kpc radius). The line is therefore indicative of the X-ray reflection from a cold neutral gas in the central regions of NGC\,3894. We discuss the implications of our findings in the context of the X-ray Baldwin effect. NGC\,3894 is the first young radio galaxy detected in HE $\gamma$-rays with the iron K$\alpha$ line., Comment: Revised version, accepted for publication in the Astrophysical Journal

The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase, validating the design choices and supporting the compelling case to extend the duration of the Fermi mission. The details provided here will be useful when designing the next generation of high-energy gamma-ray observatories., Comment: 60 pages, 28 figures. Published in ApJS

We present Chandra X-ray observations of 14 radio-loud quasars at redshifts $3 < z < 4$, selected from a well-defined sample. All quasars are detected in the 0.5-7.0 keV energy band, and resolved X-ray features are detected in five of the objects at distances of 1-12" from the quasar core. The X-ray features are spatially coincident with known radio features for four of the five quasars. This indicates that these systems contain X-ray jets. X-ray fluxes and luminosities are measured, and jet-to-core X-ray flux ratios are estimated. The flux ratios are consistent with those observed for nearby jet systems, suggesting that the observed X-ray emission mechanism is independent of redshift. For quasars with undetected jets, an upper limit on the average X-ray jet intensity is estimated using a stacked image analysis. Emission spectra of the quasar cores are extracted and modeled to obtain best-fit photon indices, and an Fe K emission line is detected from one quasar in our sample. We compare X-ray spectral properties with optical and radio emission in the context of both our sample and other quasar surveys., Comment: Accepted to ApJ, 15 pages, 7 figures, 7 tables

We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics., Comment: Published in ApJS

Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over different phases of the flare. Methods. The source was observed with a number of instruments in radio, near infrared, optical (including polarimetry and spectroscopy), ultra-violet, X-ray and gamma-ray bands. We use dedicated optical spectroscopy results to estimate the accretion disk and the dust torus luminosity. We perform spectral energy distribution modeling in the framework of combined Synchrotron-Self-Compton and External Compton scenario in which the electron energy distribution is partially determined from acceleration and cooling processes. Results. During the enhanced state the flux of both SED components drastically increased and the peaks were shifted to higher energies. Follow up observations with the MAGIC telescopes led to the detection of very-high-energy gamma-ray emission from this source, making it one of only a handful of FSRQs known in this energy range. Modeling allows us to constrain the evolution of the magnetic field and electron energy distribution in the emission region. The gamma-ray flare was accompanied by a rotation of the optical polarization vector during a low polarization state. Also, a new, superluminal radio knot contemporaneously appeared in the radio image of the jet. The optical spectroscopy shows a prominent FeII bump with flux evolving together with the continuum emission and a MgII line with varying equivalent width., Comment: 20 pages, 15 figures, 7 tables, accepted for publication in A&A

We report the detection of extended X-ray emission from two high-redshift radio quasars. These quasars, J1405+0415 at $z$=3.208 and J1610+1811 at $z$=3.118, were observed in a Chandra snapshot survey selected from a complete sample of the radio-brightest quasars in the overlap area of the VLA-FIRST radio survey and the Sloan Digital Sky Survey. The extended X-ray emission is located along the line connecting the core to a radio knot or hotspot, favoring the interpretation of X-ray jets. The inferred rest frame jet X-ray luminosities from 2--30 keV would be of order 10$^{45}$ erg~s$^{-1}$ if emitted isotropically and without relativistic beaming. In the scenario of inverse Compton scattering of the cosmic microwave background (CMB), X-ray jets without a coincident radio counterpart may be common, and should be readily detectable to redshifts even beyond 3.2 due to the (1+$z$)$^4$ increase of the CMB energy density compensating for the (1+$z$)$^{-4}$ cosmological diminution of surface brightness. If these can be X-ray confirmed, they would be the second and third examples of quasar X-ray jets without detection of underlying continuous radio jets., Comment: 10 pages, 4 figures. Accepted for publication in The Astrophysical Journal

We describe Glowbug, a gamma-ray telescope for bursts and other transients in the 30 keV to 2 MeV band. It was recently selected for funding by the NASA Astrophysics Research and Analysis program, with an expected launch in the early 2020s. Similar in concept to the Fermi Gamma Burst Monitor (GBM) and with similar sensitivity, Glowbug will join and enhance future networks of burst telescopes to increase sky coverage to short Gamma-Ray Bursts (SGRBs) from binary neutron star (BNS) mergers, including possible SGRBs from NS-black hole mergers. With the recent discovery of the SGRB coincident with the gravitational wave transient GW170817, we know such events occur with reasonable frequency. Expanded sky coverage in gamma rays is essential, as more detections of gravitational waves are expected with the improved sensitivity of the upgraded ground-based interferometers in the coming years., Comment: 3 pages, 2 figures, Proceeding of the Yamada Conference LXXI: Gamma-ray Bursts in the Gravitational Wave Era 2019

A significant fraction of all $\gamma$-ray sources detected by the Large Area Telescope aboard the \fer\ satellite is still lacking a low-energy counterpart. In addition, there is still a large population of $\gamma$-ray sources with associated low-energy counterparts that lack firm classifications. In the last 10 years we have undertaken an optical spectroscopic campaign to address the problem of unassociated/unidentified $\gamma$-ray sources (UGSs), mainly devoted to observing blazars and blazar candidates because they are the largest population of $\gamma$-ray sources associated to date. Here we describe the overall impact of our optical spectroscopic campaign on sources associated in \fer-LAT catalogs, coupled with objects found in the literature. In the literature search, we kept track of efforts by different teams that presented optical spectra of counterparts or potential counterparts of \fer-LAT catalog sources. Our summary includes an analysis of an additional 30 newly-collected optical spectra of counterparts or potential counterparts of \fer-LAT sources of previously unknown nature.New spectra were acquired at the Blanco 4-m and OAN-SPM 2.1-m telescopes, and those available in the Sloan Digital Sky Survey (data release 15) archive. All new sources with optical spectra analyzed here are classified as blazars. Thanks to our campaign, we altogether discovered and classified 394 targets with an additional 123 objects collected from a literature search. We began our optical spectroscopic campaign between the release of the second and third \fer-LAT source catalogs (2FGL and 3FGL, respectively), and classified about 25\% of the sources with uncertain nature and discovered a blazar-like potential counterpart for $\sim$10\% of UGSs listed therein. In the 4FGL catalog, about 350 \fer-LAT sources are classified to date thanks to our campaign. [incomplete abstract], Comment: 31 pages, 12 tables, 32 figures, 2 appendices, accepted for publication on A&A (pre-proofs version)

We present a comprehensive analysis of Chandra X-ray observations of 15 young radio quasars at redshifts $4.5 < z < 5.0$. All sources are detected in the $0.5-7.0$ keV energy band. Emission spectra are extracted, and the average photon index for the sample is measured to be $1.5\pm0.1$. Unabsorbed rest-frame $2-10$ keV luminosities are found to range between $(0.5-23.2) \times 10^{45}$ erg s$^{-1}$. The optical-X-ray power-law spectral index $\alpha_{ox}$ is calculated for each source using optical/UV data available in the literature. The $\alpha_{ox}$-UV relationship is compared with other quasar surveys, and an anticorrelation is observed that agrees with independent estimates. Rest-frame radio and X-ray luminosities are established for the sample, and a correlation between the luminosities is detected. These multiwavelength results reinforce a lack of spectral evolution for quasars over a broad redshift range. We additionally identify three quasars from our multiwavelength analysis that are statistically significant outliers, with one source being a Compton-thick candidate in the early universe, and discuss each in detail., Comment: Accepted to ApJ, 13 pages, 5 figures, 4 tables

We present an analysis of the new, deep (94\,ksec) {\it Chandra} ACIS-S observation of radio-loud active galaxy CGCG\,292$-$057, characterized by a LINER-type nucleus and a complex radio structure that indicates intermittent jet activity. On the scale of the host galaxy bulge, we detected excess X-ray emission with a spectrum best fit by a thermal plasma model with a temperature of $\sim 0.8$\,keV. We argue that this excess emission results from compression and heating of the hot diffuse fraction of the interstellar medium displaced by the expanding inner, $\sim 20$\, kpc-scale lobes observed in this restarted radio galaxy. The nuclear X-ray spectrum of the target clearly displays an ionized iron line at $\sim 6.7$\,keV, and is best fitted with a phenomenological model consisting of a power-law (photon index $\simeq 1.8$) continuum absorbed by a relatively large amount of cold matter (hydrogen column density $\simeq 0.7 \times 10^{23}$\,cm$^{-2}$), and partly scattered (fraction $\sim 3\%$) by ionized gas, giving rise to a soft excess component and K$\alpha$ line from iron ions. We demonstrate that the observed X-ray spectrum, particularly the equivalent width of Fe\,\texttt{XXV} K$\alpha$ (of order $0.3$\,keV) can in principle, be explained in a scenario involving a Compton-thin gas located at the scale of the broad-lined region in this source and photoionized by nuclear illumination. We compare the general spectral properties of the CGCG\,292$-$057 nucleus, with those of other nearby LINERs studied in X-rays., Comment: revised version, accepted for publication in the Astrophysical Journal

Here we present an analysis of the X-ray morphology and flux variability of the particularly bright and extended Western hotspot in the nearest powerful (FR II-type) radio galaxy, Pictor A, based on data obtained with the Chandra X-ray Observatory. The hotspot marks the position where the relativistic jet, that originates in the active nucleus of the system, interacts with the intergalactic medium, at hundreds-of-kiloparsec distances from the host galaxy, forming a termination shock that converts jet bulk kinetic energy to internal energy of the plasma. The hotspot is bright in X-rays due to the synchrotron emission of electrons accelerated to ultra-relativistic energies at the shock front. In our analysis, we make use of several Chandra observations targeting the hotspot over the last decades with various exposures and off-axis angles. For each pointing, we study in detail the PSF, which allows us to perform the image deconvolution, and to resolve the hotspot structure. In particular, the brightest segment of the X-ray hotspot is observed to be extended in the direction perpendicular to the jet, forming a thin, $\sim$3 kpc-long, feature that we identify with the front of the reverse shock. The position of this feature agrees well with the position of the optical intensity peak of the hotspot, but is clearly off-set from the position of the radio intensity peak, located $\sim$1 kpc further downstream. In addition, we measure the net count rate on the deconvolved images, finding a gradual flux decrease by about 30% over the 15-year timescale of the monitoring., Comment: revised version, accepted for publication in the Astrophysical Journal

In 2010/3, the Large Area Telescope on board Fermi revealed a transient gamma-ray source, positionally coincident with the optical nova in the symbiotic binary, V407Cyg. This event marked the first discovery of gamma-ray emission from a nova. We aimed to obtain resolved radio imaging of the material involved in the nova event; to determine the ejecta geometry and advance velocity directly in the image plane; to constrain the physical conditions of the system. We observed the source with the EVN and the VLBA over 16 epochs, between 20 days and 6 months after the optical discovery. The source is initially very dim but it later shows a substantial increase in brightness and a resolved shell-like structure 40 to 90 days after the optical event. The shell has a projected elliptical shape and is asymmetric in brightness and spectral index, being brighter and characterised by a rising spectrum at the S-E edge. We determine a projected velocity of ~3500 km/s in the initial phase, and ~2100 km/s between day 20 and 91. We also found an emitting feature about 350 mas (940 AU) to the N-W, advancing at a projected velocity of ~700 km/s along the polar axis of the binary. The total flux density in the VLBI images is significantly lower than that previously reported at similar epochs and over much wider angular scales with the VLA. Optical spectra demonstrated that in 2010 we were viewing V407Cyg along the equatorial plane and from behind the Mira. Our radio observations image the bipolar flow of the ejecta perpendicular to the orbital plane, where deceleration is much lower than through the equatorial plane probed by the truncated profile of optical emission lines. The separated polar knot at 350 mas and the bipolar flow strictly resemble the similar arrangement seen in Hen 2-104. The observed ~700 km/s expansion constrains the launch-date of the polar knot around 2004. [Abridged], Comment: Accepted for publication in A&A; 18 pages, 11 figures

Context. According to radiative models, radio galaxies may produce gamma-ray emission from the first stages of their evolution. However, very few such galaxies have been detected by the Fermi Large Area Telescope (LAT) so far. Aims. NGC 3894 is a nearby (z = 0.0108) object that belongs to the class of compact symmetric objects (CSOs, i.e., the most compact and youngest radio galaxies), which is associated with a gamma-ray counterpart in the Fourth Fermi-LAT source catalog. Here we present a study of the source in the gamma-ray and radio bands aimed at investigating its high-energy emission and assess its young nature. Methods. We analyzed 10.8 years of Fermi-LAT data between 100 MeV and 300 GeV and determined the spectral and variability characteristics of the source. Multi-epoch very long baseline array (VLBA) observations between 5 and 15 GHz over a period of 35 years were used to study the radio morphology of NGC 3894 and its evolution. Results. NGC 3894 is detected in gamma-rays with a significance >9 sigma over the full period, and no significant variability has been observed in the gamma-ray flux on a yearly time-scale. The spectrum is modeled with a flat power law ($\Gamma$ = 2.0$\pm$0.1) and a flux on the order of 2.2 $\times$ 10$^{-9}$ ph cm$^{-2}$ s$^{-1}$. For the first time, the VLBA data allow us to constrain with high precision the apparent velocity of the jet and counter-jet side to be $\beta_{\mathrm{app,NW}}$ = 0.132$\pm$0.004 and $\beta_{\mathrm{app,SE}}$ = 0.065$\pm$0.003, respectively. Conclusions. Fermi-LAT and VLBA results favor the youth scenario for the inner structure of this object, with an estimated dynamical age of 59$\pm$5 years. The estimated range of viewing angle (10{\deg} < $\theta$ < 21{\deg}) does not exclude a possible jet-like origin of the gamma-ray emission., Comment: Accepted by A&A on February 28, 2020. 8 pages, 7 figures

We present a systematic analysis of the mid-infrared (MIR) properties of the youngest radio galaxies, based on low-resolution data provided by the {\it WISE} and {\it IRAS} satellites. We restrict our analysis to sources with available X-ray data that constitute the earliest phase of radio galaxy evolution, i.e. those classified as Gigahertz Peaked Spectrum (GPS) and/or Compact Symmetric Objects (CSOs). In our sample of 29 objects, we find that the host galaxies are predominantly red/yellow ellipticals, with some of them displaying distorted morphology. We find a variety of MIR colors, and observe that the sources in which the MIR emission is dominated by the ISM component uniformly populate the region occupied by galaxies with a wide range of pronounced ($\geq 0.5 M_{\odot}$\,yr$^{-1}$) star formation activity. We compare the MIR color distribution in our sample to that in the general population of local AGN, in the population of evolved FR\,II radio galaxies, and also in the population of radio galaxies with recurrent jet activity. We conclude that the triggering of radio jets in AGN does not differentiate between elliptical hosts with substantially different fractions of young stars; instead there is a relationship between the jet duty cycle and the ongoing star formation. The distribution of the sub-sample of our sources with $z<0.4$ on the low-resolution MIR vs. absorption-corrected X-ray luminosity plane is consistent with the distribution of a sample of local AGN. Finally, we comment on the star formation rates of the two $\gamma$-ray detected sources in our sample, 1146+596 \& 1718--649., Comment: Revised version, accepted by the Astrophysical Journal

Here we discuss the jet production efficiency in a sample of 17 young radio galaxies with measured redshifts, kinematic ages, and nuclear X-ray fluxes, for which the observed luminosities of compact jets/lobes and accretion disks correspond to the same episode of the AGN activity. For the targets, we analyze the available optical data, estimating the bolometric luminosities of the accretion disks $L_{\rm bol}$, and the black hole masses; we also derive the minimum jet kinetic luminosities, $P_{\rm j}$. With such, we investigate the distribution of our sample in the three-dimensional space of the accretion rate $\lambda_{\rm Edd} \equiv L_{\rm bol}/L_{\rm Edd}$, the nuclear X-ray luminosity $L_{\rm X}$ considered here as a limit for the emission of the disk coronae, and $P_{\rm j}$, expressing the latter two parameters either in the Eddington units, or in the units of the disk luminosity. We find that (i) the accretion rate $\lambda_{\rm Edd}$ in our sample is distributed within a narrow range $\lambda_{\rm Edd} \sim 0.01 - 0.2$; (ii) the normalized jet power $P_{\rm j}/L_{\rm Edd}$ formally correlates with the accretion rate $\lambda_{\rm Edd}$, with some saturation at the largest values $\lambda_{\rm Edd}> 0.05$; (iii) the jet production efficiency $\eta_{\rm jet} \equiv P_{\rm j}/\dot{M}_{\rm acc} c^2$ spans a range from $\eta_{\rm jet} \lesssim 10^{-3}$ up to $\sim 0.2$ at maximum, which is below the level expected for magnetically arrested disks around maximally spinning black holes; and (iv) there is a diversification in $\eta_{\rm jet}$ on the hardness--intensity diagram $L_{\rm X}/L_{\rm bol} - \lambda_{\rm Edd}$, with the jets being produced most efficiently during the high/hard states, and suppressed during the soft states., Comment: the Astrophysical Journal, in press

Chandra HRC observations are investigated for evidence of proper motion and brightness changes in the X-ray jet of the nearby radio galaxy M87. Using images spanning 5 yr, proper motion is measured in the X-ray knot HST-1, with a superluminal apparent speed of $6.3 \pm 0.4 c$, or $24.1 \pm 1.6\rm\ mas\ yr^{-1}$, and in Knot D, with a speed of $2.4\pm 0.6c$. Upper limits are placed on the speeds of the remaining jet features. The X-ray knot speeds are in excellent agreement with existing measurements in the radio, optical, and ultraviolet. Comparing the X-ray results with images from the Hubble Space Telescope indicates that the X-ray and optical/UV emitting regions co-move. The X-ray knots also vary by up to 73% in brightness, whereas there is no evidence of brightness changes in the optical/UV. Using the synchrotron cooling models, we determine lower limits on magnetic field strengths of $\sim 420~\mu \rm G$ and $\sim 230~\mu \rm G$ for HST-1 and Knot A, respectively, consistent with estimates of the equipartition fields. Together, these results lend strong support to the synchrotron cooling model for Knot HST-1, which requires that its superluminal motion reflects the speed of the relativistic bulk flow in the jet., Comment: Accepted to ApJ, 12 pages, 4 figures, 4 tables

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources., Comment: White Paper submitted to the Astro2020 Astronomy and Astrophysics Decadal Survey organized by the Board of Physics and Astronomy and the Space Studies Board of the National Academies

We present a low-frequency imaging study of an extended sample of X-shaped radio sources using the Giant Metrewave radio telescope (GMRT) at two frequencies (610 and 240 MHz). The sources were drawn from a Very Large Array FIRST-selected sample and extends an initial GMRT study at the same frequencies, of 12 X-shaped radio galaxies predominantly from the 3CR catalog (Lal & Rao 2007). Both the intensity maps and spectral index maps of the 16 newly observed sources are presented. With the combined sample of 28 X-shaped radio sources, we found no systematic differences in the spectral properties of the higher surface brightness, active lobes versus the lower surface brightness, off-axis emission. The properties of the combined sample are discussed, including the possible role of a twin active galactic nuclei model in the formation of such objects., Comment: 11 pages, 3 Figures (plus 1 Figure in Appendix; note that Figure 2 is of lower quality), 3 Tables (plus 1 Table in Appendix)

We present a radio and X-ray study of PKS 1353-341, the brightest cluster galaxy radio source at the center of a recent Chandra-discovered X-ray cluster. Our multi-frequency VLA images reveal an edge-brightened (FR-II), double-lobed structure with total ~50 kpc extent and 1.5 GHz power of $1.2\times10^{25}$ W Hz$^{-1}$, separated from the bright, arcsecond-scale core. We reanalyzed the Chandra data and found the X-ray emitting AGN is offset by ~9 kpc from a compact X-ray cool-core with temperature, $kT=3.1\pm0.5$ keV, and a radius of ~22 kpc, surrounded by a hotter $kT=6.3\pm0.7$ keV gas out to ~50 kpc. The offset suggests sloshing inside the cool-core induced by a minor merger or a past outburst of the AGN that produced the large-scale radio lobes. The comparable spatial scales of the lobes with the interface between the different temperature X-ray plasma indicate the lobes are actively heating the outer layers of what is now a remnant compact cool-core. Our dual-frequency VLBA images reveal substructure in the central radio source, consisting of a radio core with double-sided pc-scale jets pointing toward the kpc-scale structures. The northern jet is detected only at 8.4 GHz, indicating its emission is behind an absorbing torus or disk. We also measured faster apparent motions in the southern jet up to $1.9\pm1.1c$ than in the northern jet ($0.8\pm0.5c$). While the VLBA observations indicate the southern jet is aligned slightly closer to our line of sight, the asymmetries are overall modest and imply minimal projection effects in the large-scale radio structures., Comment: 16 pages, ApJ accepted

High-redshift ($z>2$) blazars are the most powerful members of the blazar family. Yet, only a handful of them have both X-ray and $\gamma$-ray detection, thereby making it difficult to characterize the energetics of the most luminous jets. Here, we report, for the first time, the Fermi-Large Area Telescope detection of the significant $\gamma$-ray emission from the high-redshift blazar DA 193 ($z=2.363$). Its time-averaged $\gamma$-ray spectrum is soft ($\gamma$-ray photon index = $2.9\pm0.1$) and together with a relatively flat hard X-ray spectrum (14$-$195 keV photon index = $1.5\pm0.4$), DA 193 presents a case to study a typical high-redshift blazar with inverse Compton peak being located at MeV energies. An intense GeV flare was observed from this object in the first week of 2018 January, a phenomenon rarely observed from high-redshift sources. What makes this event a rare one is the observation of an extremely hard $\gamma$-ray spectrum (photon index = $1.7\pm0.2$), which is somewhat unexpected since high-redshift blazars typically exhibit a steep falling spectrum at GeV energies. The results of our multi-frequency campaign, including both space- (Fermi, NuSTAR, and Swift) and ground-based (Steward and Nordic Optical Telescope) observatories, are presented and this peculiar $\gamma$-ray flare is studied within the framework of a single-zone leptonic emission scenario., Comment: 14 pages, 9 figures, 2 tables, to appear in the Astrophysical Journal

We present the first optical spectroscopy of five confirmed (or strong candidate) redback millisecond pulsar binaries, obtaining complete radial velocity curves for each companion star. The properties of these millisecond pulsar binaries with low-mass, hydrogen-rich companions are discussed in the context of the 14 confirmed and 10 candidate field redbacks. We find that the neutron stars in redbacks have a median mass of 1.78 +/- 0.09 M_sun with a dispersion of sigma = 0.21 +/- 0.09. Neutron stars with masses in excess of 2 M_sun are consistent with, but not firmly demanded by, current observations. Redback companions have median masses of 0.36 +/- 0.04 M_sun with a scatter of sigma = 0.15 +/- 0.04, and a tail possibly extending up to 0.7-0.9 M_sun. Candidate redbacks tend to have higher companion masses than confirmed redbacks, suggesting a possible selection bias against the detection of radio pulsations in these more massive candidate systems. The distribution of companion masses between redbacks and the less massive black widows continues to be strongly bimodal, which is an important constraint on evolutionary models for these systems. Among redbacks, the median efficiency of converting the pulsar spindown energy to gamma-ray luminosity is ~10%., Comment: Apj in press. Long/large tables 6 and 7 also included as ancillary machine-readable files; numbering different from published version for ease of use

The Fermi $\gamma$-ray source 1FGL J1417.7--4407 (J1417) is a compact X-ray binary with a neutron star primary and a red giant companion in a $\sim$5.4 day orbit. This initial conclusion, based on optical and X-ray data, was confirmed when a 2.66 ms radio pulsar was found at the same location (and with the same orbital properties) as the optical/X-ray source. However, these initial studies found conflicting evidence about the accretion state and other properties of the binary. We present new optical, radio, and X-ray observations of J1417 that allow us to better understand this unusual system. We show that one of the main pieces of evidence previously put forward for an accretion disk---the complex morphology of the persistent H$\alpha$ emission line---can be better explained by the presence of a strong, magnetically driven stellar wind from the secondary and its interaction with the pulsar wind. The radio spectral index derived from VLA/ATCA observations is broadly consistent with that expected from a millisecond pulsar, further disfavoring an accretion disk scenario. X-ray observations show evidence for a double-peaked orbital light curve, similar to that observed in some redback millisecond pulsar binaries and likely due to an intrabinary shock. Refined optical light curve fitting gives a distance of 3.1$\pm$0.6 kpc, confirmed by a Gaia DR2 parallax measurement. At this distance the X-ray luminosity of J1417 is (1.0$^{+0.4}_{-0.3}$) $\times 10^{33}$ erg s$^{-1}$, which is more luminous than all known redback systems in the rotational-powered pulsar state, perhaps due to the wind from the giant companion. The unusual phenomenology of this system and its differing evolutionary path from redback millisecond pulsar binaries points to a new eclipsing pulsar "spider" subclass that is a possible progenitor of normal field millisecond pulsar binaries., Comment: 14 pages, 7 figures, 6 tables, Accepted for publication in ApJ

Centaurus A (Cen A) is the nearest radio galaxy discovered as a very-high-energy (VHE; 100 GeV-100 TeV) $\gamma$-ray source by the High Energy Stereoscopic System (H.E.S.S.). It is a faint VHE $\gamma$-ray emitter, though its VHE flux exceeds both the extrapolation from early Fermi-LAT observations as well as expectations from a (misaligned) single-zone synchrotron-self Compton (SSC) description. The latter satisfactorily reproduces the emission from Cen A at lower energies up to a few GeV. New observations with H.E.S.S., comparable in exposure time to those previously reported, were performed and eight years of Fermi-LAT data were accumulated to clarify the spectral characteristics of the $\gamma$-ray emission from the core of Cen A. The results allow us for the first time to achieve the goal of constructing a representative, contemporaneous $\gamma$-ray core spectrum of Cen A over almost five orders of magnitude in energy. Advanced analysis methods, including the template fitting method, allow detection in the VHE range of the core with a statistical significance of 12$\sigma$ on the basis of 213 hours of total exposure time. The spectrum in the energy range of 250 GeV-6 TeV is compatible with a power-law function with a photon index $\Gamma=2.52\pm0.13_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$. An updated Fermi-LAT analysis provides evidence for spectral hardening by $\Delta\Gamma\simeq0.4\pm0.1$ at $\gamma$-ray energies above $2.8^{+1.0}_{-0.6}$ GeV at a level of $4.0\sigma$. The fact that the spectrum hardens at GeV energies and extends into the VHE regime disfavour a single-zone SSC interpretation for the overall spectral energy distribution (SED) of the core and is suggestive of a new $\gamma$-ray emitting component connecting the high-energy emission above the break energy to the one observed at VHE energies., Comment: Accepted for publication in A&A, Abstract abridged for arXiv submission

Three periastron passages of the PSR B1259$-$63/LS 2883 binary system, consisting of a 48 ms rotation-powered pulsar and a $\sim30$ M$_{\odot}$ Be star, have been observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, in 2010, 2014, and 2017. During the most-recent periastron passage, sustained low-level gamma-ray emission was observed over a $\sim3$-week long interval immediately after periastron, which was followed by an interval of no emission. Sporadic flares were detected starting 40 days post-periastron and lasted approximately 50 days, during which the emission displayed significant spectral curvature, variability on timescales as short as 1.5 minutes, and peak flux levels well in excess of the pulsar spin-down power. By contrast, during the 2010 and 2014 periastron passages, significant gamma-ray emission was not observed with the LAT until 30 and 32 days post-periastron, respectively. The previous flares did not exhibit spectral curvature, showed no short term variability, and did not exceed the pulsar spin-down power. The high flux and short timescales observed in 2017 suggest significant beaming of the emission is required and constrain the size of the emission region. The flares occur long enough after periastron that the neutron star should already have passed through the extended disk-like outflow, thus constraining options for target material and seed photon sources for inverse Compton models., Comment: Accepted for publication in the Astrophysical Journal (2 July 2018), v2 matches accepted version

We present an analysis of 8 years of Fermi-LAT ( > 0.1 GeV) gamma-ray data obtained for the radio galaxy NGC 1275. The gamma-ray flux from NGC 1275 is highly variable on short (~ days to weeks) timescales, and has steadily increased over this 8-year timespan. By examining the changes in its flux and spectral shape in the LAT energy band over the entire dataset, we found that its spectral behavior changed around 2011 February (~ MJD 55600). The gamma-ray spectra at the early times evolve largely at high energies, while the photon indices were unchanged in the latter times despite rather large flux variations. To explain these observations, we suggest that the flux changes in the early times were caused by injection of high-energy electrons into the jet, while later, the gamma-ray flares were caused by a changing Doppler factor owing to variations in the jet Lorentz factor and/or changes in the angle to our line of sight. To demonstrate the viability of these scenarios, we fit the broad-band spectral energy distribution data with a one-zone synchrotron self-Compton (SSC) model for flaring and quiescent intervals before and after 2011 February. To explain the gamma-ray spectral behavior in the context of the SSC model, the maximum electron Lorentz factor would have changed in the early times, while a modest change in the Doppler factor adequately fits the quiescent and flaring state gamma-ray spectra in the later times., Comment: 13 pages, 9 figures; accepted for publication in ApJ

We describe the MERger-event Gamma-Ray (MERGR) Telescope intended for deployment by ~2021. MERGR will cover from 20 keV to 2 MeV with a wide field of view (6 sr) using nineteen gamma-ray detectors arranged on a section of a sphere. The telescope will work as a standalone system or as part of a network of sensors, to increase by ~50% the current sky coverage to detect short Gamma-Ray Burst (SGRB) counterparts to neutron-star binary mergers within the ~200 Mpc range of gravitational wave detectors in the early 2020's. Inflight software will provide realtime burst detections with mean localization uncertainties of 6 deg for a photon fluence of 5 ph cm^-2 (the mean fluence of Fermi-GBM SGRBs) and <3 deg for the brightest ~5% of SGRBs to enable rapid multi-wavelength follow-up to identify a host galaxy and its redshift. To minimize cost and time to first light, MERGR is directly derived from demonstrators designed and built at NRL for the DoD Space Test Program (STP). We argue that the deployment of a network that provides all-sky coverage for SGRB detection is of immediate urgency to the multi-messenger astrophysics community., Comment: 5 pages, 5 figures, Presented at the 7th Fermi Symposium 2017 October 15-20, Garmisch-Partenkirchen, Germany. Submitted to Proceedings of Science, Eds. J. Greiner & O. Reimer

We present optical photometric and spectroscopic observations of the likely stellar counterpart to the unassociated \emph{Fermi}-Large Area Telescope (LAT) $\gamma$-ray source 2FGL J0846.0+2820, selected for study based on positional coincidences of optical variables with unassociated LAT sources. Using optical spectroscopy from the SOAR telescope, we have identified a late-G giant in an eccentric ($e$ = 0.06) 8.133 day orbit with an invisible primary. Modeling the spectroscopy and photometry together lead us to infer a heavy neutron star primary of $\sim 2 M_{\odot}$ and a partially stripped giant secondary of $\sim 0.8 M_{\odot}$. H$\alpha$ emission is observed in some of the spectra, perhaps consistent with the presence of a faint accretion disk. We find the $\gamma$-ray flux of 2FGL J0846.0+2820 dropped substantially in mid-2009, accompanied by an increased variation in the optical brightness, and since then it has not been detected by \emph{Fermi}. The long period and giant secondary are reminiscent of the $\gamma$-ray bright binary 1FGL J1417.7--4407, which hosts a millisecond pulsar apparently in the final stages of the pulsar recycling process. The discovery of 2FGL J0846.0+2820 suggests the identification of a new subclass of millisecond pulsar binaries that are the likely progenitors of typical field millisecond pulsars., Comment: 12 pages, 8 figures, 3 tables. Accepted for publication in ApJ

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA., Comment: Published on Journal of High Energy Astrophysics (Elsevier)

Context. A number of novae have been found to emit high-energy gamma rays (>100 MeV). However, the origin of this emission is not yet understood. We report on the search for gamma-ray emission from 75 optically-detected Galactic novae in the first 7.4 years of operation of the Fermi Large Area Telescope using the Pass 8 data set. Aims. We compile an optical nova catalog including light curves from various resources and estimate the optical peak time and optical peak magnitude in order to search for gamma-ray emission to test if all novae are gamma-ray emitters. Methods. We repeat the analysis of the six novae previously identified as gamma-ray sources and develop a unified analysis strategy which we then apply to all novae in our catalog. We search for emission in a 15-day time window in two-day steps ranging from 20 days before to 20 days after the optical peak time. We perform a population study with Monte Carlo simulations to set constraints on the properties of the gamma-ray emission of novae. Results. Two new novae candidates have been found at 2sigma global significance. Although these two novae candidates were not detected at a significant level individually, taking them together with the other non-detected novae, we found a sub-threshold nova population with a cumulative 3sigma significance. We report the measured gamma-ray flux for detected sources and flux upper limits for novae without significant detection. Our results can be reproduced by several gamma-ray emissivity models (e.g. a power-law distribution with a slope of 2), while a constant emissivity model (i.e. assuming novae are standard candles) can be rejected., Comment: accepted for publication in A&A

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise., Comment: ApJ accepted

The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain., Comment: 33 pages, 22 figures & 3 tables. Published by The Astrophysical Journal. Available on the Fermi Science Support Center (FSSC) together with the 3FHL catalog

It has recently been discovered that some, if not all, classical novae emit GeV gamma rays during outburst, but the mechanisms involved in the production of the gamma rays are still not well understood. We present here a comprehensive multi-wavelength dataset---from radio to X-rays---for the most gamma-ray luminous classical nova to-date, V1324 Sco. Using this dataset, we show that V1324 Sco is a canonical dusty Fe-II type nova, with a maximum ejecta velocity of 2600 km s$^{-1}$ and an ejecta mass of few $\times 10^{-5}$ M$_{\odot}$. There is also evidence for complex shock interactions, including a double-peaked radio light curve which shows high brightness temperatures at early times. To explore why V1324~Sco was so gamma-ray luminous, we present a model of the nova ejecta featuring strong internal shocks, and find that higher gamma-ray luminosities result from higher ejecta velocities and/or mass-loss rates. Comparison of V1324~Sco with other gamma-ray detected novae does not show clear signatures of either, and we conclude that a larger sample of similarly well-observed novae is needed to understand the origin and variation of gamma rays in novae., Comment: 26 pages, 13 figures

We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6\sigma$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2., Comment: Accepted for publication by The Astrophysical Journal Supplement. Online analysis results available at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/FAVA/index.php . Corresponding authors: M. Giomi, R. Buehler, D. Kocevski, and M. Ajello

Blazars, one of the most extreme class of active galaxies, constitute so far the largest known population of $\gamma$-ray sources and their number is continuously growing in the Fermi catalogs. However in the latest release of the Fermi catalog there is still a large fraction of sources that are classified as blazar candidate of uncertain type (BCUs) for which optical spectroscopic observations are necessary to confirm their nature and their associations. In addition about 1/3 of the $\gamma$-ray point sources listed in the Third Fermi-LAT Source Catalog (3FGL) are still unassociated and lacking an assigned lower energy counterpart. Since 2012 we have been carrying out an optical spectroscopic campaign to observe blazar candidates to confirm their nature. In this paper, the sixth of the series, we present optical spectroscopic observations for 30 $\gamma$-ray blazar candidates from different observing programs we carried out with the TNG, WHT, OAN, SOAR and Magellan telescopes. We found that 21 out of 30 sources investigated are BL Lac objects while the remaining targets are classified as flat spectrum radio quasars showing the typical broad emission lines of normal quasi stellar objects. We conclude that our selection of $\gamma$-ray blazars candidates based on their multifrequency properties continues to be a successful way to discover potential low-energy counterparts of the FermiUnidentified Gamma-ray Sources (UGSs) and to confirm the nature of BCUs.

We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present in this sample. We apply a Bayesian multi-scale image reconstruction method to detect and measure the X-ray emission from the jets. We compute deviations from a baseline model that does not include the jet, and compare observed X-ray images with those computed with simulated images where no jet features exist. This allows us to compute p-value upper bounds on the significance that an X- ray jet is detected in a pre-determined region of interest. We detected 12 of the features unambiguously, and an additional 6 marginally. We also find residual emission in the cores of 3 quasars and in the background of 1 quasar that suggest the existence of unresolved X-ray jets. The dependence of the X-ray to radio luminosity ratio on redshift is a potential diagnostic of the emission mechanism, since the inverse Compton scattering of cosmic microwave background photons (IC/CMB) is thought to be redshift dependent, whereas in synchrotron models no clear redshift dependence is expected. We find that the high-redshift jets have X-ray to radio flux ratios that are marginally inconsistent with those from lower redshifts, suggesting that either the X-ray emissions is due to the IC/CMB rather than the synchrotron process, or that high redshift jets are qualitatively different., Comment: 42 pages, 14 figures, submitted to ApJ

We investigate the nature of the compact, and possibly variable nuclear radio source in the centre of WISE J0716-19, the proposed host galaxy of fast radio burst, FRB 150418. We observed WISE J0716-19 at 5.0 GHz with the European VLBI Network four times between 2016 March 16 and June 2. At three epochs, we simultaneously observed the source with e-MERLIN at the same frequency. We detected a compact source in the EVN data in each epoch with a significance up to ~8sigma. The four epochs yielded consistent results within their uncertainties, for both peak surface intensity and positions. The mean values for these quantities are I_peak = (115+-9) {\mu}Jy/beam and r.a. = 07:16:34.55496(7), dec. = -19:00:39.4754(8), respectively. The e-MERLIN data provided ~3-5sigma detections, at a position consistent with those of the EVN data. The presence of emission on angular scales intermediate between the EVN and e-MERLIN is consistent with being null. The brightness temperature of the EVN core is Tb~10^8.5K, close to the value required by Akiyama & Johnson (2016) to explain the radio properties of WISE J0716-19 in terms of interstellar induced short-term variability. Our observations provide direct, independent evidence of the existence of a nuclear compact source in WISE J0716-19, a physical scenario with no evident connection with FRB 150418. However, the EVN data do not show indication of the variability observed with the VLA., Comment: 4 pages, accepted for publication in A&A Letters to the Editor

We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) gamma-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii 2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of gamma-ray onsets beginning ~2 days after their first optical peaks. Significant gamma-ray emission was found extending to 39-55 days after their initial LAT detections, with systematically fainter and longer duration emission compared to previous gamma-ray detected classical novae. These novae were distinguished by multiple bright optical peaks that encompassed the timespans of the observed gamma rays. The gamma-ray light curves and spectra of the two novae are presented along with representative hadronic and leptonic models, and comparisons to other novae detected by the LAT are discussed., Comment: 13 pages, 6 figures, 4 tables, ApJ accepted