For most elements, the isotopic ratios seen in cosmic rays (CRs) are similar to those in the solar wind. The most important exception to this is $^{22}$Ne/$^{20}$Ne where the CR value is $\sim 5$ times that of the solar wind. According to most recent models of nucleosynthesis, a large amount of $^{22}\mathrm{Ne}$ is generated in Wolf-Rayet (WR) stars. In the winds of carbon sequence WR stars, i.e., WC stars, the isotopic ratio $^{22}\mathrm{Ne}$/$^{20}\mathrm{Ne}$ can be much larger than in the solar wind. Here, we consider CRs produced by $^{22}\mathrm{Ne}$-enriched WR winds in compact clusters of young massive stars like Westerlund~1. We assume that efficient CR acceleration in clusters occurs over the WR lifetime in an ensemble of shock waves from multiple massive star winds. We estimate the fraction of all Galactic CRs such sources may produce for a given set of parameters.

Northern line-of-sight extinction within Gale Crater during the 2018 global dust storm was monitored daily using MSL's Navcam. Additional observations with Mastcam (north) and Navcam (all directions) were obtained at a lower cadence. Using feature identification and geo-referencing, extinction was estimated in all possible directions. Peak extinction of $>1.1$ km$^{-1}$ was measured between sols 2086 and 2090, an order of magnitude higher than previous maxima. Northern and western directions show an initial decrease, followed by a secondary peak in extinction, not seen in column opacity measurements. Due to foreground topography, eastern direction results are provided only as limits, and southern results were indeterminable. Mastcam red and green filter results agree well, but blue filter results show higher extinctions, likely due to low signal-to-noise. Morning results are systematically higher than afternoon results, potentially indicative of atmospheric mixing.

It has been shown that supernova blast waves interacting with winds from massive stars in compact star clusters may be capable of producing cosmic-ray (CR) protons to above $10^{17}$ eV. We give a brief description of the colliding-shock-flows mechanism and look at generalizations of the diffusion of ~ 100 PeV CRs in the turbulent galactic magnetic field present in the galactic disk. We calculate the temporal evolution of the CR anisotropy from a possible distribution of young compact massive star clusters assuming the sources are intermittent on time scales of a few million years, i.e., comparable to their residence time in the Milky Way. Within the confines of our model, we determine the galactic/extra-galactic fraction of high-energy CRs resulting in anisotropies consistent with observed values. We find that galactic star clusters may contribute a substantial fraction of ~ 100 PeV CRs without producing anisotropies above observed limits.

Relativistic particle acceleration in collisionless shocks of supernova remnants is accompanied by magnetic field amplification from cosmic ray (CR) driven plasma instabilities. Bell's fast CR-current instability is predicted to produce turbulence with a non-zero mean electric field in the shock precursor. We present a Monte Carlo model of Fermi shock acceleration explicitly taking into account an effective mean upstream electric field. Our model is nonlinear and includes the backreaction effects of efficient Fermi acceleration on the shock structure.

This study constrains the lower bound of the scattering phase function of Martian water ice clouds (WICs) through the implementation of a new observation aboard the Mars Science Laboratory (MSL). The Phase Function Sky Survey (PFSS) was a multiple pointing all-sky observation taken with the navigation cameras (Navcam) aboard MSL. The PFSS was executed 35 times during the Aphelion Cloud Belt (ACB) season of Mars Year 34 over a solar longitude range of L_s=61.4{\deg}-156.5{\deg}. Twenty observations occurred in the morning hours between 06:00 and 09:30 LTST, and 15 runs occurred in the evening hours between 14:30 and 18:00 LTST, with an operationally required 2.5 hour gap on either side of local noon due the sun being located near zenith. The resultant WIC phase function was derived over an observed scattering angle range of 18.3{\deg} to 152.61{\deg}, normalized, and compared with 9 modeled phase functions: seven ice crystal habits and two Martian WIC phase functions currently being implemented in models. Through statistical chi-squared probability tests, the five most probable ice crystal geometries observed in the ACB WICs were aggregates, hexagonal solid columns, hollow columns, plates, and bullet rosettes with p-values greater than or equal to 0.60, 0.57,0.56,0.56, and 0.55, respectively. Droxtals and spheres had p-values of 0.35, and 0.2, making them less probable components of Martian WICs, but still statistically possible ones. Having a better understanding of the ice crystal habit and phase function of Martian water ice clouds directly benefits Martian climate models which currently assume spherical and cylindrical particles.

We give a brief review of the origin and acceleration of cosmic rays (CRs), emphasizing the production of CRs at different stages of supernova evolution by the first-order Fermi shock acceleration mechanism. We suggest that supernovae with trans-relativistic outflows, despite being rather rare, may accelerate CRs to energies above 10$^{18}$ eV over the first year of their evolution. Supernovae in young compact clusters of massive stars, and interaction powered superluminous supernovae, may accelerate CRs well above the PeV regime. We discuss the acceleration of the bulk of the galactic CRs in isolated supernova remnants and re-acceleration of escaped CRs by the multiple shocks present in superbubbles produced by associations of OB stars. The effects of magnetic field amplification by CR driven instabilities, as well as superdiffusive CR transport, are discussed for nonthermal radiation produced by nonlinear shocks of all speeds including trans-relativistic ones.

There is now substantial evidence that the progenitors of some core-collapse supernovae undergo enhanced or extreme mass loss prior to explosion. The imprint of this mass loss is observed in the spectra and dynamics of the expanding blastwave on timescales of days to years after core-collapse, and the effects on the spectral and dynamical evolution may linger long after the supernova has evolved into the remnant stage. In this paper, we present for the first time, largely self-consistent end-to-end simulations for the evolution of a massive star from the pre-main sequence, up to and through core collapse, and into the remnant phase. We present three models and compare and contrast how the progenitor mass loss history impacts the dynamics and spectral evolution of the supernovae and supernova remnants. We study a model which only includes steady mass loss, a model with enhanced mass loss over a period of $\sim$ 5000 years prior to core-collapse, and a model with extreme mass loss over a period of $\sim$ 500 years prior to core collapse. The models are not meant to address any particular supernova or supernova remnant, but rather to highlight the important role that the progenitor evolution plays in the observable qualities of supernovae and supernova remnants. Through comparisons of these three different progenitor evolution scenarios, we find that the mass loss in late stages (during and after core carbon burning) can have a profound impact on the dynamics and spectral evolution of the supernova remnant centuries after core-collapse.

We discuss a specific population of galactic PeVatrons which may be the main source of the galactic cosmic-ray (CR) component well above PeV energies. Supernovae in compact clusters of massive stars are proposed as powerful sources of CRs, neutrinos, and gamma-ray emission. Numerical simulations of non-linear Fermi acceleration at converging shock flows have revealed that these accelerators can provide very hard spectra of protons up to $10^{16}-10^{17}$ eV which is well above the "knee" in the all-particle CR spectrum at about $3\times10^{15}$ eV. We suggest that known supernova remnants interacting with stellar winds in the compact clusters of young massive stars Westerlund I and Cl*1806-20 can be associated with the sources of the TeV gamma-ray emission detected by H.E.S.S. and may be responsible for a fraction of the high-energy neutrinos detected with the IceCube observatory. A recent CR composition measurement with the LOFAR array has revealed a light-mass component possibly dominating the all-particle spectrum at energies around $10^{17}$ eV. Such a strong light component (mainly protons and helium) may require specific galactic CR sources such as supernovae interacting with compact clusters of massive stars in addition to isolated supernova remnants.

Relativistic astrophysical collisionless shocks represent outstanding dissipation agents of the huge power of relativistic outflows produced by accreting black holes, core collapsed supernovae and other objects into multi-messenger radiation (cosmic rays, neutrinos, electromagnetic radiation). This article provides a theoretical discussion of the fundamental physical ingredients of these extreme phenomena. In the context of weakly magnetized shocks, in particular, it is shown how the filamentation type instabilities, which develop in the precursor of pair dominated or electron-ion shocks, provide the seeds for the scattering of high energy particles as well as the agent which preheats and slows down the incoming precursor plasma. This analytical discussion is completed with a mesoscopic, non-linear model of particle acceleration in relativistic shocks based on Monte Carlo techniques. This Monte Carlo model uses a semi-phenomenological description of particle scattering which allows it to calculate the back-reaction of accelerated particles on the shock structure on length and momentum scales which are currently beyond the range of microscopic particle-in-cell (PIC) simulations.

Fast collisionless shocks in cosmic plasmas convert their kinetic energy flow into the hot downstream thermal plasma with a substantial fraction of energy going into a broad spectrum of superthermal charged particles and magnetic fluctuations. The superthermal particles can penetrate into the shock upstream region producing an extended shock precursor. The cold upstream plasma flow is decelerated by the force provided by the superthermal particle pressure gradient. In high Mach number collisionless shocks, efficient particle acceleration is likely coupled with turbulent magnetic field amplification (MFA) generated by the anisotropic distribution of accelerated particles. This anisotropy is determined by the fast particle transport making the problem strongly nonlinear and multi-scale. Here, we present a nonlinear Monte Carlo model of collisionless shock structure with super-diffusive propagation of high-energy Fermi accelerated particles coupled to particle acceleration and MFA which affords a consistent description of strong shocks. A distinctive feature of the Monte Carlo technique is that it includes the full angular anisotropy of the particle distribution at all precursor positions. The model reveals that the super-diffusive transport of energetic particles (i.e., Levy-walk propagation) generates a strong quadruple anisotropy in the precursor particle distribution. The resultant pressure anisotropy of the high-energy particles produces a non-resonant mirror-type instability which amplifies compressible wave modes with wavelengths longer than the gyroradii of the highest energy protons produced by the shock.

The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, $N(E) \propto E^{-p}$ for $p \gtrsim 2$. However, in standard Fermi shock acceleration a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient then the nonlinear backreaction of accelerated particles on the shock structure modifies the shape of the non-thermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well-fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future $\gamma$-ray observatories, such as CTA, very likely.

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.

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.

Starburst regions with multiple powerful winds of young massive stars and supernova remnants are favorable sites for high-energy cosmic ray acceleration. A supernova shock colliding with a fast wind from a compact cluster of young stars allows the acceleration of protons to energies well above the standard limits of diffusive shock acceleration in an isolated SN. The proton spectrum in such a wind-supernova PeV accelerator is hard with a large flux in the high-energy-end of the spectrum producing copious gamma-rays and neutrinos in inelastic nuclear collisions. We argue that SN shocks in the Westerlund 1 cluster in the Milky Way may accelerate protons to about 40 PeV. Once accelerated, these CRs will diffuse into surrounding dense clouds and produce neutrinos with fluxes sufficient to explain a fraction of the events detected by IceCube Observatory from the inner Galaxy.

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.

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 $\gamma$-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.

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.

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.

The giant molecular clouds (MCs) found in the Milky Way and similar galaxies play a crucial role in the evolution of these systems. The supernova explosions that mark the death of massive stars in these regions often lead to interactions between the supernova remnants (SNRs) and the clouds. These interactions have a profound effect on our understanding of SNRs. Shocks in SNRs should be capable of accelerating particles to cosmic ray (CR) energies with efficiencies high enough to power Galactic CRs. X-ray and gamma-ray studies have established the presence of relativistic electrons and protons is some SNRs and provided strong evidence for diffusive shock acceleration as the primary acceleration mechanism, including strongly amplified magnetic fields, temperature and ionization effects on the shock-heated plasmas, and modifications to the dynamical evolution of some systems. Because protons dominate the overall energetics of the CRs, it is crucial to understand this hadronic component even though electrons are much more efficient radiators and it can be difficult to identify the hadronic component. However, near MCs the densities are sufficiently high to allow the gamma-ray emission to be dominated by protons. Thus, these interaction sites provide some of our best opportunities to constrain the overall energetics of these particle accelerators. Here we summarize some key properties of interactions between SNRs and MCs, with an emphasis on recent X-ray and gamma-ray studies that are providing important constraints on our understanding of cosmic rays in our Galaxy.

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.

Tycho's supernova remnant (SNR) is well-established as a source of particle acceleration to very high energies. Constraints from numerous studies indicate that the observed gamma-ray emission results primarily from hadronic processes, providing direct evidence of highly relativistic ions that have been accelerated by the SNR. Here we present an investigation of the dynamical and spectral evolution of Tycho's SNR by carrying out hydrodynamical simulations that include diffusive shock acceleration of particles in the amplified magnetic field at the forward shock of the SNR. Our simulations provide a consistent view of the shock positions, the nonthermal emission, the thermal X-ray emission from the forward shock, and the brightness profiles of the radio and X-ray emission. We compare these with the observed properties of Tycho to determine the density of the ambient material, the particle acceleration efficiency and maximum energy, the accelerated electron to-proton ratio, and the properties of the shocked gas downstream of the expanding SNR shell. We find that evolution of a typical Type Ia supernova in a low ambient density (n_0 ~ 0.3 cm^{-3}), with an upstream magnetic field of ~5\ microGauss, and with ~16% of the SNR kinetic energy being converted into relativistic electrons and ions through diffusive shock acceleration, reproduces the observed properties of Tycho. Under such a scenario, the bulk of observed gamma-ray emission at high energies is produced by pi^0-decay resulting from the collisions of energetic hadrons, while inverse-Compton emission is significant at lower energies, comprising roughly half of the flux between 1 and 10 GeV.

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.

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.

We discuss processes in galactic cosmic ray (GCR) acceleration sites - supernova remnants, compact associations of young massive stars, and superbubbles. Mechanisms of efficient conversion of the mechanical power of the outflows driven by supernova shocks and fast stellar winds of young stars into magnetic fields and relativistic particles are discussed. The high efficiency of particle acceleration in the sources implies the importance of nonlinear feedback effects in a symbiotic relationship where the magnetic turbulence required to accelerate the CRs is created by the accelerated CRs themselves. Non-thermal emission produced by relativistic particles (both those confined in and those that escape from the cosmic accelerators) can be used to constrain the basic physical models of the GCR sources. High resolution X-ray synchrotron imaging, combined with GeV-TeV gamma ray spectra, is a powerful tool to probe the maximum energies of accelerated particles. Future MeV regime spectroscopy will provide unique information on the composition of accelerated particles.

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 {\gamma}-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 {\gamma}-ray emission is produced primarily by leptons accelerated at the SNR forward shock and the other where {\gamma}-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 {\gamma}-ray flux observed.

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.

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.

High-resolution Chandra observations of Tycho's SNR have revealed several sets of quasi-steady, high-emissivity, nearly-parallel X-ray stripes in some localized regions of the SNR. These stripes are most likely the result of cosmic-ray (CR) generated magnetic turbulence at the SNR blast wave. However, for the amazingly regular pattern of these stripes to appear requires the simultaneous action of a number of shock-plasma phenomena and is not predicted by most models of magnetic field amplification. A consistent explanation of these stripes yields information on the complex nonlinear plasma processes connecting efficient CR acceleration and magnetic field fluctuations in strong collisionless shocks. The nonlinear diffusive shock acceleration (NL-DSA) model described here, which includes magnetic field amplification from a cosmic-ray current driven instability, does predict stripes consistent with the synchrotron emission observations of Tycho's SNR. We argue that the local ambient mean magnetic field geometry determines the orientation of the stripes and therefore it can be reconstructed with the high resolution X-ray imaging. The estimated maximum energy of the CR protons responsible for the stripes is $\sim 10^{15}$\,eV. Furthermore, the model predicts that a specific X-ray polarization pattern, with a polarized fraction ~ 50%, accompanies the stripes, which can be tested with future X-ray polarimeter missions.

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.

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.

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.

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.

In order for diffusive shock acceleration (DSA) to accelerate particles to high energies, the energetic particles must be able to interact with magnetic turbulence over a broad wavelength range. The weakly anisotropic distribution of accelerated particles, i.e., cosmic rays (CRs), is believed capable of producing this turbulence in a symbiotic relationship where the magnetic turbulence required to accelerate the CRs is created by the accelerated CRs themselves. In efficient DSA, this wave-particle interaction can be strongly nonlinear where CRs modify the plasma flow and the specific mechanisms of magnetic field amplification. Resonant interactions have long been known to amplify magnetic fluctuations on the scale of the CR gyroradius, and Bell (2004) showed that the CR current can efficiently amplify magnetic fluctuations with scales smaller than the CR gyroradius. Here, we show with a multi-scale, quasi-linear analysis that the presence of turbulence with scales shorter than the CR gyroradius enhances the growth of modes with scales longer than the gyroradius, at least for particular polarizations. We use a mean-field approach to average the equation of motion and the induction equation over the ensemble of magnetic field oscillations accounting for the anisotropy of relativistic particles on the background plasma. We derive the response of the magnetized CR current on magnetic field fluctuations and show that, in the presence of short-scale, Bell-type turbulence, long wavelength modes are amplified. The polarization, helicity, and angular dependence of the growth rates are calculated for obliquely propagating modes for wavelengths both below and above the CR mean free path. The long wavelength growth rates we estimate for typical supernova remnant parameters are sufficiently fast to suggest a fundamental increase in the maximum CR energy a given shock can produce.

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.

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.

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.

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.

This is a report on the findings of the SNR/cosmic-ray working group for the white paper on the status and future of ground-based gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

The highly amplified magnetic fields suggested by observations of some supernova remnant (SNR) shells are most likely an intrinsic part of efficient particle acceleration by shocks. This strong turbulence, which may result from cosmic ray driven instabilities, both resonant and non-resonant, in the shock precursor, is certain to play a critical role in self-consistent, nonlinear models of strong, cosmic ray modified shocks. Here we present a Monte Carlo model of nonlinear diffusive shock acceleration (DSA) accounting for magnetic field amplification through resonant instabilities induced by accelerated particles, and including the effects of dissipation of turbulence upstream of a shock and the subsequent precursor plasma heating. Feedback effects between the plasma heating due to turbulence dissipation and particle injection are strong, adding to the nonlinear nature of efficient DSA. Describing the turbulence damping in a parameterized way, we reach two important results: first, for conditions typical of supernova remnant shocks, even a small amount of dissipated turbulence energy (~10%) is sufficient to significantly heat the precursor plasma, and second, the heating upstream of the shock leads to an increase in the injection of thermal particles at the subshock by a factor of several. In our results, the response of the non-linear shock structure to the boost in particle injection prevented the efficiency of particle acceleration and magnetic field amplification from increasing. We argue, however, that more advanced (possibly, non-resonant) models of turbulence generation and dissipation may lead to a scenario in which particle injection boost due to turbulence dissipation results in more efficient acceleration and even stronger amplified magnetic fields than without the dissipation.

We present a 3-dimensional model of supernova remnants (SNRs) where the hydrodynamical evolution of the remnant is modeled consistently with nonlinear diffusive shock acceleration occuring at the outer blast wave. The model includes particle escape and diffusion outside of the forward shock, and particle interactions with arbitrary distributions of external ambient material, such as molecular clouds. We include synchrotron emission and cooling, bremsstrahlung radiation, neutral pion production, inverse-Compton (IC), and Coulomb energy-loss. Boardband spectra have been calculated for typical parameters including dense regions of gas external to a 1000 year old SNR. In this paper, we describe the details of our model but do not attempt a detailed fit to any specific remnant. We also do not include magnetic field amplification (MFA), even though this effect may be important in some young remnants. In this first presentation of the model we don't attempt a detailed fit to any specific remnant. Our aim is to develop a flexible platform, which can be generalized to include effects such as MFA, and which can be easily adapted to various SNR environments, including Type Ia SNRs, which explode in a constant density medium, and Type II SNRs, which explode in a pre-supernova wind. When applied to a specific SNR, our model will predict cosmic-ray spectra and multi-wavelength morphology in projected images for instruments with varying spatial and spectral resolutions. We show examples of these spectra and images and emphasize the importance of measurements in the hard X-ray, GeV, and TeV gamma-ray bands for investigating key ingredients in the acceleration mechanism, and for deducing whether or not TeV emission is produced by IC from electrons or neutral pions from protons.

Evidence is accumulating suggesting that collisionless shocks in supernova remnants (SNRs) can amplify the interstellar magnetic field to hundreds of microgauss or even milli-gauss levels, as recently claimed for SNR RX J1713.7-3946. If these fields exist, they are almost certainly created by magnetic field amplification (MFA) associated with the efficient production of cosmic rays by diffusive shock acceleration (DSA) and their existence strengthens the case for SNRs being the primary source of galactic cosmic ray ions to the `knee' and beyond. However, the high magnetic field values in SNRs are obtained exclusively from the interpretation of observations of radiation from relativistic electrons and if MFA via nonlinear DSA produces these fields the magnetic field that determines the maximum ion energy will be substantially less than the field that determines the maximum electron energy. We use results of a steady-state Monte Carlo simulation to show how nonlinear effects from efficient cosmic ray production and MFA reduce the maximum energy of protons relative to what would be expected from test-particle acceleration.

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.

This paper has been withdrawn by the authors due to missing references to earlier work and an error in the interpretation of Figures 2 and 3. A corrected version is in preparation.

We introduce a Monte Carlo model of nonlinear diffusive shock acceleration allowing for the generation of large-amplitude magnetic turbulence. The model is the first to include strong wave generation, efficient particle acceleration to relativistic energies in nonrelativistic shocks, and thermal particle injection in an internally self-consistent manner. We find that the upstream magnetic field can be amplified by large factors and show that this amplification depends strongly on the ambient Alfven Mach number. We also show that in the nonlinear model large increases in the magnetic field do not necessarily translate into a large increase in the maximum particle momentum a particular shock can produce, a consequence of high momentum particles diffusing in the shock precursor where the large amplified field converges to the low ambient value. To deal with the field growth rate in the regime of strong fluctuations, we extend to strong turbulence a parameterization that is consistent with the resonant quasi-linear growth rate in the weak turbulence limit. We believe our parameterization spans the maximum and minimum range of the fluctuation growth and, within these limits, we show that the nonlinear shock structure, acceleration efficiency, and thermal particle injection rates depend strongly on the yet to be determined details of wave growth in strongly turbulent fields. The most direct application of our results will be to estimate magnetic fields amplified by strong cosmic-ray modified shocks in supernova remnants.

In the framework of test-particle and cosmic-ray modified hydrodynamics, we calculate synchrotron emission radial profiles in young ejecta-dominated supernova remnants (SNRs) evolving in an ambient medium which is uniform in density and magnetic field. We find that, even without any magnetic field amplification by Raleigh-Taylor instabilities, the radio synchrotron emission peaks at the contact discontinuity because the magnetic field is compressed and is larger there than at the forward shock. The X-ray synchrotron emission sharply drops behind the forward shock as the highest energy electrons suffer severe radiative losses.

We present results from a fully relativistic Monte Carlo simulation of diffusive shock acceleration (DSA) in unmodified (i.e., test-particle) shocks. The computer code uses a single algorithmic sequence to smoothly span the range from nonrelativistic speeds to fully relativistic shocks of arbitrary obliquity, providing a powerful consistency check. We show the dependence of the particle spectrum, which can differ strongly from the canonical f(p) proportional to p^{-4.23} result, on the obliquity and on the strength and ``fineness'' of scattering for a range of shock Lorentz factors. The Monte Carlo results are also shown to be consistent with a simple relation for the spectral index given by Keshet and Waxman (2005) in parallel shocks when the diffusion is ``fine'' and isotropic.

Particle acceleration in collisionless astrophysical shocks, i.e., diffusive shock acceleration (DSA), is the most likely mechanism for producing cosmic rays, at least below 10^{15} eV. Despite the success of this theory, several key elements, including the injection of thermal particles, remains poorly understood. We investigate injection in strongly nonlinear shocks by comparing a semi-analytic model of DSA with a Monte Carlo model. These two models treat injection quite differently and we show, for a particular set of parameters, how these differences influence the overall acceleration efficiency and the shape of the broad-band distribution function.

The strong shocks in young supernova remnants (SNRs) should accelerate cosmic rays (CRs) and no doubt exists that relativistic electrons are produced in SNRs. However, direct evidence that SNRs produce CR nuclei depends on seeing an unambiguous pion-decay feature and this has not yet been obtained. Nevertheless, the lack of an observed pion-decay feature does not necessarily mean that CR ions are not abundantly produced since ions do not radiate efficiently. If CR ions are produced efficiently by diffusive shock acceleration (DSA), their presence will modify the hydrodynamics of the SNR and produce morphological effects which can be clearly seen in radiation produced by electrons.

The strong shocks in young supernova remnants (SNRs) should accelerate cosmic rays (CRs) and no doubt exists that \rel electrons are produced in SNRs. However, direct and convincing evidence that SNRs produce CR nuclei has not yet been obtained and may, in fact, be long in coming if current gamma-ray observatories do not see an unambiguous pion-decay feature. Nevertheless, the lack of an observed pion-decay feature does not necessarily mean that CR ions are not abundantly produced since ions do not radiate efficiently. If CR ions are produced efficiently by diffusive shock acceleration (DSA), their presence will modify the hydrodynamics of the SNR and produce morphological effects which can be clearly seen in radiation produced by electrons. We describe in some detail our CR-hydro model, which couples DSA with the remnant hydrodynamics, and the synchrotron emission expected for two distinct parameter sets representing type Ia and type II supernovae. Several morphological features emerge in radial profiles, including the forward shock precursor, which are observable with current X-ray observatories and which should definitively show if young SNRs produce CR ions efficiently or not. For the specific case of SN 1006, we conclude, as have others, that the extremely short X-ray scale heights observed near the outer shock argue convincingly for the efficient production of CR ions.

The remnant G347.3-0.5 exhibits strong shell emission in the radio and X-ray bands, and has a purported detection in the TeV gamma-ray band by the CANGAROO-II telescope. The CANGAROO results were touted as evidence for the production of cosmic ray ions, a claim that has proven controversial due to constraining fluxes associated with a proximate unidentified EGRET source 3EG J1714-3857. HESS has now seen this source in the TeV band. The complex environment of the remnant renders modeling of its broadband spectrum sensitive to assumptions concerning the nature and parameters of the circumremnant medium. This paper explores a sampling of reasonable possibilities for multiwavelength spectral predictions from this source, using a non-linear model of diffusive particle acceleration at the shocked shell. The magnetic field strength, shell size and degree of particle cross-field diffusion act as variables to which the radio to X-ray to gamma-ray signal is sensitive. The modeling of the extant data constrains these variables, and the potential impact of the recent HESS detection on such parameters is addressed. Putative pion decay signals in hard gamma-rays resulting from hadronic interactions in dense molecular clouds are briefly discussed; the requisite suppression of the GeV component needed to accommodate the 3EG J1714-3857 EGRET data provides potential bounds on the diffusive distance from the shell to the proximate clouds.

Without amplification, magnetic fields in expanding ejecta of young supernova remnants (SNRs) will be orders of magnitude below those required to shock accelerate thermal electrons, or ions, to relativistic energies or to produce radio synchrotron emission at the reverse shock. The reported observations of such emission give support to the idea that diffusive shock acceleration (DSA) can amplify magnetic fields by large factors. Furthermore, the uncertain character of the amplification process leaves open the possibility that ejecta fields, while large enough to support radio emission and DSA, may be much lower than typical interstellar medium values. We show that DSA in such low reverse shock fields is extremely nonlinear and efficient in the production of cosmic-ray (CR) ions, although CRs greatly in excess of mc^2 are not produced. These nonlinear effects, which occur at the forward shock as well, are manifested most importantly in shock compression ratios much greater than four and cause the interaction region between the forward and reverse shocks to become narrower, denser, and cooler than would be the case if efficient cosmic-ray production did not occur. The changes in the SNR structure and evolution should be clearly observable, if present, and they convey important information on the nature of DSA and magnetic field amplification with broad astrophysical implications.

We present results from a fully relativistic Monte Carlo simulation of diffusive shock acceleration (DSA) in unmodified shocks. The computer code uses a single algorithmic sequence to smoothly span the range from nonrelativistic speeds to fully relativistic shocks of arbitrary obliquity, providing a powerful consistency check. While known results are obtained for nonrelativistic and ultra-relativistic parallel shocks, new results are presented for the less explored trans- relativistic regime and for oblique, fully relativistic shocks. We find, for a wide trans-relativistic range extending to shock Lorentz factors >30, that the particle spectrum produced by DSA varies strongly from the canonical f(p) proportional to p^{-4.23} spectrum known to result in ultra-relativistic shocks. Trans- relativistic shocks may play an important role in gamma-ray bursts and other sources and most relativistic shocks will be highly oblique.

A number of supernova remnants (SNRs) show nonthermal X-rays assumed to be synchrotron emission from shock accelerated TeV electrons. The existence of these TeV electrons strongly suggests that the shocks in SNRs are sources of galactic cosmic rays (CRs). In addition, there is convincing evidence from broad-band studies of individual SNRs and elsewhere that the particle acceleration process in SNRs can be efficient and nonlinear. If SNR shocks are efficient particle accelerators, the production of CRs impacts the thermal properties of the shock heated, X-ray emitting gas and the SNR evolution. We report on a technique that couples nonlinear diffusive shock acceleration, including the backreaction of the accelerated particles on the structure of the forward and reverse shocks, with a hydrodynamic simulation of SNR evolution. Compared to models which ignore CRs, the most important hydrodynamical effects of placing a significant fraction of shock energy into CRs are larger shock compression ratios and lower temperatures in the shocked gas. We compare our results, which use an approximate description of the acceleration process, with a more complete model where the full CR transport equations are solved (i.e., Berezhko et al., 2002), and find excellent agreement for the CR spectrum summed over the SNR lifetime and the evolving shock compression ratio. The importance of the coupling between particle acceleration and SNR dynamics for the interpretation of broad-band continuum and thermal X-ray observations is discussed.

Shock jump conditions are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically without additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., about 20%) produce significant changes in the shock compression ratio, r, at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r=3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, 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.

The supernova remnant G347.3--0.5 emits a featureless power-law in X-rays, thought to indicate shock-acceleration of electrons to high energies. We here produce a broad-band spectrum of the bright NW limb of this source by combining radio observations from the Australia Telescope Compact Array (ATCA), X-ray observations from the Advanced Satellite for Cosmology and Astrophysics (ASCA), and TeV gamma-ray observations from the CANGAROO imaging Cerenkov telescope. We assume this emission is produced by an electron population generated by diffusive shock acceleration at the remnant forward shock. The nonlinear aspects of the particle acceleration force a connection between the widely different wavelength bands and between the electrons and the unseen ions, presumably accelerated simultaneously with the electrons. This allows us to infer the relativistic proton spectrum and estimate ambient parameters such as the supernova explosion energy, magnetic field, matter density in the emission region, and efficiency of the shock acceleration process. We find convincing evidence that the shock acceleration is efficient, placing >25% of the shock kinetic energy flux into relativistic ions. Despite this high efficiency, the maximum electron and proton energies, while depending somewhat on assumptions for the compression of the magnetic field in the shock, are well below the observed `knee' at about 10^{15} eV in the Galactic cosmic-ray spectrum.

We employ hydrodynamic simulations to study the effects of high shock compression ratios, as expected for fast shocks with efficient particle acceleration, on the convective instability of driven waves in supernova remnants. We find that the instability itself does not depend significantly on the compression ratio, but because the width of the interaction region between the forward and reverse shocks can shrink significantly with increasing shock compression, we find that convective instabilities can reach all the way to the forward shock front if compression ratios are high enough.

This report discusses the possibilities for observational tests of current models and ideas regarding the origin of the Galactic Cosmic Rays using facilities which are already available.

We have developed a simple model to investigate the modifications of the hydrodynamics and non-equilibrium ionization X-ray emission in young supernova remnants due to nonlinear particle acceleration. In nonlinear, diffusive shock acceleration, the heating of the gas to X-ray emitting temperatures is strongly coupled to the acceleration of cosmic ray ions. If the acceleration is efficient and a significant fraction of the shock ram energy ends up in cosmic rays, compression ratios will be higher and the shocked temperature lower than test-particle, Rankine-Hugoniot relations predict. We illustrate how particle acceleration impacts the interpretation of X-ray data using the X-ray spectra of Kepler's remnant, observed by ASCA and RXTE. The thermal X-ray emission provides important constraints on the efficiency of particle acceleration, in complement to nonthermal emission. X-ray data from Chandra and XMM Newton, plus radio observations, will be essential to quantify nonlinear effects.

In evolved supernova remnants (SNRs) interacting with molecular clouds, such as IC 443, W44, and 3C391, a highly inhomogeneous structure consisting of a forward shock of moderate Mach number, a cooling layer, a dense radiative shell and an interior region filled with hot tenuous plasma is expected. We present a kinetic model of nonthermal electron injection, acceleration and propagation in that environment and find that these SNRs are efficient electron accelerators and sources of hard X- and gamma-ray emission. The energy spectrum of the nonthermal electrons is shaped by the joint action of first and second order Fermi acceleration in a turbulent plasma with substantial Coulomb losses. Bremsstrahlung, synchrotron, and inverse Compton radiation of the nonthermal electrons produce multiwavelength photon spectra in quantitative agreement with the radio and the hard emission observed by ASCA and EGRET from IC 443. We distinguish interclump shock wave emission from molecular clump shock wave emission accounting for a complex structure of molecular cloud. Spatially resolved X- and gamma- ray spectra from the supernova remnants IC 443, W44, and 3C391 as might be observed with BeppoSAX, Chandra XRO, XMM, INTEGRAL and GLAST would distinguish the contribution of the energetic lepton component to the gamma-rays observed by EGRET.

Cosmic-ray production in young supernova remnant (SNR) shocks is expected to be efficient and strongly nonlinear. In nonlinear, diffusive shock acceleration, compression ratios will be higher and the shocked temperature lower than test-particle, Rankine-Hugoniot relations predict. Furthermore, the heating of the gas to X-ray emitting temperatures is strongly coupled to the acceleration of cosmic-ray electrons and ions, thus nonlinear processes which modify the shock, influence the emission over the entire band from radio to gamma-rays and may have a strong impact on X-ray line models. Here we apply an algebraic model of nonlinear acceleration, combined with SNR evolution, to model the radio and X-ray continuum of Kepler's SNR.

We have extended a simple model of nonlinear diffusive shock acceleration (Berezhko and Ellison 1999) to include the injection and acceleration of electrons and the production of photons from bremsstrahlung, synchrotron, inverse-Compton, and pion-decay processes. We argue that the results of this model, which is simpler to use than more elaborate ones, offer a significant improvement over test-particle, power-law spectra which are often used in astrophysical applications of diffusive shock acceleration. With an evolutionary supernova remnant (SNR) model to obtain shock parameters as functions of ambient interstellar medium parameters and time, we predict broad-band continuum photon emission from supernova remnants in general, and SN1006 in particular, showing that our results compare well with the more complete time-dependent and spherically symmetric nonlinear model of Berezhko, Ksenofontov, and Petukhov (1999). We discuss the implications nonlinear shock acceleration has for X-ray line emission, and use our model to describe how ambient conditions determine the TeV/radio flux ratio, an important parameter for gamma-ray observations of radio SNRs.

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 moderately strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He$^{+}$/H$^+$ or O$^{+}$/H$^+$ to within a factor ~ 3. 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.

We calculate particle spectra and continuum photon emission from the Cassiopeia A supernova remnant (SNR). The particle spectra, ion and electron, result from diffusive shock acceleration at the forward SNR shock and are determined with a nonlinear Monte Carlo calculation. The calculation self-consistently determines the shock structure under the influence of ion pressure, and includes a simple parameterized treatment of electron injection and acceleration. Our results are compared to photon observations, concentrating on the connection between the Radio and GeV-TeV gamma-ray range, and to cosmic ray ion observations. We include new upper limits from the Cherenkov Array at Themis (CAT) imaging Cherenkov telescope and the Whipple 10m gamma-ray telescope at > 400 GeV. These new limits support the suggestion (e.g. Cowsik & Sarkar 1980; Allen et. al. 1997) that energetic electrons are emitting synchrotron radiation in an extremely high magnetic field (~ 1000 microGauss), far greater than values routinely assigned to the ISM, and help to constrain our model. The large magnetic field allows acceleration of cosmic ray ions to well above $10^{15}$ eV per nucleon in the ~ 300 yr lifetime of Cas A.

There has recently been interest in the role of inverse bremsstrahlung, the emission of photons by fast suprathermal ions in collisions with ambient electrons possessing relatively low velocities, in tenuous plasmas in various astrophysical contexts. This follows a long hiatus in the application of suprathermal ion bremsstrahlung to astrophysical models since the early 1970s. The potential importance of inverse bremsstrahlung relative to normal bremsstrahlung, i.e. where ions are at rest, hinges upon the underlying velocity distributions of the interacting species. In this paper, we identify the conditions under which the inverse bremsstrahlung emissivity is significant relative to that for normal bremsstrahlung in shocked astrophysical plasmas. We determine that, since both observational and theoretical evidence favors electron temperatures almost comparable to, and certainly not very deficient relative to proton temperatures in shocked plasmas, these environments generally render inverse bremsstrahlung at best a minor contributor to the overall emission. Hence inverse bremsstrahlung can be safely neglected in most models invoking shock acceleration in discrete sources such as supernova remnants. However, on scales > 100pc distant from these sources, Coulomb collisional losses can deplete the cosmic ray electrons, rendering inverse bremsstrahlung, and perhaps bremsstrahlung from knock-on electrons, possibly detectable.

The recent detection of a hard X-ray component in the supernova remnant Cassiopeia A is interpreted as synchrotron emission from electrons accelerated to energies up to 40 TeV (Allen et al., 1997). It is therefore tempting to consider TeV gamma-ray emission from this object through : i) bremsstrahlung and inverse Compton radiation from electrons and/or ii) pi0 production from an associated high energy cosmic ray component hitting surrounding material. Cas A was observed by the CAT imaging Cherenkov telescope during the observing season Aug-Nov 1998. An upper limit to the integral flux above 400 GeV of 0.74x10^-11 gamma.cm^-2.s^-1 is derived. This result is used to constrain shock-acceleration models for production of VHE gamma-rays in SNRs.

Assuming that refractory elements in cosmic rays originate in dust grains, we examine the viability of cosmic ray origin models wherein the bulk of present day cosmic rays are accelerated out of fresh supernova ejecta material before it mixes with the average interstellar medium (ISM). We conclude that the fresh ejecta scenarios that have been proposed thus far have serious flaws, and are unable to account for known properties of present day cosmic rays. These flaws include: (1) the small fraction of ejecta to ISM mass processed by the forward supernova remnant (SNR) shock; (2) the difficulty fresh ejecta grains have in reaching the forward shock in isolated SNRs, and the small expected sputtering yield, especially ahead of the shock, even if grains do reach the forward shock; (3) the implausibility that fresh ejecta material can dominate cosmic ray production in diffuse superbubbles; and (4) the lack of a connection in fresh ejecta models between the production of cosmic ray refractory and volatile elements. We conclude that the near linear increase in Be abundance with metallicity observed in old, halo stars cannot imply that a significant fraction of the cosmic rays seen today come from fresh supernova ejecta. This conclusion is supported by the analysis of the present day cosmic ray composition, as shown in Meyer & Ellison, this volume.

The composition of present day cosmic rays (CRs) is inconsistent with a significant acceleration of SN ejecta material (even with a preferential acceleration of ejecta grain material). Current CRs must result mainly from the acceleration of ISM ("solar mix") and circumstellar (22Ne- and 12C-rich WR wind) gas and grain material. The CR source composition, indeed, shows no anomaly related to SN nucleosynthesis: (i) The CR source FeNi/MgSiCa ratios have precisely solar mix values, while FeNi are mainly synthesized in SN Ia's, and MgSiCa in SN II's. If CR's originate in SN ejecta, this requires tight conditions on the acceleration efficiencies of the various SN Ia's and SN II's. (ii) The non-deficiency of the main-s-process elements, not made in any SN, relative to all elements made in SNae. (iii) All CR source isotope ratios are consistent with solar mix (except for the 22Ne and 12C excesses, WR wind). (iv) The absence of 59Ni in CRs implies that the time delay between the SN nucleosynthesis and their acceleration is > 10^5 yr. (v) The physics of SNR's and of CR shock acceleration implies that the acceleration of interior ejecta material is comparatively insignificant (Ellison & Meyer, this volume). Predominant acceleration of current CRs out of superbubble material also seems implausible. These conclusions regarding current CRs do not necessarily conflict with the linear evolution of Be/H in the early Galaxy. With the near absence of heavy elements in the early Galactic ISM, indeed, the acceleration of even a minute amount of freshly processed material must have then played a dominant role for the generation of Be from C and O. The "Be indicator" is blind to a possibly dominant early Galactic CR component originating in the ISM then composed of virtually pure H and He.

Detailed composition measurements can be a very powerful means of tracing origins, a fact used regularly by forensic scientists and art historians. One of the main motivating factors for making detailed observations of cosmic rays was always the hope that a unique compositional signature could be found which pointed unambiguously to a particular source. This has proven much harder than expected, but we have now reached a point where it appears possible to begin to decipher the information contained in the compositional data; the key, we have discovered, is to read the data not in isolation, but in the context provided by our general astronomical knowledge and by recent developments in shock acceleration theory (Meyer, Drury and Ellison, 1997, 1998; Ellison, Drury and Meyer, 1997). In our view (not, it is only fair to warn the reader, yet universally accepted) the data show clearly that the Galactic cosmic ray particles originate predominantly from the gas and dust of the general interstellar medium.

Supernova remnants (SNRs) are widely believed to be the principal source of galactic cosmic rays. Such energetic particles can produce gamma-rays and lower energy photons via interactions with the ambient plasma. In this paper, we present results from a Monte Carlo simulation of non-linear shock structure and acceleration coupled with photon emission in shell-like SNRs. These non-linearities are a by-product of the dynamical influence of the accelerated cosmic rays on the shocked plasma and result in distributions of cosmic rays which deviate from pure power-laws. Such deviations are crucial to acceleration efficiency and spectral considerations, producing GeV/TeV intensity ratios that are quite different from test particle predictions. The Sedov scaling solution for SNR expansions is used to estimate important shock parameters for input into the Monte Carlo simulation. We calculate ion and electron distributions that spawn neutral pion decay, bremsstrahlung, inverse Compton, and synchrotron emission, yielding complete photon spectra from radio frequencies to gamma-ray energies. The cessation of acceleration caused by the spatial and temporal limitations of the expanding SNR shell in moderately dense interstellar regions can yield spectral cutoffs in the TeV energy range; these are consistent with Whipple's TeV upper limits on unidentified EGRET sources. Supernova remnants in lower density environments generate higher energy cosmic rays that produce predominantly inverse Compton emission at super-TeV energies; such sources will generally be gamma-ray dim at GeV energies.

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.

Diffusive shock acceleration in the environs of a remnant's expanding shell is a popular candidate for the origin of SNR gamma-rays. In this paper, results from our study of non-linear effects in shock acceleration theory and their impact on the gamma-ray spectra of SNRs are presented. These effects describe the dynamical influence of the accelerated cosmic rays on the shocked plasma at the same time as addressing how the non-uniformities in the fluid flow force the distribution of the cosmic rays to deviate from pure power-laws. Such deviations are crucial to gamma-ray spectral determination. Our self-consistent Monte Carlo approach to shock acceleration is used to predict ion and electron distributions that spawn neutral pion decay, bremsstrahlung and inverse Compton emission components for SNRs. We demonstrate how the spatial and temporal limitations imposed by the expanding SNR shell quench acceleration above critical energies in the 500 GeV - 10 TeV range, thereby spawning gamma-ray spectral cutoffs that are quite consistent with Whipple's TeV upper limits to the EGRET unidentified sources that have SNR associations. We also discuss the role of electron injection in shocks and its impact on the significance of electromagnetic components to GeV--TeV spectral formation.

This is the second paper (the first was astro-ph/9704267) of a series analysing the Galactic Cosmic Ray (GCR) composition and origin. In this we present a quantitative model of GCR origin and acceleration based on the acceleration of a mixture of interstellar and/or circumstellar gas and dust by supernova remnant blast waves. We present results from a nonlinear shock model which includes (i) the direct acceleration of interstellar gas-phase ions, (ii) a simplified model for the direct acceleration of weakly charged dust grains to energies of order 100keV/amu simultaneously with the gas ions, (iii) frictional energy losses of the grains colliding with the gas, (iv) sputtering of ions of refractory elements from the accelerated grains and (v) the further shock acceleration of the sputtered ions to cosmic ray energies. The calculated GCR composition and spectra are in good agreement with observations.

This is the first of a series of papers analysing the Galactic Cosmic Ray composition and origin. We show that the Galactic Cosmic Ray source (GCRS) composition is best described in terms of (i) a general enhancement of the refractory elements relative to the volatile ones, and (ii) among the volatile elements, an enhancement of the heavier elements relative to the lighter ones; this mass dependence most likely reflects a mass-to-charge (A/Q) dependence of the acceleration efficiency; among the refractory elements, there is NO such enhancement of heavier species, or only a much weaker one. We regard as coincidental the similarity between the GCRS composition and that of the solar corona, which is biased according to first ionization potential. In a companion paper, this GCRS composition is interpreted in terms of an acceleration by supernova shock waves of interstellar and/or circumstellar (eg Ne22 rich Wolf-Rayet wind) gas-phase and especially dust material.

While the definitive detection of gamma-rays from known supernova remnants (SNRs) remains elusive, the collection of unidentified EGRET sources that may be associated with SNRs has motivated recent modelling of TeV emission from these sources. Current theoretical models use power-law shock-accelerated protons and electrons in their predictions of expected gamma-ray TeV fluxes from those unidentified EGRET sources with remnant associations. In this paper, we explore a more detailed non-linear shock acceleration model, which generates non-thermal proton distributions and includes a self-consistent determination of shock hydrodynamics. We obtain gamma-ray spectra for SNRs allowing for the cessation of acceleration to high energies that is due to the finite ages and sizes of remnants. Gamma-ray spectral cutoffs can be observed in the TeV range for reasonable remnant parameters, and deviations from power-law behaviour are found at all energies ranging from 1 MeV up to the cutoff. Correlated observations by INTEGRAL, Whipple and other instruments may provide stringent constraints to our understanding of supernova remnants.

The most stringent test of theoretical models of the first-order Fermi mechanism at collisionless astrophysical shocks is a comparison of the theoretical predictions with observational data on particle populations. Such comparisons have yielded good agreement between observations at the quasi-parallel portion of the Earth's bow shock and three theoretical approaches, including Monte Carlo kinetic simulations. This paper extends such model testing to the realm of oblique interplanetary shocks: here observations of proton and alpha particle distributions made by the SWICS ion mass spectrometer on Ulysses at nearby interplanetary shocks are compared with test particle Monte Carlo simulation predictions of accelerated populations. The plasma parameters used in the simulation are obtained from measurements of solar wind particles and the magnetic field upstream of individual shocks. Good agreement between downstream spectral measurements and the simulation predictions are obtained for two shocks by allowing the the ratio of the mean-free scattering length to the ionic gyroradius, to vary in an optimization of the fit to the data. Generally small values of this ratio are obtained, corresponding to the case of strong scattering. The acceleration process appears to be roughly independent of the mass or charge of the species.

We have developed a Monte Carlo technique for self-consistently calculating the hydrodynamic structure of oblique, steady-state shocks, together with the first-order Fermi acceleration process and associated non-thermal particle distributions. This is the first internally consistent treatment of modified shocks that includes cross-field diffusion of particles. Our method overcomes the injection problem faced by analytic descriptions of shock acceleration, and the lack of adequate dynamic range and artificial suppression of cross-field diffusion faced by plasma simulations; it currently provides the most broad and versatile description of collisionless shocks undergoing efficient particle acceleration. We present solutions for plasma quantities and particle distributions upstream and downstream of shocks, illustrating the strong differences observed between non-linear and test-particle cases. It is found that there are only marginal differences in the injection efficiency and resultant spectra for two extreme scattering modes, namely large-angle scattering and pitch-angle diffusion, for a wide range of shock parameters, i.e., for subluminal shocks with field obliquities less than or equal to 75 degrees and de Hoffmann-Teller frame speeds much less than the speed of light.

The rate at which particles are accelerated by the first-order Fermi mechanism in shocks depends on the angle, \teq{\Tbone}, that the upstream magnetic field makes with the shock normal. The greater the obliquity the greater the rate, and in quasi-perpendicular shocks rates can be hundreds of times higher than those seen in parallel shocks. In many circumstances pertaining to evolving shocks (\eg, supernova blast waves and interplanetary traveling shocks), high acceleration rates imply high maximum particle energies and obliquity effects may have important astrophysical consequences. However, as is demonstrated here, the efficiency for injecting thermal particles into the acceleration mechanism also depends strongly on obliquity and, in general, varies inversely with \teq{\Tbone}. The degree of turbulence and the resulting cross-field diffusion strongly influences both injection efficiency and acceleration rates. The test particle \mc simulation of shock acceleration used here assumes large-angle scattering, computes particle orbits exactly in shocked, laminar, non-relativistic flows, and calculates the injection efficiency as a function of obliquity, Mach number, and degree of turbulence. We find that turbulence must be quite strong for high Mach number, highly oblique shocks to inject significant numbers of thermal particles and that only modest gains in acceleration rates can be expected for strong oblique shocks over parallel ones if the only source of seed particles is the thermal background.