Online material
Appendix A: Cosmicray interaction model
We consider a model in which lowenergy cosmic rays (LECRs) are produced in an unspecified acceleration region and penetrate a nearby cloud of neutral gas at a constant rate (see Fig. A.1). The energetic particles can produce nonthermal Xrays by atomic collisions while they slow down by ionization and radiative energy losses in the dense cloud. We further assume that the LECRs that penetrate the cloud can escape from it after an energyindependent path length Λ, which is a free parameter of the model. The differential equilibrium number of primary CRs of type i (electrons, protons, or α particles) in the cloud is then given by (A.1)where (dQ_{i} / dt) is the differential rate of LECRs injected in the cloud, [dE / dt(E)] _{i} is the CR energy loss rate, and the maximum energy is related to the escape path length Λ (expressed in units of H atoms cm^{2}) by (A.2)where (A.3)Here, v_{i} is the particle velocity, n_{H} the mean number density of H atoms in the cloud, m_{p} the proton mass, a_{He} = 0.0964 the cosmic abundance of He relative to H (Lodders 2003), and (dE / dx)_{i,H} and (dE / dx)_{i,He} the CR stopping powers (in units of MeV g^{1} cm^{2}) in ambient H and He, respectively. We used for electrons the stoppingpower tables of Berger & Seltzer (1982) below 1 GeV and the relativistic formulae given by Schlickeiser (2002) above this energy. The stopping powers for protons and αparticles were extracted from the online databases PSTAR and ASTAR, respectively (Berger et al. 2005).
The process of CR transport in the cloud, which does not need to be specified in the above formalism, is nevertheless relevant to estimate the escape path length Λ from the cloud size. It is clear that if the cloud medium is not diffusive, because of, e.g., efficient ionneutral damping of MHD waves, Λ ~ n_{H}L_{C}, where L_{C} is the characteristic size of the cloud. But otherwise, the escape path length, which can then be estimated as (A.4)can be much greater than the characteristic column density , depending on the diffusion coefficient D. For example, with the typical mean diffusion coefficient for the propagation of Galactic CR nuclei in the local interstellar magnetic field B (Berezinsky et al. 1990), (A.5)where β = v_{i} / c and R_{i} is the particle rigidity, one gets from Eq. (A.4) for nonrelativistic protons: (A.6)where , n_{H} and B are scaled to typical values for massive molecular clouds in the GC region.
For nonrelativistic particles diffusing in the cloud with a diffusion coefficient typically with 0.3 < s_{D} < 0.5, the escape path length estimated from Eq. (A.4) depends only mildly on energy as Λ ∝ E^{ − sD / 2}. However, we have adopted here a simple slab model with an energyindependent escape path length in order to limit the number of free parameters as much as possible.
Fig. A.1
Schematic illustration of the cosmicray interaction model: fast particles produced in a lowdensity acceleration region can diffusively penetrate a denser cloud (if their kinetic energy is higher than a threshold energy E_{min}) and then produce nonthermal Xrays by atomic collisions. 

Open with DEXTER 
The process of CR penetration into molecular clouds is not well known (see, e.g., Gabici et al. 2007, and references therein). The theoretical predictions range from almostfree penetration (e.g. Cesarsky & Völk 1978) to exclusion of CRs of kinetic energies up to tens of GeV (e.g. Skilling & Strong 1976). Here, for simplicity, we assume that CRs can freely penetrate the clouds if their kinetic energy is higher than a threshold energy E_{min}, which is another free parameter of the model. We further consider the differential rate of primary CRs that penetrate the nonthermal Xray production region to be a power law in kinetic energy above E_{min}: (A.7)The model finally has four free parameters that can be studied from spectral fitting of Xray data (see Sect. 5): Λ, E_{min}, the powerlaw spectral index s, and the metallicity of the Xray emitting cloud, Z. The Xray spectral analysis also provides the CR spectrum normalization C_{i}, which allows one to estimate the power injected by the primary LECRs into the nonthermal Xray production region: (A.8)In the following, the integration in the above equation is limited to E_{max} = 1 GeV. Due to CR escape, the power continuously deposited by the fast particles inside the cloud should generally be lower than dW_{i} / dt.
Appendix B: Xrays from accelerated electron interactions
In the framework of the adopted steadystate, slab interaction model, the differential Xray production rate from collisions of accelerated electrons with the cloud constituents can be written as (B.1)where a_{j} is the abundance of element j relative to H in the Xray emitting cloud, (dσ_{ej} / dE_{X}) is the differential Xray production cross section for electron interaction with atoms j, and N_{e,p} and N_{e,s} are the differential equilibrium numbers of primary and secondary LECR electrons in the ambient medium, respectively.
B.1. Secondary electron production
Primary LECR electrons injected into an interstellar molecular cloud produce secondary electrons mainly from ionization of ambient H_{2} molecules and He atoms. The corresponding differential production rate of knockon electrons is given by (B.2)where (dσ_{H2} / dE_{s}) and (dσ_{He} / dE_{s}) are the H_{2} and He differential ionization cross sections for the production of a secondary electron of energy E_{s} by impact of a primary electron of energy E_{p}. The lower limit of the integral is 2E_{s}, because the primary electron is by convention the faster of the two electrons emerging from the collision. The maximum possible energy transfer is therefore E_{s} = 0.5(E_{p} − B_{j}) ≃ 0.5E_{p}, where B_{H2} = 15.43 eV and B_{He} = 24.59 eV are the electron binding energies of H_{2} and He, respectively. This convention is consistent with the definition of the stopping powers used throughout this paper (see Eq. (A.3)), which also pertain to the outgoing electron of higher energy.
The differential ionization cross sections are calculated from the relativistic binary encounter dipole (RBED) theory (Kim et al. 1994; 2000b). This successful model combines the binaryencounter theory for hard collisions with the dipole interaction of the Bethe theory for fast incident electrons. For the differential oscillator strengths, we use the analytic fits provided by Kim et al. (1994) for H_{2} and Kim et al. (2000a) for He. For the average orbital kinetic energy of the target electrons, we take U_{H2} = 15.98 eV and U_{He} = 39.51 eV.
By inserting Eq. (A.1) into Eq. (B.2) and using for the electron energy loss rate the expression given in Eq. (A.3), we see that the secondary electron production rate does not depend on the absolute density of H atoms in the ambient medium (n_{H}). This comment also applies to the Xray production rate, which only depends on the relative abundances a_{j} (see Eq. (B.1)). This important property of the adopted steadystate, slab model will allow us to estimate unambiguously the cosmicray power dW_{i} / dt (Eq. (A.8)) from the measured Xray flux.
Fig. B.1
Calculated differential equilibrium electron numbers (N_{e}; solid lines) for two differential injection rates of primary electrons (dQ_{e,p} / dt; dotted lines): a) s = 3, E_{min} = 10 keV; b) s = 2, E_{min} = 1 MeV. Also shown are the differential production rates of secondary, knockon electrons (dQ_{e,s} / dt; dashed lines). The H density in the nonthermal Xray production region, which intervenes in the calculation of N_{e,p} and N_{e,s}, is n_{H} = 10^{4} cm^{3} and the path length of the primary electrons in this region is Λ = 10^{24} cm^{2}. The calculations are normalized to a total power of 1 erg s^{1} injected by the primary LECR electrons in the Xray production region. 

Open with DEXTER 
Calculated differential production rates of primary and knockon electrons are shown in Fig. B.1. Also shown is the corresponding steadystate differential number of secondary electrons in the ambient medium, N_{e,s}. We calculated the latter from Eqs. (A.1) and (A.2), assuming the characteristic escape path length of the secondary particles to be Λ/2. Although this assumption is uncertain, it has no significant effect on the total Xray production.
We see in Fig. B.1 that the effect of H_{2} and He ionization on the electron energy distribution is to redistribute the total kinetic energy of the injected particles to a larger number of lowerenergy electrons. Thus, for hard enough primary electron spectrum (i.e. low s and high E_{min}, see Fig. B.1b), secondary electrons of energies E_{s} ≳ 10 keV could potentially make a significant contribution to the total nonthermal Xray emission. On the other hand, one can easily check that the successive production of knockon electrons by the secondary electrons themselves can be safely neglected for the Xray emission.
B.2. Xray continuum emission
The Xray continuum emission is due to the bremsstrahlung of both primary and secondary electrons. We take electron bremsstrahlung into account only in ambient H and He and calculate the differential cross sections from the work of Strong et al. (2000, Appendix A), which is largely based on Koch & Motz (1959). We use the scattering functions from Blumenthal & Gould (1970) to take into account the arbitrary screening of the H and He nuclei by the bound electrons.
B.3. Xray line emission
Fig. B.2
Calculated Xray emission produced by LECR electrons with the source spectra shown in Fig. B.1 interacting in a gas cloud of solar metallicity. PEB: primary electron bremsstrahlung; SEB: secondary electron bremsstrahlung. Photoelectric absorption is taken into account with a H column density of 10^{22} cm^{2}. 

Open with DEXTER 
The Xray line emission results from the filling of innershell vacancies produced by fast electrons in ambient atoms. We consider the Kα and Kβ lines (2p → 1s and 3p → 1s transitions in the Siegbahn notation) from ambient C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni. The corresponding cross sections can be written as (B.3)where E_{Kk} is the energy of line Kk (Kα or Kβ), δ(E_{X} − E_{Kk}) is Dirac’s delta function, the total cross section for the Kshell ionization of atom j by an electron of energy E, and the Kk fluorescence yield for atom j (Kaastra & Mewe 1993). Note that for element j with atomic number ≤ 12 (i.e. Mg), since these atoms do not have 3p electrons in their ground level.
For the Kshell ionization cross sections, we adopted the semiempirical formula of Quarles (1976), which agrees well with the RBED cross sections for Ni and lighter elements (see Kim et al. 2000b) and is simpler to use. We checked that the Quarles’s formula correctly reproduces the data compiled in Long et al. (1990), in particular at relativistic energies.
The width of the Xray lines produced by electron impact can be estimated from the sum of the natural widths of the atomic levels involved in the transition. Indeed, broadening effects caused by multiple simultaneous ionizations can be safely neglected for LECR electrons. Thus, the Kα_{1} and Kα_{2} components of the Fe Kα line have experimental full widths at halfmaximum (FWHM) of only 2.5 and 3.2 eV, respectively (Salem & Lee 1976). However, the energy separation of the two finestructure components is 13 eV, which is much less than the energy resolution at 6.4 keV of the Xray cameras aboard XMMNewton and Chandra, but larger than the expected resolution of the ASTROH Xray Calorimeter Spectrometer (7 eV FWHM; Takahashi et al. 2010). Here, we neglect the finestructure splitting of the K lines and for simplicity adopt the same width for all the lines: ΔE_{X} = 10 eV.
Figure B.2 shows calculated nonthermal Xray spectra (L_{X} = E_{X} × dQ_{X} / dt) produced by LECR electrons injected with the differential rates shown in Fig. B.1 into a cloud of solar metallicity. We took the photoelectric absorption of Xrays into account using a H column density N_{H} = 10^{22} cm^{2} and the cross sections of Morrison & McCammon (1983). We see in Fig. B.2 that the most prominent line is that of Fe at 6.40 keV. This is because this element has the highest product of Kα fluorescence yield (, Kaastra & Mewe 1993) and cosmic abundance. The EW of the Fe Kα line is equal to 293 and 394 eV in the spectra shown in panels a and b, respectively. The second strongest line in these spectra is the Si Kα line at 1.74 keV; its EW is equal to 80 and 90 eV in panels a and b, respectively. We also see in this figure that (i) the shape of the continuum emission reflects the hardness of the primary electron injection spectrum; and (ii) the total Xray emission is dominated by the contribution of the primary electrons. The emission from the secondary electrons is negligible in panel a and accounts for 10–20% of the total emission below 10 keV in panel b.
Appendix C: Xrays from accelerated ion interactions
The differential Xray production rate from accelerated ion interactions can be written with a slight modification of Eq. (B.1), as follows: (C.1)where the index i runs over the constituents of the nonthermal ion population. The first term in the integral represents the Xray production by the primary LECR ions and the second term the contribution of the secondary electrons. As a starting point, we assume in the present work that the LECR ion population is mainly composed of protons and α particles and that the contributions of accelerated metals to the total Xray emission can be neglected. We therefore do not consider the broad Xray line emission that can arise from atomic transitions in fast C and heavier species following electron captures and excitations (Tatischeff et al. 1998), except in Sect. 7.2. However, for typical compositions of accelerated cosmic particles, the fast metals significantly contribute neither to the production of the Xray lines from the ambient atoms nor to the bremsstrahlung continuum radiation (see Tatischeff et al. 1998). We further assume that the accelerated protons and α particles are in solar proportion, that is, C_{α} / C_{p} = a_{He} (see Eq. (A.7)).
In the calculations of the equilibrium spectra (N_{p} and N_{α}), we neglect the nuclear destruction and catastrophic energy loss (e.g. interaction involving pion production) of the fast ions in the cloud. Indeed these processes are not important in comparison with the ionization losses below ~300 MeV nucleon^{1} kinetic energy (see, e.g., Schlickeiser 2002) and most of the Xray emission below 10 keV, which is the prime focus of the present work, is produced by ions in this low energy range (see Fig. 5b).
C.1. Secondary electron production
Fig. C.1
Calculated differential equilibrium numbers of fast particles (N_{i}; solid lines) for the differential injection rate of primary protons given by s = 2 and E_{min} = 100 keV (dQ_{p} / dt; dotted lines). Also shown are the differential production rates of secondary knockon electrons (dQ_{e,s} / dt; dashed lines). a) Λ = 10^{21} cm^{2}; b) Λ = 10^{24} cm^{2}. The H density in the nonthermal Xray production region is n_{H} = 10^{4} cm^{3}. The calculations are normalized to a total power of 1 erg s^{1} injected by the primary LECR protons in this region. 

Open with DEXTER 
We calculate the production of secondary electrons associated to the ionization of ambient H_{2} molecules and He atoms. The corresponding differential ionization cross sections are obtained as in Tatischeff et al. (1998) from the work of Chu et al. (1981). We neglect the production of secondary electrons and positrons that follows the production of charged pions in hadronic collisions. In fact, the corresponding electron and positron source functions can dominate the one of knockon electrons only at energies >10 MeV (Schlickeiser 2002), and these highenergy leptons are not important for the production of Xrays <10 keV (see Fig. 5a).
Differential production rates of knockon electrons are shown in Fig. C.1, together with the corresponding equilibrium spectra of primary protons and secondary electrons. This figure illustrates the effects of changing the CR escape path length from Λ = 10^{21} cm^{2} (panel a) to 10^{24} cm^{2} (panel b). In the first case, protons of energies up to 1.4 MeV are stopped in the cloud, whereas in the second case the transition energy between proton stopping and escape is at 71 MeV. We see that above this transition energy the equilibrium spectrum has a similar slope than the source spectrum, whereas at lower energies the equilibrium proton distribution is harder due to the ionization losses. We can anticipate that the total Xray production rate will be much higher for the case Λ = 10^{24} cm^{2}, as a result of the higher proton number at equilibrium above a few MeV.
C.2. Xray continuum emission
The Xray continuum emission is due to inverse bremsstrahlung from the fast ions (the radiation of a single photon in the collision of a highspeed ion with an electron effectively at rest) and classical bremsstrahlung from the secondary knockon electrons. In the nonrelativistic domain, the bremsstrahlung produced by a proton of kinetic energy E in a collision with a H atom at rest has the same cross section as that of an electron of kinetic energy (m_{e} / m_{p})E in a collision with a stationary proton (m_{e} and m_{p} are the electron and proton masses, respectively). We calculate this cross section as in Sect. 3.2, but without taking the screening of the H nucleus by the bound electron into account. The cross section for interaction of a proton with a H atom is then multiplied by (1 + 2a_{He}) to take the ambient He into account. For α particles, we replace the proton energy E by the energy per nucleon of the projectile and multiply the proton cross section by 4 to account for the nuclear charge dependence of the bremsstrahlung cross section.
In the relativistic case, the cross section for proton inverse bremsstrahlung is different from the one for classical electron bremsstrahlung, owing to the appearance of angular and energy abberations in the transformation between the two rest frames of the interacting particles (Haug 2003). We checked that these effects can be neglected in good approximation in the present work.
In Fig. C.2 we show two Xray spectra corresponding to the particle equilibrium spectra presented in Fig. C.1. We see that the continuum emission is dominated by inverse bremsstrahlung, which is a general rule independent of the model parameters (see Tatischeff et al. 1998). We also see that, as expected, the Xray production rate is much higher for Λ = 10^{24} cm^{2} than for Λ = 10^{21} cm^{2}, the difference being a factor of 22, 337 and 1054 at 1, 10, and 100 keV, respectively.
Fig. C.2
Calculated Xray emission produced by LECR protons and αparticles interacting in a gas cloud of solar metallicity, for the differential injection rate of primary protons shown in Fig. 9. The contribution of accelerated αparticles is included as explained in the text, assuming in particular the solar abundance C_{α} / C_{p} = 0.0964. a) Λ = 10^{21} cm^{2}; b) Λ = 10^{24} cm^{2}. IB: inverse bremsstrahlung; SEB: secondary electron bremsstrahlung. Photoelectric absorption is taken into account with a H column density of 10^{22} cm^{2}. 

Open with DEXTER 
C.3. Xray line emission
For producing Xray lines from the ambient atoms, we take both the contribution from secondary electrons (see Eq. (B.3)) and that from primary ions into account. The cross sections for Kshell ionization by proton and αparticle impacts are extracted from the data library implemented by Pia et al. (2009) in the Geant4 toolkit for the simulation of particle induced Xray emission (PIXE). We use the cross sections calculated in the ECPSSR theory with highvelocity corrections (Lapicki et al. 2008). These cross sections are more accurate for mildly relativistic projectiles than those previously employed by Tatischeff et al (1998).
Proton and αparticle collisions with target atoms do not lead to significant line broadening effects caused by multiple simultaneous ionizations. We thus adopt as before a width of 10 eV for all the lines (see Sect. 3.3). We note, however, that the Xray lines produced by collisions of ions heavier than ^{4}He can be
shifted by several tens of eV, significantly broadened and split up into several components (Garcia et al. 1973). For example, the Fe Kα line produced by impacts of O ions of 1.9 MeV nucleon^{1} is blueshifted by ~50 eV in comparison with the one produced by 5MeV proton impacts, and has a FWHM of ~100 eV (see Garcia et al. 1973, Fig. 3.55).
The most intense line produced by LECR protons and αparticles is also the neutral Fe Kα line at 6.40 keV (Fig. C.2). This line has an EW of 2.31 and 0.80 keV in the spectra shown in Figs. C.2a and b, respectively. The second strongest line in these spectra is the Si Kα line at 1.74 keV; its EW is equal to 309 eV in panel a and 152 eV in panel b.
© ESO, 2012