Showing 1–50 of 357 results for author: Pohl, M

Poynting flux transport channels formed in polar cap regions of neutron star magnetospheres

Authors: Jan Benáček, Andrey Timokhin, Patricio A. Muñoz, Axel Jessner, Tatiana Rievajová, Martin Pohl, Jörg Büchner

Abstract: Pair cascades in polar cap regions of neutron stars are considered to be an essential process in various models of coherent radio emissions of pulsars. The cascades produce pair plasma bunch discharges in quasi-periodic spark events. The cascade properties, and therefore also the coherent radiation, depend strongly on the magnetospheric plasma properties and vary significantly across and along the… ▽ More Pair cascades in polar cap regions of neutron stars are considered to be an essential process in various models of coherent radio emissions of pulsars. The cascades produce pair plasma bunch discharges in quasi-periodic spark events. The cascade properties, and therefore also the coherent radiation, depend strongly on the magnetospheric plasma properties and vary significantly across and along the polar cap. It is furthermore still uncertain from where the radio emission emanates in polar cap region. We investigate the generation of electromagnetic waves by pair cascades and their propagation in the polar cap for three representative inclination angles of a magnetic dipole, $0^\circ$, $45^\circ$, and $90^\circ$. 2D particle-in-cell simulations that include quantum-electrodynamic pair cascades are used in a charge limited flow from the star surface. We found that the discharge properties are strongly dependent on the magnetospheric current profile in the polar cap and that transport channels for high intensity Poynting flux are formed along magnetic field lines where the magnetospheric currents approach zero and where the plasma cannot carry the magnetospheric currents. There, the parallel Poynting flux component is efficiently transported away from the star and may eventually escape the magnetosphere as coherent radio waves. The Poynting flux decreases faster with the distance from the star in regions of high magnetospheric currents. Our model shows that no process of energy conversion from particles to waves is necessary for the coherent radio wave emission. Moreover, the pulsar radio beam does not have a cone structure, but rather the radiation generated by the oscillating electric gap fields directly escapes along open magnetic field lines in which no pair creation occurs. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: Submitted to Astronomy & Astrophysics

Supernova remnants of red supergiants: from barrels to Cygnus loops

Authors: D M A Meyer, P F Velazquez, M Pohl, K Egberts, M Petrov, M A Villagran, D F Torres, R Batzofin

Abstract: Core-collapse supernova remnants are the nebular leftover of defunct massive stars which have died during a supernova explosion, mostly while undergoing the red supergiant phase of their evolution. The morphology and emission properties of those remnants are a function of the distribution of circumstellar material at the moment of the supernova, the intrisic properties of the explosion, as well as… ▽ More Core-collapse supernova remnants are the nebular leftover of defunct massive stars which have died during a supernova explosion, mostly while undergoing the red supergiant phase of their evolution. The morphology and emission properties of those remnants are a function of the distribution of circumstellar material at the moment of the supernova, the intrisic properties of the explosion, as well as those of the ambient medium. By means of 2.5 dimensional numerical magnetohydrodynamics simulations, we model the long term evolution of supernova remnants generated by runaway rotating massive stars moving into a magnetised interstellar medium. Radiative transfer calculations reveal that the projected non-thermal emission of the supernova remnants decreases with time, i.e. older remnants are fainter than younger ones. Older (80 kyr) supernova remnants whose progenitors were moving with space velocity corresponding to a Mach number M = 1 (v_star = 20 km/s ) in the Galactic plane of the ISM (nISM = 1/cm3 ) are brighter in synchrotron than when moving with a Mach number M = 2 (v_star = 40 km/s ). We show that runaway red supergiant progenitors first induce an asymmetric non thermal 1.4 GHz barrel like synchrotron supernova remnants (at the age of about 8 kyr), before further evolving to adopt a Cygnus loop like shape (at about 80 kyr). It is conjectured that a significative fraction of supernova remnants are currently in this bilateral-to-Cygnus-loop evolutionary sequence, and that this should be taken into account in the data interpretation of the forthcoming Cherenkov Telescope Array (CTA) observatory. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted at Astronomie and Astrophysics

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

Robust inference of the Galactic centre gamma-ray excess spatial properties

Authors: Deheng Song, Christopher Eckner, Chris Gordon, Francesca Calore, Oscar Macias, Kevork N. Abazajian, Shunsaku Horiuchi, Manoj Kaplinghat, Martin Pohl

Abstract: The gamma-ray Fermi-LAT Galactic centre excess (GCE) has puzzled scientists for over 15 years. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and poi… ▽ More The gamma-ray Fermi-LAT Galactic centre excess (GCE) has puzzled scientists for over 15 years. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology - potentially arising from the annihilation of dark matter (DM) particles - with a boxy morphology - expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VVV survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a GALPROP-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology of the GCE, regardless of the mask applied to the Galactic plane. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 21 pages, 11 figures

Report number: LAPTH-009/24

Journal ref: 2024MNRAS.530.4395S

Nonlinear feedback of the electrostatic instability on the blazar-induced pair beam and GeV cascade

Authors: Mahmoud Alawashra, Martin Pohl

Abstract: Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade that is missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields or, alternatively, an energy loss of the beam due to the beam-plasma instability. Here, we study how t… ▽ More Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade that is missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields or, alternatively, an energy loss of the beam due to the beam-plasma instability. Here, we study how the beam-plasma instability feeds back on the beam, using a realistic two-dimensional beam distribution. We find that the instability broadens the beam opening angles significantly without any significant energy loss, thus confirming a recent feedback study on a simplified one-dimensional beam distribution. However, narrowing diffusion feedback of the beam particles with Lorentz factors less than $10^6$ might become relevant even though initially it is negligible. Finally, when considering the continuous creation of TeV pairs, we find that the beam distribution and the wave spectrum reach a new quasi-steady state, in which the scattering of beam particles persists and the beam opening angle may increase by a factor of hundreds. Understanding the implications on the GeV cascade emission requires accounting for inverse-Compton cooling. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: Accepted for publication in the Astrophysical Journal. 16 pages and 12 figures

An Angular Diameter Measurement of $β$ UMa via Stellar Intensity Interferometry with the VERITAS Observatory

Authors: A. Acharyya, J. P. Aufdenberg, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. J. Chromey, J. D. Davis, Q. Feng, G. M. Foote, A. Furniss, W. Hanlon, C. E. Hinrichs, J. Holder, W. Jin, P. Kaaret, M. Kertzman, D. Kieda, T. K. Kleiner, N. Korzoun, T. LeBohec, M. A. Lisa, M. Lundy, N. Matthews, C. E McGrath , et al. (22 additional authors not shown)

Abstract: We use the VERITAS imaging air Cherenkov Telescope (IACT) array to obtain the first measured angular diameter of $β$ UMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak ($β$ UMa), an A1-type subgiant… ▽ More We use the VERITAS imaging air Cherenkov Telescope (IACT) array to obtain the first measured angular diameter of $β$ UMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak ($β$ UMa), an A1-type subgiant, by comparing effective temperature and luminosity constraints to model stellar evolution tracks. Previous interferometric limb-darkened angular-diameter measurements of $β$ UMa in the near-infrared (CHARA Array, $1.149 \pm 0.014$ mas) and mid-infrared (Keck Nuller, $1.08 \pm 0.07$ mas), together with the measured parallax and bolometric flux, have constrained the effective temperature. This paper presents current VERITAS-SII observation and analysis procedures to derive squared visibilities from correlation functions. We fit the resulting squared visibilities to find a limb-darkened angular diameter of $1.07 \pm 0.04 {\rm (stat)} \pm 0.05$ (sys) mas, using synthetic visibilities from a stellar atmosphere model that provides a good match to the spectrum of $β$ UMa in the optical wave band. The VERITAS-SII limb-darkened angular diameter yields an effective temperature of $9700\pm200\pm 200$ K, consistent with ultraviolet spectrophotometry, and an age of $390\pm 29 \pm 32 $ Myr, using MESA Isochrones and Stellar Tracks (MIST). This age is consistent with $408 \pm 6$ Myr from the CHARA Array angular diameter. △ Less

Submitted 3 January, 2024; originally announced January 2024.

AM$^3$: An Open-Source Tool for Time-Dependent Lepto-Hadronic Modeling of Astrophysical Sources

Authors: Marc Klinger, Annika Rudolph, Xavier Rodrigues, Chengchao Yuan, Gaëtan Fichet de Clairfontaine, Anatoli Fedynitch, Walter Winter, Martin Pohl, Shan Gao

Abstract: We present the AM$^3$ (``Astrophysical Multi-Messenger Modeling'') software, which has been successfully used in the past to simulate the multi-messenger emission, including neutrinos, from active galactic nuclei, including the blazar sub-class, gamma-ray bursts, and tidal disruption events. AM$^3$ is a documented state-of-the-art open source software that efficiently solves the coupled integro-di… ▽ More We present the AM$^3$ (``Astrophysical Multi-Messenger Modeling'') software, which has been successfully used in the past to simulate the multi-messenger emission, including neutrinos, from active galactic nuclei, including the blazar sub-class, gamma-ray bursts, and tidal disruption events. AM$^3$ is a documented state-of-the-art open source software that efficiently solves the coupled integro-differential equations for the spectral and temporal evolution of the relevant particle densities (photons, electrons, positrons, protons, neutrons, pions, muons, and neutrinos). AM$^3$ includes all relevant non-thermal processes (synchrotron, inverse Compton scattering, photon-photon annihilation, proton-proton and proton-photon pion production, and photo-pair production). The software self-consistently calculates the full cascade of primary and secondary particles, outperforming simple test-particle approaches, and allows for non-linear feedback and predictions in the time domain. It also allows to track separately the contributions of different radiative processes to the overall photon and neutrino spectra, including the different hadronic interaction channels. With its efficient hybrid solver combining analytical and numerical techniques, AM$^3$ combines efficiency and accuracy at a user-adjustable level. We describe the technical details of the numerical framework and present examples of applications to various astrophysical environments. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: The BSD-licensed source code can be found under https://gitlab.desy.de/am3/am3, for the user guide and documentation see https://am3.readthedocs.io/en/latest/

VERITAS contributions to the 38th International Cosmic Ray Conference

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, Z. Hughes, T. B. Humensky, W. Jin, M. N. Johnson, M. Kertzman , et al. (39 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: html page. ICRC 2023, Nagoya, Japan

Kinetic simulations of non-relativistic high-Mach-number perpendicular shocks propagating in a turbulent medium

Authors: Karol Fulat, Artem Bohdan, Gabriel Torralba Paz, Martin Pohl

Abstract: Strong non-relativistic shocks are known to accelerate particles up to relativistic energies. However, for Diffusive Shock Acceleration electrons must have a highly suprathermal energy, implying a need for very efficient pre-acceleration. Most published studies consider shocks propagating through homogeneous plasma, which is an unrealistic assumption for astrophysical environments. Using 2D3V part… ▽ More Strong non-relativistic shocks are known to accelerate particles up to relativistic energies. However, for Diffusive Shock Acceleration electrons must have a highly suprathermal energy, implying a need for very efficient pre-acceleration. Most published studies consider shocks propagating through homogeneous plasma, which is an unrealistic assumption for astrophysical environments. Using 2D3V particle-in-cell simulations, we investigate electron acceleration and heating processes at non-relativistic high-Mach-number shocks in electron-ion plasma with a turbulent upstream medium. For this purpose slabs of plasma with compressive turbulence are separately simulated and then inserted into shock simulations, which requires matching of the plasma slabs at the interface. Using a novel procedure of matching electromagnetic fields and currents, we perform simulations of perpendicular shocks setting different intensities of density fluctuations ($\lesssim 10\%$) in the upstream. The new simulation technique provides a framework for studying shocks propagating in turbulent media. We explore the impact of the fluctuations on electron heating, the dynamics of upstream electrons, and the driving of plasma instabilities. Our results indicate that while the presence of the turbulence enhances variations in the upstream magnetic field, their levels remain too low to influence significantly the behavior of electrons at perpendicular shocks. △ Less

Submitted 22 October, 2023; originally announced October 2023.

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

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

A multi-wavelength investigation of PSR J2229+6114 and its pulsar wind nebula in the radio, X-ray, and gamma-ray bands

Authors: I. Pope, K. Mori, M. Abdelmaguid, J. D. Gelfand, S. P. Reynolds, S. Safi-Harb, C. J. Hailey, H. An, VERITAS Collaboration, :, P. Bangale, P. Batista, W. Benbow, J. H. Buckley, M. Capasso, J. L. Christiansen, A. J. Chromey, A. Falcone, Q. Feng, J. P. Finley, G. M Foote, G. Gallagher, W. F Hanlon, D. Hanna, O. Hervet , et al. (35 additional authors not shown)

Abstract: G106.3$+$2.7, commonly considered a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct ("head" & "tail") regions in the radio band. A discovery of very-high-energy (VHE) gamma-ray emission ($E_γ> 100$ GeV) followed by the recent detection of ultra-high-energy (UHE) gamma-ray emission ($E_γ> 100$ TeV) from the tail region suggests tha… ▽ More G106.3$+$2.7, commonly considered a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct ("head" & "tail") regions in the radio band. A discovery of very-high-energy (VHE) gamma-ray emission ($E_γ> 100$ GeV) followed by the recent detection of ultra-high-energy (UHE) gamma-ray emission ($E_γ> 100$ TeV) from the tail region suggests that G106.3$+$2.7 is a PeVatron candidate. We present a comprehensive multi-wavelength study of the Boomerang PWN (100" around PSR J2229+6114) using archival radio and Chandra data obtained from two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with $Γ= 1.52\pm0.06$ up to 20 keV. Contrary to the previous radio study by Kothes et al. 2006, we prefer a much lower PWN B-field ($B\sim3$ $μ$G) and larger distance ($d \sim 8$ kpc) based on (1) the non-varying X-ray flux over the last two decades, (2) the energy-dependent X-ray PWN size resulting from synchrotron burn-off and (3) the multi-wavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock $\sim 1000$ years ago. In this case, the head region should be formed by GeV--TeV electrons injected earlier by the pulsar propagating into the low density environment. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Multiwavelength Observations of the Blazar PKS~0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: M. Pohl, K. Mori

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2 degrees away from the best-fit position of the IceCube neutrino event 211208A. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2 degrees away from the best-fit position of the IceCube neutrino event 211208A. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution, and the gamma-ray data from Fermi-LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both measurements provide strong constraints on leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed spectral energy distribution (SED). The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models, and it would allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. A numerical lepto-hadronic model with external target photons reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: Proceedings paper presented at the 38th International Cosmic Ray Conference (ICRC2023), held 26 July - 3 August, 2023, in Nagoya, Japan

Multiwavelength Observations of the Blazar PKS 0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, T. B. Humensky , et al. (185 additional authors not shown)

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ra… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution (SED), and the gamma-ray data from Fermi -LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both X-ray and gamma-ray measurements provide strong constraints on the leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed SED. The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models and allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. We show a numerical lepto-hadronic model with external target photons that reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 21 pages, 3 figures, accepted by ApJ

A VERITAS/Breakthrough Listen Search for Optical Technosignatures

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Pedro Batista, Wystan Benbow, Aryeh Brill, M Capasso, Manel Errando, Abraham Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Sean Griffin, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, T. Humensky, Weidong Jin, Philip Kaaret , et al. (43 additional authors not shown)

Abstract: The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectabl… ▽ More The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectable over interstellar distances. We report here on the analysis and results of dedicated VERITAS observations of Breakthrough Listen targets conducted in 2019 and 2020 and of archival VERITAS data collected since 2012. Thirty hours of dedicated observations of 136 targets and 249 archival observations of 140 targets were analyzed and did not reveal any signals consistent with a technosignature. The results are used to place limits on the fraction of stars hosting transmitting civilizations. We also discuss the minimum-pulse sensitivity of our observations and present VERITAS observations of CALIOP: a space-based pulsed laser onboard the CALIPSO satellite. The detection of these pulses with VERITAS, using the analysis techniques developed for our technosignature search, allows a test of our analysis efficiency and serves as an important proof-of-principle. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 15 pages, 7 figures

Evolution of Relativistic Pair Beams: Implications for Laboratory and TeV Astrophysics

Authors: Marvin Beck, Oindrila Ghosh, Florian Grüner, Martin Pohl, Carl B. Schroeder, Günter Sigl, Ryan D. Stark, Benno Zeitler

Abstract: Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam thr… ▽ More Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam through diffusion in momentum space. We chalk out this evolution using a Fokker-Planck approach in which the drift and the diffusion terms respectively describe these phenomena in a compact form. We present particle-in-cell simulations to trace the complete evolution of the unstable beam-plasma system for a generic narrow Gaussian pair beam for which the growth rate is reactive. We show that the instability leads to an energetic broadening of the pair beam, slowing down the instability growth in the linear phase, in line with the analytical and numerical solutions of the Fokker-Planck equation. Whereas in a laboratory experiment the change in the momentum distribution is an easily measured observable as a feedback of the instability, the consequence of diffusive broadening in an astrophysical scenario can be translated to an increase in the opening angle of the pair beam. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 20 pages, 10 figures

VERITAS discovery of very high energy gamma-ray emission from S3 1227+25 and multiwavelength observations

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Wystan Benbow, Aryeh Brill, Jodi Christiansen, Alisha Chromey, Manel Errando, Abe Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Greg Gallagher, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, Weidong Jin, Madalyn Johnson, Philip Kaaret , et al. (46 additional authors not shown)

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

Submitted 4 May, 2023; originally announced May 2023.

Comments: 18 pages, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Electron-Ion Temperature Ratio in Astrophysical Shocks

Authors: John C. Raymond, Parviz Ghavamian, Artem Bohdan, Dongsu Ryu, Jacek Niemiec, Lorenzo Sironi, Aaron Tran, Elena Amato, Masahiro Hoshino, Martin Pohl, Takanobu Amano, Federico Fiuza

Abstract: Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Her… ▽ More Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T$_e$/T$_p$ for supernova remnant shocks using H alpha profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the H alpha determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence and varying preshock conditions. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Search for Ultraheavy Dark Matter from Observations of Dwarf Spheroidal Galaxies with VERITAS

Authors: A. Acharyya, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, M. Baumgart, W. Benbow, J. H. Buckley, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. F. Hanlon, O. Hervet, J. Hoang, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (29 additional authors not shown)

Abstract: Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has bee… ▽ More Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has been suggested as an under-explored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS gamma-ray observatory. With 216 hours of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a $γ$-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle. △ Less

Submitted 24 February, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 10 pages, 7 figures

VERITAS and Fermi-LAT constraints on the Gamma-ray Emission from Superluminous Supernovae SN2015bn and SN2017egm

Authors: A. Acharyya, C. B. Adams, P. Bangale, W. Benbow, J. H. Buckley, M. Capasso, V. V. Dwarkadas, M. Errando, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, W. F Hanlon, O. Hervet, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (34 additional authors not shown)

Abstract: Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnet… ▽ More Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar's spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV - 100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 20 pages, 7 figures, 2 tables

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. H. Buckley, A. Brill, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, K. A Farrell, Q. Feng, J. P. Finley, G. M Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, C. Giuri, O. Gueta , et al. (64 additional authors not shown)

Abstract: The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four… ▽ More The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hours of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star forming galaxies, dwarf-spheroidal galaxies, and highly-variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications. △ Less

Submitted 13 January, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: Version with corrected author list

Journal ref: Research Notes of the AAS, Volume 7, Number 1, 2023

Pre-acceleration in the Electron Foreshock II: Oblique Whistler Waves

Authors: Paul J. Morris, Artem Bohdan, Martin S. Weidl, Michelle Tsirou, Karol Fulat, Martin Pohl

Abstract: Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large scal… ▽ More Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large scale particle-in-cell simulation that generates a single shock with parameters motivated from supernova remnants. Within the foreshock, reflected electrons excite the oblique whistler instability and produce electromagnetic whistler waves, which co-move with the upstream flow and as non-linear structures eventually reach radii of up to 5 ion-gyroradii. We show that the inner electromagnetic configuration of the whistlers evolves into complex non-linear structures bound by a strong magnetic field around 4 times the upstream value. Although these non-linear structures do not in general interact with co-spatial upstream electrons, they resonate with electrons that have been reflected at the shock. We show that they can scatter, or even trap, reflected electrons, confining around $0.8\%$ of the total upstream electron population to the region close to the shock where they can undergo substantial pre-acceleration. This acceleration process is similar to, yet approximately 3 times more efficient than, stochastic shock drift acceleration. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: 16 pages, 11 figures, accepted for publication in ApJ

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS expe… ▽ More MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Nonthermal Electron Acceleration at Collisionless Quasi-perpendicular Shocks

Authors: T. Amano, Y. Matsumoto, A. Bohdan, O. Kobzar, S. Matsukiyo, M. Oka, J. Niemiec, M. Pohl, M. Hoshino

Abstract: Shock waves propagating in collisionless heliospheric and astrophysical plasmas have been studied extensively over the decades. One prime motivation is to understand the nonthermal particle acceleration at shocks. Although the theory of diffusive shock acceleration (DSA) has long been the standard for cosmic-ray acceleration at shocks, plasma physical understanding of particle acceleration remains… ▽ More Shock waves propagating in collisionless heliospheric and astrophysical plasmas have been studied extensively over the decades. One prime motivation is to understand the nonthermal particle acceleration at shocks. Although the theory of diffusive shock acceleration (DSA) has long been the standard for cosmic-ray acceleration at shocks, plasma physical understanding of particle acceleration remains elusive. In this review, we discuss nonthermal electron acceleration mechanisms at quasi-perpendicular shocks, for which substantial progress has been made in recent years. The discussion presented in this review is restricted to the following three specific topics. The first is stochastic shock drift acceleration (SSDA), which is a relatively new mechanism for electron injection into DSA. The basic mechanism, related in-situ observations and kinetic simulations results, and how it is connected with DSA will be discussed. Second, we discuss shock surfing acceleration (SSA) at very high Mach number shocks relevant to young supernova remnants (SNRs). While the original proposal under the one-dimensional assumption is unrealistic, SSA has now been proven efficient by a fully three-dimensional kinetic simulation. Finally, we discuss the current understanding of the magnetized Weibel-dominated shock. Spontaneous magnetic reconnection of self-generated current sheets within the shock structure is an interesting consequence of Weibel-generated strong magnetic turbulence. We argue that high Mach number shocks with both Alfven and sound Mach numbers exceeding 20-40 will likely behave as a Weibel-dominated shock. Despite a number of interesting recent findings, the relative roles of SSDA, SSA, and magnetic reconnection for electron acceleration at collisionless shocks and how the dominant particle acceleration mechanisms change depending on shock parameters remain to be answered. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: To appear in Reviews of Modern Plasma Physics as an invited review

Modeling kilonova afterglows: Effects of the thermal electron population and interaction with GRB outflows

Authors: Vsevolod Nedora, Tim Dietrich, Masaru Shibata, Martin Pohl, Ludovica Crosato Menegazzi

Abstract: Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. We employ semi-analytic afterglow model,… ▽ More Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. We employ semi-analytic afterglow model, $\texttt{PyBlastAfterglow}$. We consider kilonova ejecta profiles from ab-initio numerical relativity binary neutron star merger simulations, targeted to GW170817. We do not perform model selection. We find that the emission from thermal electrons dominates at early times. If the interstellar medium density is high (${\simeq}0.1\,\ccm$) it adds an early time peak to the light curve. As ejecta decelerates the spectral and temporal indexes change in a characteristic way that, if observed, can be used to reconstruct the ejecta velocity distribution. For the low interstellar medium density, inferred for GRB 170817A, the emission from the non-thermal electron population generally dominates. We also assess how kilonova afterglow light curves change if the interstellar medium has been partially removed and pre-accelerated by laterally expanding gamma-ray burst ejecta. We find that the main effect is the emission suppression at early time ${\lesssim}10^{3}\,$days, and at its maximum it reaches ${\sim}40\%$ when the fast tail of the kilonova ejecta moves subsonically through the wake of laterally spreading gamma-ray burst ejecta. The subsequent rebrightening, when these ejecta break through and shocks form, is very mild (${\lesssim}10\%$), and may not be observable. △ Less

Submitted 13 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 20 pages, 14 figures

Rectangular core-collapse supernova remnants: application to Puppis A

Authors: D. M. -A. Meyer, P. F. Velazquez, O. Petruk, A. Chiotellis, M. Pohl, F. Camps-Farina, M. Petrov, E. M. Reynoso, J. C. Toledo-Roy, E. M. Schneiter, A. Castellanos-Ramirez, A. Esquivel

Abstract: Core-collapse supernova remnants are the gaseous nebulae of galactic interstellar media (ISM) formed after the explosive death of massive stars. Their morphology and emission properties depend both on the surrounding circumstellar structure shaped by the stellar wind-ISM interaction of the progenitor star and on the local conditions of the ambient medium. In the warm phase of the Galactic plane (n… ▽ More Core-collapse supernova remnants are the gaseous nebulae of galactic interstellar media (ISM) formed after the explosive death of massive stars. Their morphology and emission properties depend both on the surrounding circumstellar structure shaped by the stellar wind-ISM interaction of the progenitor star and on the local conditions of the ambient medium. In the warm phase of the Galactic plane (n = 1/cm3, T = 8000 K), an organised magnetic field of strength 7 microG has profound consequences on the morphology of the wind bubble of massive stars at rest. In this paper we show through 2.5D magneto-hydrodynamical simulations, in the context of a Wolf-Rayet-evolving 35 Mo star, that it affects the development of its supernova remnant. When the supernova remnant reaches its middle age (15 to 20 kyr), it adopts a tubular shape that results from the interaction between the isotropic supernova ejecta and the anisotropic, magnetised, shocked stellar progenitor bubble into which the supernova blast wave expands. Our calculations for non-thermal emission, i.e. radio synchrotron and inverse Compton radiation, reveal that such supernova remnants can, due to projection effects, appear as rectangular objects in certain cases. This mechanism for shaping a supernova remnant is similar to the bipolar and elliptical planetary nebula production by wind-wind interaction in the low-mass regime of stellar evolution. If such a rectangular core-collapse supernova remnant is created, the progenitor star must not have been a runaway star. We propose that such a mechanism is at work in the shaping of the asymmetric core-collapse supernova remnant Puppis A. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: Accepted at MNRAS

Pre-acceleration in the Electron Foreshock I: Electron Acoustic Waves

Authors: Paul J. Morris, Artem Bohdan, Martin S. Weidl, Martin Pohl

Abstract: To undergo diffusive shock acceleration, electrons need to be pre-accelerated to increase their energies by several orders of magnitude, else their gyro-radii are smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel or perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot t… ▽ More To undergo diffusive shock acceleration, electrons need to be pre-accelerated to increase their energies by several orders of magnitude, else their gyro-radii are smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel or perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot to a region called the electron foreshock. To determine the pre-acceleration in this region, we undertake PIC simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for $θ_{\rm Bn} = 74.3^\circ$, it increases to $R \sim 5\%$ for $θ_{\rm Bn} = 30^\circ$, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to $R^{0.6}$ and a wavelength $\approx 2 λ_{\rm{se}}$, where $λ_{\rm{se}}$ is the electron skin length. While the initial reflection mechanism is typically a combination of shock surfing acceleration and magnetic mirroring, we show that once the electrostatic waves have been generated upstream they themselves can increase the momenta of upstream electrons parallel to the magnetic field. In $\lesssim 1\%$ of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction are re-scattered back towards the shock after reflection, re-interact with the shock with energies much greater than thermal, and cross into the downstream. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ. 15 pages, 15 figures

The electron foreshock at high-Mach-number nonrelativistic oblique shocks

Authors: Artem Bohdan, Martin S. Weidl, Paul J. Morris, Martin Pohl

Abstract: In the Universe matter outside of stars and compact objects is mostly composed of collisionless plasma. The interaction of a supersonic plasma flow with an obstacle results in collisionless shocks that are often associated with intense nonthermal radiation and the production of cosmic ray particles. Motivated by simulations of non-relativistic high-Mach-number shocks in supernova remnants, we inve… ▽ More In the Universe matter outside of stars and compact objects is mostly composed of collisionless plasma. The interaction of a supersonic plasma flow with an obstacle results in collisionless shocks that are often associated with intense nonthermal radiation and the production of cosmic ray particles. Motivated by simulations of non-relativistic high-Mach-number shocks in supernova remnants, we investigate the instabilities excited by relativistic electron beams in the extended foreshock of oblique shocks. The phase-space distributions in the inner and outer foreshock regions are derived with a Particle-in-Cell simulation of the shock and used as initial conditions for simulations with periodic boundary conditions to study their relaxation towards equilibrium. We find that the observed electron-beam instabilities agree very well with the predictions of a linear dispersion analysis: the electrostatic electron-acoustic instability dominates in the outer region of the foreshock, while the denser electron beams in the inner foreshock drive the gyroresonant oblique-whistler instability. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: accepted for publication in Physics of Plasmas. 14 pages, 12 figures

Assessing the Impact of Hydrogen Absorption on the Characteristics of the Galactic Center Excess

Authors: Martin Pohl, Oscar Macias, Phaedra Coleman, Chris Gordon

Abstract: We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals… ▽ More We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals seen with the H$I$4PI survey and, b), the accounting for gas that primarily manifests itself through absorption. We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter. We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H$I$ model is used to estimate the cosmic-ray induced diffuse gamma-ray emission. The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence for a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favoured by the data. This finding is robust under the variation of various parameters, for example the excitation temperature of atomic hydrogen, and a number of tests for systematic issues. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: ApJ in press, 19 pages

Spectral softening in core-collapse supernova remnant expanding inside wind-blown bubble

Authors: Samata Das, Robert Brose, Dominique M. -A. Meyer, Martin Pohl, Iurii Sushch, Pavlo Plotko

Abstract: Context. Galactic cosmic rays are widely assumed to arise from diffusive shock acceleration, specifically at shocks in supernova remnants (SNRs). These shocks expand in a complex environment, particularly in the core-collapse scenario as these SNRs evolve inside the wind-blown bubbles created by their progenitor stars. The cosmic rays (CRs) at core-collapse SNRs may carry spectral signatures of th… ▽ More Context. Galactic cosmic rays are widely assumed to arise from diffusive shock acceleration, specifically at shocks in supernova remnants (SNRs). These shocks expand in a complex environment, particularly in the core-collapse scenario as these SNRs evolve inside the wind-blown bubbles created by their progenitor stars. The cosmic rays (CRs) at core-collapse SNRs may carry spectral signatures of that complexity. Aims. We study particle acceleration in the core-collapse SNR of a progenitor with initial mass 60 $M_\odot$ and realistic stellar evolution. The SNR shock interacts with discontinuities inside the wind-blown bubble and generates several transmitted and reflected shocks. We analyse their impact on particle spectra and the resulting emission from the remnant. Methods. The hydrodynamic equations for the evolution of SNR inside the pre-supernova circumstellar medium have been solved simultaneously with the transport equation for cosmic rays in test-particle approximation and with the induction equation for the magnetohydrodynamics (MHD) in 1-D spherical symmetry. Results. The evolution of core-collapse SNRs inside complex wind-blown bubbles modifies the spectra of both the particles and their emission. We have found softer particle spectra with spectral indices close to 2.5 during shock propagation inside the shocked wind, and this softness persists at later evolutionary stages. Further, our calculated total production spectrum released into the interstellar medium demonstrates spectral consistency at high energy with the galactic CRs injection spectrum, required in propagation models. The magnetic field structure effectively influences the emission morphology of SNR as it governs the transportation of particles and the synchrotron emissivity. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 14 pages, 13 figures, Accepted in Astronomy & Astrophysics

Journal ref: A&A 661, A128 (2022)

Suppression of the TeV pair-beam plasma instability by a tangled weak intergalactic magnetic field

Authors: Mahmoud Alawashra, Martin Pohl

Abstract: We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which ope… ▽ More We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which operates the best in the absence of IGMF. Considering IGMF with correlation lengths smaller than a kpc, we find that weak magnetic fields increase the transverse momentum of the pair beam particles, which dramatically reduces the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. We show that the beam-plasma instability is eliminated as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 11 pages, 5 figures

Diffusive shock acceleration at oblique high Mach number shocks

Authors: Allard Jan van Marle, Artem Bohdan, Paul J. Morris, Martin Pohl, Alexandre Marcowith

Abstract: The current paradigm of cosmic ray (CR) origin states that the most part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) t… ▽ More The current paradigm of cosmic ray (CR) origin states that the most part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. Then, with this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate the upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely triggered at high Mach number shocks. △ Less

Submitted 4 March, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

MSC Class: 85-08 ACM Class: J.2

Leptonic non-thermal emission from supernova remnants evolving in the circumstellar magnetic field

Authors: Iurii Sushch, Robert Brose, Martin Pohl, Pavlo Plotko, Samata Das

Abstract: The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of Supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse… ▽ More The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of Supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse Compton scattering on ambient photon fields) are responsible for this emission. For a successful diagnostic, a good understanding of the spatial and energy distribution of the underlying particle population is crucial. Most SNRs are created in core-collapse explosions and expand into the wind bubble of their progenitor stars. This circumstellar medium features a complex spatial distribution of gas and magnetic field which naturally strongly affects the resulting particle population. In this work, we conduct a detailed study of the spectro-spatial evolution of the electrons accelerated at the forward shock of core-collapse SNRs and their non-thermal radiation, using the RATPaC code that is designed for the time- and spatially dependent treatment of particle acceleration at SNR shocks. We focus on the impact of the spatially inhomogeneous magnetic field through the efficiency of diffusion and synchrotron cooling. It is demonstrated that the structure of the circumstellar magnetic field can leave strong signatures in the spectrum and morphology of the resulting non-thermal emission. △ Less

Submitted 23 November, 2021; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJ

The throughput calibration of the VERITAS telescopes

Authors: C. B. Adams, W. Benbow, A. Brill, J. H. Buckley, J. L. Christiansen, A. Falcone, Q. Feng, J. P. Finley, G. M Foote, L. Fortson, A. Furniss, C. Giuri, D. Hanna, T. Hassan, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin, P. Kaaret, T. K Kleiner, S. Kumar, M. J. Lang, M. Lundy, G. Maier , et al. (20 additional authors not shown)

Abstract: Context. The response of imaging atmospheric Cherenkov telescopes to incident γ-ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These aging processes affect the reconstructed energies of the events and γ-ray fluxes. Aims. This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescop… ▽ More Context. The response of imaging atmospheric Cherenkov telescopes to incident γ-ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These aging processes affect the reconstructed energies of the events and γ-ray fluxes. Aims. This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) to account for changes in the optical throughput and detector performance over time. Methods. The total throughput of a Cherenkov telescope is the product of camera-dependent factors, such as the photomultiplier tube gains and their quantum efficiencies, and the mirror reflectivity and Winston cone response to incoming radiation. This document summarizes different methods to determine how the camera gains and mirror reflectivity have evolved over time and how we can calibrate this changing throughput in reconstruction pipelines for imaging atmospheric Cherenkov telescopes. The implementation is validated against seven years of observations with the VERITAS telescopes of the Crab Nebula, which is a reference object in very-high-energy astronomy. Results. Regular optical throughput monitoring and the corresponding signal calibrations are found to be critical for the reconstruction of extensive air shower images. The proposed implementation is applied as a correction to the signals of the photomultiplier tubes in the telescope simulation to produce fine-tuned instrument response functions. This method is shown to be effective for calibrating the acquired γ-ray data and for recovering the correct energy of the events and photon fluxes. At the same time, it keeps the computational effort of generating Monte Carlo simulations for instrument response functions affordably low. △ Less

Submitted 15 November, 2021; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: 19 pages, 20 figures, 2 tables Received: 22 September 2021 / Accepted: 03 November 2021

Journal ref: A&A 658, A83 (2022)

Variability and Spectral Characteristics of Three Flaring Gamma-ray Quasars Observed by VERITAS and Fermi-LAT

Authors: C. B. Adams, J. Batshoun, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, B. Cavins, J. L. Christiansen, P. Coppi, M. Errando, K. A Farrell, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, A. Gent, C. Giuri, D. Hanna, T. Hassan, O. Hervet, J. Holder, M. Houck, T. B. Humensky, W. Jin , et al. (41 additional authors not shown)

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

Submitted 25 October, 2021; originally announced October 2021.

Comments: 34 pages, 13 figures. Accepted for publication in the Astrophysical Journal

Multi-Wavelength Observation Campaign of the TeV Gamma-Ray Binary HESS J0632+057 with NuSTAR, VERITAS, MDM, and Swift

Authors: Y. M. Tokayer, H. An, J. P. Halpern, J. Kim, K. Mori, C. J. Hailey, C. B. Adams, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, M. Errando, A. Falcone, K. A Farrell, G. M Foote, L. Fortson, A. Furniss, A. Gent, C. Giuri, D. Hanna, T. Hassan, O. Hervet, J. Holder, B. Hona, T. B. Humensky , et al. (31 additional authors not shown)

Abstract: HESS J0632+057 belongs to a rare subclass of binary systems which emits gamma-rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp ``primary'' peak and broader ``secondary'' peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between Dec. 2019 and Feb. 2020, when the orbital phase ($φ$) is bet… ▽ More HESS J0632+057 belongs to a rare subclass of binary systems which emits gamma-rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp ``primary'' peak and broader ``secondary'' peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between Dec. 2019 and Feb. 2020, when the orbital phase ($φ$) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multi-wavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from Oct. 2019 to Mar. 2020, the MDM observatory did not detect significant variation in H$α$ and H$β$ line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light curve data with an intra-binary shock model constrained the orbital parameters, suggesting two orbital phases (at $φ_D = 0.13$ and 0.37) where the pulsar crosses the Be-disk, as well as phases for the periastron ($φ_0 = 0.30$) and inferior conjunction ($φ_{\text{IFC}} = 0.75$). The broad-band X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases. △ Less

Submitted 3 October, 2021; originally announced October 2021.

Comments: 21 pages, 9 figures, accepted to ApJ

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

Authors: C. B. Adams, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, A. J. Chromey, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, J P. Finley, G. Foote, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin , et al. (387 additional authors not shown)

Abstract: The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these obs… ▽ More The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: Accepted for publication by The Astrophysical Journal

VERITAS contributions to the 37th International Cosmic Ray Conference

Authors: C. B. Adams, A. Archer, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin , et al. (36 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany. Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: html page. 2021 ICRC, Berlin, Germany

Time-dependent treatment of cosmic-ray spectral steepening due to turbulence driving

Authors: Martin Pohl

Abstract: Cosmic-ray acceleration at non-relativistic shocks relies on scattering by turbulence that the cosmic rays drive upstream of the shock. We explore the rate of energy transfer from cosmic rays to non-resonant Bell modes and the spectral softening it implies. Accounting for the finite time available for turbulence driving at supernova-remnant shocks yields a smaller spectral impact than found earlie… ▽ More Cosmic-ray acceleration at non-relativistic shocks relies on scattering by turbulence that the cosmic rays drive upstream of the shock. We explore the rate of energy transfer from cosmic rays to non-resonant Bell modes and the spectral softening it implies. Accounting for the finite time available for turbulence driving at supernova-remnant shocks yields a smaller spectral impact than found earlier with steady-state considerations. Generally, for diffusion scaling with the Bohm rate by a factor $η$, the change in spectral index is at most $η$ divided by the Alfvénic Mach number of the thermal sub-shock. For $M_\mathrm{A}\lesssim 50$ it is well below this limit. Only for very fast shocks and very efficient cosmic-ray acceleration the change in spectral index may reach $0.1$. For standard SNR parameters it is negligible. Independent confirmation is derived by considering the synchrotron energy losses of electrons: if intense nonthermal multi-keV emission is produced, the energy loss, and hence the spectral steepening, is very small for hadronic cosmic rays that produce TeV-band gamma-ray emission. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: ApJ in press, 8 pages

Morphology of supernova remnants and their halos

Authors: Robert Brose, Martin Pohl, Iurii Sushch

Abstract: Supernova remnants (SNRs) are known to accelerate particles to relativistic energies, on account of their nonthermal emission. The observational progress from radio to gamma-ray observations reveals more and more morphological features that need to be accounted for when modeling the emission from those objects. We use our time-dependent acceleration code RATPaC to study the formation of extended… ▽ More Supernova remnants (SNRs) are known to accelerate particles to relativistic energies, on account of their nonthermal emission. The observational progress from radio to gamma-ray observations reveals more and more morphological features that need to be accounted for when modeling the emission from those objects. We use our time-dependent acceleration code RATPaC to study the formation of extended gamma-ray halos around supernova remnants and the morphological implications that arise when the high-energetic particles start to escape from the SNRs. We performed spherically symmetric 1D simulations in which we simultaneously solved the transport equations for cosmic rays, magnetic turbulence, and the hydrodynamical flow of the thermal plasma. Our simulations span 25,000 years, thus covering the free-expansion and the Sedov-Taylor phase of the SNR's evolution. We find a strong difference in the morphology of the gamma-ray emission from SNRs at later stages dependent on the emission process. At early times, both the inverse-Compton and the Pion-decay morphology are shell-like. However, as soon as the maximum-energy of the freshly accelerated particles starts to fall, the inverse-Compton morphology starts to become center-filled, whereas the Pion-decay morphology keeps its shell-like structure. Escaping high-energy electrons start to form an emission halo around the SNR at this time. There are good prospects for detecting this spectrally hard emission with the future Cerenkov Telescope Array, as there are for detecting variations in the gamma-ray spectral index across the interior of the SNR. Further, we find a constantly decreasing nonthermal X-ray flux that makes a detection of X-ray unlikely after the first few thousand years of the SNR's evolution. The radio flux is increasing throughout the SNR's lifetime and changes from a shell-like to a more center-filled morphology later on. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: 10 pages, 7 figures, to be published in Astronomy and Astrophysics

Electron Acceleration at Rippled Low-Mach-number Shocks in High-beta Collisionless Cosmic Plasmas

Authors: Oleh Kobzar, Jacek Niemiec, Takanobu Amano, Masahiro Hoshino, Shuichi Matsukiyo, Yosuke Matsumoto, Martin Pohl

Abstract: Using large-scale fully-kinetic two-dimensional particle-in-cell simulations, we investigate the effects of shock rippling on electron acceleration at low-Mach-number shocks propagating in high-$β$ plasmas, in application to merger shocks in galaxy clusters. We find that the electron acceleration rate increases considerably when the rippling modes appear. The main acceleration mechanism is stochas… ▽ More Using large-scale fully-kinetic two-dimensional particle-in-cell simulations, we investigate the effects of shock rippling on electron acceleration at low-Mach-number shocks propagating in high-$β$ plasmas, in application to merger shocks in galaxy clusters. We find that the electron acceleration rate increases considerably when the rippling modes appear. The main acceleration mechanism is stochastic shock-drift acceleration, in which electrons are confined at the shock by pitch-angle scattering off turbulence and gain energy from the motional electric field. The presence of multi-scale magnetic turbulence at the shock transition and the region immediately behind the main shock overshoot is essential for electron energization. Wide-energy non-thermal electron distributions are formed both upstream and downstream of the shock. The maximum energy of the electrons is sufficient for their injection into diffusive shock acceleration. We show for the first time that the downstream electron spectrum has a~power-law form with index $p\approx 2.5$, in agreement with observations. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: 15 pages, 14 figures, to be published in ApJ

Origin and role of relativistic cosmic particles

Authors: A. Araudo, G. Morlino, B. Olmi, F. Acero, I. Agudo, R. Adam, R. Alves Batista, E. Amato, E. O. Anguner, L. A. Antonelli, Y. Ascasibar, C. Balazs, J. Becker Tjus, C. Bigongiari, E. Bissaldi, J. Bolmont, C. Boisson, P. Bordas, Ž. Bošnjak, A. M. Brown, M. Burton, N. Bucciantini, F. Cangemi, P. Caraveo, M. Cardillo , et al. (99 additional authors not shown)

Abstract: This white paper briefly summarizes the importance of the study of relativistic cosmic rays, both as a constituent of our Universe, and through their impact on stellar and galactic evolution. The focus is on what can be learned over the coming decade through ground-based gamma-ray observations over the 20 GeV to 300 TeV range. The majority of the material is drawn directly from "Science with the C… ▽ More This white paper briefly summarizes the importance of the study of relativistic cosmic rays, both as a constituent of our Universe, and through their impact on stellar and galactic evolution. The focus is on what can be learned over the coming decade through ground-based gamma-ray observations over the 20 GeV to 300 TeV range. The majority of the material is drawn directly from "Science with the Cherenkov Telescope Array", which describes the overall science case for CTA. We request that authors wishing to cite results contained in this white paper cite the original work. △ Less

Submitted 15 June, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted as input to ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium. arXiv admin note: text overlap with arXiv:1709.07997

An Archival Search for Neutron-Star Mergers in Gravitational Waves and Very-High-Energy Gamma Rays

Authors: C. B. Adams, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, C. Giuri, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, W. Jin, P. Kaaret , et al. (37 additional authors not shown)

Abstract: The recent discovery of electromagnetic signals in coincidence with neutron-star mergers has solidified the importance of multimessenger campaigns in studying the most energetic astrophysical events. Pioneering multimessenger observatories, such as LIGO/Virgo and IceCube, record many candidate signals below the detection significance threshold. These sub-threshold event candidates are promising ta… ▽ More The recent discovery of electromagnetic signals in coincidence with neutron-star mergers has solidified the importance of multimessenger campaigns in studying the most energetic astrophysical events. Pioneering multimessenger observatories, such as LIGO/Virgo and IceCube, record many candidate signals below the detection significance threshold. These sub-threshold event candidates are promising targets for multimessenger studies, as the information provided by them may, when combined with contemporaneous gamma-ray observations, lead to significant detections. Here we describe a new method that uses such candidates to search for transient events using archival very-high-energy gamma-ray data from imaging atmospheric Cherenkov telescopes (IACTs). We demonstrate the application of this method to sub-threshold binary neutron star (BNS) merger candidates identified in Advanced LIGO's first observing run. We identify eight hours of archival VERITAS observations coincident with seven BNS merger candidates and search them for TeV emission. No gamma-ray emission is detected; we calculate upper limits on the integral flux and compare them to a short gamma-ray burst model. We anticipate this search method to serve as a starting point for IACT searches with future LIGO/Virgo data releases as well as in other sub-threshold studies for multimessenger transients, such as IceCube neutrinos. Furthermore, it can be deployed immediately with other current-generation IACTs, and has the potential for real-time use that places minimal burden on experimental operations. Lastly, this method may serve as a pilot for studies with the Cherenkov Telescope Array, which has the potential to observe even larger fields of view in its divergent pointing mode. △ Less

Submitted 2 June, 2021; originally announced June 2021.

A Search for TeV Gamma-ray Emission from Pulsar Tails by VERITAS

Authors: Wystan Benbow, A. Brill, James Buckley, M Capasso, A Chromey, M. Errando, Abraham Falcone, K. A. Farrell, Qi Feng, J Finley, G. M. Foote, Lucy Fortson, Amy Furniss, Alasdair Gent, C Giuri, David Hanna, Tarek Hassan, Olivier Hervet, Jamie Holder, G Hughes, T. B. Humensky, Weidong Jin, Philip Kaaret, Oleg Kargaltsev, Mary P. Kertzman , et al. (34 additional authors not shown)

Abstract: We report on the search for very-high-energy gamma-ray emission from the regions around three nearby supersonic pulsars (PSR B0355+54, PSR J0357+3205 and PSR J1740+1000) that exhibit long X-ray tails. To date there is no clear detection of TeV emission from any pulsar tail that is prominent in X-ray or radio. We provide upper limits on the TeV flux, and luminosity, and also compare these limits wi… ▽ More We report on the search for very-high-energy gamma-ray emission from the regions around three nearby supersonic pulsars (PSR B0355+54, PSR J0357+3205 and PSR J1740+1000) that exhibit long X-ray tails. To date there is no clear detection of TeV emission from any pulsar tail that is prominent in X-ray or radio. We provide upper limits on the TeV flux, and luminosity, and also compare these limits with other pulsar wind nebulae detected in X-rays and the tail emission model predictions. We find that at least one of the three tails is likely to be detected in observations that are a factor of 2-3 more sensitive. The analysis presented here also has implications for deriving the properties of pulsar tails, for those pulsars whose tails could be detected in TeV. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Comments: 21 pages, 6 figures

Magneto-hydrodynamic simulations of young supernova remnants and their energy-conversion phase

Authors: O. Petruk, T. Kuzyo, S. Orlando, M. Pohl, R. Brose

Abstract: Supernova remnants (SNRs) can be rich sources of information on the parent SN explosion. Thus investigating the transition from the phase of SN to that of SNR can be crucial to link these two phases of evolution. Here we aim to study the early development of SNR in more details, paying the major attention to the transition from the early-expansion stage to the Sedov stage and the role played by ma… ▽ More Supernova remnants (SNRs) can be rich sources of information on the parent SN explosion. Thus investigating the transition from the phase of SN to that of SNR can be crucial to link these two phases of evolution. Here we aim to study the early development of SNR in more details, paying the major attention to the transition from the early-expansion stage to the Sedov stage and the role played by magnetic field in this transition. To this end, spherical magneto-hydrodynamic simulations of SNRs have been performed to study the evolution of magnetic field in young SNRs and explore a sequence of the SNR evolutionary stages in the pre-radiative epoch. Remnants of three supernova types are considered, namely, SNIa, SNIc and SNIIP, that covers a wide space of parameters relevant for SNRs. Changes in global characteristics and development of spatial distributions are analysed. It is shown that the radial component of magnetic field rapidly drops downstream of the forward shock. Therefore, the radially-aligned polarization patterns observed in few young SNRs cannot be reproduced in the one-dimensional MHD simulations. The period SNR takes for the transition from the earliest ejecta-driven phase to the Sedov phase is long enough, with its distinctive physical features, headed by the energy conversion from mostly kinetic one to a fixed ratio between the thermal and kinetic components. This transition worth to be distinguished as a phase in SNR evolutionary scheme. The updated sequence of stages in SNR evolution could be the free expansion (of gas) -- energy-conversion -- Sedov-Taylor -- post-adiabatic -- radiative. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Comments: Accepted by MNRAS

VERITAS Observations of the Galactic Center Region at Multi-TeV Gamma-Ray Energies

Authors: C. B. Adams, W. Benbow, A. Brill, R. Brose, M. Buchovecky, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, W. Jin, P. Kaaret , et al. (34 additional authors not shown)

Abstract: The Galactic Center (GC) region hosts a variety of powerful astronomical sources and rare astrophysical processes that emit a large flux of non-thermal radiation. The inner 375 pc x 600 pc region, called the Central Molecular Zone, is home to the supermassive black hole Sagittarius A*, massive cloud complexes, and particle accelerators such as supernova remnants. We present the results of our impr… ▽ More The Galactic Center (GC) region hosts a variety of powerful astronomical sources and rare astrophysical processes that emit a large flux of non-thermal radiation. The inner 375 pc x 600 pc region, called the Central Molecular Zone, is home to the supermassive black hole Sagittarius A*, massive cloud complexes, and particle accelerators such as supernova remnants. We present the results of our improved analysis of the very-high-energy (VHE) gamma-ray emission above 2 TeV from the GC using 125 hours of data taken with the VERITAS imaging-atmospheric Cherenkov telescope between 2010 and 2018. The central source VER J1745-290, consistent with the position of Sagittarius A*, is detected at a significance of 38 standard deviations above the background level $(38σ)$, and we report its spectrum and light curve. Its differential spectrum is consistent with a power law with exponential cutoff, with a spectral index of $2.12^{+0.22}_{-0.17}$, a flux normalization at 5.3 TeV of $1.27^{+0.22}_{-0.23}\times 10^{-13}$ TeV-1 cm-2 s-1, and cutoff energy of $10.0^{+4.0}_{-2.0}$ TeV. We also present results on the diffuse emission near the GC, obtained by combining data from multiple regions along the GC ridge which yield a cumulative significance of $9.5σ$. The diffuse GC ridge spectrum is best fit by a power law with a hard index of 2.19 $\pm$ 0.20, showing no evidence of a cutoff up to 40 TeV. This strengthens the evidence for a potential accelerator of PeV cosmic rays being present in the GC. We also provide spectra of the other sources in our field of view with significant detections, composite supernova remnant G0.9+0.1 and HESS J1746-285. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: 19 pages, 8 figures, Accepted for publication in Astrophysical Journal

Non-thermal radio supernova remnants of exiled Wolf-Rayet stars

Authors: D. M. -A. Meyer, M. Pohl, M. Petrov, L. Oskinova

Abstract: A signification fraction of Galactic massive stars (> 8Mo) are ejected from their parent cluster and supersonically sail away through the interstellar medium (ISM). The winds of these fast-moving stars blow asymmetric bubbles thus creating a circumstellar environment in which stars eventually die with a supernova explosion. The morphology of the resulting remnant is largely governed by the circums… ▽ More A signification fraction of Galactic massive stars (> 8Mo) are ejected from their parent cluster and supersonically sail away through the interstellar medium (ISM). The winds of these fast-moving stars blow asymmetric bubbles thus creating a circumstellar environment in which stars eventually die with a supernova explosion. The morphology of the resulting remnant is largely governed by the circumstellar medium of the defunct progenitor star. In this paper, we present 2D magneto-hydrodynamical simulations investigating the effect of the ISM magnetic field on the shape of the supernova remnants of a 35Mo star evolving through a Wolf-Rayet phase and running with velocity 20 and 40 km/s, respectively. A 7 microG ambient magnetic field is sufficient to modify the properties of the expanding supernova shock front and in particular to prevent the formation of filamentary structures. Prior to the supernova explosion, the compressed magnetic field in the circumstellar medium stabilises the wind/ISM contact discontinuity in the tail of the wind bubble. A consequence is a reduced mixing efficiency of ejecta and wind materials in the inner region of the remnant, where the supernova shock wave propagates. Radiative transfer calculations for synchrotron emission reveal that the non-thermal radio emission has characteristic features reflecting the asymmetry of exiled core-collapse supernova remnants from Wolf-Rayet progenitors. Our models are qualitatively consistent with the radio appearance of several remnants of high-mass progenitors, namely the bilateral G296.5+10.0 and the shell-type remnants CTB109 and Kes 17, respectively. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 15 pages, 10 figures. Accepted at MNRAS

Magnetic field amplification by the Weibel instability at planetary and astrophysical high-Mach-number shocks

Authors: Artem Bohdan, Martin Pohl, Jacek Niemiec, Paul J. Morris, Yosuke Matsumoto, Takanobu Amano, Masahiro Hoshino, Ali Sulaiman

Abstract: Collisionless shocks are ubiquitous in the Universe and often associated with strong magnetic field. Here we use large-scale particle-in-cell simulations of non-relativistic perpendicular shocks in the high-Mach-number regime to study the amplification of magnetic field within shocks. The magnetic field is amplified at the shock transition due to the ion-ion two-stream Weibel instability. The norm… ▽ More Collisionless shocks are ubiquitous in the Universe and often associated with strong magnetic field. Here we use large-scale particle-in-cell simulations of non-relativistic perpendicular shocks in the high-Mach-number regime to study the amplification of magnetic field within shocks. The magnetic field is amplified at the shock transition due to the ion-ion two-stream Weibel instability. The normalized magnetic-field strength strongly correlates with the Alfvénic Mach number. Mock spacecraft measurements derived from PIC simulations are fully consistent with those taken in-situ at Saturn's bow shock by the Cassini spacecraft. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: Accepted to PRL. 7 pages, 4 figure

Mildly relativistic magnetized shocks in electron-ion plasmas -- II. Particle acceleration and heating

Authors: Arianna Ligorini, Jacek Niemiec, Oleh Kobzar, Masanori Iwamoto, Artem Bohdan, Martin Pohl, Yosuke Matsumoto, Takanobu Amano, Shuichi Matsukiyo, Masahiro Hoshino

Abstract: Particle acceleration and heating at mildly relativistic magnetized shocks in electron-ion plasma are investigated with unprecedentedly high-resolution two-dimensional particle-in-cell simulations that include ion-scale shock rippling. Electrons are super-adiabatically heated at the shock, and most of the energy transfer from protons to electrons takes place at or downstream of the shock. We are t… ▽ More Particle acceleration and heating at mildly relativistic magnetized shocks in electron-ion plasma are investigated with unprecedentedly high-resolution two-dimensional particle-in-cell simulations that include ion-scale shock rippling. Electrons are super-adiabatically heated at the shock, and most of the energy transfer from protons to electrons takes place at or downstream of the shock. We are the first to demonstrate that shock rippling is crucial for the energization of electrons at the shock. They remain well below equipartition with the protons. The downstream electron spectra are approximately thermal with a limited supra-thermal power-law component. Our results are discussed in the context of wakefield acceleration and the modelling of electromagnetic radiation from blazar cores. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: 11 pages, 11 figures; This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review