Showing 1–50 of 69 results for author: Goyal, A

Peas-in-a-Pod Across the Radius Valley: Rocky Systems are Less Uniform in Mass but More Uniform in Size and Spacing

Authors: Armaan V. Goyal, Songhu Wang

Abstract: The ubiquity of "peas-in-a-pod" architectural patterns and the existence of the radius valley each present a striking population-level trend for planets with $R_{p} \leq 4 R_{\oplus}$ that serves to place powerful constraints on the formation and evolution of these subgiant worlds. As it has yet to be determined whether the strength of this peas-in-a-pod uniformity differs on either side of the ra… ▽ More The ubiquity of "peas-in-a-pod" architectural patterns and the existence of the radius valley each present a striking population-level trend for planets with $R_{p} \leq 4 R_{\oplus}$ that serves to place powerful constraints on the formation and evolution of these subgiant worlds. As it has yet to be determined whether the strength of this peas-in-a-pod uniformity differs on either side of the radius valley, we separately assess the architectures of systems containing only small ($R_{p} \leq 1.6 R_{\oplus}$), rocky planets from those harboring only intermediate-size ($1.6 R_{\oplus} < R_{p} \leq 4 R_{\oplus}$), volatile-rich worlds to perform a novel statistical comparison of intra-system planetary uniformity across compositionally distinct regimes. We find that, compared to their volatile-rich counterparts, rocky systems are less uniform in mass ($2.6σ$), but more uniform in size ($4.0σ$) and spacing ($3.0σ$). We provide further statistical validation for these results, demonstrating that they are not substantially influenced by the presence of mean motion resonances, low-mass host stars, alternative bulk compositional assumptions, sample size effects, or detection biases. We also obtain tentative evidence ($>2 σ$ significance) that the enhanced size uniformity of rocky systems is dominated by the presence of super-Earths ($1 R_{\oplus} \leq R_{p} \leq 1.6 R_{\oplus}$), while their enhanced mass diversity is driven by the presence of sub-Earth ($R_{p} < 1 R_{\oplus}$) worlds. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Accepted to ApJ Letters (May 2024). 17 pages (including 3 for Appendix), 4 figures, 3 tables

Radio-only and Radio-to-far-ultraviolet Spectral Energy Distribution Modeling of 14 ULIRGs: Insights into the Global Properties of Infrared Bright Galaxies

Authors: Subhrata Dey, Arti Goyal, Katarzyna Małek, Tanio Díaz-Santos

Abstract: We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ultraviolet--infrared (FIR) and radio bands ($\sim$50 MHz to $\sim$30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet--FIR--radio SED. For the radio-only SED modeling, we use the UltraNest packag… ▽ More We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ultraviolet--infrared (FIR) and radio bands ($\sim$50 MHz to $\sim$30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet--FIR--radio SED. For the radio-only SED modeling, we use the UltraNest package, leveraging its nested sampling algorithm. Combining the results from our previous study on 11 luminous infrared galaxies (LIRGs), we discuss the global astrophysical properties of a sample of 25 starburst galaxies (z $<$ 0.5). Their radio spectra are frequently characterized by bends and turnovers, with no indication of ULIRGs exhibiting more complicated SEDs than LIRGs despite showing more signs of interactions. Including radio measurements in the CIGALE modeling constrained the dust luminosity and star formation rate (SFR) estimates by more than 1 order of magnitude better than previously reported for starburst galaxies. We show that total and nonthermal radio luminosity at 1.4 and 4.8 GHz frequencies can be good estimators of recent SFRs for all LIRGs and those ULIRGS with an insignificant influence of active galactic nuclei. A weaker but still significant correlation is observed between radio SFRs at 1.4 GHz and old (averaged over 100 Myr) SFRs based on SED modeling, indicative of multiple episodes of starburst activity during their lifetime. The thermal radio luminosity at 4.8 GHz is a better tracer of recent star formation than the thermal luminosity at 1.4 GHz. Statistically, our modeled nonthermal radio spectral indices do not significantly correlate with redshift, stellar mass, SFR, specific SFR, and dust mass. △ Less

Submitted 7 May, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: Published in the ApJ

Journal ref: ApJ 966 61 (2024)

Solar flare catalog from 3 years of Chandrayaan-2 XSM observations

Authors: Aravind Bharathi Valluvan, Ashwin Goyal, Devansh Jain, Abhinna Sundar Samantaray, Abhilash Sarwade, Kasiviswanathan Sankarasubramanian

Abstract: We present a catalog of 6266 solar flares detected by the X-Ray Solar Monitor onboard the Chandrayaan-2 lunar orbiter between 1.55 and 12.4 keV (1 and 8 Å) from 2019 September 12 to 2022 November 4, including 1469 type A flares. The catalog represents the first large sample, including both type A, hot thermal flares, and type B, impulsive flares, with a sub-A class sensitive instrument. We also de… ▽ More We present a catalog of 6266 solar flares detected by the X-Ray Solar Monitor onboard the Chandrayaan-2 lunar orbiter between 1.55 and 12.4 keV (1 and 8 Å) from 2019 September 12 to 2022 November 4, including 1469 type A flares. The catalog represents the first large sample, including both type A, hot thermal flares, and type B, impulsive flares, with a sub-A class sensitive instrument. We also detect 213 sub-A and 1330 A class flares. Individual flares are fit with an exponentially-modified Gaussian function and multi-flare groups are decomposed into individual flares. We validate our findings with flare catalogs made using visual inspection as well as automatic pipelines on Geostationary Operational Environmental Satellite and Solar Dynamics Observatory data. We find a clear bimodality in the ratio of the width to decay time between type A and B flares. We infer a power-law index of $α_F = 1.92 \pm 0.09$ for the background-subtracted peak flux distribution of XSM flares, which is consistent with the value $\sim 2$ reported in the literature. We also infer $α_F = 1.90 \pm 0.09$ for type B, and $α_F = 1.94 \pm 0.08$ for type A flares, which has previously not been reported in the literature. These comparable values hint at a similarity in their generative processes. △ Less

Submitted 8 January, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 29 pages, 15 figures, 5 tables

The LOFAR Two-Metre Sky Survey (LoTSS): VI. Optical identifications for the second data release

Authors: M. J. Hardcastle, M. A. Horton, W. L. Williams, K. J. Duncan, L. Alegre, B. Barkus, J. H. Croston, H. Dickinson, E. Osinga, H. J. A. Röttgering, J. Sabater, T. W. Shimwell, D. J. B. Smith, P. N. Best, A. Botteon, M. Brüggen, A. Drabent, F. de Gasperin, G. Gürkan, M. Hajduk, C. L. Hale, M. Hoeft, M. Jamrozy, M. Kunert-Bajraszewska, R. Kondapally , et al. (27 additional authors not shown)

Abstract: The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of $\sim 5,700$ deg$^2$. The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many ext… ▽ More The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of $\sim 5,700$ deg$^2$. The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many extended radio sources found in LOFAR images as a result of its excellent sensitivity to extended structure. In this paper we present source associations and identifications for sources in the second data release based on optical and near-infrared data, using a combination of a likelihood-ratio cross-match method developed for our first data release, our citizen science project Radio Galaxy Zoo: LOFAR, and new approaches to algorithmic optical identification, together with extensive visual inspection by astronomers. We also present spectroscopic or photometric redshifts for a large fraction of the optical identifications. In total 4,116,934 radio sources lie in the area with good optical data, of which 85% have an optical or infrared identification and 58% have a good redshift estimate. We demonstrate the quality of the dataset by comparing it with earlier optically identified radio surveys. This is by far the largest ever optically identified radio catalogue, and will permit robust statistical studies of star-forming and radio-loud active galaxies. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Comments: 29 pages. Accepted by A&A; data products available at https://lofar-surveys.org/dr2_release.html

Journal ref: A&A 678, A151 (2023)

Enhanced Size Uniformity for Near-resonant Planets

Authors: Armaan V. Goyal, Fei Dai, Songhu Wang

Abstract: Super-Earths within the same close-in, compact planetary system tend to exhibit a striking degree of uniformity in their radius, mass, and orbital spacing, and this 'peas-in-a-pod' phenomenon itself serves to provide one of the strongest constrains on planet formation at large. While it has been recently demonstrated from independent samples that such planetary uniformity occurs for both configura… ▽ More Super-Earths within the same close-in, compact planetary system tend to exhibit a striking degree of uniformity in their radius, mass, and orbital spacing, and this 'peas-in-a-pod' phenomenon itself serves to provide one of the strongest constrains on planet formation at large. While it has been recently demonstrated from independent samples that such planetary uniformity occurs for both configurations near and distant from mean motion resonance, the question thus remains if the strength of this uniformity itself differs between near-resonant and nonresonant configurations such that the two modes may be astrophysically distinct in their evolution. We thus provide in this work a novel comparative size uniformity analysis for 48 near-resonant and 251 nonresonant multi-planet systems from the California Kepler Survey catalog, evaluating uniformity both across systems and between planetary pairs within the same system. We find that while multiplanet configurations exhibit strong peas-in-a-pod size uniformity regardless of their proximity to resonance, near-resonant configurations display enhanced intra-system size uniformity as compared to their analogous nonresonant counterparts at the level of both entire systems and subsystem planetary pairs and chains. These results are broadly consistent with a variety of formation paradigms for multiplanet systems, such as convergent migration within a turbulent protoplanetary disk or planet-planet interactions incited by postnebular dynamical instabilities. Nevertheless, further investigation is necessary to ascertain whether the nonresonant and near-resonant planetary configurations respectively evolve via a singular process or mechanisms that are dynamically distinct. △ Less

Submitted 15 August, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 15 pages, 6 figures. Accepted to ApJ July 2023

Optical variability of eight FRII-type quasars with 13-yr photometric light curves

Authors: Agnieszka Kuźmicz, Arti Goyal, Stanisław Zola, Marek Jamrozy, Marek Dróżdż, Waldemar Ogłoza, Michał Siwak, Daniel E. Reichart, Vladimir V. Kouprianov, Daniel B. Caton

Abstract: We characterize the optical variability properties of eight lobe-dominated radio quasars (QSOs): B2 0709$+$37, FBQS J095206.3$+$235245, PG 1004$+$130, [HB89] 1156$+$631, [HB89] 1425$+$267, [HB89] 1503$+$691, [HB89] 1721$+$343, 4C $+$74.26, systematically monitored for a duration of 13 years since 2009. The quasars are radio-loud objects with extended radio lobes that indicate their orientation clo… ▽ More We characterize the optical variability properties of eight lobe-dominated radio quasars (QSOs): B2 0709$+$37, FBQS J095206.3$+$235245, PG 1004$+$130, [HB89] 1156$+$631, [HB89] 1425$+$267, [HB89] 1503$+$691, [HB89] 1721$+$343, 4C $+$74.26, systematically monitored for a duration of 13 years since 2009. The quasars are radio-loud objects with extended radio lobes that indicate their orientation close to the sky plane. Five of the eight QSOs are classified as giant radio quasars. All quasars showed variability during our monitoring, with magnitude variations between 0.3 and 1 mag for the least variable and the most variable QSO, respectively. We performed both structure function (SF) analysis and power spectrum density (PSD) analysis for the variability characterization and search for characteristic timescales and periodicities. As a result of our analysis, we obtained relatively steep SF slopes ($α$ ranging from 0.49 to 0.75) that are consistent with the derived PSD slopes ($\sim$2--3). All the PSDs show a good fit to single power law forms, indicating a red-noise character of variability between $\sim$13 years and weeks timescales. We did not measure reliable characteristic timescales of variability from the SF analysis which indicates that the duration of the gathered data is too short to reveal them. The absence of bends in the PSDs (change of slope from $\geq$1 to $\sim$0) on longer timescales indicates that optical variations are most likely caused by thermal instabilities in the accretion disk. △ Less

Submitted 8 October, 2022; originally announced October 2022.

Comments: Accepted for publication in ApJS; 17 pages, 5 figures, 5 tables

Low frequency radio continuum imaging and SED modeling of 11 LIRGs: radio-only and FUV to radio bands

Authors: Subhrata Dey, Arti Goyal, Katarzyna Małek, Timothy J. Galvin, Nicholas Seymour, Tanio Díaz Santos, Julia Piotrowska, Vassilis Charmandaris

Abstract: We present the detailed analysis of 11 local luminous infrared galaxies (LIRGs) from ultraviolet through far-infrared to radio ($\sim$70 MHz to $\sim$15 GHz) bands. We derive the astrophysical properties through spectral energy distribution (SED) modeling using the Code Investigating GALaxy Emission (CIGALE) and UltraNest codes. The radio SEDs include our new observations at 325 and 610 MHz from t… ▽ More We present the detailed analysis of 11 local luminous infrared galaxies (LIRGs) from ultraviolet through far-infrared to radio ($\sim$70 MHz to $\sim$15 GHz) bands. We derive the astrophysical properties through spectral energy distribution (SED) modeling using the Code Investigating GALaxy Emission (CIGALE) and UltraNest codes. The radio SEDs include our new observations at 325 and 610 MHz from the GMRT and the measurements from public archives. Our main results are (1) radio SEDs show turnovers and bends, (2) the synchrotron spectral index of the fitted radio spectra ranges between $-$0.5 and $-$1.7, and (3) the infrared luminosity, dust mass, dust temperature, stellar mass, star-formation rates (SFRs) and AGN fraction obtained from CIGALE falls in the range exhibited by galaxies of the same class. The ratio of 60$μ$m infrared and 1.4 GHz radio luminosity, the 1.4 GHz thermal fraction, and emission measure range between 2.1 and 2.9, 0.1% and 10%, 0.02 and 269.5$\times$10$^{6}$ cm$^{-6}$ pc, respectively. We conclude that the turnovers seen in the radio SEDs are due to free-free absorption; this is supported by the low AGN fraction derived from the CIGALE analysis. The decomposed 1.4 GHz thermal and nonthermal radio luminosities allowed us to compute the star formation rate (SFR) using scaling relations. A positive correlation is observed between the SFR$_{IR}$ obtained 10 Myr ago (compared to 100 Myr ago) and 1.4 GHz radio (total and nonthermal) because similar synchrotron lifetimes are expected for typical magnetic field strengths observed in these galaxies ($\approx$50$μ$G). △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: ApJ accepted. Comments are welcome

Generalized Peas-in-a-Pod: Extending Intra-System Mass Uniformity to Non-TTV Systems via the Gini Index

Authors: Armaan Goyal, Songhu Wang

Abstract: It has been demonstrated that planets belonging to the same close-in, compact multiple-planet system tend to exhibit a striking degree of uniformity in their sizes. A similar trend has also been found to hold for the masses of such planets, but considerations of such intra-system mass uniformity have generally been limited to statistical samples wherein a majority of systems have constituent plane… ▽ More It has been demonstrated that planets belonging to the same close-in, compact multiple-planet system tend to exhibit a striking degree of uniformity in their sizes. A similar trend has also been found to hold for the masses of such planets, but considerations of such intra-system mass uniformity have generally been limited to statistical samples wherein a majority of systems have constituent planetary mass measurements obtained via analysis of transit timing variations (TTVs). Since systems with strong TTV signals typically lie in or near mean motion resonance, it remains to be seen whether intra-system mass uniformity is still readily emergent for non-resonant systems with non-TTV mass provenance. We thus present in this work a mass uniformity analysis of 17 non-TTV systems with masses measured via radial velocity (RV) measurements. Using the Gini index, a common statistic for economic inequality, as our primary metric for uniformity, we find that our sample of 17 non-TTV systems displays intra-system mass uniformity at a level of $\sim 2.5 σ$ confidence. We provide additional discussion of possible statistical and astrophysical underpinnings for this result. We also demonstrate the existence of a correlation ($r = 0.25$) between characteristic solid surface density ($Σ_0$) of the minimum mass extrasolar nebula (MMEN) and system mass Gini index, suggesting that more massive disks may generally form systems with more unequal planetary masses. △ Less

Submitted 7 June, 2022; v1 submitted 31 May, 2022; originally announced June 2022.

Comments: 9 pages, 5 figures. Accepted to ApJ

A Radio-selected Population of Dark, Long Gamma-ray Bursts: Comparison to the Long Gamma-ray Burst Population and Implications for Host Dust Distributions

Authors: Genevieve Schroeder, Tanmoy Laskar, Wen-fai Fong, Anya E. Nugent, Edo Berger, Ryan Chornock, Kate D. Alexander, Jennifer Andrews, R. Shane Bussmann, Alberto J. Castro-Tirado, Armaan V. Goyal, Charles D. Kilpatrick, Maura Lally, Adam Miller, Peter Milne, Kerry Paterson, Alicia Rouco Escorial, Michael C. Stroh, Giacomo Terreran, Bevin Ashley Zauderer

Abstract: We present cm-band and mm-band afterglow observations of five long-duration $γ$-ray bursts (GRBs; GRB 130131A, 130420B, 130609A, 131229A, 140713A) with dust-obscured optical afterglow emission, known as "dark" GRBs. We detect the radio afterglow of two of the dark GRBs (GRB 130131A and 140713A), along with a tentative detection of a third (GRB 131229A) with the Karl G. Jansky Very Large Array (VLA… ▽ More We present cm-band and mm-band afterglow observations of five long-duration $γ$-ray bursts (GRBs; GRB 130131A, 130420B, 130609A, 131229A, 140713A) with dust-obscured optical afterglow emission, known as "dark" GRBs. We detect the radio afterglow of two of the dark GRBs (GRB 130131A and 140713A), along with a tentative detection of a third (GRB 131229A) with the Karl G. Jansky Very Large Array (VLA). Supplemented by three additional VLA-detected dark GRBs from the literature, we present uniform modeling of their broadband afterglows. We derive high line-of-sight dust extinctions of $A_{V, \rm GRB} \gtrsim 2.2 - 10.6~{\rm mag}$. Additionally, we model the host galaxies of the six bursts in our sample, and derive host galaxy dust extinctions of $A_{V, \rm Host} \approx 0.3-4.7~{\rm mag}$. Across all tested $γ$-ray (fluence and duration) and afterglow properties (energy scales, geometries and circumburst densities), we find dark GRBs to be representative of more typical unobscured long GRBs, except in fluence, for which observational biases and inconsistent classification may influence the dark GRB distribution. Additionally, we find that $A_{V, \rm GRB}$ is not related to a uniform distribution of dust throughout the host, nor to the extremely local environment of the burst, indicating that a larger scale patchy dust distribution is the cause of the high line-of-sight extinction. Since radio observations are invaluable to revealing heavily dust-obscured GRBs, we make predictions for the detection of radio emission from host star formation with the next generation VLA. △ Less

Submitted 6 May, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Comments: 40 pages, 15 figures, submitted to ApJ

Multiwavelength variability power spectrum analysis of the blazars 3C 279 and PKS 1510-089 on multiple timescales

Authors: Arti Goyal, Marian Soida, Lukasz Stawarz, Paul J. Wiita, Kari Nilsson, Svetlana Jorstad, Alan P. Marscher, Margo F. Aller, Hugh D. Aller, Anne Lahteenmaki, Talvikki Hovatta, Staszek Zola, Krzysztof Nalewajko, Merja Tornikoski, Joni Tammi, Mark Hodges, Sebastian Kiehlmann, Anthony C. S. Readhead, Walter Max-Moerbeck, Elina Lindfors, Vandad Fallah Ramazani, D. E. Reichart, D. B. Caton, Janeth Valverde, Deirdre Horan , et al. (2 additional authors not shown)

Abstract: We present the results of variability power spectral density (PSD) analysis using multiwavelength radio to GeV\,$γ$-ray light curves covering decades/years to days/minutes timescales for the blazars 3C 279 and PKS 1510-089. The PSDs are modeled as single power-laws, and the best-fit spectral shape is derived using the `power spectral response' method. With more than ten years of data obtained with… ▽ More We present the results of variability power spectral density (PSD) analysis using multiwavelength radio to GeV\,$γ$-ray light curves covering decades/years to days/minutes timescales for the blazars 3C 279 and PKS 1510-089. The PSDs are modeled as single power-laws, and the best-fit spectral shape is derived using the `power spectral response' method. With more than ten years of data obtained with weekly/daily sampling intervals, most of the PSDs cover ~2-4 decades in temporal frequency; moreover, in the optical band, the PSDs cover ~6 decades for 3C 279 due to the availability of intranight light curves. Our main results are the following: (1) on timescales ranging from decades to days, the synchrotron and the inverse Compton spectral components, in general, exhibit red-noise (slope ~2) and flicker-noise (slope ~1) type variability, respectively; (2) the slopes of $γ$-ray variability PSDs obtained using a 3-hr integration bin and a 3-weeks total duration exhibit a range between ~1.4 and ~2.0 (mean slope = 1.60$\pm$0.70), consistent within errors with the slope on longer timescales; (3) comparisons of fractional variability indicate more power on timescales $\leq$100\, days at $γ$-ray frequencies as compared to longer wavelengths, in general (except between $γ$-ray and optical frequencies for PKS 1510$-$089); (4) the normalization of intranight optical PSDs for 3C\,279 appears to be a simple extrapolation from longer timescales, indicating a continuous (single) process driving the variability at optical wavelengths; (5) the emission at optical/infrared wavelengths may involve a combination of disk and jet processes for PKS\,1510$-$089. △ Less

Submitted 20 January, 2022; v1 submitted 8 December, 2021; originally announced December 2021.

Comments: ApJ accepted. Comments are welcome

Optical variability power spectrum analysis of blazar sources on intranight timescales

Authors: Arti Goyal

Abstract: We report the first results of a systematic investigation to characterize blazar variability power spectral densities (PSDs) at optical frequencies using densely sampled (5--15 minutes integration time), high photometric accuracy ($\lesssim$0.2--0.5\%) R-band intranight light curves, covering timescales ranging from several hours to $\sim$15\,minutes. Our sample consists of 14 optically bright bla… ▽ More We report the first results of a systematic investigation to characterize blazar variability power spectral densities (PSDs) at optical frequencies using densely sampled (5--15 minutes integration time), high photometric accuracy ($\lesssim$0.2--0.5\%) R-band intranight light curves, covering timescales ranging from several hours to $\sim$15\,minutes. Our sample consists of 14 optically bright blazars, including nine BL Lacertae objects (BL Lacs) and five flat-spectrum radio quasars (FSRQs) which have shown statistically significant variability during 29 monitoring sessions. We model the intranight PSDs as simple power--laws and derive the best-fit slope along with uncertainty using the `power spectral response' method. Our main results are the following: (1) on 19 out of 29 monitoring sessions, the intranight PSDs show an acceptable fit to simple power-laws at the rejection confidence $\leq$ 90\%; (2) for these 19 instances, the PSD slopes show a large range from 1.4 to 4.0, consistent with statistical characters of red (slope$\sim$2) and black (slope$\geq$3) noise stochastic processes; (3) the average PSD slopes for the BL Lacs and FSRQs are indistinguishable from one another; (4) the normalization of intranight PSDs for individual blazar sources which were monitored on more than one occasion turns out to be consistent with one another with a few exceptions. The average PSD slope, 2.9$\pm$0.3 (1$σ$ uncertainty) is steeper than the red-noise type character of variability found on longer timescales (many decades to days), indicative of a cutoff in the variability spectrum on timescales around a few days at the synchrotron frequencies of the emission spectrum. △ Less

Submitted 2 January, 2021; originally announced January 2021.

Comments: ApJ accepted

Identifying radio active galactic nuclei among radio-emitting galaxies

Authors: Dorota Kozieł-Wierzbowska, Natalia Vale Asari, Grażyna Stasińska, Fabio R. Herpich, Marek Sikora, Natalia Żywucka, Arti Goyal

Abstract: Basing our analysis on ROGUE I, a catalog of over 32,000 radio sources associated with optical galaxies, we provide two diagnostics to select the galaxies where the radio emission is due to an active galactic nucleus (AGN). Each of these diagnostics can be applied independently. The first one, dubbed MIRAD, compares the flux $F_{W3}$ in the $W3$ mid-infrared band of the WISE telescope, with the ra… ▽ More Basing our analysis on ROGUE I, a catalog of over 32,000 radio sources associated with optical galaxies, we provide two diagnostics to select the galaxies where the radio emission is due to an active galactic nucleus (AGN). Each of these diagnostics can be applied independently. The first one, dubbed MIRAD, compares the flux $F_{W3}$ in the $W3$ mid-infrared band of the WISE telescope, with the radio flux at 1.4 GHz, $\Frad$. MIRAD requires no optical spectra. The second diagnostic, dubbed DLM, relates the 4000 Å break strength, $D_{\rm n}(4000)$, with the radio luminosity per unit stellar mass. The DLM diagram has already been used in the past, but not as standalone. For these two diagrams, we propose simple, empirical dividing lines that result in the same classification for the objects in common. These lines correctly classify as radio-AGN 99.5 percent of the extended radio sources in the ROGUE~I catalog, and as star-forming (SF) galaxies 98--99 percent of the galaxies identified as such by their emission line ratios. Both diagrams clearly show that radio AGNs are preferentially found among elliptical galaxies and among galaxies hosting the most massive black holes. Most of the radio sources classified as radio-AGNs in the MIRAD or DLM diagrams are either optically weak AGNs or retired galaxies. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 14 pages, 14 figures, submitted to ApJ

The LOFAR Two Metre Sky Survey: Deep Fields Data Release 1 -- III. Host-galaxy identifications and value added catalogues

Authors: R. Kondapally, P. N. Best, M. J. Hardcastle, D. Nisbet, M. Bonato, J. Sabater, K. J. Duncan, I. McCheyne, R. K. Cochrane, R. A. A. Bowler, W. L. Williams, T. W. Shimwell, C. Tasse, J. H. Croston, A. Goyal, M. Jamrozy, M. J. Jarvis, V. H. Mahatma, H. J. A. R\öttgering, D. J. B. Smith, A. Wo\łowska, M. Bondi, M. Brienza, M. J. I. Brown, M. Br\üggen , et al. (18 additional authors not shown)

Abstract: We present the source associations, cross-identifications, and multi-wavelength properties of the faint radio source population detected in the deep tier of the LOFAR Two Metre Sky Survey (LoTSS): the LoTSS Deep Fields. The first LoTSS Deep Fields data release consists of deep radio imaging at 150~MHz of the ELAIS-N1, Lockman Hole, and Boötes fields, down to RMS sensitives of around 20, 22, and 32… ▽ More We present the source associations, cross-identifications, and multi-wavelength properties of the faint radio source population detected in the deep tier of the LOFAR Two Metre Sky Survey (LoTSS): the LoTSS Deep Fields. The first LoTSS Deep Fields data release consists of deep radio imaging at 150~MHz of the ELAIS-N1, Lockman Hole, and Boötes fields, down to RMS sensitives of around 20, 22, and 32$~μ$Jy\,beam$^{-1}$, respectively. These fields are some of the best studied extra-galactic fields in the northern sky, with existing deep, wide-area panchromatic photometry from X-ray to infrared wavelengths, covering a total of $\approx$~26~\mbox{deg$^{2}$}. We first generated improved multi-wavelength catalogues in ELAIS-N1 and Lockman Hole; combined with the existing catalogue for Boötes, we present forced, matched aperture photometry for over 7.2 million sources across the three fields. We identified multi-wavelength counterparts to the radio detected sources, using a combination of the Likelihood Ratio method and visual classification, which greatly enhances the scientific potential of radio surveys and allows for the characterisation of the photometric redshifts and the physical properties of the host galaxies. The final radio-optical cross-matched catalogue consists of 81\,951 radio-detected sources, with counterparts identified and multi-wavelength properties presented for 79\,820 ($>$97\%) sources. We also examine the properties of the host galaxies, and through stacking analysis find that the radio population with no identified counterpart is likely dominated by AGN at $z\sim3-4$. This dataset contains one of the largest samples of radio-selected star-forming galaxies and active galactic nuclei (AGN) at these depths, making it ideal for studying the history of star-formation, and the evolution of galaxies and AGN across cosmic time. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 34 pages, 14 figures and 14 tables. Accepted for publication in A&A special issue. This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication

Journal ref: A&A 648, A3 (2021)

The Broad-band Counterpart of the Short GRB 200522A at $z=0.5536$: A Luminous Kilonova or a Collimated Outflow with a Reverse Shock?

Authors: W. Fong, T. Laskar, J. Rastinejad, A. Rouco Escorial, G. Schroeder, J. Barnes, C. D. Kilpatrick, K. Paterson, E. Berger, B. D. Metzger, Y. Dong, A. E. Nugent, R. Strausbaugh, P. K. Blanchard, A. Goyal, A. Cucchiara, G. Terreran, K. D. Alexander, T. Eftekhari, C. Fryer, B. Margalit, R. Margutti, M. Nicholl

Abstract: We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short GRB 200522A, located at a small projected offset of $\approx 1$ kpc from the center of a young, star-forming host galaxy at $z=0.5536$. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our HST observat… ▽ More We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short GRB 200522A, located at a small projected offset of $\approx 1$ kpc from the center of a young, star-forming host galaxy at $z=0.5536$. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our HST observations at a rest-frame time of $\approx2.3$ days, has a luminosity of $\approx (1.3-1.7) \times 10^{42}$ erg s$^{-1}$. This is substantially lower than on-axis short GRB afterglow detections, but is a factor of $\approx 8$-$17$ more luminous than the kilonova of GW170817, and significantly more luminous than any kilonova candidate for which comparable observations exist. The combination of the counterpart's color ($i-y = -0.08\pm 0.21$; rest-frame) and luminosity cannot be explained by standard radioactive heating alone. We present two scenarios to interpret the broad-band behavior of GRB 200522A: a synchrotron forward shock with a luminous kilonova (potentially boosted by magnetar energy deposition), or forward and reverse shocks from a $\approx14^{\circ}$, relativistic ($Γ_0 \gtrsim 80$) jet. Models which include a combination of enhanced radioactive heating rates, low-lanthanide mass fractions, or additional sources of heating from late-time central engine activity may provide viable alternate explanations. If a stable magnetar was indeed produced in GRB 200522A, we predict that late-time radio emission will be detectable starting $\approx 0.3$-$6$ years after the burst for a deposited energy of $\approx 10^{53}$ erg. Counterparts of similar luminosity to GRB 200522A associated with gravitational wave events will be detectable with current optical searches to $\approx\!250$ Mpc. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 33 pages, 13 figures, 5 tables. Submitted to AAS Journals

A Late-time Radio Survey of Short GRBs at $z<0.5$: New Constraints on the Remnants of Neutron Star Mergers

Authors: Genevieve Schroeder, Ben Margalit, Wen-fai Fong, Brian D. Metzger, Peter K. G. Williams, Kerry Paterson, Kate D. Alexander, Tanmoy Laskar, Armaan V. Goyal, Edo Berger

Abstract: Massive, rapidly-spinning magnetar remnants produced as a result of binary neutron star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. We present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (SGRBs; $z<0.5$) on r… ▽ More Massive, rapidly-spinning magnetar remnants produced as a result of binary neutron star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. We present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (SGRBs; $z<0.5$) on rest-frame timescales of $\approx2.4-13.9$ yr following the bursts. We place $3σ$ limits on radio continuum emission of $F_ν\lesssim6-20\,μ$Jy at the burst positions, or $L_ν\lesssim(0.6-8.3)\times10^{28}$erg s$^{-1}$Hz$^{-1}$. Comparing these limits with new light curve modeling which properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of $E_{ej}\lesssim(0.6-6.7)\times 10^{52}$erg ($E_{ej}\lesssim(1.8-17.6)\times10^{52}$erg) for an ejecta mass of $0.03\,M_{\odot}$ ($0.1\,M_{\odot}$). We present a uniform re-analysis of 27 SGRBs with $5.5-6.0$ GHz observations, and find that $\gtrsim50\%$ of SGRBs did not form stable magnetar remnants in their mergers. Assuming SGRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction $f_{GW190425}$ of cases, we place constraints on the maximum mass of a non-rotating neutron star (NS) ($M_{TOV}$), finding $M_{TOV}\lesssim2.23\,M_{\odot}$ for $f_{GW190425}=0.4$; this limit increases for larger values of $f_{GW190425}$. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If $\gtrsim30-300$ radio remnants are discovered in VLASS, this suggests that SGRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce. △ Less

Submitted 1 September, 2020; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: 21 pages, 7 figures, submitted to ApJ

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155-304

Authors: Arti Goyal

Abstract: We present the results of the power spectral density (PSD) analysis for the blazars Mrk\,421 and PKS\,2155$-$304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability timescales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from OVRO, optical and… ▽ More We present the results of the power spectral density (PSD) analysis for the blazars Mrk\,421 and PKS\,2155$-$304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability timescales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from OVRO, optical and infrared (B, V, R, I, J, H, and K-bands), X-rays from the {\it Swift} and the {\it Rossi} X-ray Timing Explorer, high and very high energy $γ-$rays from the {\it Fermi} and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infra-red and optical frequencies have slopes $\sim$1.8, indicative of random-walk type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to very high energy \,$γ$-rays, however, have slopes $\sim$1.2, suggesting a flicker noise type process; (3) there is significantly more variability power at X-rays, high and very high energy $γ$-rays on timescales $\lesssim$ 100 days, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the timescales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: MNRAS accepted

Catalogue with visual morphological classification of 32,616 radio galaxies with optical hosts

Authors: Natalia Żywucka, Dorota Kozieł-Wierzbowska, Arti Goyal

Abstract: We present the catalogue of Radio sources associated with Optical Galaxies and having Unresolved or Extended morphologies I (ROGUE I). It was generated by cross-matching galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR 7) as well as radio sources from the First Images of Radio Sky at Twenty Centimetre (FIRST) and the National Radio Astronomical Observatory VLA Sky Survey (NVSS) c… ▽ More We present the catalogue of Radio sources associated with Optical Galaxies and having Unresolved or Extended morphologies I (ROGUE I). It was generated by cross-matching galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR 7) as well as radio sources from the First Images of Radio Sky at Twenty Centimetre (FIRST) and the National Radio Astronomical Observatory VLA Sky Survey (NVSS) catalogues. We created the largest handmade catalogue of visually classified radio objects and associated with them optical host galaxies, containing 32,616 galaxies with a FIRST core within 3 arcsec of the optical position. All listed objects possess the good quality SDSS DR 7 spectra with the signal-to-noise ratio $>$10 and spectroscopic redshifts up to $z=0.6$. The radio morphology classification was performed by a visual examination of the FIRST and the NVSS contour maps overlaid on a DSS image, while an optical morphology classification was based on the 120 arcsec snapshot images from SDSS DR 7. The majority of radio galaxies in ROGUE I, i.e. $\sim$ 93%, are unresolved (compact or elongated), while the rest of them exhibit extended morphologies, such as Fanaroff-Riley (FR) type I, II, and hybrid, wide-angle tail, narrow-angle tail, head-tail sources, and sources with intermittent or reoriented jet activity, i.e. double-double, X-shaped, and Z-shaped. Most of FR IIs have low radio luminosities, comparable to the luminosities of FR Is. Moreover, due to visual check of all radio maps and optical images, we were able to discover or reclassify a number of radio objects as giant, double-double, X-shaped, and Z-shaped radio galaxies. The presented sample can serve as a database for training automatic methods of identification and classification of optical and radio galaxies. △ Less

Submitted 23 March, 2020; originally announced March 2020.

Comments: Proceedings paper under the IAU symposium "Nuclear Activity in Galaxies Across Cosmic Time", eds. M. Povic, P. Marziani, J. Masegosa, H. Netzer, S. H. Negu, and S. B. Tessema

Radio sources associated with Optical Galaxies and having Unresolved or Extended morphologies (ROGUE). I. A catalog of SDSS galaxies with FIRST core identifications

Authors: Dorota Kozieł-Wierzbowska, Arti Goyal, Natalia Żywucka

Abstract: We present the catalog of Radio sources associated with Optical Galaxies and having Unresolved or Extended morphologies I (ROGUE~I), consisting of 32,616 spectroscopically selected galaxies. It is the largest handmade catalog of this kind, obtained by cross-matching galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 7 and radio sources from both the First Images of Radio Sky at Twenty… ▽ More We present the catalog of Radio sources associated with Optical Galaxies and having Unresolved or Extended morphologies I (ROGUE~I), consisting of 32,616 spectroscopically selected galaxies. It is the largest handmade catalog of this kind, obtained by cross-matching galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 7 and radio sources from both the First Images of Radio Sky at Twenty Centimetre (FIRST) survey and the NRAO VLA Sky Survey \textit{without imposing a limit to the radio flux densities}. The catalog provides a \textit{visual} classification of radio and optical morphologies of galaxies presenting a FIRST core within 3\arcsec\ of the optical position. The radio morphological classification is performed by examining the radio-optical overlays of linear sizes equal to 1 Mpc at the source distance, while the 120\arcsec\ image snapshots from the SDSS database are used for optical classification. The results of our search are: (i) single-component unresolved and elongated, radio sources constitute the major group in the ROGUE I catalog ($\sim$90%), and $\sim$8% exhibiting {\it extended} morphologies, (ii) samples of 269, 730, and 115 Fanaroff-Riley (FR) type I, II, and hybrid galaxies, respectively, are presented (iii) we report 55 newly discovered giant/possible giant, 16 double-double, 9 X-shaped, and 25 Z-shaped radio sources, (iv) on the optical front, most galaxies have elliptical morphologies ($\sim$62%) while spirals form the second major category ($\sim$17%) followed by distorted ($\sim$ 12%) and lenticular ($\sim$7%) morphologies, (v) division between the FR I and the FR~II sources in the radio-optical luminosity plane is blurred, in tune with recent studies. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Comments: 43 pages, 14 tables, 11 figures, accepted for publication in ApJS

Multiband optical flux density and polarization microvariability study of optically bright blazars

Authors: Magdalena Pasierb, Arti Goyal, Michał Ostrowski, Łukasz Stawarz, Paul J. Wiita, Gopal-Krishna, Valeri M. Larionov, Daria A. Morozova, Ryosuke Itoh, Fahri Alicavus, Ahmet Erdem, Santosh Joshi, Staszek Zola, Georgy A. Borman, Tatiana S. Grishina, Evgenia N. Kopatskaya, Elena G. Larionova, Sergey S. Savchenko, Anna A. Nikiforova, Yulia V. Troitskaya, Ivan S. Troitsky, Hiroshi Akitaya, Miho Kawabata, Tatsuya Nakaoka

Abstract: We present the results of flux density, spectral index, and polarization intra-night monitoring studies of a sample of eight optically bright blazars, carried out by employing several small to moderate aperture (0.4\,m to 1.5\,m diameter) telescopes fitted with CCDs and polarimeters located in Europe, India, and Japan. The duty cycle of flux variability for the targets is found to be $\sim 45$ per… ▽ More We present the results of flux density, spectral index, and polarization intra-night monitoring studies of a sample of eight optically bright blazars, carried out by employing several small to moderate aperture (0.4\,m to 1.5\,m diameter) telescopes fitted with CCDs and polarimeters located in Europe, India, and Japan. The duty cycle of flux variability for the targets is found to be $\sim 45$ percent, similar to that reported in earlier studies. The computed two-point spectral indices are found to be between 0.65 to 1.87 for our sample, comprised of low- and intermediate frequency peaked blazars, with one exception; they are also found to be statistically variable for about half the instances where `confirmed' variability is detected in flux density. In the analysis of the spectral evolution of the targets on hourly timescale, a counter-clockwise loop (soft-lagging) is noted in the flux-spectral index plane on two occasions, and in one case a clear spectral flattening with the decreasing flux is observed. In our data set, we also observe a variety of flux-polarization degree variability patterns, including instances with a relatively straightforward anti-correlation, correlation, or counter-clockwise looping. These changes are typically reflected in the flux-polarization angle plane: the anti-correlation between the flux and polarization degree is accompanied by an anti-correlation between the polarization angle and flux, while the counter-clockwise flux-PD looping behaviour is accompanied by a clockwise looping in the flux-polarization angle representation. We discuss our findings in the framework of the internal shock scenario for blazar sources. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: MNRAS accepted

The LOFAR Two-metre Sky Survey (LoTSS) III. First Data Release: optical/IR identifications and value-added catalogue

Authors: W. L. Williams, M. J. Hardcastle, P. N. Best, J. Sabater, J. H. Croston, K. J. Duncan, T. W. Shimwell, H. J. A. Röttgering, D. Nisbet, G. Gürkan, L. Alegre, R. K. Cochrane, A. Goyal, C. L. Hale, N. Jackson, M. Jamrozy, R. Kondapally, M. Kunert-Bajraszewska, V. H. Mahatma, B. Mingo, L. K. Morabito, I. Prandoni, C. Roskowinski, A. Shulevski, D. J. B. Smith , et al. (16 additional authors not shown)

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of th… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 8 figures, 5 tables and 23 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available at https://lofar-surveys.org/ upon publication in the journal

Journal ref: A&A 622, A2 (2019)

The LOFAR Two-metre Sky Survey - II. First data release

Authors: T. W. Shimwell, C. Tasse, M. J. Hardcastle, A. P. Mechev, W. L. Williams, P. N. Best, H. J. A. Röttgering, J. R. Callingham, T. J. Dijkema, F. de Gasperin, D. N. Hoang, B. Hugo, M. Mirmont, J. B. R. Oonk, I. Prandoni, D. Rafferty, J. Sabater, O. Smirnov, R. J. van Weeren, G. J. White, M. Atemkeng, L. Bester, E. Bonnassieux, M. Brüggen, G. Brunetti , et al. (82 additional authors not shown)

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,μ$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 16 figures, 1 table and 22 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Identification of blazar candidates behind Small and Large Magellanic Clouds

Authors: Natalia Żywucka, Arti Goyal, Marek Jamrozy, Łukasz Stawarz, Michał Ostrowski, Szymon Kozłowski, Andrzej Udalski

Abstract: We report the identification of blazar candidates behind the Magellanic Clouds. The objects were selected from the Magellanic Quasars Survey (MQS), which targeted the entire Large Magellanic Cloud (LMC) and 70\% of the Small Magellanic Cloud (SMC). Among the 758 MQS quasars and 898 of unidentified (featureless spectra) objects, we identified a sample of 44 blazar candidates, including 27 flat spec… ▽ More We report the identification of blazar candidates behind the Magellanic Clouds. The objects were selected from the Magellanic Quasars Survey (MQS), which targeted the entire Large Magellanic Cloud (LMC) and 70\% of the Small Magellanic Cloud (SMC). Among the 758 MQS quasars and 898 of unidentified (featureless spectra) objects, we identified a sample of 44 blazar candidates, including 27 flat spectrum radio quasars and 17 BL Lacertae objects, respectively. All the blazar candidates from our sample were identified with respect to their radio, optical, and mid-infrared properties. The newly selected blazar candidates possess the long-term, multi-colour photometric data from the Optical Gravitational Lensing Experiment, multi-colour mid-infrared observations, and archival radio data for one frequency at least. In addition, for nine of them the radio polarization data are available. With such data, these objects can be used to study the physics behind the blazar variability detected in the optical and mid-infrared bands, as a tool to investigate magnetic field geometry of the LMC and SMC, and as an exemplary sample of point like sources most likely detectable in $γ$-ray range with the newly emerging Cherenkov Telescope Array. △ Less

Submitted 10 September, 2018; originally announced September 2018.

Comments: 17 pages, 10 figures, accepted for publication in ApJ

Study of short term enhanced TeV $γ$- ray emission from Mrk 421 observed with TACTIC on December 28, 2014

Authors: K K Singh, K K Yadav, K Chanchalani, P Chandra, B Ghosal, A K Tickoo, R C Rannot, P Marandi, N K Agarwal, M Kothari, K K Gour, H C Goyal, A Goyal, N Kumar, C Borwankar, S R Kaul, V K Dhar, M K Koul, R Koul, K Venugopal, C K Bhat, N Chouhan

Abstract: In this work, we report on the detection of enhanced TeV gamma- ray emission from the high synchrotron-peaked blazar Mrk 421 with the TACTIC telescope on the night of December 28, 2014 (MJD 57019). We use data from the TACTIC observations of Mrk 421 for one week during December 25-31, 2014 (MJD 57016-57022) in this study. The TACTIC observation on December 28, 2014 (MJD 57019) alone results in the… ▽ More In this work, we report on the detection of enhanced TeV gamma- ray emission from the high synchrotron-peaked blazar Mrk 421 with the TACTIC telescope on the night of December 28, 2014 (MJD 57019). We use data from the TACTIC observations of Mrk 421 for one week during December 25-31, 2014 (MJD 57016-57022) in this study. The TACTIC observation on December 28, 2014 (MJD 57019) alone results in the detection of 86+-17 gamma- ray like events from Mrk 421 with a statistical significance of 5.17sigma in a livetime of ~ 2.2 hours above an energy threshold of 0.85 TeV. The high statistics (higher than three Crab Units) of TeV photons enables us to study the very high energy (VHE) gamma- ray emission from the source at shorter timescales. A minimum variability timescale of ~ 0.72 days is obtained for the TeV gamma- ray emission from Mrk 421 during the above flaring activity of the source. We have also analyzed the contemporaneous data from Fermi-LAT to study the high energy (HE) gamma- ray emission during the high activity state of the source. The HE gamma- ray emission is observed to be increasing after the TeV flaring activity detected with the TACTIC. We also use other near simultaneous archival data available from the Swift-BAT in hard X-rays and from SPOL at Steward Observatory in optical V and R bands to characterize the multi-wavelength emission of Mrk 421 during the high activity state observed at TeV energies. The TeV gamma- ray emission observed on December 28, 2014 (MJD 57019) is found to be more prominant than the emissions in lower energy bands during the same period. The TeV gamma- ray observation of Mrk 421 in high activity state with the TACTIC telescope is also used to understand the physical mechanism for blazar emission under the frame work of the leptonic single zone synchrotron self Compton process. △ Less

Submitted 8 August, 2018; originally announced August 2018.

Comments: 21 pages, 5 figures, Accepted for publication in Astroparticle Physics

Characterizing optical variability of OJ 287 in 2016 - 2017

Authors: Alok C. Gupta, Haritma Gaur, Paul J. Wiita, A. Pandey, P. Kushwaha, S. M. Hu, O. M. Kurtanidze, E. Semkov, G. Damljanovic, A. Goyal, M. Uemura, A. Darriba, Xu Chen, O. Vince, M. F. Gu, Z. Zhang, R. Bachev, R. Chanishvili, R. Itoh, M. Kawabata, S. O. Kurtanidze, T. Nakaoka, M. G. Nikolashvili, L. Stawarz, A. Strigachev

Abstract: We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 1… ▽ More We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability. △ Less

Submitted 13 January, 2019; v1 submitted 11 March, 2018; originally announced March 2018.

Comments: 16 pages, 8 figures, 7 tables; Accepted for Publication to AJ

Journal ref: 2019, AJ, 157, 95

H.E.S.S. discovery of very high energy gamma-ray emission from PKS 0625-354

Authors: H. E. S. S. Collaboration, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun , et al. (231 additional authors not shown)

Abstract: PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012 during 5.5 hours. The source was detected above an energy threshold of 200 GeV at a significance level of 6.1$σ$. No significant variability is found in these observations. The source is well described with a power-law spectrum with photon index $Γ =2.84 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 Te… ▽ More PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012 during 5.5 hours. The source was detected above an energy threshold of 200 GeV at a significance level of 6.1$σ$. No significant variability is found in these observations. The source is well described with a power-law spectrum with photon index $Γ =2.84 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 TeV) $N_0(E_0)=(0.58 \pm 0.22_{stat} \pm 0.12_{syst})\times10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$. Multi-wavelength data collected with Fermi-LAT, Swift-XRT, Swift-UVOT, ATOM and WISE are also analysed. Significant variability is observed only in the Fermi-LAT $γ$-ray and Swift-XRT X-ray energy bands. Having a good multi-wavelength coverage from radio to very high energy, we performed a broadband modelling from two types of emission scenarios. The results from a one zone lepto-hadronic, and a multi-zone leptonic models are compared and discussed. On the grounds of energetics, our analysis favours a leptonic multi-zone model. Models associated to the X-ray variability constraint supports previous results suggesting a BL Lac nature of PKS 0625-354, with, however, a large-scale jet structure typical of a radio galaxy. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 13 pages, 6 figures, accepted for publication in MNRAS

A search for new supernova remnant shells in the Galactic plane with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy , et al. (241 additional authors not shown)

Abstract: A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A… ▽ More A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912+101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account. △ Less

Submitted 27 April, 2018; v1 submitted 18 January, 2018; originally announced January 2018.

Comments: Accepted for publication in A&A, 24 pages, 13 figures. Minor corrections in v2

Journal ref: A&A 612, A8 (2018)

Contributions of the High Energy Stereoscopic System (H.E.S.S.) to the 35th International Cosmic Ray Conference (ICRC), Busan, Korea

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjaniany, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy , et al. (234 additional authors not shown)

Abstract: Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: Index of H.E.S.S. conference proceedings to the 35th ICRC

Stochastic modeling of multiwavelength variability of the classical BL Lac object OJ 287 on timescales ranging from decades to hours

Authors: A. Goyal, L. Stawarz, S. Zola, V. Marchenko, M. Soida, K. Nilsson, S. Ciprini, A. Baran, M. Ostrowski, P. J. Wiita, Gopal-Krishna, A. Siemiginowska, M. Sobolewska, S. Jorstad, A. Marscher, M. F. Aller H. D. Aller T. Hovatta, D. B. Caton, D. Reichart, K. Matsumoto, K. Sadakane, K. Gazeas, M. Kidger, V. Piirola, H. Jermak, F. Alicavus , et al. (87 additional authors not shown)

Abstract: 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 $γ$-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) pr… ▽ More 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 $γ$-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 $γ$-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. △ Less

Submitted 10 July, 2018; v1 submitted 13 September, 2017; originally announced September 2017.

Comments: ApJ accepted

Measurement of the EBL spectral energy distribution using the VHE gamma-ray spectra of H.E.S.S. blazars

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjaniany, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy , et al. (234 additional authors not shown)

Abstract: Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extrag… ▽ More Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult to determine. An EBL-model independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon. △ Less

Submitted 20 July, 2017; v1 submitted 19 July, 2017; originally announced July 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 606, A59 (2017)

Constraints on particle acceleration in SS433/W50 from MAGIC and H.E.S.S. observations

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo , et al. (386 additional authors not shown)

Abstract: The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capabili… ▽ More The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses and precession of the circumstellar disk periodically covering the central binary system is expected to be at its minimum. The eastern and western SS433/W50 interaction regions are also examined. We aim to constrain some theoretical models previously developed for this system. We made use of dedicated observations from MAGIC and H.E.S.S. from 2006 to 2011 which were combined for the first time and accounted for a total effective observation time of 16.5 h. Gamma-ray attenuation does not affect the jet/medium interaction regions. The analysis of a larger data set amounting to 40-80 h, depending on the region, was employed. No evidence of VHE gamma-ray emission was found. Upper limits were computed for the combined data set. We place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Our findings suggest that the fraction of the jet kinetic power transferred to relativistic protons must be relatively small to explain the lack of TeV and neutrino emission from the central system. At the SS433/W50 interface, the presence of magnetic fields greater 10$μ$G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with energies up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions. △ Less

Submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in A&A, 9 pages, 2 figures

Journal ref: A&A 612, A14 (2018)

Characterising the VHE diffuse emission in the central 200 parsecs of our Galaxy with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjaniany, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy , et al. (234 additional authors not shown)

Abstract: The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density pro… ▽ More The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain for the first time a detailed characterisation of the correlation with matter and to search for additional features and individual gamma-ray sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and improved analysis techniques we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of the full ridge diffuse emission is compatible with the one previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE gamma-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13-0.11. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Journal ref: A&A 612, A9 (2018)

Systematic search for very-high-energy gamma-ray emission from bow shocks of runaway stars

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas , et al. (238 additional authors not shown)

Abstract: Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E $\geq 0.1$TeV) gamma rays. Subs… ▽ More Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E $\geq 0.1$TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the \textit{Fermi}-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1 % of the kinetic wind energy. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: 7 pages, 2 figures, 1 table. Accepted for publication in A&A

Journal ref: A&A 612, A12 (2018)

The population of TeV pulsar wind nebulae in the H.E.S.S. Galactic Plane Survey

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon , et al. (233 additional authors not shown)

Abstract: The nine-year H.E.S.S. Galactic Plane Survey (HGPS) yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also pre… ▽ More The nine-year H.E.S.S. Galactic Plane Survey (HGPS) yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spin-down power $\dot{E}$. This seems to be caused both by an increase of extension with decreasing $\dot{E}$, and hence with time, compatible with a power law $R_\mathrm{PWN}(\dot{E}) \sim \dot{E}^{-0.65 \pm 0.20}$, and by a mild decrease of TeV gamma-ray luminosity with decreasing $\dot{E}$, compatible with $L_{1-10\,\mathrm{TeV}} \sim \dot{E}^{0.59 \pm 0.21}$. We also find that the offsets of pulsars with respect to the wind nebula centres with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar offset is correlated with a high apparent TeV efficiency $L_{1-10\,\mathrm{TeV}}/\dot{E}$. In addition to 14 HGPS sources considered as firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that form likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theory. △ Less

Submitted 10 March, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: 28 pages, 12 figures, 8 tables. Accepted for publication in A&A

Journal ref: A&A 612, A2 (2018)

Multi-wavelength variability study of the classical BL Lac object PKS 0735+178 on timescales ranging from decades to minutes

Authors: Arti Goyal, Lukasz Stawarz, Michal Ostrowski, Valeri Larionov, Gopal-Krishna, Paul J. Wiita, Santosh Joshi, Marian Soida, Ivan Agudo

Abstract: We present the results of our power spectral analysis for the BL Lac object PKS 0735+178 utilizing the Fermi-LAT survey at high-energy $γ$-rays, several ground-based optical telescopes, and single-dish radio telescopes operating at GHz frequencies. The novelty of our approach is that, by combining long-term and densely sampled intra-night light curves in the optical regime, we were able to constru… ▽ More We present the results of our power spectral analysis for the BL Lac object PKS 0735+178 utilizing the Fermi-LAT survey at high-energy $γ$-rays, several ground-based optical telescopes, and single-dish radio telescopes operating at GHz frequencies. The novelty of our approach is that, by combining long-term and densely sampled intra-night light curves in the optical regime, we were able to construct for the first time the optical power spectrum of the blazar for a time domain extending from 23 years down to minutes. Our analysis reveals that: (i) the optical variability is consistent with a pure red noise, for which the power spectral density can well be approximated by a single power-law throughout the entire time domain probed; (ii) the slope of power spectral density at high-energy $γ$-rays ($\sim 1$), is significantly flatter than that found at radio and optical frequencies ($\sim 2$) within the corresponding time variability range; (iii) for the derived power spectra we did not detect any low-frequency flattening, nor do we see any evidence for cut-offs at the highest frequencies down to the noise floor levels due to measurement uncertainties. We interpret our findings in terms of a model where the blazar variability is generated by the underlying single stochastic process (at radio and optical frequencies), or a linear superposition of such processes (in the $γ$-ray regime). Along with the detailed PSD analysis, we also present the results of our extended (1998-2015) intra-night optical monitoring program and newly acquired optical photo-polarimetric data for the source. △ Less

Submitted 5 March, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: The Astrophysical Journal, in press

Multiwavelength study of VHE emission from Markarian 501 using TACTIC observations during April-May, 2012

Authors: P Chandra, K K Singh, R C Rannot, K K Yadav, H Bhatt, A K Tickoo, B Ghosal, M Kothari, K K Gaur, A Goyal, H C Goyal, N Kumar, P Marandi, N Chouhan, S Sahayanathan, K Chanchalani, N K Agarwal, V K Dhar, S R Kaul, M K Koul, R Koul, K Venugopal, C K Bhat, C Borwankar, J Bhagwan , et al. (1 additional authors not shown)

Abstract: We have observed Markarian 501 in Very High Energy (VHE) gamma-ray wavelength band for 70.6 hours from 15 April to 30 May, 2012 using TACTIC telescope. Detailed analysis of $\sim$66.3 hours of clean data revealed the presence of a TeV $γ$-ray signal (686$\pm$77 $γ$-ray events) from the source direction with a statistical significance of 8.89$σ$ above 850 GeV. Further, a total of 375 $\pm$ 47 $γ$-r… ▽ More We have observed Markarian 501 in Very High Energy (VHE) gamma-ray wavelength band for 70.6 hours from 15 April to 30 May, 2012 using TACTIC telescope. Detailed analysis of $\sim$66.3 hours of clean data revealed the presence of a TeV $γ$-ray signal (686$\pm$77 $γ$-ray events) from the source direction with a statistical significance of 8.89$σ$ above 850 GeV. Further, a total of 375 $\pm$ 47 $γ$-ray like events were detected in 25.2 hours of observation from 22 - 27 May, 2012 with a statistical significance of 8.05$σ$ indicating that the source has possibly switched over to a relatively high gamma-ray emission state. We have derived time-averaged differential energy spectrum of the state in the energy range 850 GeV - 17.24 TeV which fits well with a power law function of the form $dF/dE=f_{0}E^{-Γ}$ with $f_{0}= (2.27 \pm 0.38) \times 10 ^{-11} $ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ and $Γ=2.57 \pm 0.15$. In order to investigate the source state, we have also used almost simultaneous multiwavelength observations viz: high energy data collected by $\it{Fermi}$-LAT, X-ray data collected by $\it{Swift}$-XRT and MAXI, optical and UV data collected by $\it{Swift}$-UVOT, and radio data collected by OVRO, and reconstructed broad-band Spectral Energy Distribution (SED). The obtained SED supports leptonic model (homogeneous single zone) for VHE gamma-ray emission involving synchrotron and synchrotron self Compton (SSC) processes. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: 10

Journal ref: New Astronomy 54C (2017) pp. 42-51

A time dependent approach to model X-ray and $γ$--ray light curves of Mrk 421 observed during the flare in February 2010

Authors: K K Singh, S Sahayanathan, A Sinha, N Bhatt, A K Tickoo, K K Yadav, R C Rannot, P Chandra, K Venugopal, P Marandi, N Kumar, H C Goyal, A Goyal, N K Agarwal, M Kothari, K Chanchalani, V K Dhar, N Chouhan, C K Bhat, M K Koul, R Koul

Abstract: Strong X-ray and $γ$--ray flares have been detected in February 2010 from the high synchrotron peaked blazar Mrk 421 (z=0.031). With the motivation of understanding the physics involved in this flaring activity, we study the variability of the source in X-ray and $γ$--ray energy bands during the period February 10-23, 2010 (MJD 55237-55250). We use near simultaneous X-ray data collected by \emph{M… ▽ More Strong X-ray and $γ$--ray flares have been detected in February 2010 from the high synchrotron peaked blazar Mrk 421 (z=0.031). With the motivation of understanding the physics involved in this flaring activity, we study the variability of the source in X-ray and $γ$--ray energy bands during the period February 10-23, 2010 (MJD 55237-55250). We use near simultaneous X-ray data collected by \emph{MAXI}, \emph{Swift}-XRT and $γ$--ray data collected by \emph{Fermi}-LAT and \emph{TACTIC} along with the optical V-band observations by \emph{SPOL} at Steward Observatory. We observe that the variation in the one day averaged flux from the source during the flare is characterized by fast rise and slow decay. Besides, the TeV $γ$--ray flux shows a strong correlation with the X-ray flux, suggesting the former to be an outcome of synchrotron self Compton emission process. To model the observed X-ray and $γ$--ray light curves, we numerically solve the kinetic equation describing the evolution of particle distribution in the emission region. The injection of particle distribution into the emission region, from the putative acceleration region, is assumed to be a time dependent power law. The synchrotron and synchrotron self Compton emission from the evolving particle distribution in the emission region are used to reproduce the X-ray and $γ$--ray flares successfully. Our study suggests that the flaring activity of Mrk 421 can be an outcome of an efficient acceleration process associated with the increase in underlying non-thermal particle distribution. △ Less

Submitted 6 January, 2017; originally announced January 2017.

Comments: Accepted for publication in New Astronomy

Journal ref: New Astronomy 54 (2017) 24-29

Gamma-ray blazar spectra with H.E.S.S. II mono analysis: the case of PKS 2155-304 and PG 1553+113

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, F. Brun , et al. (311 additional authors not shown)

Abstract: The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft sp… ▽ More The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 13 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 600, A89 (2017)

First limits on the very-high energy gamma-ray afterglow emission of a fast radio burst: H.E.S.S. observations of FRB 150418

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas , et al. (244 additional authors not shown)

Abstract: Aims: Following the detection of the fast radio burst FRB150418 by the SUPERB project at the Parkes radio telescope, we aim to search for very-high energy gamma-ray afterglow emission. Methods: Follow-up observations in the very-high energy gamma-ray domain were obtained with the H.E.S.S. imaging atmospheric Cherenkov telescope system within 14.5 hours of the radio burst. Results: The obtained 1.4… ▽ More Aims: Following the detection of the fast radio burst FRB150418 by the SUPERB project at the Parkes radio telescope, we aim to search for very-high energy gamma-ray afterglow emission. Methods: Follow-up observations in the very-high energy gamma-ray domain were obtained with the H.E.S.S. imaging atmospheric Cherenkov telescope system within 14.5 hours of the radio burst. Results: The obtained 1.4 hours of gamma-ray observations are presented and discussed. At the 99 % C.L. we obtained an integral upper limit on the gamma-ray flux of (E>350 GeV) < 1.33 x 10^-8 m^-2s^-1. Differential flux upper limits as function of the photon energy were derived and used to constrain the intrinsic high-energy afterglow emission of FRB 150418. Conclusions: No hints for high-energy afterglow emission of FRB 150418 were found. Taking absorption on the extragalactic background light into account and assuming a distance of z = 0.492 based on radio and optical counterpart studies and consistent with the FRB dispersion, we constrain the gamma-ray luminosity at 1 TeV to L < 5.1 x 10^47 erg/s at 99% C.L. △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: accepted for publication in A&A

Journal ref: A&A 597, A115 (2017)

Deeper H.E.S.S. observations of Vela Junior (RX J0852.0$-$4622): Morphology studies and resolved spectroscopy

Authors: H. E. S. S. Collaboration, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arakawa, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon , et al. (238 additional authors not shown)

Abstract: Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and sp… ▽ More Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$ TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644. △ Less

Submitted 26 April, 2018; v1 submitted 6 November, 2016; originally announced November 2016.

Comments: 14 pages, 5 figures, 8 tables. Published in A&A

Journal ref: A&A 612, A7 (2018)

Characterizing the gamma-ray long-term variability of PKS 2155-304 with H.E.S.S. and Fermi-LAT

Authors: H. E. S. S. Collaboration, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian T. Andersson, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun, P. Brun, M. Bryan, T. Bulik , et al. (222 additional authors not shown)

Abstract: Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain.… ▽ More Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain. Over the course of ~9 yr of H.E.S.S observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index ß_VHE = 1.10 +0.10 -0.13) on time scales larger than one day. An analysis of 5.5 yr of HE Fermi LAT data gives consistent results (ß_HE = 1.20 +0.21 -0.23, on time scales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior (ß ~ 2) seen on shorter time scales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection. △ Less

Submitted 11 October, 2016; originally announced October 2016.

Comments: 11 pages, 16 figures

Journal ref: A&A 598, A39 (2017)

H.E.S.S. observations of RX J1713.7-3946 with improved angular and spectral resolution; evidence for gamma-ray emission extending beyond the X-ray emitting shell

Authors: H. E. S. S. Collaboration, :, H. Abdalla, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas , et al. (231 additional authors not shown)

Abstract: Supernova remnants exhibit shock fronts (shells) that can accelerate charged particles up to very high energies. In the past decade, measurements of a handful of shell-type supernova remnants in very-high-energy gamma rays have provided unique insights into the acceleration process. Among those objects, RX$~$J1713.7-3946 (also known as G347.3-0.5) has the largest surface brightness, allowing us in… ▽ More Supernova remnants exhibit shock fronts (shells) that can accelerate charged particles up to very high energies. In the past decade, measurements of a handful of shell-type supernova remnants in very-high-energy gamma rays have provided unique insights into the acceleration process. Among those objects, RX$~$J1713.7-3946 (also known as G347.3-0.5) has the largest surface brightness, allowing us in the past to perform the most comprehensive study of morphology and spatially resolved spectra of any such very-high-energy gamma-ray source. Here we present extensive new H.E.S.S. measurements of RX$~$J1713.7-3946, almost doubling the observation time compared to our previous publication. Combined with new improved analysis tools, the previous sensitivity is more than doubled. The H.E.S.S. angular resolution of $0.048^\circ$ ($0.036^\circ$ above 2 TeV) is unprecedented in gamma-ray astronomy and probes physical scales of 0.8 (0.6) parsec at the remnant's location. The new H.E.S.S. image of RX$~$J1713.7-3946 allows us to reveal clear morphological differences between X-rays and gamma rays. In particular, for the outer edge of the brightest shell region, we find the first ever indication for particles in the process of leaving the acceleration shock region. By studying the broadband energy spectrum, we furthermore extract properties of the parent particle populations, providing new input to the discussion of the leptonic or hadronic nature of the gamma-ray emission mechanism. △ Less

Submitted 9 February, 2017; v1 submitted 27 September, 2016; originally announced September 2016.

Comments: 32 pages, 12 figures, accepted for publication in A&A, v2 is final accepted version

Journal ref: A&A 612, A6 (2018)

H.E.S.S. limits on line-like dark matter signatures in the 100 GeV to 2 TeV energy range close to the Galactic Centre

Authors: H. E. S. S. Collaboration, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon , et al. (224 additional authors not shown)

Abstract: A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected wi… ▽ More A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% CL the presence of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross-section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given. △ Less

Submitted 29 September, 2016; v1 submitted 26 September, 2016; originally announced September 2016.

Comments: Accepted by PRL. 7 pages, 3 figures, 1 table

The supernova remnant W49B as seen with H.E.S.S. and Fermi-LAT

Authors: H. E. S. S. Collaboration, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun , et al. (231 additional authors not shown)

Abstract: The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). $γ$-ray observations of SNR/MC associations are a powerful tool to constrain the origin of Galactic cosmic-rays, as they can probe the acceleration of hadrons thro… ▽ More The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). $γ$-ray observations of SNR/MC associations are a powerful tool to constrain the origin of Galactic cosmic-rays, as they can probe the acceleration of hadrons through their interaction with the surrounding medium and subsequent emission of non-thermal photons. The detection of a $γ$-ray source coincident with W49B at very high energies (VHE; E > 100 GeV) with the H.E.S.S. Cherenkov telescopes is reported together with a study of the source with 5 years of Fermi-LAT high energy $γ$-ray (0.06 - 300 GeV) data. The smoothly-connected combined source spectrum, measured from 60 MeV to multi-TeV energies, shows two significant spectral breaks at $304\pm20$ MeV and $8.4_{-2.5}^{+2.2}$ GeV, the latter being constrained by the joint fit from the two instruments. The detected spectral features are similar to those observed in several other SNR/MC associations and are found to be indicative of $γ$-ray emission produced through neutral-pion decay. △ Less

Submitted 2 September, 2016; originally announced September 2016.

Comments: 10 pages, 6 figures. Accepted for publication in A&A

Journal ref: A&A 612, A5 (2018)

Search for dark matter annihilations towards the inner Galactic halo from 10 years of observations with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdallah, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, E. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon , et al. (223 additional authors not shown)

Abstract: The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using $γ$-ray observations towards the inner 300 parsecs of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likel… ▽ More The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using $γ$-ray observations towards the inner 300 parsecs of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likelihood method using Galactic center (GC) data accumulated by H.E.S.S. over the last 10 years (2004-2014), and does not show any significant $γ$-ray signal above background. Assuming Einasto and Navarro-Frenk-White DM density profiles at the GC, we derive upper limits on the annihilation cross section $\langle σv\rangle$. These constraints are the strongest obtained so far in the TeV DM mass range and improve upon previous limits by a factor 5. For the Einasto profile, the constraints reach $\langle σv\rangle$ values of $\rm 6\times10^{-26} cm^3s^{-1}$ in the $W^+W^-$ channel for a DM particle mass of 1.5 TeV, and $\rm 2\times10^{-26} cm^3s^{-1}$ in the $τ^+τ^-$ channel for 1 TeV mass. For the first time, ground-based $γ$-ray observations have reached sufficient sensitivity to probe $\langle σv\rangle$ values expected from the thermal relic density for TeV DM particles. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Comments: Accepted by PRL. Includes Supplemental Material. 7+4 pages, 2+4 figures, 1 table

Journal ref: Phys. Rev. Lett. 117, 111301 (2016)

A search for very high-energy flares from the microquasars GRS 1915+105, Circinus X-1, and V4641 Sgr using contemporaneous H.E.S.S. and RXTE observations

Authors: H. E. S. S. Collaboration, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun , et al. (222 additional authors not shown)

Abstract: Microquasars are potential $γ$-ray emitters. Indications of transient episodes of $γ$-ray emission were recently reported in at least two systems: Cyg X-1 and Cyg X-3. The identification of additional $γ$-ray-emitting microquasars is required to better understand how $γ$-ray emission can be produced in these systems. Theoretical models have predicted very high-energy (VHE) $γ$-ray emission from mi… ▽ More Microquasars are potential $γ$-ray emitters. Indications of transient episodes of $γ$-ray emission were recently reported in at least two systems: Cyg X-1 and Cyg X-3. The identification of additional $γ$-ray-emitting microquasars is required to better understand how $γ$-ray emission can be produced in these systems. Theoretical models have predicted very high-energy (VHE) $γ$-ray emission from microquasars during periods of transient outburst. Observations reported herein were undertaken with the objective of observing a broadband flaring event in the $γ$-ray and X-ray bands. Contemporaneous observations of three microquasars, GRS 1915+105, Circinus X-1, and V4641 Sgr, were obtained using the High Energy Spectroscopic System (H.E.S.S.) telescope array and the Rossi X-ray Timing Explorer (RXTE) satellite. X-ray analyses for each microquasar were performed and VHE $γ$-ray upper limits from contemporaneous H.E.S.S. observations were derived. No significant $γ$-ray signal has been detected in any of the three systems. The integral $γ$-ray photon flux at the observational epochs is constrained to be $I(> 560 {\rm\ GeV}) < 7.3\times10^{-13}$ cm$^{-2}$ s$^{-1}$, $I(> 560 {\rm\ GeV}) < 1.2\times10^{-12}$ cm$^{-2}$ s$^{-1}$, and $I(> 240 {\rm\ GeV}) < 4.5\times10^{-12}$ cm$^{-2}$ s$^{-1}$ for GRS 1915+105, Circinus X-1, and V4641 Sgr, respectively. The $γ$-ray upper limits obtained using H.E.S.S. are examined in the context of previous Cherenkov telescope observations of microquasars. The effect of intrinsic absorption is modelled for each target and found to have negligible impact on the flux of escaping $γ$-rays. When combined with the X-ray behaviour observed using RXTE, the derived results indicate that if detectable VHE $γ$-ray emission from microquasars is commonplace, then it is likely to be highly transient. △ Less

Submitted 18 July, 2016; v1 submitted 15 July, 2016; originally announced July 2016.

Comments: 23 pages, 12 figures, accepted for publication by A&A

Journal ref: A&A 612, A10 (2018)

Extended VHE gamma-ray emission towards SGR1806-20, LBV1806-20, and stellar cluster Cl*1806-20

Authors: H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun, P. Brun , et al. (219 additional authors not shown)

Abstract: Using the High Energy Spectroscopic System (H.E.S.S.) telescopes we have discovered a steady and extended very high-energy (VHE) gamma-ray source towards the luminous blue variable candidate LBV1806-20, massive stellar cluster Cl*1806-20, and magnetar SGR1806-20. The new VHE source, HESSJ1808-204, was detected at a statistical significance of >6sigma (post-trial) with a photon flux normalisation… ▽ More Using the High Energy Spectroscopic System (H.E.S.S.) telescopes we have discovered a steady and extended very high-energy (VHE) gamma-ray source towards the luminous blue variable candidate LBV1806-20, massive stellar cluster Cl*1806-20, and magnetar SGR1806-20. The new VHE source, HESSJ1808-204, was detected at a statistical significance of >6sigma (post-trial) with a photon flux normalisation $(2.9 \pm 0.4_{\rm stat} \pm 0.5_{\rm sys})\times 10^{-13}$\,ph\,cm$^{-2}$\,s$^{-1}$\,TeV$^{-1}$ at 1TeV and a power-law photon index of $2.3\pm0.2_{\rm stat}\pm 0.3_{\rm sys}$. The luminosity of this source (0.2 to 10 TeV; scaled to distance d=8.7kpc) is $L_{\rm VHE}\sim1.6 \times 10^{34}(d/{8.7kpc})^2$\,erg\,s$^{-1}$. The VHE gamma-ray emission is extended and is well fit by a single Gaussian with statistical standard deviation of $0.095^\circ \pm 0.015^\circ$. This extension is similar to that of the synchrotron radio nebula G10.0-0.3, which is thought to be powered by LBV1806-20. The VHE gamma-ray luminosity could be provided by the stellar wind luminosity of LBV1806-20 by itself and/or the massive star members of Cl*1806-20. Alternatively, magnetic dissipation (e.g. via reconnection) from SGR1806-20 can potentially account for the VHE luminosity. The origin and hadronic and/or leptonic nature of the accelerated particles responsible for HESSJ1808-204 is not yet clear. If associated with SGR1806-20, the potentially young age of the magnetar (650yr) can be used to infer the transport limits of these particles to match the VHE source size. This discovery provides new interest in the potential for high-energy particle acceleration from magnetars, massive stars, and/or stellar clusters. △ Less

Submitted 27 June, 2016; v1 submitted 16 June, 2016; originally announced June 2016.

Comments: Accepted for publication in A&A (This version reduced some latex in the title and abstract)

Journal ref: A&A 612, A11 (2018)

Acceleration of petaelectronvolt protons in the Galactic Centre

Authors: H. E. S. S. collaboration, :, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, E. O. Angüner, M. Backes, A. Balzer, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun, P. Brun, M. Bryan, T. Bulik , et al. (207 additional authors not shown)

Abstract: Galactic cosmic rays reach energies of at least a few Peta-electronvolts (1 PeV =$10^\mathbf{15}$ electron volts). This implies our Galaxy contains PeV accelerators (PeVatrons), but all proposed models of Galactic cosmic-ray accelerators encounter non-trivial difficulties at exactly these energies. Tens of Galactic accelerators capable of accelerating particle to tens of TeV (1 TeV =… ▽ More Galactic cosmic rays reach energies of at least a few Peta-electronvolts (1 PeV =$10^\mathbf{15}$ electron volts). This implies our Galaxy contains PeV accelerators (PeVatrons), but all proposed models of Galactic cosmic-ray accelerators encounter non-trivial difficulties at exactly these energies. Tens of Galactic accelerators capable of accelerating particle to tens of TeV (1 TeV =$10^\mathbf{12}$ electron volts) energies were inferred from recent gamma-ray observations. None of the currently known accelerators, however, not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays, have shown the characteristic tracers of PeV particles: power-law spectra of gamma rays extending without a cutoff or a spectral break to tens of TeV. Here we report deep gamma-ray observations with arcminute angular resolution of the Galactic Centre regions, which show the expected tracer of the presence of PeV particles within the central 10~parsec of the Galaxy. We argue that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts and an outflow from the Galactic Centre. Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last $\gtrsim 10^{6-7}$ years, and therefore should be considered as a viable alternative to supernova remnants as a source of PeV Galactic cosmic rays. △ Less

Submitted 24 March, 2016; originally announced March 2016.

Comments: 29 pages, 7 figures, 3 tables, in press, for official published article see http://www.nature.com/nature/journal/v531/n7595/full/nature17147.html

Journal ref: Nature 531, 476 (2016)

Detailed spectral and morphological analysis of the shell type SNR RCW 86

Authors: H. E. S. S. Collaboration, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, E. O. Angüner, M. Backes, A. Balzer, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun, P. Brun, M. Bryan, T. Bulik, J. Carr , et al. (210 additional authors not shown)

Abstract: Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies… ▽ More Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1~keV X-ray band, the 2-5~keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results:We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5~keV X-ray data reveals a correlation with the 0.4-50~TeV gamma-ray emission.The spectrum of RCW~86 is best described by a power law with an exponential cutoff at $E_{cut}=(3.5\pm 1.2_{stat})$ TeV and a spectral index of $Γ$~$1.6\pm 0.2$. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW~86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to $\sim$0.1\% of the initial kinetic energy of a Type I a supernova explosion. When using a hadronic model, a magnetic field of $B$~100$μ$G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E$^{-2}$ spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of $Γ_p$~1.7 would be required, which implies that ~$7\times 10^{49}/n_{cm^{-3}}$erg has been transferred into high-energy protons with the effective density $n_{cm^{-3}}=n/ 1$ cm^-3. This is about 10\% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1~cm^-3. △ Less

Submitted 18 January, 2016; originally announced January 2016.

Comments: accepted for publication by A&A

Discovery of a Highly Polarized Optical Microflare in the Blazar S5 0716+714 During 2014 WEBT Campaign

Authors: Gopal Bhatta, Arti Goyal, Michal Ostrowski, Łukasz Stawarz, A. A. Arkharov, Hiroshi Akitaya, Rumen Bachev, Erika Benítez, G. A. Borman, Daniele Carosati, Andy Cason, Goran Damljanovic, Sarah Dhalla, Antonio Frasca, David Hiriart, Shao Ming Hu, Ryosuke Itoh, Damian Jableka, Svetlana Jorstad, Koji Kawabata, Sergey Klimanov, Omar Kurtanidze, Valeri Larionov, Douglas Laurence, Giuseppe Leto , et al. (20 additional authors not shown)

Abstract: The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we pr… ▽ More The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings. △ Less

Submitted 31 July, 2015; v1 submitted 30 July, 2015; originally announced July 2015.

Comments: 9 pages, 4 Figures, ApJ Letter accepted

Discovery of variable VHE gamma-ray emission from the binary system 1FGL J1018.6-5856

Authors: H. E. S. S. Collaboration, :, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, E. O. Angüner, M. Backes, A. Balzer, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, E. Birsin, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, J. Bregeon, F. Brun, P. Brun, M. Bryan, T. Bulik , et al. (208 additional authors not shown)

Abstract: Re-observations with the H.E.S.S. telescope array of the very-high-energy (VHE) source HESS J1018-589 A coincident with the Fermi-LAT $γ$-ray binary 1FGL J1018.6-5856 have resulted in a source detection significance of more than 9$σ$, and the detection of variability ($χ^2$/$ν$ of 238.3/155) in the emitted $γ$-ray flux. This variability confirms the association of HESS J1018-589 A with the high-en… ▽ More Re-observations with the H.E.S.S. telescope array of the very-high-energy (VHE) source HESS J1018-589 A coincident with the Fermi-LAT $γ$-ray binary 1FGL J1018.6-5856 have resulted in a source detection significance of more than 9$σ$, and the detection of variability ($χ^2$/$ν$ of 238.3/155) in the emitted $γ$-ray flux. This variability confirms the association of HESS J1018-589 A with the high-energy $γ$-ray binary detected by Fermi-LAT, and also confirms the point-like source as a new very-high-energy binary system. The spectrum of HESS J1018-589 A is best fit with a power-law function with photon index $Γ= 2.20 \pm 0.14_{\rm stat} \pm 0.2_{\rm sys}$. Emission is detected up to ~20 TeV. The mean differential flux level is $(2.9 \pm 0.4)\times10^{-13}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV, equivalent to ~1% of the flux from the Crab Nebula at the same energy. Variability is clearly detected in the night-by-night lightcurve. When folded on the orbital period of 16.58 days, the rebinned lightcurve peaks in phase with the observed X-ray and high-energy phaseograms. The fit of the H.E.S.S. phaseogram to a constant flux provides evidence of periodicity at the level of 3$σ$. The shape of the VHE phaseogram and measured spectrum suggest a low inclination, low eccentricity system with a modest impact from VHE $γ$-ray absorption due to pair production ($τ$ < 1 at 300 GeV). △ Less

Submitted 9 March, 2015; originally announced March 2015.

Comments: Accepted for publication in A&A

Journal ref: A&A 577, A131 (2015)