Advances in Astronomy

Volume 2013 (2013), Article ID 487606, 4 pages

http://dx.doi.org/10.1155/2013/487606

Research Article

Incorporation of ³⁶Cl into Calcium-Aluminum-Rich Inclusions in the Solar Wind Implantation Model

Glynn E. Bricker¹ and Marc W. Caffee²

¹Department of Physics, Purdue University North Central, Westville, IN 46391, USA

²PRIMELab, Department of Physics, West Lafayette, IN 47907, USA

Received 11 September 2013; Accepted 22 October 2013

Academic Editor: Alberto J. Castro-Tirado

Copyright © 2013 Glynn E. Bricker and Marc W. Caffee. 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.

Abstract.

We consider the short-lived radionuclide (SLR) ³⁶Cl in calcium-aluminum-rich inclusions (CAIs) found in primitive meteorites with the solar wind implantation model. In this model, SLRs are produced via nuclear reaction with solar energetic particles (SEPs) interacting with gaseous targets in the protosolar atmosphere in T-Tauri stars. These SLRs are captured by the solar wind and then implanted in CAI precursor materials, which have dropped from the funnel flow leading onto the protostar. This method of incorporating SLRs into solar system materials is currently active in our solar system and has been measured with SLRs from the solar wind being implanted on the lunar surface. T-Tauri stars are capable of SEP fluxes ~10⁵ greater than contemporary SEP fluxes. Here we scale the production rate of ³⁶Cl to the ancient SEP activity. From the enhanced production rates and the refractory mass inflow rate at 0.06 AU from the protosun, we model the ancient ³⁶Cl content in CAIs. We find the initial isotopic ratio of ³⁶Cl/³⁵Cl to range from about 1 × 10⁻⁵ to 5 × 10⁻⁵ and the concentration of ³⁶Cl to range from about 3 × 10¹³ to 1.5 × 10¹⁴ atoms g⁻¹.

1. Introduction

Studies report evidence for the incorporation of the short-lived radionuclide (SLR) (T1/2 = 0.3 Myr) ³⁶Cl into early solar system materials, including calcium-aluminum-rich inclusions (CAIs). Lin et al. [1] infer an initial ³⁶Cl/³⁵Cl ratio of ≥ in sodalite from the carbonaceous chondrite Ningqiang based on Al-Mg systematics. Also, relying upon Al-Mg systematics, Jacobsen et al. [2] report that initial ratio of ³⁶Cl/³⁵Cl in wadalite from Allende would have been > had the initial ²⁶Al/²⁷Al ratio been the canonical value of ~.

Clues to the source of the ³⁶Cl can be gleaned from this initial ratio. If the initial ³⁶Cl/³⁵Cl had been ~, the source of ³⁶Cl most likely would not have been from a nearby super nova or AGB star [3]. A remaining mechanism for the production of this radionuclide is local irradiation from solar energeticparticle (SEP) born in the protosolar atmosphere. The initial SLR ratios can be used to constrain early solar system evolution and conditions. This is especially true in terms of solar luminosity and flaring characterization.

Bricker and Caffee [4] proposed a solar wind implantation model for incorporation of ¹⁰Be in CAI precursor materials. In this model, ¹⁰Be and possibly other SLRs are produced by SEP reactions in the protosolar atmosphere of a more energetic T-Tauri sun, characterized by SEP fluxes many orders of magnitude greater than contemporary particle fluxes. Studies of the Orion Nebulae indicate that premain sequence (PMS) stars exhibit X-ray luminosity and hence SEP fluxes on the order of ~10⁵ over contemporary SEP flux levels [5]. The SLRs are produced through nuclear reactions with these SEPs, and the irradiation produced SLRs are then entrained in the solar wind and subsequently implanted into CAI precursor material. This production mechanism is operational in the contemporary solar system and is responsible for implantation of solar wind nuclei, including ¹⁰Be [6] and ¹⁴C [7], in lunar material. In this work, we consider ³⁶Cl found in CAIs in primitive carbonaceous meteorites in accordance with a solar wind implantation model. Table 1 characterizes ³⁶Cl found in CAIs.

Table 1: Characteristics of ³⁶Cl found in CAIs.

Stable nuclide
Half-life
Initial isotopic ratio
Stable nuclide (ppm)
Radionuclide (g⁻¹)

35Cl
0.3 Myr
1.6 × 10⁻⁴
228
1.1 × 10¹⁵

References Lin et al. [1] and Leya et al. [8].

2. Solar Wind Implantation Model

2.1. Overview

In the solar wind implantation model, SLRs are produced in the solar nebula ~4.6 Gyr ago by the bombardment of target material in the solar atmosphere by solar energetic particles. These SLRs escape the solar atmosphere entrained in the solar wind. Some fraction of these outward flowing SLRs are incorporated into inflowing material which has fallen from the main accretion flow from the protoplanetary accretion disk. In the region in which the inflowing material and outflowing solar wind intersect, SLRs may be incorporated into the precursor CAI material. The fluctuating -wind model of Shu et al. [9–11] provides the basic framework for incorporation of SLRs into CAI-precursor materials and the subsequent transportation of these implanted refractory materials to asteroidal distances. Figure 1 is a cartoon of the basic magnetic field geometry, main accretion flow, and SLRs, including ³⁶Cl outflow.

Figure 1: Solar wind implantation model of magnetic field geometry for SLRs production via nuclear reactions. The gray area represents accretion flow onto “hot spots” on the PMS star. SLRs produced close to the protosolar surface are implanted into CAI material which has fallen from the accretion flow (figure after Shu et al. [10]).

2.2. Refractory Mass Inflow Rate

The refractory mass inflow rate, that is, the mass that falls from the funnel flow onto the star at the -region, is given by

where is disk mass accretion rate, is the cosmic mass fraction, and is the fraction of material that enters the -region [12]. For , we adopt solar masses year⁻¹. The mass accretion rates can vary from ~10⁻⁷ to ~10⁻¹⁰ solar masses year⁻¹ for T-Tauri stars from 1 to 3 Myr [13]. Embedded class 0 and class I PMS stars can have accretion rates from ~10⁻⁵ to ~10⁻⁶ solar masses year⁻¹ [14]. The value we adopt corresponds to class II or III PMS stars. Following Lee et al. [12], we adopt a cosmic mass fraction, , and fraction of refractory material fraction , of and 0.01, respectively. represents the fraction of material that is refractory and represents the fraction of mass that does not accrete onto the protosun. The choice 0.01 is a maximum for and corresponds to all the mass which comprises the planets falling from the accretion flow. If some of the rocky material accreting towards the proto-Sun never reached the edge of the accretion disk, would be smaller by a factor of 20 at most, but this scenario is unlikely. The choice of is the preferred value of Lee et al. [12] in their model. (See Lee et al. [12] for a detailed discussion of and .) From (1) and the values described above, we find the rate at which this refractory material is carried into the -region, called here the refractory mass inflow rate, , is  g s⁻¹. On the extremes, could be two orders of magnitude greater if the PMS star was classes 0 or 1; could also be four orders of magnitude less if the mass accretion rate was to solar masses year⁻¹ and .

2.3. Ancient Effective Production Rate

The effective ancient ³⁶Cl outflow rate, in units of s⁻¹, is given by

where is the ancient production rate and is the fraction of the solar wind ³⁶Cl that was captured into the CAI-forming region; . (See Bricker and Caffee [4] for a discussion of factor .) We calculate the ³⁶Cl production rates assuming that solar energetic particles are characterized by a power law relationship:

where ranges from 2.5 to 4. For impulsive flares, that is, , we use ³He/H = 0.1 and ³He/H = 0.3, and for gradual flares, that is, , we use ³He/H = 0. Contemporary SEP fluxes at 1 AU are ~100 protons cm⁻²s⁻¹ for  MeV [15]. We assume an increase in ancient particle fluxes over the current particle flux of , yielding an energetic particle flux of protons cm⁻²s⁻¹ for  MeV at the surface of the protosun.

The production rates for cosmogenic nuclides can be calculated via

where represents the target elements for the production of the considered nuclide, is the abundance of the target element (g g⁻¹), indicates the energetic particles that cause the reaction, is the cross-section for the production of the nuclide from the interaction of particle with energy from target for the considered reaction (cm²), and is the differential energetic particle flux of particle at energy (cm⁻² s⁻¹) [15]. We assume gaseous targets, Cl, K, S, and Ca, of solar composition [16].

The cross-sections we use are from Gounelle et al. [17], with the exception of (³He,x)³⁶Cl, which  is a new experimental cross-section from Herzog et al. [18]. The Gounelle et al. [17] cross-sections are a combination of experimental data, fragmentation, and Hauser-Feshbach codes. The uncertainty associated with model codes is at best a factor of two. We have used experimental cross-sections, whenever possible, in order to limit uncertainties associated with the calculations. The reactions we have considered here are the primary production pathways. This takes into account both abundance of target material and cross-sections. Any other reaction would add little to the overall ³⁶Cl production rate. Table 2 shows the cross-sections used in the calculations.

Table 2: Nuclear reactions considered in this paper.

36S(p, n)³⁶Cl

37Cl(p, pn)³⁶Cl

34S(³He, p)³⁶Cl

35Cl(³He, 2p)³⁶Cl

39K(³He, x)³⁶Cl

natCa(³He, x)³⁶Cl

3. Results

The concentration of ³⁶Cl found in refractory rock predicted by our model is given by

where is given atoms s⁻¹ and is given in g s⁻¹.

From the value of given above from (1) and calculations of from (4), we calculate the concentration of ³⁶Cl in CAIs using (5) and find the associated isotopic ratio for different flare parameters given in Table 3, and we plot the predicted ³⁶Cl content for spallation production from energetic protons in Figure 2.

Table 3: Predicted ³⁶Cl content in CAIs.

Flare parameter
Atoms g⁻¹ (in CAIs)
Isotopic ratio

, ³He/H = 0
3.35 × 10¹³
1.14 × 10⁻⁵

, ³He/H = 0.1
5.26 × 10¹³
1.79 × 10⁻⁵

, ³He/H = 0.3
1.54 × 10¹⁴
5.25 × 10⁻⁵

Figure 2: Predicted ³⁶Cl content in CAIs from energetic protons as a function of solar flare parameter.

4. Discussion

From the results above, the isotopic ratio predicted by the solar wind implantation model is from a factor of 4 to a factor of 10 below the inferred initial ³⁶Cl/³⁵Cl ratio of . Given the uncertainties in the parameters, that is, , , and , the model is viable for ³⁶Cl/³⁵Cl initial ratio of . Jacobsen et al. [2] report that, based on Al-Mg systematics, the initial ³⁶Cl/³⁵Cl ratio may have been >. The model calculations underproduce this value by several orders of magnitude for the flare parameters given here. It is not possible to produce an order of magnitude correct ³⁶Cl/³⁵Cl ratio of without overproducing ¹⁰Be/⁹Be, ⁴¹Ca/⁴⁰Ca, and ⁵³Mn/⁵⁵Mn by several orders of magnitude.

Clearly, some other mechanism is needed to explain to provenance of ³⁶Cl at the levels described by Jacobsen et al. [2]. A possible scenario is that the irradiation of target material occurred in a volatile rich region away from the protosun. This region would contain much more target materials, that is, Cl, S, and K. Greater target material would lead to greater initial ³⁶Cl content in CAIs. This argument for the enhanced ³⁶Cl found in wadalite samples is also invoked by Jacobsen et al. [2].

Evidence suggests that the distribution of ²⁶Al was homogenous in the solar system, but this is not always the case [3]. Marhas and Goswami [19] found several CAIs that had the “canonical” ¹⁰Be/⁹Be ratio, but were devoid of ²⁶Al, demonstrating a decoupling of ²⁶Al from ¹⁰Be. If ³⁶Cl is also decoupled from ²⁶Al, it is difficult to infer what the initial ³⁶Cl/³⁵Cl ratios were based solely on Al-Mg systems although Jacobsen et al. [2] do report a canonical initial ²⁶Al/²⁷Al ratio for primary minerals in their sample.

The irradiation origin of ²⁶Al found in CAIs is a matter of considerable debate. The proposed region where the enhanced ³⁶Cl production may have taken place could have already been seeded with ²⁶Al from some nonirradiation origin, that is, supernova. Conversely, the origin of ¹⁰Be is likely from SEP irradiation close to the protosun [4, 17]. It would be beneficial to establish a correlation between ¹⁰Be and ³⁶Cl to find if the irradiation that most likely produced ¹⁰Be also produced ³⁶Cl. More studies are needed in this area to establish this relationship.

References

Y. Lin, Y. Guan, L. A. Leshin, Z. Ouyang, and D. Wang, “Short-lived chlorine-36 in a Ca- and Al-rich inclusion from the Ningqiang carbonaceous chondrite,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 5, pp. 1306–1311, 2005.

B. Jacobsen, J. Matzel, I. D. Hutcheon et al., “Formation of the short-lived radionuclide ³⁶Cl in the protoplanetary disk during late-stage irradiation of a volatile-rich reservoir,” Astrophysical Journal Letters, vol. 731, no. 2, pp. L28–L34, 2011.

G. J. Wasserburg, M. Busso, R. Gallino, and K. M. Nollett, “Short-lived nuclei in the early solar system: possible AGB sources,” Nuclear Physics A, vol. 777, pp. 5–69, 2006.

G. E. Bricker and M. W. Caffee, “Solar wind implantation model for ¹⁰Be in calcium-aluminum inclusions,” Astrophysical Journal, vol. 725, no. 1, pp. 443–449, 2010.

E. D. Feigelson, G. P. Garmtextre, and S. H. Pravdo, “Magnetic flaring in the pre-main-sequence sun and implications for the early solar system,” Astrophysical Journal Letters, vol. 572, no. 1, pp. 335–349, 2002.

K. Nishiizumtext and M. W. Caffee, “Beryllium-10 from the sun,” Science, vol. 294, no. 5541, pp. 352–354, 2001.

A. J. T. Jull, D. Lal, L. R. McHargue, G. S. Burr, and D. J. Donahue, “Cosmogenic and implanted radionuclides studied by selective etching of lunar soils,” Nuclear Instruments & Methods in Physics Research B, vol. 172, no. 1–4, pp. 867–872, 2000.

I. Leya, A. N. Halliday, and R. Wieler, “The predictable collateral consequences of nucleosynthesis by spallation reactions in the early solar system,” Astrophysical Journal Letters, vol. 594, no. 1, pp. 605–616, 2003.

F. H. Shu, H. Shang, and T. Lee, “Toward an astrophysical theory of chondrites,” Science, vol. 271, no. 5255, pp. 1545–1552, 1996.

F. H. Shu, H. Shang, A. E. Glassgold, and T. Lee, “X-rays and fluctuating x-winds from protostars,” Science, vol. 277, no. 5331, pp. 1475–1479, 1997.

F. H. Shu, H. Shang, M. Gounelle, A. E. Glassgold, and T. Lee, “The origin of chondrules and refractory inclusions in chondritic meteorites,” Astrophysical Journal Letters, vol. 548, no. 2, pp. 1029–1050, 2001.

T. Lee, F. H. Shu, H. Shang, A. E. Glassgold, and K. E. Rehm, “Protostellar cosmtextc rays and extinct radioactivities in meteorites,” Astrophysical Journal Letters, vol. 506, no. 2, pp. 898–912, 1998.

N. Calvet, C. Briceño, J. Hernández et al., “Disk evolution in the orion OBI association,” Astronomtextcal Journal, vol. 129, no. 2, pp. 935–946, 2005.

D. Ward-Thompson, “The formation and evolution of low mass protostars,” Astrophysics & Space Science, vol. 239, no. 1, pp. 151–170, 1996.

R. C. Reedy and K. Marti, “Solar-cosmtextc-ray fluxes during the last 10 mtextllion years,” in The Sun in Time, C. P. Sonnet, M. S. Giampapa, and M. S. Mathews, Eds., pp. 260–287, University of Arizona Press, Tucson, Ariz, USA, 1991.

K. Lodders, “Solar system abundances and condensation temperatures of the elements,” Astrophysical Journal Letters, vol. 591, no. 2, pp. 1220–1247, 2003.

M. Gounelle, F. H. Shu, H. Shang, A. E. Glassgold, K. E. Rehm, and T. Lee, “The irradiation origin of beryllium radioisotopes and other short-lived radionuclides,” Astrophysical Journal Letters, vol. 640, no. 2, pp. 1163–1170, 2006.

G. F. Herzog, M. W. Caffee, T. Faestermann et al., “Cross sections from 5 to 35 MeV for the reactions (³He,x)²⁶Al, ²⁷Al(³He,x)²⁶Al, (³He,x)⁴¹Ca, and (³He,x)³⁶Cl: implications for early irradiation in the solar system,” Meteoritics and Planetary Science, vol. 46, no. 10, pp. 1427–1446, 2011.

K. K. Marhas and J. N. Goswami, “Low energy particle production of short-lived nuclides in the early solar system,” New Astronomy Reviews, vol. 48, no. 1–4, pp. 139–144, 2004.