Mildly relativistic shocks that are embedded in colliding magnetohydrodynamic flows are prime sites for relativistic particle acceleration and the production of strongly variable, polarized ...multi-wavelength emission from relativistic jet sources such as blazars and gamma-ray bursts. The principal energization mechanisms at these shocks are diffusive shock acceleration and shock drift acceleration. In recent work, we had self-consistently coupled shock acceleration and radiation transfer simulations in blazar jets in a basic one-zone scenario. These one-zone models revealed that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the hydromagnetic turbulence in the shock layer. In this paper, we expand our previous work by including full time dependence and treating two zones, one being the site of acceleration and the other a larger emission zone. This construction is applied to multi-wavelength flares of the flat-spectrum radio quasar (FSRQ) 3C 279, fitting snapshot SEDs and generating light curves that are consistent with observed variability timescales. We also present a generic study for the typical flaring behavior of the BL Lac object Mrk 501. The model predicts correlated variability across all wavebands, but cross-band time lags depending on the type of blazar (FSRQ versus BL Lac), as well as distinctive spectral hysteresis patterns in all wavelength bands, from millimeter radio waves to gamma-rays. These evolutionary signatures serve to provide diagnostics on the competition between acceleration and radiative cooling.
Diffusive shock acceleration (DSA) at relativistic shocks is widely thought to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud active galactic nuclei ...such as blazars. Such acceleration can produce the non-thermal particles that emit the broad-band continuum radiation that is detected from extragalactic jets. An important recent development for blazar science is the ability of Fermi-Large Area Telescope spectroscopy to pin down the shape of the distribution of the underlying non-thermal particle population. This paper highlights how multiwavelength spectra spanning optical to X-ray to gamma-ray bands can be used to probe diffusive acceleration in relativistic, oblique, magnetohydrodynamic (MHD) shocks in blazar jets. Diagnostics on the MHD turbulence near such shocks are obtained using thermal and non-thermal particle distributions resulting from detailed Monte Carlo simulations of DSA. These probes are afforded by the characteristic property that the synchrotron vF(sub v), peak energy does not appear in the gamma-ray band above 100 MeV. We investigate self-consistently the radiative synchrotron and inverse Compton signatures of the simulated particle distributions. Important constraints on the diffusive mean free paths of electrons, and the level of electromagnetic field turbulence are identified for three different case study blazars, Mrk 501, BL Lacertae and AO 0235+164. The X-ray excess of AO 0235+164 in a flare state can be modelled as the signature of bulk Compton scattering of external radiation fields, thereby tightly constraining the energy-dependence of the diffusion coefficient for electrons. The concomitant interpretations that turbulence levels decline with remoteness from jet shocks, and the probable significant role for non-gyroresonant diffusion, are posited.
Dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter through gamma-ray emission due to their proximity, lack of astrophysical backgrounds and high dark ...matter density. They are often used to place restrictive bounds on the dark matter annihilation cross section. In this paper, we analyze six years of Fermi-LAT gamma-ray data from 19 dSphs that are satellites of the Milky Way, and derive from a stacked analysis of 15 dSphs, robust 95% confidence level lower limits on the dark matter lifetime for several decay channels and dark matter masses between ~1 GeV and 10 TeV. Our findings are based on a bin-by-bin maximum likelihood analysis treating the J-factor as a nuisance parameter using the Pass 8 event class. Our constraints from this ensemble are among the most stringent and solid in the literature, and competitive with existing ones coming from the extragalactic gamma-ray background, galaxy clusters, AMS-02 cosmic ray data, Super-K and ICECUBE neutrino data, while rather insensitive to systematic uncertainties. In particular, among gamma-ray searches, we improve existing limits for dark matter decaying into bb (mu super(+)mu super(-)) for dark matter masses below ~30(200) GeV, demonstrating that dSphs are compelling targets for constraining dark matter decay lifetimes.
We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh & Timokhin (WT19). Motivated by the power-law fluence distributions of FRB 121102, ...we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol that directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates that the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, 3 decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs, which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low-energy character depends on the population model, while the high-energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at ∼1037-1040 erg with a low-energy cutoff at ∼1037 erg. Lastly, we consider the local fluence distribution of FRBs and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal a flattening of the global FRB fluence distribution and saturation in FRB rates.
Hard, nonthermal, persistent pulsed X-ray emission extending between 10 and ∼150 keV has been observed in nearly 10 magnetars. For inner-magnetospheric models of such emission, resonant inverse ...Compton scattering of soft thermal photons by ultrarelativistic charges is the most efficient production mechanism. We present angle-dependent upscattering spectra and pulsed intensity maps for uncooled, relativistic electrons injected in inner regions of magnetar magnetospheres, calculated using collisional integrals over field loops. Our computations employ a new formulation of the QED Compton scattering cross section in strong magnetic fields that is physically correct for treating important spin-dependent effects in the cyclotron resonance, thereby producing correct photon spectra. The spectral cutoff energies are sensitive to the choices of observer viewing geometry, electron Lorentz factor, and scattering kinematics. We find that electrons with energies 15 MeV will emit most of their radiation below 250 keV, consistent with inferred turnovers for magnetar hard X-ray tails. More energetic electrons still emit mostly below 1 MeV, except for viewing perspectives sampling field-line tangents. Pulse profiles may be singly or doubly peaked dependent on viewing geometry, emission locale, and observed energy band. Magnetic pair production and photon splitting will attenuate spectra to hard X-ray energies, suppressing signals in the Fermi-LAT band. The resonant Compton spectra are strongly polarized, suggesting that hard X-ray polarimetry instruments such as X-Calibur, or a future Compton telescope, can prove central to constraining model geometry and physics.
Magnetars are neutron stars characterized by strong surface magnetic fields generally exceeding the quantum critical value of 44.1 TG. High-energy photons propagating in their magnetospheres can be ...attenuated by QED processes like photon splitting and magnetic pair creation. In this paper, we compute the opacities due to photon splitting and pair creation by photons emitted anywhere in the magnetosphere of a magnetar. Axisymmetric, twisted dipole field configurations embedded in the Schwarzschild metric are treated. The paper computes the maximum energies for photon transparency that permit propagation to infinity in curved spacetime. Special emphasis is given to cases where photons are generated along magnetic field loops and/or in polar regions; these cases directly relate to resonant inverse Compton scattering models for the hard X-ray emission from magnetars and Comptonized soft gamma-ray emission from giant flares. We find that increases in magnetospheric twists raise or lower photon opacities, depending on both the emission locale and the competition between field-line straightening and field strength enhancement. Consequently, given the implicit spectral transparency of hard X-ray bursts and persistent "tail" emission of magnetars, photon splitting considerations constrain their emission region locales and the twist angle of the magnetosphere; these constraints can be probed by future soft gamma-ray telescopes such as COSI and AMEGO. The inclusion of twists generally increases the opaque volume of pair creation by photons above its threshold, except when photons are emitted in polar regions and approximately parallel to the field.
ABSTRACT
Over the last four decades, persistent and flaring emission of magnetars observed by various telescopes has provided us with a suite of light curves and spectra in soft and hard X-rays, with ...no emission yet detected above around 1 MeV. Attenuation of such high-energy photons by magnetic pair creation and photon splitting is expected to be active in the magnetospheres of magnetars, possibly accounting for the paucity of gamma-rays in their signals. This paper explores polarization-dependent opacities for these two QED processes in static vacuum dipole magnetospheres of highly magnetized neutron stars, calculating attenuation lengths and determining escape energies, which are the maximum photon energies for transparency out to infinity. The numerical trajectory integral analysis in flat and curved space–times provides upper bounds of a few MeV or less to the visible energies for magnetars for locales proximate to the stellar surface. Photon splitting opacity alone puts constraints on the possible emission locales in their magnetospheres: regions within field loops of maximum altitudes $\, r_{{\rm max}}\sim 2\!-\!4\,$ stellar radii are not commensurate with maximum detected energies of around 250 keV. These constraints apply not only to magnetar flares but also to their quiescent hard X-ray tail emission. An exploration of photon splitting attenuation in the context of a resonant inverse Compton scattering model for the hard X-ray tails derives distinctive phase-resolved spectroscopic and polarimetric signatures, of significant interest for future MeV-band missions such as AMEGO and e-ASTROGAM.
Multiwavelength follow-up of unidentified Fermi sources has vastly expanded the number of known galactic-field "black widow" and "redback" millisecond pulsar binaries. Focusing on their ...rotation-powered state, we interpret the radio to X-ray phenomenology in a consistent framework. We advocate the existence of two distinct modes differing in their intrabinary shock orientation, distinguished by the phase centering of the double-peaked X-ray orbital modulation originating from mildly relativistic Doppler boosting. By constructing a geometric model for radio eclipses, we constrain the shock geometry as functions of binary inclination and shock standoff R0. We develop synthetic X-ray synchrotron orbital light curves and explore the model parameter space allowed by radio eclipse constraints applied on archetypal systems B1957+20 and J1023+0038. For B1957+20, from radio eclipses the standoff is R0 ∼ 0.15-0.3 fraction of binary separation from the companion center, depending on the orbit inclination. Constructed X-ray light curves for B1957+20 using these values are qualitatively consistent with those observed, and we find occultation of the shock by the companion as a minor influence, demanding significant Doppler factors to yield double peaks. For J1023+0038, radio eclipses imply R0 0.4, while X-ray light curves suggest 0.1 R0 0.3 (from the pulsar). Degeneracies in the model parameter space encourage further development to include transport considerations. Generically, the spatial variation along the shock of the underlying electron power-law index should yield energy dependence in the shape of light curves, motivating future X-ray phase-resolved spectroscopic studies to probe the unknown physics of pulsar winds and relativistic shock acceleration therein.
We report on NICER observations of the magnetar SGR 1935+2154, covering its 2020 burst storm and long-term persistent emission evolution up to ∼90 days postoutburst. During the first 1120 s taken on ...April 28 00:40:58 UTC, we detect over 217 bursts, corresponding to a burst rate of >0.2 bursts s−1. Three hours later, the rate was 0.008 bursts s−1, remaining at a comparatively low level thereafter. The T90 burst duration distribution peaks at 840 ms; the distribution of waiting times to the next burst is fit with a lognormal with an average of 2.1 s. The 1-10 keV burst spectra are well fit by a blackbody, with an average temperature and area of kT = 1.7 keV and R2 = 53 km2. The differential burst fluence distribution over ∼3 orders of magnitude is well modeled with a power-law form dN/dF ∝ F−1.5 0.1. The source persistent emission pulse profile is double-peaked hours after the burst storm. We find that the burst peak arrival times follow a uniform distribution in pulse phase, though the fast radio burst associated with the source aligns in phase with the brighter peak. We measure the source spin-down from heavy-cadence observations covering days 21-39 postoutburst, Hz s−1, a factor of 2.7 larger than the value measured after the 2014 outburst. Finally, the persistent emission flux and blackbody temperature decrease rapidly in the early stages of the outburst, reaching quiescence 40 days later, while the size of the emitting area remains unchanged.
Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active galactic ...nuclei and gamma-ray bursts. These sources remain good candidate sites for the generation of ultrahigh energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are germane to the production of relativistic electrons and ions are outlined. The technique employed to identify these characteristics is a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, magnetohydrodynamic (MHD) shocks. Using a compact prescription for diffusion of charges in MHD turbulence, this approach generates particle angular and momentum distributions at any position upstream or downstream of the shock. Simulation output is presented for both small angle and large angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffmann-Teller frame does not exist. The distribution function power-law indices compare favorably with results from other techniques. They are found to depend sensitively on the mean magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along fields in shock environs. An interesting regime of flat-spectrum generation is addressed; we provide evidence for it being due to shock drift acceleration, a phenomenon well known in heliospheric shock studies. The impact of these theoretical results on blazar science is outlined. Specifically, Fermi Large Area Telescope gamma-ray observations of these relativistic jet sources are providing significant constraints on important environmental quantities for relativistic shocks, namely, the field obliquity, the frequency of scattering, and the level of field turbulence.