Issue |
A&A
Volume 666, October 2022
|
|
---|---|---|
Article Number | A27 | |
Number of page(s) | 7 | |
Section | Astrophysical processes | |
DOI | https://doi.org/10.1051/0004-6361/202142551 | |
Published online | 30 September 2022 |
Fast methods for tracking grain coagulation and ionization
II. Extension to thermal ionization
1
Institut de Recherche en Astrophysique et Planétologie, Université Paul, Sabatier Toulouse 3, 118 Rte de Narbonne, 31062 Toulouse, France
2
Department of Astrophysics, American Museum of Natural History, Central Park West at 79th Street, New York, NY 10024, USA
e-mail: pierre.marchand.astr@gmail.com
3
Université Paris-Saclay, CNRS, Institut d’astrophysique spatiale, 91405 Orsay, France
4
Laboratoire Univers et Particules de Montpellier, Université de Montpellier, CNRS/IN2P3, CC 72, Place Eugène Bataillon, 34095 Montpellier Cedex 5, France
5
AIM, CEA, CNRS, Université Paris-Saclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 Gif-sur-Yvette, France
Received:
29
October
2021
Accepted:
18
February
2022
Thermal ionization is a critical process at temperatures T > 103 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include thermal ionization. The model can be used to inexpensively calculate the density of ions and electrons and the electric charge of each size of grains for an arbitrary size distribution. This tool should be particularly useful for the self-consistent calculation of nonideal MHD resistivities in multidimensional simulations, especially of protostellar collapse and protoplanetary disks.
Key words: dust, extinction / stars: formation / ISM: abundances
© P. Marchand et al. 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
In Marchand et al. (2021, hereafter Paper I), we presented a novel method to calculate the coagulation and ionization of grains at a low computational cost. The only ionization source in that model was cosmic rays, which are dominant for isolated dense cores in the interstellar medium. At high temperatures, however, thermal ionization plays a major role in determining the ionization equilibrium of the gas-grain mixture, with a significant impact on the nonideal magnetohydrodynamics (MHD) resistivities.
The first stage of the protostellar collapse is isothermal at ∼10 K until density reaches ∼10−13 g cm−3. At this point, the dust-gas mixture becomes opaque to its own thermal radiation and the temperature rises as a core forms and contracts slowly (Larson 1969). At 2000 K, the dissociation of H2 molecules absorbs energy, allowing a rapid second collapse that leads to the creation of the protostar when all H2 is depleted. The thermal ionization of hydrogen occurs during the second collapse, and all hydrogen becomes ionized early in the protostar’s life, provoking a drop in resistivities to virtually zero (Marchand et al. 2016). Accounting for the thermal ionization of hydrogen is therefore critical to accurately describe the transition between the first core and the protostar, and thus from the nonideal to ideal MHD regime.
Another example of how this could be applicable pertains to chondrule formation. Chondrules are molten grains found in meteorites, whose formation requires rapid heating to ∼2000 K (Ebel et al. 2012). While there is no consensus on this topic, it has been proposed that their creation takes place in magnetic current sheets in protoplanetary disks, which may reach temperatures > 1500 K (Joung et al. 2004; McNally et al. 2014). Those high temperatures would trigger the thermal ionization of K and Na, leading to a sharp decrease in resistivities. Subsequently, the downward gradient in resistivities may create an instability that would allow the magnetic field to pile up in the current sheet (Hubbard et al. 2012). That phenomenon may be the origin of thunderclaps and extremely localized heating of the grains and gas (McNally et al. 2013), which are necessary conditions for chrondrule formation. However, readers should refer to Desch & Turner (2015), who argue that insufficient alkali metals would evaporate from grains to allow for this instability to act.
In this paper, we focus on extending the ionization model of Paper I by including the thermal ionization of one gas-phase species. In Sect. 2 we analytically derive the grain charge (Sect. 2.1) and thermal ionization equilibrium (Sect. 2.2), whose numerical implementations are described in Sect. 3. In Sect. 4 we discuss applications, and in Sect. 5 we present our conclusions.
2. Analytical method
2.1. Grain charge
Let us consider the two ionic species i and s of number density ni and ns. Species i corresponds to all the ions that are exclusively ionized by cosmic rays in the same manner as in Paper I, for which we assume an average atomic mass μi = mi/mH. Species s, however, undergoes both cosmic-ray and thermal ionization. We define ns, 0 as the total abundance of species s (both neutral and charged), so that ns, 0 is an upper bound for ns. We also consider an arbitrary size distribution of dust grains.
Grain charges Ze, with e being the electron charge, fluctuate stochastically due to the collection of electrons and the recombination of ions on their surfaces. The grain charge equilibrium (Eq. (4.3) of Draine & Sutin 1987, and Eq. (24) of Paper I) can be written as
where
is the reduced temperature of a grain of radius ak and temperature T, while kB is the Boltzmann constant. Furthermore, f(Z, τk) is the distribution function of grain charges for a grain of reduced temperature τk, and Ji, Js, and Je are the fluxes of species i, s, and electrons onto the grains, respectively. The fluxes are (Draine & Sutin 1987)
with nj being the abundance of species j (for j = i, s, or e), sj being the sticking probability of species j on grains, mj and qj being the mass and charge of species j, and vj = (8kBT/πmj)1/2 being the thermal speed of j. The polarization factor of grains for species j is , which depends on the relative signs of Z and qj.
Similarly to the case without thermal ionization, presented in Paper I, we need to solve Eq. (1) for small grains and low temperatures, for which τk ≪ 1 and f(Z) is only significant for Z = −1, 0 and 1; that is, these grains hold a maximum of one charge. Therefore
Equation (1) can be rewritten for Z = 0 and Z = −1 as follows:
As in Paper I, we have
where appears in Je(0, τk), Ji(0, τk), and Js(0, τk), while
appears in Je(1, τk), Ji(−1, τk), and Js(−1, τk). We can then solve the equation system (4)–(6). With si = ss = 1, we obtain
where ,
, and
The main difference with Paper I is the appearance of the term qis. While ni + ns is the total abundance of ions, ni + qisns is an effective abundance that reflects the relative flux of ions onto grains. The average charge of grains
and the grain-ion recombination enhancement factor
thus are given by the same expressions as in Paper I,
where we neglected the recombination of ions on positively charged grains ().
For the larger grains (τk ≫ 1), the same kind of change needs to be made to the Spitzer equation (Spitzer 1949; Draine & Sutin 1987) that governs the grain’s electric potential ψ. For ψ < 0, eψ represents the repulsion of the flux of electrons by the negatively charged grains, while 1 − ψ characterizes the attraction of the flux of ions. The flux equilibrium can thus be written as
Introducing the same notations as above, we can write
which is the same equation as the one-ion model of Paper I (Eq. (34)) with the modified expression for ϵ. The average charge of large grains and the grain-ion recombination enhancement factor yield the same expressions as in Paper I,
The average charge and recombination enhancement factor for a mix of small and large grains is assumed to be the sum of the contributions from both Eqs. (15),(19) and Eqs. (16),(20) (Draine & Sutin 1987).
2.2. Ionization equilibrium
We always assume charge neutrality,
This allows us to find the ionization equilibrium for species i. In Paper I, we considered the balance between the creation of species i by cosmic-ray ionization, and the destruction of species i by recombination with electrons and with grains. In this two-ion model, we also need to consider the charge exchange reactions between species i and s. The one-ion model hides and summarizes all the charge transfer reactions between gas-phase species in the choice of μi. Here, we have to explicitly account for the creation of species i by the destruction of species s, and vice versa. The ionization equilibrium is then
where ζ is the cosmic-ray ionization rate, ks,i and ki,s are the chemical reaction rates of species s → i and i → s, respectively, and we consider the recombination rate of ions i with electrons to be cm3 s−1, based on the recombination rate of HCO+ taken from the UMIST database (McElroy et al. 2013). The terms of the form ka, bnanc are the transformation of species c to species b, through the chemical reaction with species a. Hence the term (nH − ns, 0 − ni) represents the total abundance of neutral species that can be ionized into ion i, as (ns, 0 − ns) is the total abundance of neutral species s that can be ionized to ion s.
The ionization equilibrium for species s is similar, with the addition of a thermal ionization term (Pneuman & Mitchell 1965, also see Sect. 4)
with the values of β and T0 depending on the species. Table 1 summarizes the values of constants specific to species s for the cases of sodium, potassium, and hydrogen. The ionization equilibrium equation is then
Constants depending on the choice of species s for sodium, potassium, and hydrogen.
Those equations are valid if the Saha equation is valid as well, meaning that there should be a large number of particles within a Debye length of each other. The validity condition is then
with
3. Numerical implementation and tests
Equations (18), (21), (22), and (24) need to be solved for ψ, ϵ, ni, and ns. In this section, we discuss the solution for this system with four equations and four unknowns. The system could be reduced to three equations, as Eq. (18) is an explicit expression of ψ as a function of ϵ. That would, however, significantly increase the analytical and numerical complexity of the calculation, and it is unclear whether this would lead to better performances or not.
3.1. Numerical convergence
Although the system of equations is valid for a wide range of physically valid parameters, we need to be cautious to ensure numerical convergence toward the solution, especially at high density and temperature. At high density, ψ converges toward zero by a negative value and becomes very small in an absolute value. At high temperature, ns overwhelmingly dominates ni due to the thermal ionization. Therefore, the numerical implementation has to be robust for the cases |ψ|≪1 and ns/ni ≫ 1.
For this purpose, the four equations must be normalized so that they can be written in the form 1 + x = 0 to avoid sums of very large or very small numbers. It is therefore necessary to include species s in the normalization of the equations to avoid convergence issues at large temperatures when ns grows much larger than ni. We therefore normalized Eq. (21) by the total number of ions ni + ns, we used both the cosmic-ray ionization rate and the chemical reaction rate s → i to normalize Eq. (22), and we included the thermal ionization term in the normalization of Eq. (24).
Another issue arises from the average grain charge (15). Before the thermal ionization starts to be relevant and ni ≫ ns, ϵ converges toward 1/Θ as density increases (see Fig. 2 of Paper I). When the difference between ϵ and 1/Θ becomes close to machine precision, the term 1 − (ϵΘ)2 reaches a lower bound1 which prevents the convergence of the charge neutrality (Eq. (21)), as the grain charge fails to decrease. A solution to avoid this issue is to replace ϵ by ϵ′+ϵ0, with ϵ0 = 1/Θ. This substitution has to be made in all the equations. The new variable to find is ϵ′ which converges toward zero instead of 1/Θ. The term 1 − (ϵΘ)2 in Eq. (15) is then mathematically equal to and can be replaced by −2ϵ′Θ − (ϵ′Θ)2, which avoids the lower-bound issue.
3.2. The normalized system of equations
The numerical solution of this system requires its normalization, as discussed in Sect. 3.1. We define the function F(ψ, ϵ′,ni, ns) = (f1, f2, f3, f4), with
where
and
Here, vi, s and μi, s stand for either vi and μi or vs and μs.
3.3. Solving the system
Similarly to Paper I, we used a Newton–Raphson method to solve the equation system. Let X = (ψ, ϵ, ni, ns). Starting from an educated guess X0, we iterated
until ||F(Xn)||< δ ≪ 1. The matrix 𝕁 is the Jacobian of the system defined by 𝕁i, j = ∂fi/∂Xj. The full analytic components of the Jacobian matrix are given in Appendix A.
For reliable convergence, this iterative solution for the system of equations is best started from as close an estimate as possible. For low density (nH < 107 cm−3) and low temperature (typically 10 K), a good starting point is
where ψ0 is the solution of
For larger densities and temperatures, we recommend solving the system for a gradual increase in those quantities, using the previous solution as a first estimate. In particular, ns increases very quickly with density and temperature once the thermal ionization starts, and large leaps may lead to convergence failure.
3.4. Tests
We solved the normalized system of Eqs. (27)–(30). For testing purposes, we used typical parameters of star-forming environments to compare with the existing literature, but the reader should keep in mind that a wide range of physically sound parameters is possible. The density spans nH = 104–1025 cm−3, starting from the lowest density and increasing nH gradually. We assumed the same barotropic equation of state as in Marchand et al. (2016) to emulate the rise in temperature during a protostellar collapse
with T0 = 10 K, n1 = 1011 cm−3, n2 = 1016 cm−3, and n3 = 1021 cm−3. We assumed a nonevolving Mathis, Rumpl, Nordsieck (MRN) grain-size distribution (Mathis et al. 1977), with a slope of −3.5 between minimum grain size amin = 5 nm and maximum grain size amax = 250 nm, sampled by 26 bins. The grain bulk density is ρ = 2.9 g cm−3 and the dust-to-gas mass ratio is 1%, which is typical of the insterstellar medium (Bohlin et al. 1978). We set ζ = 5 × 10−17 s−1 (Padovani et al. 2013), se = 0.5 (Umebayashi & Nakano 1990), and μi = 25 (close to the molecular mass of Mg, Fe, or HCO+, Marchand et al. 2016).
In this test, the temperature exceeds several 103 K, above which all grains should be quickly destroyed by evaporation or sputtering (Lenzuni et al. 1995). This is, however, not an issue since at such high temperatures, the contribution of the (computed) charge of grains is negligible compared to that from ions and electrons. It is also possible to combine our method with any grain destruction model. Marchand et al. (2016) assumed that the grain evaporation would occur between ∼800 K and ∼1600 K in the density range nH ≈ 1016 − 1019 cm−3, which coincides with the beginning of the thermal ionization of K and Na. The species s is assumed to be K, Na, or H, using the values of Table 1. Desch & Turner (2015) show that K and Na may originally be confined to grains and have to be evaporated before being available for thermal ionization. Our method is compatible with models accounting for that process since ns, 0 can be freely modified at any time. For simplicity, however, we assume here that K and Na are already present in the gas phase in quantities given by ns, 0 in the table. Here, we present the evolution of ni/nH, ns/nH, ne/nH and the average grain charge for several bins. The results are displayed in Figs. 1–3, respectively.
![]() |
Fig. 1. Test of our method for protostellar collapse conditions. Evolution of the average charge of grains (left axis) for several grain sizes (color lines), and the fractional abundance (right axis) of cosmic-ray ionized ions ni/nH (solid line), thermally ionized K+ ions ns/nH (dotted line), and electrons ne/nH (dashed line). |
![]() |
Fig. 2. Same as Fig. 1, but for Na. The only ions present are the cosmic-ray ionized ones i and Na (ion s), which is both cosmic-ray and thermally ionized. |
![]() |
Fig. 3. Same as Fig. 1, but for H. The only ions presents are the cosmic-ray ionized ones i and H (ions s), which is both cosmic-ray and thermally ionized. |
These figures can be compared to Fig. 7 of Marchand et al. (2016). Although the abundances of Na+ and K+ seem overestimated compared to ne, and that of H+ seems underestimated (with a negligible impact on the resistivities), our model reproduces the evolution of abundances of this more detailed calculation at the key points fairly well: K and Na start their ionization around nH = 1018 cm−3 and T = 1600 K, and they saturate around nH = 1022 cm−3 and T = 1.5 × 104 K at the maximum fractional abundance of their respective species. Furethermore, H has a similar behavior, starting its thermal ionization at nH = 1020 cm−3 and T = 2650 K, turning virtually all neutrals into ions by nH = 1023 cm−3 and T = 7 × 104 K. The main difference with Fig. 7 of Marchand et al. (2016) is the earlier rise of electron density at nH ≈ 1016 cm−3 due to the thermionic emission of grains included in the complete chemical calculation. This emission is associated with a drop in neutral grain density and rise in positively and negatively charged grains.
In all three cases, the large input of electrons into the gas significantly increases the grain charges as well. We note that at the hydrogen ionization fraction nH+/nH > 0.1 reached at nH > 1023 cm−3, the Saha equation is no longer valid, so our model becomes imprecise. However, the resistivities are so low in this regime that ideal MHD is a valid approximation; thus this is never an issue in practice.
Figures 4–6 display the associated nonideal MHD resistivities using the formulae provided by Marchand et al. (2016). For display purposes only, we prescribe a magnetic field (Li et al. 2011)
![]() |
Fig. 4. Evolution of the Ohmic (yellow), Hall (blue), and ambipolar (dark red) resistivities for the thermal ionization of K. The dashed line represents the Hall resistivity in negative values. The thin dotted lines are the resistivities in the absence of thermal ionization. |
which corresponds to the critical magnetic field strength of a spherical cloud. The magnetic field strength and the resistivities are therefore overestimated compared to protostellar collapse simulations with nonideal MHD for nH ≳ 1012 cm−3.
The nonideal MHD terms significantly affect the MHD evolution of protostellar collapse and protoplanetary disks for resistivities larger than ≈1018 cm2 s−1. Nonideal MHD terms are then important at all densities before thermal ionization starts (except the Ohmic diffusion at a low density). At this point, the abundance of charged species in the gas significantly increases, leading to a sharp decrease in all resistivities. The overall behavior is consistent with previous works (Kunz & Mouschovias 2010; Marchand et al. 2016; Wurster et al. 2016; Koga et al. 2019).
4. Discussion
The method presented here is applicable in a wide variety of environments, and it is particularly suited for modeling protostar formation and protoplanetary disks. At later stages of the star formation process, or in the presence of nearby massive stars, photoionization by UV or X-rays may become relevant (Getman & Feigelson 2021). In this case, their ionization rate can simply be added to the cosmic-ray ionization rate ζ in the system of Eqs. (27)–(30).
We describe our method to calculate the resistivities as fast in comparison to solving a full chemical network. We have implemented the algorithm in the 3D MHD RAMSES code (Teyssier 2002). The code previously calculated the resistivities by interpolating on the precalculated table of Marchand et al. (2016). Without the Hall effect, our thermal ionization algorithm is faster than reading the chemical table, as the calculation needs to be performed only once per cell per time-step. This is different with the Hall effect, which requires, in addition, the self-consistent calculation of resistivities on cell edges (Marchand et al. 2018). In this case, the code runs at similar speeds for both methods (it is important to note that this may vary with different implementations). However, the method presented in this paper is much more flexible than a precalculated table because the physical conditions, the chemical composition, and the grain size-distribution can be changed at any point for a self-consistent calculation.
In Table 1, we provide the thermal ionization coefficients for K, Na, and H from Pneuman & Mitchell (1965), which is a theoretical work, to match the rates used in Marchand et al. (2016). In the 1960s and 1970s, there were many discussions about the ionization rates of alkali metals. Flame experiments to measure those rates (Hollander et al. 1963; Ashton & Hayhurst 1973) resulted in much larger cross sections than theoretically predicted (Aller 1961; Hollenbach & Salpeter 1969, see Schofield 1965; Shui 1977 for review). Desch & Turner (2015) argue that the experimental value of Ashton & Hayhurst (1973) should be preferred to the theoretical value of Pneuman & Mitchell (1965). This is debatable as the conditions in flames may be difficult to control and different from astrophysical plasmas (pressure, chemical composition), where Na and K ionize by colliding with H2. The larger coefficient rates would suggest thermal ionization of K and Na at a lower temperature, typically T = 1200 K instead of T = 1600 K, which could be of importance in protoplanetary disks. For reference, we provide the experimental values of Ashton & Hayhurst (1973) for K and Na:
5. Conclusions
We detail an extension of the ionization model of Marchand et al. (2021) to include the thermal ionization of one gas species. It is possible to increase the number of ionized species by adding their contribution in the same manner, at the price of a higher numerical cost. We have presented the examples of K, Na, and H. Both the chemical abundances and the resistivities show a behavior consistent with previous work. This method is then a powerful tool to self-consistently calculate the ionization of the dust-grain mixture and the nonideal MHD resistivities in hydrodynamical simulations, faster than a complete chemical network. The method is flexible and valid in a wide variety of environments, including star formation and protoplanetary disks.
Acknowledgments
We thank the referee for their insightful comments that helped improve the manuscript. We thank Jérémy E. Cohen for his insight on the solving of linear systems. P. M. acknowledges financial support by the Kathryn W. Davis Postdoctoral Fellowship of the American Museum of Natural History. U. L. acknowledges financial support from the European Research Council (ERC) via the ERC Synergy Grant ECOGAL (grant 855130). M.-M. M. L. acknowledges partial support from NSF grant AST18-15461.
References
- Aller, L. H. 1961, The Abundances of the Elements (New York: Interscience Publ) [Google Scholar]
- Ashton, A., & Hayhurst, A. 1973, Combust. Flame, 21, 69 [CrossRef] [Google Scholar]
- Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 [Google Scholar]
- Desch, S. J., & Turner, N. J. 2015, ApJ, 811, 156 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T., & Sutin, B. 1987, ApJ, 320, 803 [NASA ADS] [CrossRef] [Google Scholar]
- Ebel, D. S., Hubbard, A., McNally, C., et al. 2012, Meteorit. Planetary Sci. Suppl., 75, 5387 [NASA ADS] [Google Scholar]
- Getman, K. V., & Feigelson, E. D. 2021, ApJ, 916, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Hollander, T., Kalff, P. J., & Alkemade, C. T. J. 1963, J. Chem. Phys., 39, 2558 [NASA ADS] [CrossRef] [Google Scholar]
- Hollenbach, D. J., & Salpeter, E. E. 1969, J. Chem. Phys., 50, 4157 [NASA ADS] [CrossRef] [Google Scholar]
- Hubbard, A., McNally, C. P., & Mac Low, M.-M. 2012, ApJ, 761, 58 [NASA ADS] [CrossRef] [Google Scholar]
- Joung, M. K. R., Mac Low, M.-M., & Ebel, D. S. 2004, ApJ, 606, 532 [NASA ADS] [CrossRef] [Google Scholar]
- Koga, S., Tsukamoto, Y., Okuzumi, S., & Machida, M. N. 2019, MNRAS, 484, 2119 [Google Scholar]
- Kunz, M. W., & Mouschovias, T. C. 2010, MNRAS, 408, 322 [NASA ADS] [CrossRef] [Google Scholar]
- Larson, R. B. 1969, MNRAS, 145, 271 [Google Scholar]
- Lenzuni, P., Gail, H.-P., & Henning, T. 1995, ApJ, 447, 848 [NASA ADS] [CrossRef] [Google Scholar]
- Li, Z.-Y., Krasnopolsky, R., & Shang, H. 2011, ApJ, 738, 180 [Google Scholar]
- Marchand, P., Masson, J., Chabrier, G., et al. 2016, A&A, 592, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marchand, P., Commerçon, B., & Chabrier, G. 2018, A&A, 619, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marchand, P., Guillet, V., Lebreuilly, U., & Mac Low, M. M. 2021, A&A, 649, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 [Google Scholar]
- McElroy, D., Walsh, C., Markwick, A. J., et al. 2013, A&A, 550, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- McNally, C. P., Hubbard, A., Mac Low, M.-M., Ebel, D. S., & D’Alessio, P. 2013, ApJ, 767, L2 [NASA ADS] [CrossRef] [Google Scholar]
- McNally, C. P., Hubbard, A., Yang, C.-C., & Mac Low, M.-M. 2014, ApJ, 791, 62 [NASA ADS] [CrossRef] [Google Scholar]
- Padovani, M., Hennebelle, P., & Galli, D. 2013, A&A, 560, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pneuman, G. W., & Mitchell, T. P. 1965, Icarus, 4, 494 [NASA ADS] [CrossRef] [Google Scholar]
- Schofield, K., & Sugden, S. T. 1965, Some Observations on the Ionization of Alkali and Alkaline-earth Elements in Hydrogen Flames (Cambridge: University of Cambridge), 10 [Google Scholar]
- Shui, V. H. 1977, Phys. Fluids, 20, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Spitzer, L., Jr 1949, Leaflet Astron. Soc. Pac., 5, 336 [Google Scholar]
- Teyssier, R. 2002, A&A, 385, 337 [CrossRef] [EDP Sciences] [Google Scholar]
- Umebayashi, T., & Nakano, T. 1990, MNRAS, 243, 103 [NASA ADS] [CrossRef] [Google Scholar]
- Wurster, J., Price, D. J., & Bate, M. R. 2016, MNRAS, 457, 1037 [Google Scholar]
Appendix A: Jacobian components
The components of the normalized Jacobian matrix required for the solution of equation (36) are
with
All Tables
Constants depending on the choice of species s for sodium, potassium, and hydrogen.
All Figures
![]() |
Fig. 1. Test of our method for protostellar collapse conditions. Evolution of the average charge of grains (left axis) for several grain sizes (color lines), and the fractional abundance (right axis) of cosmic-ray ionized ions ni/nH (solid line), thermally ionized K+ ions ns/nH (dotted line), and electrons ne/nH (dashed line). |
In the text |
![]() |
Fig. 2. Same as Fig. 1, but for Na. The only ions present are the cosmic-ray ionized ones i and Na (ion s), which is both cosmic-ray and thermally ionized. |
In the text |
![]() |
Fig. 3. Same as Fig. 1, but for H. The only ions presents are the cosmic-ray ionized ones i and H (ions s), which is both cosmic-ray and thermally ionized. |
In the text |
![]() |
Fig. 4. Evolution of the Ohmic (yellow), Hall (blue), and ambipolar (dark red) resistivities for the thermal ionization of K. The dashed line represents the Hall resistivity in negative values. The thin dotted lines are the resistivities in the absence of thermal ionization. |
In the text |
![]() |
Fig. 5. Same as Fig. 4, but for Na. |
In the text |
![]() |
Fig. 6. Same as Fig. 4, but for H. |
In the text |
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.