- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
ISRN Astronomy and Astrophysics
Volume 2011 (2011), Article ID 905015, 10 pages
On the Cosmic Ray-Induced Ionization Rate in Molecular Clouds
V. A. Ambartsumian Byurakan Astrophysical Observatory, Byurakan 0213, Armenia
Received 30 January 2011; Accepted 13 March 2011
Academic Editors: C. Meegan and T. Wilson
Copyright © 2011 Ararat G. Yeghikyan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The transformation of the energy dependence of the cosmic ray proton flux in the keV to GeV region is investigated theoretically when penetrating inside molecular clouds ( mag). The computations suggest that energy losses of the cosmic ray particles by interaction with the matter of the molecular cloud are principally caused by the inelastic (electronic) interaction potential; the transformed energy distribution of energetic protons is determined mainly by the column density of the absorbing medium. A cutoff of the cosmic ray spectrum inside clouds by their magnetic fields is also phenomenologically taken into account. This procedure allows a determination of environment-dependent ionization rates of molecular clouds. The theoretically predicted ionization rates are in good agreement with those derived from astronomical observations of absorption lines in the spectrum of the cloud connected with the Herbig Be star LkH 101.
The physical characteristics of interstellar clouds such as temperature, pressure, and radiation balance as well as energy deposition from cosmic ray (CR) particles and ultraviolet photons (UV) depend strongly on the state of hydrogen (atomic versus molecular, neutral versus ionized) . Interstellar clouds can be categorized thoroughly into three classes based on their size, gas densities, temperature, and ionization fraction [4–7]: (i) diffuse, (ii) translucent, and (iii) molecular clouds. Diffuse cold clouds are characterized by typical number densities of less than 102 cm−3 and temperatures of about 100 K. Those diffuse regions in which hydrogen is present mainly in its neutral, atomic form are called HI clouds . Interstellar regions in which the gas exists predominantly in molecular form, that is primarily molecular hydrogen (H2), are referred to as molecular clouds. These areas hold average number densities of cm−3 and temperatures of about 10 K . As intermediate objects between the low-extinction diffuse clouds (Av < 1) and the optically thick molecular clouds (Av > 5), van Dishoeck and Black  introduced the class of translucent clouds with visual extinction spanning the range Av~1–5 mag. The interior of these clouds are entirely blocked from the external UV field at 13.6 eV (the ionization energy of atomic hydrogen) and beyond; recall that the column density along the line of sight, (in cm−2), is related to AV (in mag) by Av . As a matter of fact, each molecular cloud has a translucent layer. Here, the edges of each molecular cloud with average concentration of ~103 cm−3 and, in the ultimate sense, the edge of every clump with cm−3 must be considered as a translucent region . However, although UV does not penetrate into the deeper layers, energetic CR particles can induce an ionization and destruction of molecular hydrogen in the interior of molecular clouds (1); here, presents a proton of the CR with a kinetic energy before and after the ionization process. Further, the singly ionized hydrogen molecule () can react with H2 to form via(2) Therefore, the measurement of the column density can be used as a direct probe of the CR ionization rate () in dense clouds . This ionization rate presents a critical and fundamental quantity to physical and chemical models of molecular clouds, for which no direct measurement exists, and usually it is chosen as a free parameter . Furthermore, the cation is expected to be responsible for initiating ion-molecule reactions that lead, for instance, to the formation of simple hydrides (water, ammonia, methane) and possibly to higher order, saturated hydrocarbons ; unsaturated hydrocarbons and their radicals, on the other hand, can be synthesized via rapid neutral-neutral reactions . Prior to the direct identification of in infrared absorption against continuum radiation emanating from background protostars (W33A and GL2136) , molecular ions like HCO+ and N2H+ have been utilized to estimate the ionization rate in molecular clouds . On the other hand, the number density of in the steady-state, , is directly related to the CR ionization rate, (in s−1), via (3), with the number densities of molecular hydrogen, , and of carbon monoxide,  Here, it is assumed that is produced solely via (1-2) followed by reaction with any gas-phase constituent of the cloud, , with the exceptions of He, O2, and Holding fractional abundances between 10−4 and 10−5, CO is the dominant reactant in the molecular cloud; the ion-molecule reaction rate of CO with has no entrance barrier and is very fast with gas-kinetics limits ( cm3s−1 . Considering that the ratio of the number densities, , is approximately constant in dark clouds with [17, 18], Equation (3) reduces to (5). Assuming a constant value of s−1, Geballe  estimated number densities of () cm−3 However, it is clear that can depend on the cloud’s column density, . According to observations, the ionization rate varies between 10−18 and 10−16 s−1 in dense and diffuse clouds, respectively . Comparing chemical models of molecular clouds with actual astronomical observations, Caselli  derived an empirical relationship between the number density dependent ionization rate and the number density of H2 with an empirical parameter (6) Because the cross-section, (cm2), of the ionization of H2 via reaction (1) is well known , then it is possible to relate the CR flux, (cm−2s−1eV−1), in an energy interval (eV) to the ionization fraction and the cross-section via (7)  For instance, considering , cm2 (Figure 1, ), this yields . Hence, a detailed, theoretical determination of a function demands a thorough incorporation of self-consistent cloud models  to include depth-dependent ionization rates.
In this paper, I present a numerical approach to compute , in a molecular cloud provided that the energetic spectral flux of the CR particles outside of the given cloud and a radial distribution of the number density of H2 inside the cloud are known. It was shown long ago that molecular clouds can magnetically reflect low-energy CR particles by a screening mechanism involving the generation of Alfven waves [21, 22]. This means that the CR fluxes may be cut off at low energies and one should be able to choose a corresponding limiting value of energy (Section 2.1). Further, the column density of , , can be calculated as the obvious consequence of 5, via (8) and can then be compared with observational column densities of . This enables one to derive the depth-dependent ionization rate where is the cloud size
To compute the environment-dependent ionization rates, it is crucial to extract first the energy-dependent fluxes of the galactic CR particles outside the molecular clouds taking into account a cutting off of the fluxes by magnetic fields (Section 2.1). We will investigate then the energy loss processes of these charged particles upon their penetration inside the molecular cloud structure (Section 2.2) and calculate the change in flux of the CR particles inside the molecular cloud (Section 2.3). There were many earlier works devoted to such calculations (see [23–26] and references therein) where only diffuse HI clouds and protostars were considered. These authors used galactic CR fluxes that gave overestimated ionization rates at clouds edges, and clouds magnetic fields have not been taken into account. The aim of this work is to show that the ionization rate inside contracting molecular clouds is strongly dependent on environment and the limiting value of CR energy. The ionization rates in appropriate astrochemical models should be chosen taking into account both CR particles energy loss processes and the cutting off their spectral fluxes by magnetic fields of the clouds.
2. Penetration of Cosmic Ray Particles Inside the Clouds
2.1. The Non-Attenuated Cosmic Ray Spectrum
CR particles—predominantly protons (90%) and α-particles (10%)—are accelerated by magnetic fields in hydrodynamically active regions in which the kinetic energy of ejected stellar mass motions can be transferred magnetically into the energy of the particle, thus reaching up to the TeV region . In the diffuse interstellar medium (ISM), these energetic particles can propagate almost freely without losing a significant fraction of their kinetic energy to the matter in the diffuse cloud . Only protons will be considered here. Very recently, Cooper et al.  compiled observational and model data of the CR proton flux in the very local ISM (Figure 2), which is the power-law extension of the CR spectrum to an intersection with the high-energy side of the thermal plasma distribution. It should be noted that the non attenuated particle flux in  is different from the usually used distribution  (Figure 2). The total proton flux in the energy range 1 eV−10 GeV calculated from  is particlescm−2s−1 with an energy deposition of eV cm−2s−1, while on the basis of data from  one got 9.6 particles cm−2s−1 and eV cm−2s−1, respectively. In our model, we assume that the CR particle flux in the solar neighborhood is identical to that of in the ISM.
Most importantly, we have to consider that even the non attenuated spectrum of the CR particles can be influenced by Alfven waves (see [21, 22] and references therein). Only particles with energies greater than a minimum energy of about 100 MeV will enter the cloud. A more refined investigation depicted that the value of ( MeV) is given by (9) where (in cm−3), (in pc), and (in μG) are normalization parameters of the cloud [21, 22]. Here , , and are a number density, the radius, and the magnetic field strength, respectively. Observations reported measurements of molecular cloud magnetic field strengths ranging from 0.1 μG to about 6 μG . The typical value of the field strength in the ISM is 4-5 μG ; thus, varies between about 10 to 100 MeV for typical molecular clouds It is obvious that in order to estimate the environment-dependent ionization rates in molecular clouds systematically, a parameter study is crucial. Here, we compute the ionization rates for two flux distributions of the CR ( versus ); the influence of each of the distributions with the minimum energy induced by the magnetic shielding is investigated for MeV, 50 MeV, and 100 MeV. Finally, these models are applied at a typical cloud scenario: a molecular cloud with H2 number densities of cm−3 and an average size of 10 pc.
2.2. Interaction of Energetic Particles with Matter
The interaction and hence energy loss processes of each charged particle from the cosmic radiation field with the matter of the molecular cloud depend primarily on the kinetic energy of the particles [6, 28] and the physical characteristics of the ISM (density, spatial size, magnetic field). The elementary composition of CR in the keV to GeV energy range for the most important elements is similar to the cosmic elementary abundance. Here, protons (p) and α-particles dominate with an α/p ratio of less than 0.1; the abundance of heavier ions accounts for less than a few percent .
When an energetic particle (projectile) progresses through a molecular cloud, it can interact by two principal energy loss mechanisms with the surrounding matter . These are: (i) an interaction of the projectile with the electronic potential (electrons) of the atoms/molecules and (ii) with the nuclei of the atoms/molecules in the interstellar cloud (nuclear interaction/collision; elastic interaction). Inelastic energy loss processes change the kinetic energy of the projectile without modifying noticeably the propagation direction of CR particle. On the other hand, each elastic collision with an atomic nucleus changes the direction of the trajectory of the projectile . These energy loss processes—the dominance of inelastic versus elastic energy loss processes—depend strongly on the kinetic energy of the projectile. A common feature for each CR particle is that the nuclear energy loss, also called nuclear stopping power (), with a maximum at a relatively low energy of the order of 1 keVamu−1; as the kinetic energy increases, the nuclear stopping of the ion decreases. This means that the nuclear stopping power is important only for low ion holding very low velocities. For projectiles with kinetic energies larger than 1 keVamu−1, the energy loss processes are dictated by electronic interaction (electronic stopping power, ). The latter leads predominantly to electronic excitation and ionization of the atoms/molecules in the molecular cloud (molecules can also be excited rotationally and vibrationally). It is important to stress that at the beginning of the energy transfer processes from the implant to the surrounding matter at high energies, the ion decelerates mainly by electronic stopping at high energies, therefore almost propagating in a straight path. When the ion has slowed down sufficiently to the keVamu−1 region, the nuclear interaction potential becomes increasingly important. Therefore, to quantify the energy loss of an ion while crossing a molecular cloud, one has to investigate the energy-dependence of the nuclear, (), and electronic, (), stopping power. The total stopping power is simply the sum of both contributions (10). Here, we utilize the SRIM (Stopping and Range of Ions in Matter) code which is based on the modified Bethe-Bloch approximation  and is more convenient to use. Note that at kinetic energies close to 1 GeV, the energy loss of each charged particle is also determined by Bremsstrahlung, Cherenkov radiation, and nuclear reactions . Here, we neglect these processes since the differential flux of the CR particles at these energies is at least 4 orders of magnitude lower than, for instance, in the MeV region 
2.3. Energy Loss and Ranges of Charged Particles in Molecular Clouds
It is possible to quantify now the energy loss and the penetration (range) of charged particles inside molecular clouds. To describe a particle penetration depth into a cloud, the projected mean range, , is often used in collision dynamics. Compared to the actual path length, , of a projectile, the projected range, that is, the straight-line distance from implantation origin to the end point of the trajectory projected in the initial velocity vector, is always smaller. Here, significant angular deflections of the projectile via nuclear interaction happen to be eminent at energies of less than 1 keVamu−1. For example, considering 1 MeV protons, we can make use of (e.g., ) a crucial relationship between the mean free path, , that is, the average distance between two collisions, of a proton from the CR radiation field, the ionization cross-section of molecular hydrogen (1 MeV), and the number density of H2 (11). For cm2  and cm−3, a 1 MeV proton has a free mean path of about 31013 cm We computed the projected ranges and electronic, nuclear, and total stopping powers of protons in molecular hydrogen gas with the SRIM code for a density of gcm−3 . SRIM does not allow a computation of the stopping powers at the low density of interstellar clouds directly. Since the range is inversely proportional to the density of a medium, we can utilize (11) to scale the range of a proton of an energy in a medium of a densityto compute the range of the same proton with an identical energy in a medium of a density via (12). Applying this procedure for molecular clouds with number density of 103 cm−3, we find that . Since H2 is the dominating component of these environments, our scaling procedure is well justified. For instance, applying this procedure to 1 MeV protons gives a computed range in molecular clouds with number densities of 103 cm−3 of cm. Thus, one can conclude from such computations that (i) CR protons with energies cannot penetrate inside molecular clouds with ~103 cm−3 deeper than 0.1 pc, (ii) for energies , the CR particles travel essentially straight paths through the cloud since deflections from its trajectory via nuclear interactions are negligible, and (iii) the non attenuated flux of the CR particles is changed due to the energy loss; this in turn slows the particles down and shifts the flux distribution from high values to lower In a one-dimensional steady-state approximation, neglecting diffusion and internal sources, the proton flux inside a cloud after shielding by a homogeneous layer, with a depth and a density can be numerically computed from the stopping power for H2 (Figure 3) and the initial, non attenuated proton flux, (Figure 2). Accounting for (10), the energy-dependent energy loss of a particle with a kinetic energy , , can be calculated via(13) Considering energy conservation, we have to bear in mind that the external energy flux —in units of eV cm−2 s−1—of the protons outside the molecular cloud must be equal to the energy flux from the decelerated protons inside the cloud, , plus the energy flux given to the target atoms, Here, the external energy flux can be computed as the product of the energy and the non attenuated flux (Figure 4) to Likewise, the internal energy flux presents the product of the energy and the internal flux —the quantity we want to compute—to To calculate the energy flux transferred from particles of an energy given to the target molecules (the H2 gas inside the cloud), we have to account for the following scenario. For particles at a given energy , the energy flux contributing from these particles is being reduced (they are actually slowed down from an energy to lower energies) via interaction with the target gas by (c.f. (13) for the definition of ); likewise, particles of an energy , which is higher than , can be decelerated to the energy ; this process actually increases the energy flux for a given energy by Thus, a particle with an energy transverses a column density and exits with an energy Combining these relationships with (14)–(16) and inserting these into (13) gives the method to calculate via (18) from which we can extract then One can now utilizes these attenuated fluxes to investigate quantitatively to what extent the CR protons ionize molecular hydrogen via (1). Accounting for the ionization cross-section (Figure 1) and the attenuated CR proton fluxes (Figure 4), the depth-dependent ionization rate, (s−1), can be derived via (19). Here, the integration limits range from (the lowest kinetic energy of the ions penetrating to the depth, ) to the maximum energy of 1010 eV (Figures 1 and 2).
Table 2 shows the results of our computations for six molecular cloud models; the parameters of these models are compiled in Table 1. The corresponding internal fluxes and depth-dependent ionization rates are depicted in Figure 4. Based on these data, we can derive the following principal conclusions. (1)The energy absorbed by the molecular hydrogen inside the molecular cloud at a distinct depth is very small ( eV cm−2s−1 at the most) compared to the non attenuated energy flux of eV cm−2s−1  and eV cm−2s−1 . This translates to an energy loss of the CR protons of less than 2% irrespective on the model; therefore, the charged particles transverse the cloud almost in straight trajectories nearly unperturbed. Further, as expected from the computed stopping power (Figure 3), electronic energy loss presents the dominating energy transfer mode and is three orders of magnitude more important than the nuclear energy transfer. Summarized, the energy loss of the charged particles is irrespective on the choice of the non attenuated flux, on the number density of molecular hydrogen (at least in the 103 to 104 cm−3 regime), and on the selection of (10 MeV, 50 MeV, and 100 MeV). (2)Even though in , total non attenuated proton flux of cm−2 s−1 is four orders of magnitude larger than the value of  (9.3 cm−2 s−1), the majority of flux distribution in  are low-energy protons. These species are magnetically deflected and hence do not penetrate inside the cloud. This yields to a very narrow range of proton fluxes of—depending on the depth-6 protons cm−2s−1 to 21 protons cm−2s−1  and 6 protons cm−2s−1 to 9 protons cm−2s−1 . (3)Our computed ionization rates range between s−1 and s−1. It is interesting comparing these data with ionization rates claimed by observations, that is, 10−18 s−1–10−16 s−1. The present computations predict that—depending on the parameters utilized in models M1–M6—the ionization rate can be in certain areas of the molecular cloud higher than previously postulated. Likewise, our calculations suggest minimal ionization rates of s−1; this value is in strong contrast to estimated data of about 10−18 s−1. Also, the resulting ionization rates decreases slightly toward the center of the molecular cloud as predicted by .
It should be stressed here that these results depend essentially on the cut-off procedure of the CR spectra at low energies (10 to 100 MeV); otherwise the contradiction between calculated and observed ionization rates may be even stronger. For example, the CR spectra when protons with energies 1 MeV ≤ ≤ 10 GeV penetrate inside clouds cause ionization rates close to 10−15 s−1 under typical clouds conditions (~103 cm−3) (Figure 5) or to 10−16 s−1 ( cm−3) Figure 6).
4. Astrophysical Implications
It is interesting to compare our calculations with those data derived from astronomical observations. We would like to stress that our computed values of the ionization rate are in good agreement with an (indirectly) observed data point. Brittain et al.  measured absorption lines in the spectrum of the Herbig Be star LkHα 101 suggesting a column density of cm−2 . A more detailed investigation of the data revealed that the absorption lines must be correlated with a dense cloud with number densities of ~104 cm−3, the length along the line of sight of ~0.4 pc, and an averaged kinetic temperature of ~30 K. To illustrate a validity of the calculated value of the fractional ionization (Figure 4), we can apply (3)–(5) to observations . Here, utilizing a ratio of the column density of CO and H2 in this object of , we derive and, hence, a theoretically predicted column density of cm−2; this compares very well with the observed value of . Moreover, the calculated value of the ionization rate throughout the cloud, s−1, while the observed value is of s−1. Here, it is also worth mentioning that the CR flux—assumed similar in the solar neighborhood and at the edge of the non attenuated cloud—is adequately transformed after shielding corrections inside the dense cloud toward LkHα 101. This model has been shown to correctly reproduce the astronomical observations of column densities of . However, one should point out that the calculated fluxes are sensitive to the density distribution in the cloud; in approximation used one assumed a homogeneously mixed cloud; non homogeneously mixed clouds should demand on more refined procedure.
Further studies have the goal to expand this model to investigate how CR particles heavier than protons influence the ionization rate . In addition, one will examine the effect of the chemical composition of the molecular cloud. Recall that in the present model it is assumed the cloud consists entirely of molecular hydrogen; therefore, the incorporation of atomic hydrogen, helium, and carbon monoxide is worth investigating. Finally, the sensitivity of the model to dust particles inside molecular clouds—although having very low number densities of about grains cm−3—will be probed. These extensions will also yield valuable information on the energy deposition inside grain particles which in turn can help to understand the energy dependent formation of molecules on icy grains deep inside molecular clouds . This would allow an objective, unbiased comparison of the energy deposition inside ice-coated grains by CR versus UV photons and would expand Shen et al.’s  preliminary study significantly.
In summary, it is shown that the distribution of the CR particle flux inside a molecular cloud can be calculated if the energy loss of energetic protons—the principal constituent of the cosmic radiation field—in molecular hydrogen is known. The latter can be computed utilizing the SRIM program package and rescaling the computed data to the actual number densities of the cold molecular clouds. Integrating over the energy-dependent ionization cross-section of molecular hydrogen, a simple numerical procedure finally yields the ionization rates of molecular clouds. A phenomenological cutting-off procedure of the CR spectra at low energies permits to bring together the calculated and observed ionization rates. Based on astronomically observed absorption lines, calculated values of the ionization rates are well within those derived from astronomical observations. These environment-dependent ionization rates are in turn critical parameters to establish updated physical and chemical models of molecular clouds.
Partially this work was supported by the Fulbright Scholar Program (USA). The author thanks Ralf Kaiser for helpful discussions. Useful comments from two anonymous referees helped to improve the presentation of this paper.
- L. Nagy and L. Végh, “Ionization of molecular hydrogen by proton impact. I. Single ionization,” Physical Review A, vol. 46, no. 1, pp. 284–289, 1992.
- J. F. Cooper, E. R. Christian, J. D. Richardson, and C. Wang, “Proton irradiation of Centaur, Kuiper Belt, and Oort Cloud objects at plasma to cosmic ray energy,” Earth, Moon and Planets, vol. 92, no. 1-4, pp. 261–277, 2003.
- W. R. Webber and S. M. Yushak, “A measurement of the energy spectra and relative abundance of the cosmic-ray H and He isotopes over a broad energy,” The Astrophysical Journal, vol. 275, pp. 391–404, 1983.
- M. A. Dopita and R. S. Sutherland, Astrophysics of the Diffuse Universe, Springer, New York, NY, USA, 2003.
- C. F. McKee, “Multiphase gas in quasar absorption systems,” in The Physics of the Interstellar Medium and Intergalactic Medium, A. Ferrara, C. McKee, C. Heiles, and P. Shapiro, Eds., vol. 80 of ASP Conference Series, pp. 468–477, ASP Press, San Francisco, Calif, USA, 1995.
- R. I. Kaiser, “Experimental investigation on the formation of carbon-bearing molecules in the interstellar medium via neutral-neutral reactions,” Chemical Reviews, vol. 102, no. 5, pp. 1309–1358, 2002.
- A. G. G. M. Tielens, The Physics and Chemistry of the Interstellar Medium, CUP, 2005.
- L. Blitz and J. P. Williams, “Molecular clouds,” in The Origins of Stars and Planetary Systems, C. Lada and N. D. Kilafis, Eds., pp. 3–28, Kluwer Academic Press, 1998.
- E. F. van Dishoeck and J. H. Black, “The photodissociation and chemistry of interstellar CO,” The Astrophysical Journal, vol. 334, pp. 771–802, 1988.
- J. E. Dyson and D. A. Williams, Physics of the Interstellar Medium, IoP, 1997.
- R. Stark and K. Reif, “Surface brightness measurements of extended galactic nebulae,” Publications Astronomical Society of Australia, vol. 15, no. 1, pp. 86–90, 1998.
- E. Herbst, “The astrochemistry of ,” Philosophical Transactions of the Royal Society A, vol. 358, no. 1774, pp. 2523–2534, 2000.
- P. Caselli, “Constraining chemical-physical properties of pre-stellar cores,” Astrophysics and Space Science, vol. 285, no. 3-4, pp. 619–631, 2003.
- P. Ehrenfreund and S. B. Charnley, “Organic molecules in the interstellar medium, comets, and meteorites: a voyage from dark clouds to the early Earth,” Annual Review of Astronomy and Astrophysics, vol. 38, no. 1, pp. 427–488, 2000.
- T. R. Geballe and T. Oka, “Detection of interstellar space,” Nature, vol. 384, no. 6607, pp. 334–335, 1996.
- P. Caselli, C. M. Walmsley, A. Zucconi, M. Tafalla, L. Dore, and P. C. Myers, “Molecular ions in L1544. II. The ionization degree,” The Astrophysical Journal, vol. 565, no. 1, pp. 344–358, 2002.
- T. R. Geballe, “ between the stars,” Philosophical Transactions of the Royal Society A, vol. 358, no. 1774, pp. 2503–2513, 2000.
- H. H. Lee, R. P. A. Bettens, and E. Herbst, “Fractional abundances of molecules in dense interstellar clouds: a compendium of recent model results,” Astronomy and Astrophysics Supplement Series, vol. 119, no. 1, pp. 111–114, 1996.
- P. Jenniskens, G. A. Baratta, A. Kouchi, M. S. de Groot, J. M. Greenberg, and G. Strazzulla, “Carbon dust formation on interstellar grains,” Astronomy and Astrophysics, vol. 273, p. 583, 1993.
- J. P. Williams, E. A. Bergin, P. Caselli, P. C. Myers, and R. Plume, “The ionization fraction in dense molecular gas. I. Low-mass cores,” The Astrophysical Journal, vol. 503, no. 2, pp. 689–699, 1998.
- J. Skilling and A. W. Strong, “Cosmic ray exclusion from dense molecular clouds,” Astronomy and Astrophysics, vol. 53, no. 2, pp. 253–258, 1976.
- C. J. Cesarsky and H. J. Volk, “Cosmic ray penetration into molecular clouds,” Astronomy and Astrophysics, vol. 70, p. 367, 2001.
- L. Spitzer Jr. and M. G. Tomasko, “Heating of H i regions by energetic particles,” The Astrophysical Journal, vol. 152, p. 971, 1968.
- P. M. Solomon and M. W. Werner, “Low-energy cosmic rays and the abundance of atomic hydrogen in dark clouds,” The Astrophysical Journal, vol. 165, p. 41, 1971.
- T. Cravens and A. Dalgarno, “Ionization, dissociation, and heating efficiencies of cosmic rays in a gas of molecular hydrogen,” The Astrophysical Journal, vol. 219, no. 15, pp. 750–752, 1978.
- T. Umebayashi and T. Nakano, “Fluxes of energetic particles and the ionization rate in very dense interstellar clouds,” Publications of the Astronomical Society of Japan, vol. 33, p. 617, 1981.
- R. L. Curran and A. Chrysostomou, “Magnetic fields in massive star-forming regions,” Monthly Notices of the Royal Astronomical Society, vol. 382, no. 2, pp. 699–716, 2007.
- E. G. Zweibel and C. Heiles, “Magnetic fields in galaxies and beyond,” Nature, vol. 385, no. 6612, pp. 131–136, 1997.
- J. P. Meyer, L. O. Drury, and D. C. Ellison, “Galactic cosmic rays from supernova remnants. I. A cosmic-ray composition controlled by volatility and mass-to-charge ratio,” The Astrophysical Journal, vol. 487, no. 1, pp. 182–196, 1997.
- J. F. Ziegler, J. P. Biersack, and U. Littmark, The Stopping and Range of Ions in Matter, Pergamon, 1985.
- J. F. Ziegler and J. P. Biersack, 2010, http://www.srim.org/.
- S. D. Brittain, T. Simon, C. Kulesa, and T. W. Rettig, “Interstellar line absorption toward LkHα 101,” The Astrophysical Journal, vol. 606, no. 2, pp. 911–916, 2004.
- C. J. Shen, J. M. Greenberg, W. A. Schutte, and E. F. van Dishoeck, “Cosmic ray induced explosive chemical desorption in dense clouds,” Astronomy and Astrophysics, vol. 415, no. 1, pp. 203–215, 2004.