Volume 616, August 2018
|Number of page(s)||9|
|Section||Numerical methods and codes|
|Published online||15 August 2018|
Variability of the adiabatic parameter in monoatomic thermal and non-thermal plasmas★
Department of Mathematics, University of Évora,
R. Romão Ramalho 59,
2 Zentrum für Astronomie und Astrophysik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
Accepted: 19 April 2018
Context. Numerical models of the evolution of interstellar and integalactic plasmas often assume that the adiabatic parameter γ (the ratio of the specific heats) is constant (5/3 in monoatomic plasmas). However, γ is determined by the total internal energy of the plasma, which depends on the ionic and excitation state of the plasma. Hence, the adiabatic parameter may not be constant across the range of temperatures available in the interstellar medium.
Aims. We aim to carry out detailed simulations of the thermal evolution of plasmas with Maxwell–Boltzmann and non-thermal (κ and n) electron distributions in order to determine the temperature variability of the total internal energy and of the adiabatic parameter.
Methods. The plasma, composed of H, He, C, N, O, Ne, Mg, Si, S, and Fe atoms and ions, evolves under collisional ionization equilibrium conditions, from an initial temperature of 109 K. The calculations include electron impact ionization, radiative and dielectronic recombinations and line excitation. The ionization structure was calculated solving a system of 112 linear equations using the Gauss elimination method with scaled partial pivoting. Numerical integrations used in the calculation of ionization and excitation rates are carried out using the double-exponential over a semi-finite interval method. In both methods a precision of 10−15 is adopted.
Results. The total internal energy of the plasma is mainly dominated by the ionization energy for temperatures lower than 8 × 104 K with the excitation energy having a contribution of less than one percent. In thermal and non-thermal plasmas composed of H, He, and metals, the adiabatic parameter evolution is determined by the H and He ionizations leading to a profile in general having three transitions. However, for κ distributed plasmas these three transitions are not observed for κ < 15 and for κ < 5 there are no transitions. In general, γ varies from 1.01 to 5/3. Lookup tables of the γ parameter are presented as supplementary material.
Key words: atomic processes / atomic data / hydrodynamics / methods: numerical / ISM: general / intergalactic medium
Tables of the gamma parameter variability with temperature (Fig. 5) + data for n = 25 and κ = 25 are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (126.96.36.199) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/616/A58
© ESO 2018
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.