Source author record

Donald C. Ellison

Donald C. Ellison 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

22works
4topics
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

22 published item(s)

preprint2015arXiv

Are Models for Core-Collapse Supernova Progenitors Consistent with the Properties of Supernova Remnants?

The recent discovery that the Fe-K line luminosities and energy centroids observed in nearby SNRs are a strong discriminant of both progenitor type and circumstellar environment has implications for our understanding of supernova progenitor evolution. Using models for the chemical composition of core-collapse supernova ejecta, we model the dynamics and thermal X-ray emission from shocked ejecta and circumstellar material, modeled as an $r^{-2}$ wind, to ages of 3000 years. We compare the X-ray spectra expected from these models to observations made with the Suzaku satellite. We also model the dynamics and X-ray emission from Type Ia progenitor models. We find a clear distinction in Fe-K line energy centroid between core-collapse and Type Ia models. The core-collapse supernova models predict higher Fe-K line centroid energies than the Type Ia models, in agreement with observations. We argue that the higher line centroids are a consequence of the increased densities found in the circumstellar environment created by the expansion of the slow-moving wind from the massive progenitors.

preprint2015arXiv

Electron and Ion Acceleration in Relativistic Shocks with Applications to GRB Afterglows

We have modeled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parameterizing the particle diffusion, our steady-state Monte Carlo simulation allows us to follow particles from particle injection at nonthermal thermal energies to above PeV energies, including the nonlinear smoothing of the shock structure due to cosmic-ray (CR) backpressure. We observe the mass-to-charge (A/Z) enhancement effect believed to occur in efficient Fermi acceleration in non-relativistic shocks and we parameterize the transfer of ion energy to electrons seen in particle-in-cell (PIC) simulations. For a given set of environmental and model parameters, the Monte Carlo simulation determines the absolute normalization of the particle distributions and the resulting synchrotron, inverse-Compton, and pion-decay emission in a largely self-consistent manner. The simulation is flexible and can be readily used with a wide range of parameters typical of gamma-ray burst (GRB) afterglows. We describe some preliminary results for photon emission from shocks of different Lorentz factors and outline how the Monte Carlo simulation can be generalized and coupled to hydrodynamic simulations of GRB blast waves. We assume Bohm diffusion for simplicity but emphasize that the nonlinear effects we describe stem mainly from an extended shock precursor where higher energy particles diffuse further upstream. Quantitative differences will occur with different diffusion models, particularly for the maximum CR energy and photon emission, but these nonlinear effects should be qualitatively similar as long as the scattering mean free path is an increasing function of momentum.

preprint2015arXiv

Modeling Bright Gamma-ray and Radio Emission at Fast Cloud Shocks

Recent observations by the Large Area Telescope (LAT) onboard the Fermi satellite have revealed bright gamma-ray emission from middle-aged supernova remnants (SNRs) inside our Galaxy. These remnants, which also possess bright non-thermal radio shells, are often found to be interacting directly with surrounding gas clouds. We explore the non-thermal emission mechanism at these dynamically evolved SNRs by constructing a hydrodynamical model. Two scenarios of particle acceleration, either a re-acceleration of Galactic cosmic rays (CRs) or an efficient nonlinear diffusive shock acceleration (NLDSA) of particles injected from downstream, are considered. Using parameters inferred from observations, our models are contrasted with the observed spectra of SNR W44. For the re-acceleration case, we predict a significant enhancement of radio and GeV emission as the SNR undergoes a transition into the radiative phase. If sufficiently strong magnetic turbulence is present in the molecular cloud, the re-acceleration scenario can explain the observed broadband spectral properties. The NLDSA scenario also succeeds in explaining the $γ$-ray spectrum but fails to reproduce the radio spectral index. Efficient NLDSA also results in a significant post-shock non-thermal pressure that limits the compression during cooling and prevents the formation of a prominent dense shell. Some other interesting differences between the two models in hydrodynamical behavior and resulting spectral features are illustrated.

preprint2015arXiv

Particle spectra and efficiency in nonlinear relativistic shock acceleration: survey of scattering models

We include a general form for the scattering mean free path in a nonlinear Monte Carlo model of relativistic shock formation and Fermi acceleration. Particle-in-cell (PIC) simulations, as well as analytic work, suggest that relativistic shocks tend to produce short-scale, self-generated magnetic turbulence that leads to a scattering mean free path (mfp) with a stronger momentum dependence than the mfp ~ p dependence for Bohm diffusion. In unmagnetized shocks, this turbulence is strong enough to dominate the background magnetic field so the shock can be treated as parallel regardless of the initial magnetic field orientation, making application to gamma-ray bursts (GRBs), pulsar winds, Type Ibc supernovae, and extra-galactic radio sources more straightforward and realistic. In addition to changing the scale of the shock precursor, we show that, when nonlinear effects from efficient Fermi acceleration are taken into account, the momentum dependence of the mfp has an important influence on the efficiency of cosmic-ray production as well as the accelerated particle spectral shape. These effects are absent in nonrelativistic shocks and do not appear in relativistic shock models unless nonlinear effects are self-consistently described. We show, for limited examples, how the changes in Fermi acceleration translate to changes in the intensity and spectral shape of gamma-ray emission from proton-proton interactions and pion-decay radiation.

preprint2014arXiv

Magnetic field amplification in nonlinear diffusive shock acceleration including resonant and non-resonant cosmic-ray driven instabilities

We present a nonlinear Monte Carlo model of efficient diffusive shock acceleration (DSA) where the magnetic turbulence responsible for particle diffusion is calculated self-consistently from the resonant cosmic-ray (CR) streaming instability, together with non-resonant short- and long-wavelength CR-current-driven instabilities. We include the backpressure from CRs interacting with the strongly amplified magnetic turbulence which decelerates and heats the super-alfvenic flow in the extended shock precursor. Uniquely, in our plane-parallel, steady-state, multi-scale model, the full range of particles, from thermal (~eV) injected at the viscous subshock, to the escape of the highest energy CRs (~PeV) from the shock precursor, are calculated consistently with the shock structure, precursor heating, magnetic field amplification (MFA), and scattering center drift relative to the background plasma. In addition, we show how the cascade of turbulence to shorter wavelengths influences the total shock compression, the downstream proton temperature, the magnetic fluctuation spectra, and accelerated particle spectra. A parameter survey is included where we vary shock parameters, the mode of magnetic turbulence generation, and turbulence cascading. From our survey results, we obtain scaling relations for the maximum particle momentum and amplified magnetic field as functions of shock speed, ambient density, and shock size.

preprint2014arXiv

Reverse and Forward Shock X-ray Emission in an Evolutionary Model of Supernova Remnants undergoing Efficient Diffusive Shock Acceleration

We present new models for the forward and reverse shock thermal X-ray emission from core-collapse and Type Ia supernova remnants (SNRs) which include the efficient production of cosmic rays via non-linear diffusive shock acceleration (DSA). Our CR-hydro-NEI code takes into account non-equilibrium ionization (NEI), hydrodynamic effects of efficient CR production on the SNR evolution, and collisional temperature equilibration among heavy ions and electrons in both the shocked supernova (SN) ejecta and the shocked circumstellar material. While X-ray emission is emphasized here, our code self-consistently determines both thermal and non-thermal broadband emission from radio to TeV energies. We include Doppler broadening of the spectral lines by thermal motions of the ions and by the remnant expansion. We study, in general terms, the roles which the ambient environment, progenitor models, temperature equilibration, and processes related to DSA have on the thermal and non-thermal spectra. The study of X-ray line emission from young SNRs is a powerful tool for determining specific SN elemental contributions, and for providing critical information that helps to understand the type and energetics of the explosion, the composition of the ambient medium in which the SN exploded, and the ionization and dynamics of the hot plasma in the shocked SN ejecta and interstellar medium. With the approaching launch of the next-generation X-ray satellite Astro-H, observations of spectral lines with unprecedented high resolution will become a reality. Our self-consistent calculations of the X-ray spectra from various progenitors will help interpret future observations of SNRs.

preprint2013arXiv

A CR-hydro-NEI Model of Multi-wavelength Emission from the Vela Jr. Supernova Remnant (SNR RX J0852.0-4622)

Based largely on energy budget considerations and the observed cosmic-ray (CR) ionic composition, supernova remnant (SNR) blast waves are the most likely sources of CR ions with energies at least up to the "knee" near 3 PeV. Shocks in young shell-type TeV-bright SNRs are surely producing TeV particles, but the emission could be dominated by ions producing neutral pion-decay emission or electrons producing inverse-Compton gamma-rays. Unambiguously identifying the GeV-TeV emission process in a particular SNR will not only help pin down the origin of CRs, it will add significantly to our understanding of the diffusive shock acceleration (DSA) mechanism and improve our understanding of supernovae and the impact SNRs have on the circumstellar medium. In this study, we investigate the Vela Jr. SNR, an example of TeV-bright non-thermal SNRs. We perform hydrodynamic simulations coupled with non-linear DSA and non-equilibrium ionization near the forward shock (FS) to confront currently available multi-wavelength data. We find, with an analysis similar to that used earlier for SNR RX J1713.7-3946, that self-consistently modeling the thermal X-ray line emission with the non-thermal continuum in our one-dimensional model strongly constrains the fitting parameters, and this leads convincingly to a leptonic origin for the GeV-TeV emission for Vela Jr. This conclusion is further supported by applying additional constraints from observation, including the radial brightness profiles of the SNR shell in TeV gamma-rays, and the spatial variation of the X-ray synchrotron spectral index. We will discuss implications of our models on future observations by the next-generation telescopes.

preprint2013arXiv

Monte Carlo Simulations of Nonlinear Particle Acceleration in Parallel Trans-relativistic Shocks

We present results from a Monte Carlo simulation of a parallel collisionless shock undergoing particle acceleration. Our simulation, which contains parameterized scattering and a particular thermal leakage injection model, calculates the feedback between accelerated particles ahead of the shock, which influence the shock precursor and "smooth" the shock, and thermal particle injection. We show that there is a transition between nonrelativistic shocks, where the acceleration efficiency can be extremely high and the nonlinear compression ratio can be substantially greater than the Rankine-Hugoniot value, and fully relativistic shocks, where diffusive shock acceleration is less efficient and the compression ratio remains at the Rankine-Hugoniot value. This transition occurs in the trans-relativistic regime and, for the particular parameters we use, occurs around a shock Lorentz factor = 1.5. We also find that nonlinear shock smoothing dramatically reduces the acceleration efficiency presumed to occur with large-angle scattering in ultra-relativistic shocks. Our ability to seamlessly treat the transition from ultra-relativistic to trans-relativistic to nonrelativistic shocks may be important for evolving relativistic systems, such as gamma-ray bursts and type Ibc supernovae. We expect a substantial evolution of shock accelerated spectra during this transition from soft early on to much harder when the blast-wave shock becomes nonrelativistic.

preprint2012arXiv

A Generalized Model of Nonlinear Diffusive Shock Acceleration Coupled to an Evolving Supernova Remnant

To better model the efficient production of cosmic rays (CRs) in supernova remnants (SNRs) with the associated coupling between CR production and SNR dynamics, we have generalized an existing cr-hydro-NEI code (i.e., Ellison et al. 2012) to include the following processes: (1) an explicit calculation of the upstream precursor structure including the position dependent flow speed, density, temperature, and magnetic field strength; (2) a momentum and space dependent CR diffusion coefficient; (3) an explicit calculation of magnetic field amplification (MFA); (4) calculation of the maximum CR momentum using the amplified magnetic field; (5) a finite Alfven speed for the particle scattering centers; and (6) the ability to accelerate a superthermal seed population of CRs as well as the ambient thermal plasma. While a great deal of work has been done modeling SNRs, most work has concentrated on either the continuum emission from relativistic electrons or ions, or the thermal emission from the shock heated plasma. Our generalized code combines these elements and describes the interplay between CR production and SNR evolution, including the nonlinear coupling of efficient diffusive shock acceleration (DSA), based mainly on the work of P. Blasi and co-workers, and a non-equilibrium ionization (NEI) calculation of thermal X-ray line emission. We believe our generalized model will provide a consistent modeling platform for SNRs, including those interacting with molecular clouds, and improve the interpretation of current and future observations, including the high-quality spectra expected from Astro-H. SNR RX J1713.7-3946 is modeled as an example.

preprint2012arXiv

Fermi-LAT Observations and A Broadband Study of Supernova Remnant CTB 109

CTB 109 (G109.1-1.0) is a Galactic supernova remnant (SNR) with a hemispherical shell morphology in X-rays and in the radio band. In this work we report the detection of γ-ray emission coincident with CTB 109, using 37 months of data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. We study the broadband characteristics of the remnant using a model that includes hydrodynamics, efficient cosmic ray acceleration, nonthermal emission and a self-consistent calculation of the X-ray thermal emission. We find that the observations can be successfully fit with two distinct parameter sets, one where the γ-ray emission is produced primarily by leptons accelerated at the SNR forward shock and the other where γ-rays produced by forward shock accelerated cosmic-ray ions dominate the high-energy emission. Consideration of thermal X-ray emission introduces a novel element to the broadband fitting process, and while it does not rule out either the leptonic or the hadronic scenarios, it constrains the parameter sets required by the model to fit the observations. Moreover, the model which best fits the thermal and nonthermal emission observations is an intermediate case, where both radiation from accelerated electrons and hadrons contribute almost equally to the γ-ray flux observed.

preprint2011arXiv

Core-collapse model of broadband emission from SNR RX J1713.7-3946 with thermal X-rays and Gamma-rays from escaping cosmic rays

We present a spherically symmetric, core-collapse model of SNR RX J1713.7-3946 that includes a hydrodynamic simulation of the remnant evolution coupled to the efficient production of cosmic rays (CRs) by nonlinear diffusive shock acceleration (DSA). High-energy CRs that escape from the forward shock (FS) are propagated in surrounding dense material that simulates either a swept-up, pre-supernova shell or a nearby molecular cloud. The continuum emission from trapped and escaping CRs, along with the thermal X-ray emission from the shocked heated ISM behind the FS, integrated over the remnant, is compared against broadband observations. Our results show conclusively that, overall, the GeV-TeV emission is dominated by inverse-Compton from CR electrons if the supernova is isolated regardless of its type, i.e., not interacting with a >>100 Msun shell or cloud. If the SNR is interacting with a much larger mass >10^4 Msun, pion-decay from the escaping CRs may dominate the TeV emission, although a precise fit at high energy will depend on the still uncertain details of how the highest energy CRs are accelerated by, and escape from, the FS. Based on morphological and other constraints, we consider the 10^4 Msun pion-decay scenario highly unlikely for SNR RX J1713.7-3946 regardless of the details of CR escape. Importantly, even though CR electrons dominate the GeV-TeV emission, the efficient production of CR ions is an essential part of our leptonic model.

preprint2011arXiv

Gamma-ray emission of accelerated particles escaping a supernova remnant in a molecular cloud

We present a model of gamma-ray emission from core-collapse supernovae originating from the explosions of massive young stars. The fast forward shock of the supernova remnant (SNR) can accelerate particles by diffusive shock acceleration (DSA) in a cavern blown by a strong, pre-supernova stellar wind. As a fundamental part of nonlinear DSA, some fraction of the accelerated particles escape the shock and interact with a surrounding massive dense shell producing hard photon emission. To calculate this emission, we have developed a new Monte Carlo technique for propagating the cosmic rays (CRs) produced by the forward shock of the SNR, into the dense, external material. This technique is incorporated in a hydrodynamic model of an evolving SNR which includes the nonlinear feedback of CRs on the SNR evolution, the production of escaping CRs along with those that remain trapped within the remnant, and the broad-band emission of radiation from trapped and escaping CRs. While our combined CR-hydro-escape model is quite general and applies to both core collapse and thermonuclear supernovae, the parameters we choose for our discussion here are more typical of SNRs from very massive stars whose emission spectra differ somewhat from those produced by lower mass progenitors directly interacting with a molecular cloud.

preprint2011arXiv

Magnetic fields in cosmic particle acceleration sources

We review here some magnetic phenomena in astrophysical particle accelerators associated with collisionless shocks in supernova remnants, radio galaxies and clusters of galaxies. A specific feature is that the accelerated particles can play an important role in magnetic field evolution in the objects. We discuss a number of CR-driven, magnetic field amplification processes that are likely to operate when diffusive shock acceleration (DSA) becomes efficient and nonlinear. The turbulent magnetic fields produced by these processes determine the maximum energies of accelerated particles and result in specific features in the observed photon radiation of the sources. Equally important, magnetic field amplification by the CR currents and pressure anisotropies may affect the shocked gas temperatures and compression, both in the shock precursor and in the downstream flow, if the shock is an efficient CR accelerator. Strong fluctuations of the magnetic field on scales above the radiation formation length in the shock vicinity result in intermittent structures observable in synchrotron emission images. Resonant and non-resonant CR streaming instabilities in the shock precursor can generate mesoscale magnetic fields with scale-sizes comparable to supernova remnants and even superbubbles. This opens the possibility that magnetic fields in the earliest galaxies were produced by the first generation Population III supernova remnants and by clustered supernovae in star forming regions.

preprint2011arXiv

The Impact of Efficient Particle Acceleration on the Evolution of Supernova Remnants in the Sedov-Taylor Phase

We investigate the effects of the efficient production of cosmic rays on the evolution of supernova remnants (SNRs) in the adiabatic Sedov-Taylor phase. We model the SNR by coupling the hydrodynamic evolution with nonlinear diffusive shock acceleration (DSA), and track self-consistently the ionization state of the shock-heated plasma. Using a plasma emissivity code and the results of the model, we predict the thermal X-ray emission and combine it with the non-thermal component in order to obtain the complete spectrum in this energy range. Hence, we study how the interpretation of thermal X-ray observations is affected by the efficiency of the DSA process, and find that, compared to test particle cases, the efficient DSA example yields a smaller shock radius and speed, a larger compression ratio, and lower intensity X-ray thermal emission. We also find that a model where the shock is not assumed to produce CRs can fit the X-ray observational properties of an example with efficient particle acceleration, with a different set of input parameters, and in particular a much lower explosion energy. Additionally, we model the broadband non-thermal emission, and investigate what signatures result from the acceleration of particles.

preprint2010arXiv

Efficient cosmic ray acceleration, hydrodynamics, and Self-consistent Thermal X-ray Emission applied to SNR RX J1713.7-3946

We model the broad-band emission from SNR RX J1713.7-3946 including, for the first time, a consistent calculation of thermal X-ray emission together with non-thermal emission in a nonlinear diffusive shock acceleration (DSA) model. Our model tracks the evolution of the SNR including the plasma ionization state between the forward shock and the contact discontinuity. We use a plasma emissivity code to predict the thermal X-ray emission spectrum assuming the initially cold electrons are heated either by Coulomb collisions with the shock heated protons (the slowest possible heating), or come into instant equilibration with the protons. For either electron heating model, electrons reach >10^7 K rapidly and the X-ray line emission near 1 keV is more than 10 times as luminous as the underlying thermal bremsstrahlung continuum. Since recent Suzaku observations show no detectable line emission, this places strong constraints on the unshocked ambient medium density and on the relativistic electron to proton ratio. For the uniform circumstellar medium (CSM) models we consider, the low densities and high relativistic electron to proton ratios required to match the Suzaku X-ray observations definitively rule out pion-decay as the emission process producing GeV-TeV photons. We show that leptonic models, where inverse-Compton scattering against the cosmic background radiation dominates the GeV-TeV emission, produce better fits to the broad-band thermal and non-thermal observations in a uniform CSM.

preprint2010arXiv

The Role of Diffusive Shock Acceleration on Nonequilibrium Ionization in Supernova Remnant Shocks II: Emitted Spectra

We present a grid of nonequilibrium ionization models for the X-ray spectra from supernova remnants undergoing efficient diffusive shock acceleration. The calculation follows the hydrodynamics of the blast wave as well as the time-dependent ionization of the plasma behind the shock. The ionization state is passed to a plasma emissivity code to compute the thermal X-ray emission, which is combined with the emission from nonthermal synchrotron emission to produce a self-consistent model for the thermal and nonthermal emission from cosmic-ray dominated shocks. We show how plasma diagnostics such as the G'-ratio of He-like ions, defined as the ratio of the sum of the intercombination, forbidden, and satellite lines to the resonance line, can vary with acceleration efficiency, and discuss how the thermal X-ray emission, when the time-dependent ionization is not calculated self-consistently with the hydrodynamics, can differ from the thermal X-ray emission from models which do account for the hydrodynamics. Finally we compare the thermal X-ray emission from models which show moderate acceleration (~ 35%) to the thermal X-ray emission from test-particle models.

preprint2009arXiv

Spectra of magnetic fluctuations and relativistic particles produced by a nonresonant wave instability in supernova remnant shocks

We model strong forward shocks in young supernova remnants with efficient particle acceleration where a nonresonant instability driven by the cosmic ray current amplifies magnetic turbulence in the shock precursor. Particle injection, magnetic field amplification (MFA) and the nonlinear feedback of particles and fields on the bulk flow are derived consistently. The shock structure depends critically on the efficiency of turbulence cascading. If cascading is suppressed, MFA is strong, the shock precursor is stratified, and the turbulence spectrum contains several discrete peaks. These peaks, as well as the amount of MFA, should influence synchrotron X-rays, allowing observational tests of cascading and other assumptions intrinsic to the nonlinear model of nonresonant wave growth.

preprint2009arXiv

The Role of Diffusive Shock Acceleration on Non-Equilibrium Ionization in Supernova Remnants

We present results of semi-analytic calculations which show clear evidence for changes in the non-equilibrium ionization behind a supernova remnant forward shock undergoing efficient diffusive shock acceleration (DSA). The efficient acceleration of particles (i.e., cosmic rays) lowers the shock temperature and raises the density of the shocked gas, thus altering the ionization state of the plasma in comparison to the test particle approximation where cosmic rays gain an insignificant fraction of the shock energy. The differences between the test particle and efficient acceleration cases are substantial and occur for both slow and fast temperature equilibration rates: in cases of higher acceleration efficiency, particular ion states are more populated at lower electron temperatures. We also present results which show that, in the efficient shock acceleration case, higher ionization fractions are reached noticeably closer to the shock front than in the test-particle case, clearly indicating that DSA may enhance thermal X-ray production. We attribute this to the higher postshock densities which lead to faster electron temperature equilibration and higher ionization rates. These spatial differences should be resolvable with current and future X-ray missions, and can be used as diagnostics in estimating the acceleration efficiency in cosmic-ray modified shocks.

preprint2008arXiv

Dots, clumps and filaments: the intermittent images of synchrotron emission in random magnetic fields of young supernova remnants

Non-thermal X-ray emission in some supernova remnants originates from synchrotron radiation of ultra-relativistic particles in turbulent magnetic fields. We address the effect of a random magnetic field on synchrotron emission images and spectra. A random magnetic field is simulated to construct synchrotron emission maps of a source with a steady distribution of ultra-relativistic electrons. Non-steady localized structures (dots, clumps and filaments), in which the magnetic field reaches exceptionally high values, typically arise in the random field sample. These magnetic field concentrations dominate the synchrotron emission (integrated along the line of sight) from the highest energy electrons in the cut-off regime of the distribution, resulting in an evolving, intermittent, clumpy appearance. The simulated structures resemble those observed in X-ray images of some young supernova remnants. The lifetime of X-ray clumps can be short enough to be consistent with that observed even in the case of a steady particle distribution. The efficiency of synchrotron radiation from the cut-off regime in the electron spectrum is strongly enhanced in a turbulent field compared to emission from a uniform field of the same magnitude.

preprint2007arXiv

Particle Acceleration in Supernova Remnants and the Production of Thermal and Nonthermal Radiation

If highly efficient, cosmic ray production can have a significant effect on the X-ray emission from SNRs as well as their dynamical evolution. Using hydrodynamical simulations including diffusive shock acceleration, we produce spectra for both the thermal and nonthermal forward shock emission. For a given ambient density and explosion energy, we find that the position of the forward shock at a given age is a strong function of the acceleration efficiency, providing a signature of cosmic-ray production. Using an approximate treatment for the ionization state of the plasma, we investigate the effects of slow vs. rapid heating of the postshock electrons on the ratio of thermal to nonthermal X-ray emission at the forward shock. We also investigate the effects of magnetic field strength on the observed spectrum for efficient cosmic-ray acceleration. The primary effect of a large field is a considerable flattening of the nonthermal spectrum in the soft X-ray band. Spectral index measurements from X-ray observations may thus be indicators of the postshock magnetic field strength. The predicted gamma-ray flux from inverse-Compton (IC) scattering and neutral pion decay is strongly affected by the ambient conditions and, for the particular parameters used in our examples, the IC emission at E ~ 1 TeV exceeds that from pion decay, although at both lower and higher energies this trend is reversed for cases of high ambient density. More importantly, high magnetic fields produce a steepening of the electron spectrum over a wide energy range which may make it more difficult to differentiate between IC and pion-decay emission solely by spectral shape.

preprint2002arXiv

Nonlinear Particle Acceleration in Relativistic Shocks

Monte Carlo techniques are used to model nonlinear particle acceleration in parallel collisionless shocks of various speeds, including mildly relativistic ones. When the acceleration is efficient, the backreaction of accelerated particles modifies the shock structure and causes the compression ratio, r, to increase above test-particle values. Modified shocks with Lorentz factors less than about 3 can have compression ratios considerably greater than 3 and the momentum distribution of energetic particles no longer follows a power law relation. These results may be important for the interpretation of gamma-ray bursts if mildly relativistic internal and/or afterglow shocks play an important role accelerating particles that produce the observed radiation. For shock Lorentz factors greater than about 10, r approaches 3 and the so-called `universal' test-particle result of N(E) proportional to E^{-2.3} is obtained for sufficiently energetic particles. In all cases, the absolute normalization of the particle distribution follows directly from our model assumptions and is explicitly determined.

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.