We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, A/Z, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock temperature proportional to A. When diffusive shock acceleration is efficient, ions develop a non-thermal tail whose extent scales with Z and whose normalization is enhanced as (A/Z) 2 , so that incompletely-ionized heavy ions are preferentially accelerated. We discuss how these findings can explain observed heavy-ion enhancements in Galactic cosmic rays.
PACS numbers:Introduction.-Non-relativistic shocks are well-known as sources of energetic particles. Prominent examples of such shocks are the blast waves of supernova remnants (SNRs), which are thought to be the sources of Galactic cosmic rays (GCRs) [e.g., 1, 2], and heliospheric shocks, where solar energetic particles (SEPs) are measured in situ [e.g., 3, 4]. Chemical abundances in GCRs and SEPs provide crucial information about their sources and the processes responsible for their acceleration.At trans-relativistic energies, the chemical composition of GCRs roughly resembles the composition of the solar system [5], the most evident deviation being the enhancement of secondaries produced by spallation of primary GCRs during their propagation in the Milky Way. A more careful analysis, however, reveals that the GCR composition is controlled by volatility and mass/charge ratios: refractory elements show larger enhancements than volatile ones, and heavier volatile elements are more abundant than lighter ones [6,7]. Moreover, elemetns with low first ionization potential tend to be overrepresented in GCRs [e.g.,5]. At TeV energies, where spallation is negligible, the fluxes of H, He, C-N-O, and Fe do not differ by more than one order of magnitude [e.g., 8, and references therein]. Since their typical solar number abundances relative to H are χ He = 0.0963, χ CN O = 9.54 × 10 −4 , χ F e = 8.31 × 10 −5 [9], the abundances observed in GCRs suggest that heavy ions must be preferentially injected and accelerated compared to protons.Diffusive shock acceleration (DSA) [e.g., 10, 11] at SNR shocks is likely the mechanism responsible for ion acceleration up to ∼ 10 17 eV [12]. DSA produces universal power-law momentum spectra f (p) ∝ p −3r/(r−1) , where r is the shock compression ratio; for strong shocks r → 4 and f (p) ∝ p −4 . For relativistic particles the energy spectrum is then f (E) = 4πp 2 f (p)dp/dE ∝ E −2 , while at non-relativistic energies one gets f (E) ∝ E −3/2 [13].