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.

Cosmic dust particles effectively attenuate starlight. Their absorption of starlight produces emission spectra from the near- to far-infrared, which depends on the sizes and properties of the dust grains, and spectrum of the heating radiation field. The near- to mid-infrared is dominated by the emissions by very small grains. Modeling the absorption of starlight by these particles is, however, computationally expensive and a significant bottleneck for self-consistent radiation transport codes treating the heating of dust by stars. In this paper, we summarize the formalism for computing the stochastic emissivity of cosmic dust, which was developed in earlier works, and present a new library HEATCODE implementing this formalism for the calculation for arbitrary grain properties and heating radiation fields. Our library is highly optimized for general-purpose processors with multiple cores and vector instructions, with hierarchical memory cache structure. The HEATCODE library also efficiently runs on co-processor cards implementing the Intel Many Integrated Core (Intel MIC) architecture. We discuss in detail the optimization steps that we took in order to optimize for the Intel MIC architecture, which also significantly benefited the performance of the code on general-purpose processors, and provide code samples and performance benchmarks for each step. The HEATCODE library performance on a single Intel Xeon Phi coprocessor (Intel MIC architecture) is approximately 2 times a general-purpose two-socket multicore processor system with approximately the same nominal power consumption. The library supports heterogeneous calculations employing host processors simultaneously with multiple coprocessors, and can be easily incorporated into existing radiation transport codes.

In this review we discuss some observational aspects and theoretical models of astrophysical collisionless shocks in partly ionized plasma with the presence of non-thermal components. A specific feature of fast strong collisionless shocks is their ability to accelerate energetic particles that can modify the shock upstream flow and form the shock precursors. We discuss the effects of energetic particle acceleration and associated magnetic field amplification and decay in the extended shock precursors on the line and continuum multi-wavelength emission spectra of the shocks. Both Balmer-type and radiative astrophysical shocks are discussed in connection to supernova remnants interacting with partially neutral clouds. Quantitative models described in the review predict a number of observable line-like emission features that can be used to reveal the physical state of the matter in the shock precursors and the character of nonthermal processes in the shocks. Implications of recent progress of gamma-ray observations of supernova remnants in molecular clouds are highlighted.

Recent accurate measurements of cosmic-ray (CR) protons and nuclei by ATIC-2, CREAM, and PAMELA reveal: a) unexpected spectral hardening in the spectra of CR species above a few hundred GeV per nucleon, b) a harder spectrum of He compared to protons, and c) softening of the CR spectra just below the break energy. These newly-discovered features may offer a clue to the origin of the observed high-energy Galactic CRs. We discuss possible interpretations of these spectral features and make predictions for the secondary CR fluxes and secondary to primary ratios, anisotropy of CRs, and diffuse Galactic {\gamma}-ray emission in different phenomenological scenarios. Our predictions can be tested by currently running or near-future high-energy astrophysics experiments.

GALPROP is a numerical code for calculating the galactic propagation of relativistic charged particles and the diffuse emissions produced during their propagation. The code incorporates as much realistic astrophysical input as possible together with latest theoretical developments and has become a de facto standard in astrophysics of cosmic rays. We present GALPROP WebRun, a service to the scientific community enabling easy use of the freely available GALPROP code via web browsers. In addition, we introduce the latest GALPROP version 54, available through this service.

We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study.

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.

This research was motivated by the recent observations indicating very strong magnetic fields at some supernova remnant shocks, which suggests in-situ generation of magnetic turbulence. The dissertation presents a numerical model of collisionless shocks with strong amplification of stochastic magnetic fields, self-consistently coupled to efficient shock acceleration of charged particles. Based on a Monte Carlo simulation of particle transport and acceleration in nonlinear shocks, the model describes magnetic field amplification using the state-of-the-art analytic models of instabilities in magnetized plasmas in the presence of non-thermal particle streaming. The results help one understand the complex nonlinear connections between the thermal plasma, the accelerated particles and the stochastic magnetic fields in strong collisionless shocks. Also, predictions regarding the efficiency of particle acceleration and magnetic field amplification, the impact of magnetic field amplification on the maximum energy of accelerated particles, and the compression and heating of the thermal plasma by the shocks are presented. Particle distribution functions and turbulence spectra derived with this model can be used to calculate the emission of observable nonthermal radiation.

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.

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.

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.