Source author record

Matthew G. Baring

Matthew G. Baring appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2022arXiv

Exploring the MeV Sky with a Combined Coded Mask and Compton Telescope: The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds on and improves past and currently proposed missions at such energies. We outline the motivations for combining a coded mask and a Compton telescope and we define the scientific goals of such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which enable identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

preprint2022arXiv

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals.

preprint2022arXiv

Improving the low-energy transient sensitivity of AMEGO-X using single-site events

AMEGO-X, the All-sky Medium Energy Gamma-Ray Observatory eXplorer, is a proposed instrument designed to bridge the so-called "MeV gap" by surveying the sky with unprecedented sensitivity from ~100 keV to about one GeV. This energy band is of key importance for multi-messenger and multi-wavelength studies but it is nevertheless currently under-explored. AMEGO-X addresses this situation by proposing a design capable of detecting and imaging gamma rays via both Compton interactions and pair production processes. However, some of the objects that AMEGO-X will study, such as gamma-ray bursts and magnetars, extend to energies below ~100 keV where the dominant interaction becomes photoelectric absorption. These events deposit their energy in a single pixel of the detector. In this work we show how the ~3500 cm^2 effective area of the AMEGO-X tracker to events between ~25 keV to ~100 keV will be utilized to significantly improve its sensitivity and expand the energy range for transient phenomena. Although imaging is not possible for single-site events, we show how we will localize a transient source in the sky using their aggregate signal to within a few degrees. This technique will more than double the number of cosmological gamma-ray bursts seen by AMEGO-X, allow us to detect and resolve the pulsating tails of extragalactic magnetar giant flares, and increase the number of detected less-energetic magnetar bursts -- some possibly associated with fast radio bursts. Overall, single-site events will increase the sensitive energy range, expand the science program, and promptly alert the community of fainter transient events.

preprint2022arXiv

Intensity and Polarization Characteristics of Extended Neutron Star Surface Regions

The surfaces of neutron stars are sources of strongly polarized soft X rays due to the presence of strong magnetic fields. Radiative transfer mediated by electron scattering and free-free absorption is central to defining local surface anisotropy and polarization signatures. Scattering transport is strongly influenced by the complicated interplay between linear and circular polarizations. This complexity has been captured in a sophisticated magnetic Thomson scattering simulation we recently developed to model the outer layers of fully-ionized atmospheres in such compact objects, heretofore focusing on case studies of localized surface regions. Yet, the interpretation of observed intensity pulse profiles and their efficacy in constraining key neutron star geometry parameters is critically dependent upon adding up emission from extended surface regions. In this paper, intensity, anisotropy and polarization characteristics from such extended atmospheres, spanning considerable ranges of magnetic colatitudes, are determined using our transport simulation. These constitute a convolution of varied properties of Stokes parameter information at disparate surface locales with different magnetic field strengths and directions relative to the local zenith. Our analysis includes full general relativistic propagation of light from the surface to an observer at infinity. The array of pulse profiles for intensity and polarization presented highlights how powerful probes of stellar geometry are possible. Significant phase-resolved polarization degrees in the range of 10-60% are realized when summing over a variety of surface field directions. These results provide an important background for observations to be acquired by NASA's new IXPE X-ray polarimetry mission.

preprint2020arXiv

Persistent Emission Properties of SGR J1935+2154 During Its 2020 Active Episode

We present detailed spectral and temporal characteristics of the persistent X-ray emission of SGR J1935+2154 based on our XMM-Newton and Chandra observations taken in the aftermath of its April 2020 burst storm, during which hundreds of energetic X-ray bursts were emitted, including one associated with an extraordinary fast radio burst. We clearly detect the pulsed X-ray emission in the XMM-Newton data. An average spin-down rate of 1.6$\times$10$^{-11}$ s s$^{-1}$ is obtained using our spin period measurement combined with three earlier values reported from the same active episode. Our investigations of the XMM-Newton and Chandra spectra with a variety of phenomenological and physically-motivated models, concluded that the magnetic field topology of SGR J1935+2154 is most likely highly non-dipolar. The spectral models indicate that surface field strengths in somewhat localized regions substantially exceed the polar value of 4.4$\times$10$^{14}$ G inferred from a spin-down torque associated with a rotating magnetic dipole.

preprint2020arXiv

Simultaneous magnetic polar cap heating during a flaring episode from the magnetar 1RXS J170849.0-400910

During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of $\sim$2. The flare spectrum is thermal with an average temperature of 2.1 keV. Phase resolved spectroscopy show that the two peaks possess the same temperature, but differ in size. These observational results along with simple light curve modeling indicate that two identical antipodal spots, likely the magnetic poles, are heated simultaneously at the onset of the flare and for its full duration. Hence, the origin of the flare has to be connected to the global dipolar structure of the magnetar. This might best be achieved externally, via twists to closed magnetospheric dipolar field lines seeding bombardment of polar footpoint locales with energetic pairs. Approximately 1.86 hours following the onset of the flare, a short burst with its own 3-minute thermal tail occurred. The burst tail is also pulsating at the spin period of the source and phase-aligned with the flare profile, implying an intimate connection between the two phenomena. The burst may have been caused by a magnetic reconnection event in the same twisted dipolar field lines anchored to the surface hot spots, with subsequent return currents supplying extra heat to these polar caps.

preprint2020arXiv

The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). 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 which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates 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, $\lesssim 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 $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ 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 flattening of the global FRB fluence distribution and saturation in FRB rates.

preprint2019arXiv

Time-Dependent, Multi-Wavelength Shock Acceleration Models for Active Flares of 3C 279

Jets in blazars are an excellent forum for studying acceleration at relativistic shocks using the highly-variable emission seen across the electromagnetic spectrum. Our recent work on combining multi-wavelength leptonic emission models with simulated thermal+non-thermal distributions from shock acceleration theory has resulted in new insights into plasma conditions in blazars. This has demonstrated the ability to infer the cyclotron frequency, the plasma density and thus also the Alfven speed, thereby determining the rapidity of particle energization. An important inference was that MHD turbulence levels decline with remoteness from jet shocks. This paper outlines new results from our recent extension of this program to a two-zone, time-evolving construction, modeling together both extended, enhanced emission states from larger radiative regions, and prompt flare events from compact acceleration zones. These are applied to flares in the FSRQ blazar 3C 279 monitored by Fermi-LAT in gamma-rays in late 2013. With impulsive injection episodes from the shock zone, as the acceleration first proceeds and then abates, the radiative simulations obtain a pronounced spectral hardening in the optical and gamma-ray bands as the flare grows, followed by a softening during the decay phase. For 3C 279, while model radio and X-ray synchrotron flares are temporally correlated, there is a lag in both bands relative to GeV gamma rays and optical emission on timescales of a number of hours. This delay is governed by the short cooling time associated with the bright external Compton signal.

preprint2016arXiv

Learning About the Magnetar Swift J1834.9-0846 from its Wind Nebula

The first wind nebula around a magnetar was recently discovered in X-rays around Swift~J1834.9$-$0846. We study this magnetar's global energetics and the properties of its particle wind or outflows. At a distance of $\sim4\;$kpc, Swift~J1834.9$-$0846 is located at the center of the supernova remnant (SNR) W41 whose radius is $\sim 19\;$pc, an order of magnitude larger than that of the X-ray nebula ($\sim2\;$pc). The association with SNR W41 suggests a common age of $\sim5-100\;$kyr, while its spin-down age is $4.9$~kyr. A small natal kick velocity may partly explain why a wind nebula was detected around this magnetar but not around other magnetars, most of which appear to have larger kick velocities and may have exited their birth SNR. We find that the GeV and TeV source detected by Fermi/LAT and H.E.S.S., respectively, of radius $\sim11\;$pc is most likely of hadronic origin. The dynamics and internal structure of the nebula are examined analytically to explain the nebula's current properties. Its size may naturally correspond to the diffusion-dominated cooling length of the X-ray emitting $e^+e^-$ pairs. This may also account for the spectral softening of the X-ray emission from the nebula's inner to outer parts. Analysis of the X-ray synchrotron nebula implies that (i) the nebular magnetic field is $\gtrsim 11\;μ$G (and likely $\lesssim30\;μ$G), and (ii) the nebula is not powered predominantly by the magnetar's quiescent spin-down-powered MHD wind, but by other outflows that contribute most of its energy. The latter are most likely associated with the magnetar's bursting activity, and possibly dominated by outflows associated with its past giant flares. The energy source for the required outflows cannot be the decay of the magnetar's dipole field alone, and is most likely the decay of its much stronger internal magnetic field.

preprint2016arXiv

New Limits on the Dark Matter Lifetime from Dwarf Spheroidal Galaxies using Fermi-LAT

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 {\it 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 $\sim 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 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 $\bar{b}b$, ($μ^+μ^-$) for DM masses below $\sim 30\, (200)$~GeV, demonstrating that dSphs are compelling targets for constraining dark matter decay lifetimes.

preprint2016arXiv

Probing Acceleration and Turbulence at Relativistic Shocks in Blazar Jets

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 broadband continuum radiation that is detected from extragalactic jets. An important recent development for blazar science is the ability of Fermi-LAT spectroscopy to pin down the shape of the distribution of the underlying non-thermal particle population. This paper highlights how multi-wavelength 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 $νF_ν$ 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.

preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

preprint2014arXiv

An Observed Correlation Between Thermal and Non-Thermal Emission in Gamma-Ray Bursts

Recent observations by the $Fermi$ Gamma-ray Space Telescope have confirmed the existence of thermal and non-thermal components in the prompt photon spectra of some Gamma-ray bursts (GRBs). Through an analysis of six bright Fermi GRBs, we have discovered a correlation between the observed photospheric and non-thermal $γ$-ray emission components of several GRBs using a physical model that has previously been shown to be a good fit to the Fermi data. From the spectral parameters of these fits we find that the characteristic energies, $E_{\rm p}$ and $kT$, of these two components are correlated via the relation $E_{\rm p} \propto T^α$ which varies from GRB to GRB. We present an interpretation in which the value of index $α$ indicates whether the jet is dominated by kinetic or magnetic energy. To date, this jet composition parameter has been assumed in the modeling of GRB outflows rather than derived from the data.

preprint2014arXiv

Compton scattering in strong magnetic fields: Spin-dependent influences at the cyclotron resonance

The quantum electrodynamical (QED) process of Compton scattering in strong magnetic fields is commonly invoked in atmospheric and inner magnetospheric models of x-ray and soft gamma-ray emission in high-field pulsars and magnetars. A major influence of the field is to introduce resonances at the cyclotron frequency and its harmonics, where the incoming photon accesses thresholds for the creation of virtual electrons or positrons in intermediate states with excited Landau levels. At these resonances, the effective cross section typically exceeds the classical Thomson value by over 2 orders of magnitude. Near and above the quantum critical magnetic field of 44.13 TeraGauss, relativistic corrections must be incorporated when computing this cross section. This paper presents formalism for the QED magnetic Compton differential cross section valid for both subcritical and supercritical fields, yet restricted to scattered photons that are below pair creation threshold. Calculations are developed for the particular case of photons initially propagating along the field, mathematically simple specializations that are germane to interactions involving relativistic electrons frequently found in neutron star magnetospheres. This exposition of relativistic, quantum, magnetic Compton cross sections treats electron spin dependence fully, since this is a critical feature for describing the finite decay lifetimes of the intermediate states. The formalism employs both the Johnson and Lippmann (JL) wave functions and the Sokolov and Ternov (ST) electron eigenfunctions of the magnetic Dirac equation. The ST states are formally correct for self-consistently treating spin-dependent effects that are so important in the resonances. Relatively compact analytic forms for the cross sections are presented that will prove useful for astrophysical modelers.

preprint2014arXiv

Magnetic Pair Creation Transparency in Gamma-Ray Pulsars

Magnetic pair creation $γ\to e^+e^-$ has been at the core of radio pulsar paradigms and central to polar cap models of gamma-ray pulsars for over three decades. The Fermi gamma-ray pulsar population now exceeds 140 sources and has defined an important part of Fermi's science legacy. Among the population characteristics well established is the common occurrence of exponential turnovers in their spectra in the 1--10 GeV range. These turnovers are too gradual to arise from magnetic pair creation in the strong magnetic fields of pulsar inner magnetospheres. By demanding insignificant photon attenuation precipitated by such single-photon pair creation, the energies of these turnovers for Fermi pulsars can be used to compute lower bounds for the typical altitude of GeV band emission. This paper explores such pair transparency constraints below the turnover energy, and updates earlier altitude bound determinations of that have been deployed in various Fermi pulsar papers. For low altitude emission locales, general relativistic influences are found to be important. Rotational aberration influences are also explored, and are found to be small at low altitudes, except near the magnetic pole. The altitude bounds are typically in the range of 2-7 stellar radii for the young Fermi pulsar population, and provide key information on the emission altitude in radio quiet pulsars that do not possess double-peaked pulse profiles. The bound for the Crab pulsar is at a much higher altitude, with the putative detection by MAGIC out to 350-400 GeV implying a lower bound of 310km to the emission region, i.e., approximately 20% of the light cylinder radius. These results are also extended to the super-critical field domain, where it is found that emission in magnetars originating below around 10 stellar radii will not appear in the Fermi-LAT band.

preprint2014arXiv

Petawatt laser absorption bounded

The interaction of petawatt ($10^{15}\ \mathrm{W}$) lasers with solid matter forms the basis for advanced scientific applications such as table-top particle accelerators, ultrafast imaging systems and laser fusion. Key metrics for these applications relate to absorption, yet conditions in this regime are so nonlinear that it is often impossible to know the fraction of absorbed light $f$, and even the range of $f$ is unknown. Here using a relativistic Rankine-Hugoniot-like analysis, we show for the first time that $f$ exhibits a theoretical maximum and minimum. These bounds constrain nonlinear absorption mechanisms across the petawatt regime, forbidding high absorption values at low laser power and low absorption values at high laser power. For applications needing to circumvent the absorption bounds, these results will accelerate a shift from solid targets, towards structured and multilayer targets, and lead the development of new materials.

preprint2014arXiv

Quasi-Periodic Oscillations in Short Recurring Bursts of the Soft-Gamma Repeater J1550-5418

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. The scarcity of giant flares makes a search for QPOs in the shorter, far more numerous bursts from Soft Gamma Repeaters (SGRs) desirable. In Huppenkothen et al (2013), we developed a Bayesian method for searching for QPOs in short magnetar bursts, taking into account the effects of the complicated burst structure, and have shown its feasibility on a small sample of bursts. Here, we apply the same method to a much larger sample from a burst storm of 286 bursts from SGR J1550-5418. We report a candidate signal at 260 Hz in a search of the individual bursts, which is fairly broad. We also find two QPOs at 93 Hz and one at 127 Hz, when averaging periodograms from a number of bursts in individual triggers, at frequencies close to QPOs previously observed in magnetar giant flares. Finally, for the first time, we explore the overall burst variability in the sample, and report a weak anti-correlation between the power-law index of the broadband model characterising aperiodic burst variability, and the burst duration: shorter bursts have steeper power law indices than longer bursts. This indicates that longer bursts vary over a broader range of time scales, and are not simply longer versions of the short bursts.

preprint2013arXiv

Multiwavelength Probes of the Environs of Relativistic Shocks in Blazar Jets

Diffusive shock acceleration (DSA) at relativistic shocks is likely to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud AGN. An important recent development for blazar science is the ability of Fermi-LAT data to pin down the power-law index of the high energy portion of emission in these sources, and therefore also the index of the underlying non-thermal particle population. This diagnostic potential was not possible prior to Fermi launch, when gamma-ray information was dominated by the highly-absorbed TeV band. This paper highlights how multiwavelength spectra including X-ray band and Fermi data can be used to probe diffusive acceleration in relativistic, oblique, MHD shocks in blazar jets. The spectral index of the non-thermal particle distributions resulting from Monte Carlo simulations of DSA, and the fraction of thermal particles accelerated to non-thermal energies, depend sensitively on the particles' mean free path scale, and also on the magnetic field obliquity to the shock normal. We investigate self-consistently the radiative synchrotron/Compton signatures of the resulting thermal and non-thermal particle distributions. Important constraints on the frequency of particle scattering and the level of field turbulence are identified for the blazar AO 0235+164. The possible interpretation that turbulence levels decline with remoteness from jet shocks, and a significant role for non-gyroresonant diffusion, are discussed.

preprint2012arXiv

Broadband Spectral Investigations of SGR J1550-5418 Bursts

We present the results of our broadband spectral analysis of 42 SGR J1550-5418 bursts simultaneously detected with the Swift/X-ray Telescope (XRT) and the Fermi/Gamma-ray Burst Monitor (GBM), during the 2009 January active episode of the source. The unique spectral and temporal capabilities of the XRT Windowed Timing mode have allowed us to extend the GBM spectral coverage for these events down to the X-ray domain (0.5-10 keV). Our earlier analysis of the GBM data found that the SGR J1550-5418 burst spectra were described equally well with a Comptonized model or with two blackbody functions; the two models were statistically indistinguishable. Our new broadband (0.5 - 200 keV) spectral fits show that, on average, the burst spectra are better described with two blackbody functions than with the Comptonized model. Thus, our joint XRT/GBM analysis clearly shows for the first time that the SGR J1550-5418 burst spectra might naturally be expected to exhibit a more truly thermalized character, such as a two-blackbody or even a multi-blackbody signal. Using the Swift and RXTE timing ephemeris for SGR J1550-5418 we construct the distribution of the XRT burst counts with spin phase and find that it is not correlated with the persistent X-ray emission pulse phase from SGR J1550-5418. These results indicate that the burst emitting sites on the neutron star need not be co-located with hot spots emitting the bulk of the persistent X-ray emission. Finally, we show that there is a significant pulse phase dependence of the XRT burst counts, likely demonstrating that the surface magnetic field of SGR J1550-5418 is not uniform over the emission zone, since it is anticipated that regions with stronger surface magnetic field could trigger bursts more efficiently.

preprint2012arXiv

Detection of spectral evolution in the bursts emitted during the 2008-2009 active episode of SGR J1550 - 5418

In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E 1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations \sim 170 ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to < 8 deg. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.

preprint2012arXiv

Resonant Compton Upscattering in High Field Neutron Stars

The extremely efficient process of resonant Compton upscattering by relativistic electrons in high magnetic fields is believed to be a leading emission mechanism of high field pulsars and magnetars in the production of intense X-ray radiation. New analytic developments for the Compton scattering cross section using Sokolov & Ternov (S&T) states with spin-dependent resonant widths are presented. These new results display significant numerical departures from both the traditional cross section using spin-averaged widths, and also from the spin-dependent cross section that employs the Johnson & Lippmann (J&L) basis states, thereby motivating the astrophysical deployment of this updated resonant Compton formulation. Useful approximate analytic forms for the cross section in the cyclotron resonance are developed for S&T basis states. These calculations are applied to an inner magnetospheric model of the hard X-ray spectral tails in magnetars, recently detected by RXTE and INTEGRAL. Relativistic electrons cool rapidly near the stellar surface in the presence of intense baths of thermal X-ray photons. We present resonant Compton cooling rates for electrons, and the resulting photon spectra at various magnetospheric locales, for magnetic fields above the quantum critical value. These demonstrate how this scattering mechanism has the potential to produce the characteristically flat spectral tails observed in magnetars.

preprint2012arXiv

The Black Hole Evolution and Space Time (BEST) Observatory

In this white paper, we discuss the concept of a next-generation X-ray mission called BEST (Black hole Evolution and Space Time). The mission concept uses a 3000 square centimeter effective area mirror (at 6 keV) to achieve unprecedented sensitivities for hard X-ray imaging spectrometry (5-70 keV) and for broadband X-ray polarimetry (2-70 keV). BEST can make substantial contributions to our understanding of the inner workings of accreting black holes, our knowledge about the fabric of extremely curved spacetime, and the evolution of supermassive black holes. BEST will allow for time resolved studies of accretion disks. With a more than seven times larger mirror area and a seven times wider bandpass than GEMS, BEST will take X-ray polarimetry to a new level: it will probe the time variability of the X-ray polarization from stellar mass and supermassive black holes, and it will measure the polarization properties in 30 independent energy bins. These capabilities will allow BEST to conduct tests of accretion disk models and the underlying spacetimes. With three times larger mirror area and ten times better angular resolution than NuSTAR, BEST will be able to make deep field observations with a more than 15 times better sensitivity than NuSTAR. The mission will be able to trace the evolution of obscured and unobscured black holes in the redshift range from zero to six, covering the most important epoch of supermassive black hole growth. The hard X-ray sensitivity of BEST will enable a deep census of non-thermal particle populations. BEST will give us insights into AGN feedback by measuring the particle luminosity injected by AGNs into the interstellar medium (ISM) of their hosts, and will map the emission from particles accelerated at large scale structure shocks. Finally, BEST has the potential to constrain the equation of state of neutron stars (NS).

preprint2011arXiv

Burst and Persistent Emission Properties during the Recent Active Episode of the Anomalous X-ray Pulsar 1E 1841-045

Swift/BAT detected the first burst from 1E 1841-045 in May 2010 with intermittent burst activity recorded through at least July 2011. Here we present Swift and Fermi/GBM observations of this burst activity and search for correlated changes to the persistent X-ray emission of the source. The T90 durations of the bursts range between 18-140 ms, comparable to other magnetar burst durations, while the energy released in each burst ranges between (0.8 - 25)E38 erg, which is in the low side of SGR bursts. We find that the bursting activity did not have a significant effect on the persistent flux level of the source. We argue that the mechanism leading to this sporadic burst activity in 1E 1841-045 might not involve large scale restructuring (either crustal or magnetospheric) as seen in other magnetar sources.

preprint2011arXiv

Cooling Rates for Relativistic Electrons Undergoing Compton Scattering in Strong Magnetic Fields

For inner magnetospheric models of hard X-ray and gamma-ray emission in high-field pulsars and magnetars, resonant Compton upscattering is anticipated to be the most efficient process for generating continuum radiation. This is due in part to the proximity of a hot soft photon bath from the stellar surface to putative radiation dissipation regions in the inner magnetosphere. Moreover, because the scattering process becomes resonant at the cyclotron frequency, the effective cross section exceeds the classical Thomson value by over two orders of magnitude, thereby enhancing the efficiency of continuum production and the cooling of relativistic electrons. This paper presents computations of the electron cooling rates for this process, which are needed for resonant Compton models of non-thermal radiation from such highly-magnetized pulsars. The computed rates extend previous calculations of magnetic Thomson cooling to the domain of relativistic quantum effects, sampled near and above the quantum critical magnetic field of 44.13 TeraGauss. This is the first exposition of fully relativistic, quantum magnetic Compton cooling rates for electrons, and it employs both the traditional Johnson and Lippman cross section, and a newer Sokolov and Ternov (ST) formulation of Compton scattering in strong magnetic fields. Such ST formalism is formally correct for treating spin-dependent effects that are important in the cyclotron resonance, and has not been addressed before in the context of cooling by Compton scattering. The QED effects are observed to profoundly lower the rates below extrapolations of the familiar magnetic Thomson results, as expected, when recoil and Klein-Nishina reductions become important.

preprint2011arXiv

Diffusive Acceleration of Particles at Oblique, Relativistic, Magnetohydrodynamic Shocks

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 ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are germane to 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-LAT 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.

preprint2011arXiv

Fermi/GBM Observations of SGR J0501+4516 Bursts

We present our temporal and spectral analyses of 29 bursts from SGR J0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T90 durations of the bursts can be fit with a log-normal distribution with a mean value of ~ 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T90s estimated in count space (following a log-normal distribution with a mean value of ~ 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that Epeak decreases with energy flux (and fluence) to a minimum of ~30 keV at F=8.7e-6 erg/cm2/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550-5418 and 1806-20. The isotropic luminosity corresponding to these flux values is roughly similar for all sources (0.4-1.5 e40 erg/s).

preprint2010arXiv

A Suzaku X-Ray Study of the Particle Acceleration Processes in the Relativistic Jet of Blazar Mrk 421

We report on the findings of a 364 ksec observation of the BL LAC object Mrk 421 with the X-ray observatory Suzaku. The analysis in this paper uses fluxes and hardness ratios in the broad energy range from 0.5 keV to 30 keV. During the course of the observation, the 0.5 keV - 30 keV flux decreased by a factor of $\sim$2 and was accompanied by several large flares occurring on timescales of a few hours. We find that fitting a broken power model to spectra from isolated epochs during the observation describes the data well. Different flares exhibit different spectral and hardness ratio evolutions. The cumulative observational evidence indicates that the particle acceleration mechanism in the Mrk 421 jet produces electron energy distributions with a modest range of spectral indices and maximum energies. We argue that the short-timescale X-ray spectral variability in the flares can be attributed mostly to intrinsic changes in the acceleration process, dominating other influences such as fluctuations in the Doppler beaming factor, or radiative cooling in or outside the acceleration zone.

preprint2010arXiv

Perspectives on Gamma-Ray Pulsar Emission

Pulsars are powerful sources of radiation across the electromagnetic spectrum. This paper highlights some theoretical insights into non-thermal, magnetospheric pulsar gamma-ray radiation. These advances have been driven by NASA's Fermi mission, launched in mid-2008. The Large Area Telescope (LAT) instrument on Fermi has afforded the discrimination between polar cap and slot gap/outer gap acceleration zones in young and middle-aged pulsars. Altitude discernment using the highest energy pulsar photons will be addressed, as will spectroscopic interpretation of the primary radiation mechanism in the LAT band, connecting to both polar cap/slot gap and outer gap scenarios. Focuses will mostly be on curvature radiation and magnetic pair creation, including population trends that may afford probes of the magnetospheric accelerating potential.

preprint2010arXiv

Photospheric radius expansion during magnetar bursts

On August 24th 2008 the new magnetar SGR 0501+4516 (discovered by SWIFT) emitted a bright burst with a pronounced double-peak structure in hard X-rays, reminiscent of the double-peak temporal structure seen in some bright thermonuclear bursts on accreting neutron stars. In the latter case this is due to Photospheric Radius Expansion (PRE): when the flux reaches the Eddington limit, the photosphere expands and cools so that emission becomes softer and drops temporarily out of the X-ray band, re-appearing as the photosphere settles back down. We consider the factors necessary to generate double-peaked PRE events, and show that such a mechanism could plausibly operate in magnetar bursts, despite the vastly different emission process. Identification of the magnetic Eddington limit in a magnetar would constrain magnetic field and distance and could, in principle, enable a measurement of gravitational redshift. It would also locate the emitting region at the neutron star surface, constraining the burst trigger mechanism. Conclusive confirmation of PRE events will require more detailed radiative models for bursts. However for SGR 0501+4516 the predicted critical flux (using the magnetic field strength inferred from timing and the distance suggested by its probable location in the Perseus arm of our Galaxy) is consistent with that observed in the August 24th burst.

preprint2010arXiv

Using Gamma-Ray Burst Prompt Emission to Probe Relativistic Shock Acceleration

It is widely accepted that the prompt transient signal in the 10 keV - 10 GeV band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the ultra-relativistic expansion. The detailed understanding of the dissipation and accompanying acceleration at these shocks is a currently topical subject. This paper explores the relationship between GRB prompt emission spectra and the electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The focus is on the array of possible high-energy power-law indices in accelerated populations, highlighting how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. It is emphasized that diffusive shock acceleration theory generates no canonical spectrum at relativistic MHD discontinuities. This diversity is commensurate with the significant range of spectral indices discerned in prompt burst emission. Such system diagnostics are now being enhanced by the broadband spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst Monitor (GBM) provides key diagnostics on the lower energy portions of the particle population, the focus here is on constraints in the non-thermal, power-law regime of the particle distribution that are provided by the Large Area Telescope (LAT).

preprint2009arXiv

Particle Acceleration at Relativistic Shocks in Extragalactic Systems

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 strong and interesting candidate sites for the generation of ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are salient to the issue of cosmic ray ion and electron production are outlined. Results from a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, MHD shocks are presented. Simulation output is described for both large angle and small angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffman-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, providing evidence for its origin being due to shock drift acceleration. The impact of these theoretical results on gamma-ray burst and blazar science is outlined. Specifically, Fermi gamma-ray observations of these cosmic sources are already providing significant constraints on important environmental quantities for relativistic shocks, namely the frequency of scattering and the level of field turbulence.

preprint2008arXiv

Particle Acceleration at Interplanetary Shocks

The acceleration of interstellar pick-up ions as well as solar wind species has been observed at a multitude of interplanetary (IP) shocks by different spacecraft. This paper expands upon previous work modeling the phase space distributions of accelerated ions associated with the shock event encountered on day 292 of 1991 by the Ulysses mission at 4.5 AU. A kinetic Monte Carlo simulation is employed here to model the diffusive acceleration process. This exposition presents recent developments pertaining to the incorporation into the simulation of the diffusive characteristics incurred by field line wandering (FLW), according to the work of Giacalone and Jokipii. For a pure field-line wandering construct, it is determined that the upstream spatial ramp scales are too short to accommodate the HI-SCALE flux increases for 200 keV protons, and that the distribution function for H+ somewhat underpopulates the combined SWICS/HI-SCALE spectra at the shock. This contrasts our earlier theory/data comparison where it was demonstrated that diffusive transport in highly turbulent fields according to kinetic theory can successfully account for both the proton distributions and upstream ramp scales, using a single turbulence parameter. The principal conclusion here is that, in a FLW scenario, the transport of ions across the mean magnetic field is slightly less efficient than is required to effectively trap energetic ions within a few Larmor radii of the shock layer and thereby precipitate efficient acceleration. This highlights the contrast between ion transport in highly turbulent shock environs and remote, less-disturbed interplanetary regions.

preprint2008arXiv

Search for Redshifted 2.2 MeV Neutron Capture Line From A0535+262 in Outburst

The Be/X-ray binary system A0535+262 underwent a giant outburst in May-June 2005, followed by a dimmer outburst in August-September 2005. This increased intensity provided an opportunity to search for redshifted neutron-capture lines from the surface of the neutron star. If discovered, such lines would constrain the neutron star equation of state, providing the motivation of this search. The spectrometer (SPI) on board the INTEGRAL satellite observed the dimmer outburst and provided the data for this research. We have not detected a line with enough significance, with the width-dependent upper limits on the broadened and redshifted neutron capture line in the range of (2 - 11) x 10^(-4) photons cm^(-2) s^(-1). To our knowledge, these are the strongest upper limits on the redshifted 2.2 MeV emission from an accreting neutron star. Our analysis of the transparency of the neutron star surface for 2.2 MeV photons shows that photons have a small but finite chance of leaving the atmosphere unscattered, which diminishes the possibility of detection.

preprint2006arXiv

Temporal Evolution of Pair Attenuation Signatures in Gamma-Ray Burst Spectra

The spectra obtained above 100 MeV by the EGRET experiment aboard the Compton Gamma-Ray Observatory for a handful of gamma-ray bursts has given no indication of any spectral attenuation that might preclude detection of bursts at higher energies. With the discovery of optical afterglows and counterparts to bursts in the last few years, enabling the determination of significant redshifts for these sources, it is anticipated that profound spectral attenuation will arise in the Gamma-Ray Large Area Space Telescope (GLAST) energy band of 30 MeV-300 GeV for many if not most bursts. This paper explores time-dependent expectations for burst spectral properties in the EGRET/GLAST band, focusing on how attenuation of photons by pair creation internal to the source generates distinctive spectral signatures. The energy of spectral breaks and the associated spectral indices provide valuable information that constrains the bulk Lorentz factor of the GRB outflow at a given time. Moreover, the distinct temporal behavior that is present for internal attenuation is easily distinguished from extrinsic absorption due to intervening cosmic background fields. These characteristics define palpable observational goals for both spaced-based hard gamma-ray experiments such as GLAST, and ground-based Cherenkov telescopes, and strongly impact the observability of bursts above 300 MeV.

preprint2005arXiv

Modeling Accelerated Pick-up Ion Distributions at an Interplanetary Shock

The acceleration of interstellar pick-up ions as well as solar wind species has been observed at a multitude of interplanetary (IP) shocks by different spacecraft. The efficiency of injection of the pick-up ion component differs from that of the solar wind, and is expected to be strongly enhanced at highly oblique and quasi-perpendicular shock events, in accord with inferences from {\it in situ} observations. This paper explores theoretical modeling of the phase space distributions of accelerated ions obtained by the Ulysses mission for the Day 292, 1991 shock associated with a corotating interaction region, encountered before Ulysses' fly-by of Jupiter. A Monte Carlo simulation is used to model the acceleration process, adapting a technique that has been successfully tested on earlier IP shocks possessing minimal pick-up ion presence. Phase space distributions from the simulation technique for various low mass ions are compared with SWICS and HI-SCALE data to deduce values of a ``turbulence parameter'' that controls the efficiency of injection, and the degree of cross-field diffusion. Acceptable fits are obtained for the $H^+$ and $He^+$ populations using standard prescriptions for the pick-up ion distribution; $He^{++}$ spectral data was only fit well for scenarios very close to the Bohm diffusion limit. It is also found that the simulation successfully accounts for the observation of energetic protons farther upstream of the forward shock than lower energy pick-up protons, using the same turbulence parameter that is required to achieve reasonable spectral fits.

preprint1998arXiv

Direct Acceleration of Pickup Ions at the Solar Wind Termination Shock: the Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross-field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He$^{+}$/H$^{+}$ or O$^{+}$/H$^{+}$ to within a factor $\sim 5$. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

preprint1997arXiv

The Attenuation of Gamma-Ray Emission in Strongly-Magnetized Pulsars

Gamma-rays from pulsars can be efficiently attenuated in their magnetospheres via the mechanism of single-photon pair production and also the exotic QED process of photon splitting, which become prolific in fields approaching the quantum critical value of $B_{cr}=4.41\times 10^{13}$ Gauss. Recently we have published results of our modelling of strongly-magnetized $γ$-ray pulsars, which focused on the escape or attenuation of photons emitted near the pole at the neutron star surface in dipole fields, in a Schwarzschild metric. We found that pair production and splitting totally inhibit emission above around 10--30 MeV in PSR1509-58, whose surface field is inferred to be as high as $0.7B_{cr}$. Our model pulsar spectra are consistent with the EGRET upper limits for PSR1509-58 for a wide range of polar cap sizes. Here we review the principal predictions of our attenuation analysis, and identify how its powerful observational diagnostic capabilities relate to current and future gamma-ray experiments. Diagnostics include the energy of the gamma-ray turnover and the spectral polarization, which constrain the estimated polar cap size and field strength, and can determine the relative strength of splitting and pair creation.