Search

Context. Balmer emission may be a powerful diagnostic tool for testing the paradigm of cosmic ray (CR) acceleration in young supernova remnant (SNR) shocks. The width of the broad Balmer line is a direct indicator of the downstream plasma temperature. In the case of efficient particle acceleration, an appreciable fraction of the total kinetic energy of the plasma is channeled into CRs, therefore the downstream temperature decreases and so does the broad Balmer line width. This width also depends on the level of thermal equilibration between ions and neutral hydrogen atoms in the downstream. Since generally only a few charge exchange (CE) reactions occur before ionization in young SNR shocks, equilibration between ions and neutrals is not reached, and a kinetic description of the neutrals is required to properly compute Balmer emission. Aims. We provide a method for calculating Balmer emission using a self-consistent description of the shock structure in the presence of neutrals and CRs, which also accounts for the non-Maxwellian distribution of neutrals. Methods. We use a recently developed semi-analytical approach, where neutral particles, ionized plasma, accelerated particles, and magnetic fields are all coupled together through the mass, momentum, and energy flux-conservation equations. The distribution of neutrals is obtained from the full Boltzmann equation in velocity space, coupled to Maxwellian ions through ionization and CE processes. The computation is also an improvement over previous work thanks to a better approximation of the atomic interaction rates. Results. We find that for shock speeds ≳2500 km s-1, the distribution of broad neutrals never approaches a Maxwellian and its moments differ from those of the ionized component. These differences lead to a smaller FWHM than predicted in previous calculations, where thermalization was assumed. Conclusions. The method presented here provides a realistic estimate of particle acceleration efficiency in Balmer-dominated shocks. [ABSTRACT FROM AUTHOR]

Abstract: We present here a semi-analytical solution of the problem of particle acceleration at non-linear shock waves with a free-escape boundary at some location upstream. This solution, besides allowing us to determine the spectrum of particles accelerated at the shock front, including the shape of the cutoff at some maximum momentum, also allows us to determine the spectrum of particles escaping the system from upstream. This latter aspect of the problem is crucial for establishing a connection between the accelerated particles in astrophysical sources, such as supernova remnants, and the cosmic rays observed at the Earth. An excellent approximate solution, which leads to a computationally fast calculation of the structure of shocks with an arbitrary level of cosmic ray modification, is also obtained. [Copyright &y& Elsevier]