Open Access
Volume 624, April 2019
Article Number A131
Number of page(s) 12
Section Interstellar and circumstellar matter
Published online 25 April 2019

© E. Jacquet 2019

Licence Creative Commons
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Primitive meteorites, or chondrites, are conglomerates of solids that formed and accreted in the solar protoplanetary disk. Among those, the oldest are the refractory inclusions, comprising calcium-aluminium-rich inclusions (CAIs) and amoeboid olivine aggregates (AOA) (e.g., MacPherson 2014; Krot et al. 2004). Thermodynamic calculations predict that CAIs should be the first condensates in a cooling gas of solar composition (Grossman 2010; Davis & Richter 2014). Since the required temperatures would be in the range 1500–2000 K (Wooden et al. 2007), it is widely assumed that refractory inclusions formed close to the Sun, but precisely how close and in what astrophysical setting remains unclear (Wood 2004; Jacquet 2014).

An important clue in this respect may be provided by beryllium-10, a short-lived radionuclide decaying into boron-10 with a half-life of 1.5 Ma, and which, unlike other known extinct radionuclides such as aluminum-26, is not produced by stellar nucleosynthesis but may be formed through spallation by energetic particles (Davis & McKeegan 2014; Lugaro et al. 2018). Indeed, since the original discovery by McKeegan et al. (2000), all CAIs with analytically suitable Be/B have shown 10B/11B excesses correlated therewith. The slopes of the resulting isochrons translate into initial (i.e., upon the last isotopic equilibration) 10Be/9Be ratios averaging ~6 × 10−4 in CV chondrite (type A and B) CAIs (Davis & McKeegan 2014), with less than a factor of two spread. Yet systematically lower values around 3 × 10−4 and 5 × 10−4 have been found for two CV chondrite FUN CAIs (fractionated and unknown nuclear effects; MacPherson et al. 2003; Wielandt et al. 2012) and platy hibonite crystals (PLAC) in CM chondrites (Liu et al. 2009, 2010), respectively, and, conversely, values up to 10−2 have been found for Isheyevo CAI 411 (Gounelle et al. 2013) and more recently, fine-grained group II CV chondrite CAIs (Sossi et al. 2017). The overall spread, along with correlated 50V excesses (also ascribed to spallation) in the latter objects (Sossi et al. 2017), is difficult to reconcile with simple 10Be inheritance from the (galactic cosmic ray-irradiated) protosolar cloud (Desch et al. 2004), which should be largely homogeneous, and argues in favor of local production in the disk following flares from the young Sun (Gounelle et al. 2001, 2006, 2013; Sossi et al. 2017), such as those manifested in X-rays by present-day protostars (e.g., Feigelson et al. 2002; Wolk et al. 2005; Preibisch et al. 2005; Telleschi et al. 2007; Güdel et al. 2007; Bustamante et al. 2016) or evidence for enhanced ionization in some protostellar envelopes (Ceccarelli et al. 2014; Favre et al. 2017, 2018).

Since the energetic protons would not penetrate further than ~690 kg m−2 (Umebayashi & Nakano 1981) into the gas, Sossi et al. (2017) suggested that CAIs spent a few thousand orbits at the inner edge of the protoplanetary disk, similar to earlier studies that had adopted the framework of the X-wind scenario (Lee et al. 1998; Shu et al. 2001; Gounelle et al. 2001). However, among the objections to the latter discussed by Desch et al. (2010), a general issue is that CAI-forming temperatures would be expected significantly further outward (0.1–1 AU), because of local dissipation of turbulence. If nonetheless CAIs all went somehow through this narrow region and escaped accretion to the Sun or ejection from the solar system, they would be essentially “lucky” foreign material in the chondrites that incorporated them. Yet CAIs are present in carbonaceous chondrites at about the abundances predicted by in situ condensation (Jacquet et al. 2012). Also, although (non-CI) carbonaceous chondrites are enriched in refractory elements relative to CI chondrites, this cannot be ascribed to simple CAI addition to CAI-free CI chondritic material for should the CAIs bementally subtracted, these chondrites would have subsolar refractory element abundances (Hezel et al. 2008), not even considering those CAIs that were converted into chondrules (Misawa & Nakamura 1988; Jones & Schilk 2009; Metzler & Pack 2016; Ebert & Bischoff 2016; Jacquet & Marrocchi 2017; Marrocchi et al. 2018). Isotopic systematics for elements of different volatilities also indicate a need for a non-refractory isotopically CAI-like component in chondrites (e.g. Nanne et al. 2019). So there is evidence for a genetic relationship between CAIs and (part of) their host carbonaceous chondrites, which argues against an origin at the disk inner edge, whose contribution to distant chondritic matter would likely be minor.

However, the inner edge of the disk is not the only region where CAIs or their precursors could have been exposed to energetic solar protons. CAIs floating further out in the disk would at times have reached the upper layers of the disk. While such excursions may individually incur modest proton fluences, their cumulative contributions might explain the total 10Be evidenced in CAIs. Moreover, prior to CAI condensation, the gas exposed at the surface of the disk may also have undergone spallation and passed on the then-produced 10Be to the later formed condensates. The purpose of this paper is to analytically calculate the amount of beryllium-10 produced by such channels and thence evaluate the viability of having CAIs form at 0.1–1 AU from the young Sun against their irradiation record. Figure 1 provides a sketch of the overall scenario and setting. After several generalities in Sect. 2, I will express the 10Be/9Be ratio in the disk and free-floating solid CAIs in Sect. 3. I then discuss the numerical evaluations thereof and their implications in Sect. 4 before concluding in Sect. 5.

thumbnail Fig. 1

Sketch of the scenario of 10Be production developed in this paper. Solar cosmic rays cause spallation reactions in the upper layers of the disk, in the gas and/or condensates, within a region (outlined by a dashed separatrix) where magnetic field lines originating from the proto-Sun and its outskirts hit the disk. The field and particle emission geometries are purely illustrative and do not detract from the generality of the formalism.

Open with DEXTER

2 Prolegomena

2.1 Disk model

I consider aprotoplanetary disk in a cylindrical coordinate system with heliocentric distance R. The inner regions can be described under the steady-state approximation so long as the evolution timescales of the system are long compared to their local viscous timescale (1)

where is the effective turbulent viscosity, α the dimensionless turbulence parameter (e.g., Jacquet 2013), the isothermal sound speed with kB the Boltzmann constant, T the temperature and m = 3.9 × 10−27 kg the mean molecular mass, and ΩK the Keplerian angular velocity. Since tvis(R) is also much shorter than the half-life of 10Be, I will neglect its decay during the CAI formation epoch. This is consistent with the short timescales (no more than a few hundreds of millenia) of CAI formation and (isotopically resetting) re-heating events (e.g., MacPherson 2014), which should have ended after the last production of 10Be so as to account for the isochron behavior of the Be–B system in each CAI (that is, with 10 B excesses scaling with Be rather than the target nuclides leading to its parent; see next subsection).

If infall from the parental cloud can be neglected in the inner region, the disk mass accretion rate = − 2πRΣuR (with Σ the disk surface density and uR the net, turbulence-averaged, gas radial velocity) is uniform and obeys (2)

or, equivalently, (3)

where we have assumed that the stress vanishes at the disk inner edge R0 (e.g., Balbus & Hawley 1998). When calculating 10Be/9Be in later sections, I shall show expressions for a general (R) before specializing to the case of uniform mass accretion rate. Appendix A explores the corrections of infall to Eq. (2).

If I inject the temperature due to viscous dissipation of turbulence (e.g., Appendix A of Jacquet et al. 2012) in Eq. (2), I obtain (4)

with κ the specific Rosseland mean opacity and σSB the Stefan-Boltzmann constant. This confirms that CAI-forming temperatures may be reached at a fraction of an AU for ≳ 10−7 Ma.

2.2 10Be production rate

The local production rate of 10Be by spallation normalized to 9Be may be written as (e.g., Sossi et al. 2017) (5)

with σk,i the 10Be production cross section for target (tg) nuclide k and cosmic ray species i (with energy per nucleon E) and I(i, E, Ω) the corresponding (orientation-dependent) monoenergetic cosmic ray specific intensity (defined by number and not energy, as in Gounelle et al. 2001). Here and throughout, all isotopic or elemental ratios are atomic and the Einstein convention for summation over the repeated indices k and i is adopted. I assume that the incoming cosmic rays have a solar composition (e.g. 4 He/H = 0.1; but see Mewaldt et al. 2007 for details on contemporaneous solar energetic particles), and obey a power law ∝ Ep for EE10 ≡ 10 MeV upon arrival on the disk. I will ignore the contribution of secondary neutrons, despite their comparable cross sections (e.g., Leya & Masarik 2009) and larger attenuation columns. This is because a dilute medium such as the surface of the disk will let free decay thwart significant accumulation of neutron flux (Umebayashi & Nakano 1981). However, based on this latter work, this is only marginally true (see Eq. (C.5)), so the production rates given here should be strictly viewed as a lower bounds.

In the following sections, we will be interested in the density (ρ)-weighted vertical average of the production rate at a given heliocentric distance, (6)

where Σ′ is the column density integrated vertically from the upper surface and the +subscript in the final equality denotes restriction to cosmic rays from the upper side1, assuming Σ is much larger than the attenuation column Σp,i(E) of the cosmic rays, defined here as (7)

where Σ″ is the actually traversed column density. If I envision a cosmic ray beam (with the particles gyrating along a magnetic field line) penetrating the disk with a grazing angle ϕ and ignore secondary particles, Σ″ = Σ′∕(μsinϕ) with μ the cosine of the pitch angle of the cosmic ray with respect to the magnetic field line (Padovani et al. 2018). I thus obtain (8)

where n is the downward-directed unit vector normal to the disk surface and d FH ∕dE the incoming differential proton number vector flux. If I scale the latter to the incoming > 10 MeV nucleon−1 energy vector flux F10 and extract O as a proxy for all 10Be-producing target nuclides, Eq. (8) can be further manipulated to finally yield (9)

with (10)

which is plotted in Fig. 2 (see Appendix B for further details on its calculation and data sources). Here, f(a, b) is the fraction of the total >10 MeV nucleon−1 energy luminosity L10 that reaches the disk between heliocentric distances a and b. The latter depends on the geometry of energetic particle emission and the magnetic field configuration around the young Sun, which are largely unknown (protosolar corona, star-disk fields, or disk fields; Feigelson et al. 2002; Donati & Landstreet 2009; Padovani et al. 2016). For reference, in the case of a geometrically flared disk, an isotropic point source with straight trajectories would yield f(R0, R) = HspallR with Hspall the height of the spallation layer above the midplane (see Appendix C for an estimate). For illustration in the plots of the next section, I will use a form of f assuming that the outgoing energy flux density of the energetic particles is uniform over the Sun-centered sphere of radius R0 and that the magnetic field is approximately dipolar (with a magnetic moment along the rotation axis of the disk) for field lines interior to a maximum irradiation distance Rmax. This yields2

However, in the main text, I will keep general expressions as functions of f. This parameterization of our ignorance will turn out to be quite useful in the calculation.

A means of comparison with previous models invoking irradiation of bare solids (e.g., Lee et al. 1998; Gounelle et al. 2006; Sossi et al. 2017) is to divide the rate given by Eq. (9) by that in a target of the same composition exposed to the same unattenuated proton flux. Assuming isotropic distribution of cosmic rays (in each half-space along the field line), this gives (12)

with (13)

which (for a refractory target) is 90 kg m−2 for p = 2.5 typical of gradual flares favored by Gounelle et al. (2013) and Sossi et al. (2017), and 5 kg m−2 for p = 4 commensurate with impulsive flares (with subdominant fluences in the present-day Sun; Desai & Giacalone 2016). Obviously, for Σ ≫ 104 kg m−2 (the Minimum Mass Solar Nebula at 1 AU; Hayashi 1981) this calls for (intermittent) irradiation timescales 2 ± 1 orders of magnitude above the centuries calculated by Sossi et al. (2017). As alluded to in the Introduction, the transport timescale tvis may, however, provide the correct order of magnitude. The following section undertakes to calculate more precisely the net 10Be abundances produced in gas and solids and during radial transport.

thumbnail Fig. 2

Kp as a function of the power-law exponent of the incoming cosmic ray energy distribution, for a solar (without possible implantation contributions) and a rocky target (see Sects. 3.1 and 3.2, respectively).

Open with DEXTER

3 10Be abundances

3.1 Bulk 10Be/9Be in the disk

In this subsection, I calculate the average 10Be/9Be of the disk as a function of heliocentric distance, irrespective of its (temperature-dependent) physical state (condensation fraction) to which spallation reactions are insensitive. I assume solar abundances throughout the inner disk owing to tight coupling of the condensates with the gas (e.g., Jacquet et al. 2012).

A steady-state gradient of 10Be/9Be should arise in the disk because production of 10Be is balanced by loss to the Sun by accretion. The transport equation reads (14)

with cBe, O the Be,O concentration in number per unit mass of solar gas and DR the turbulent diffusion coefficient. The Kp,disk factor here includes contributions from k = 12C, 16 O, and 14 N (even though the last one is unimportant, with N/O = 0.14, compared to C/O = 0.5; Lodders 2003) and i = 1H, 4 He. In addition, it includes contributions from indirect reactions (Desch et al. 2004), i.e. those where the heavy nuclei originate from the cosmic rays instead of the target, since, although the hereby produced 10 Be will retain the momentum of the incoming particles, it should be stopped further downstream in the disk. Since Galilean invariance dictates σi,k = σk,i, this amounts to multiplying each “direct” reaction contribution in the right-hand side of Eq. (10) by , where the second expression uses the assumption that the disk has a (solar) composition, identical to that of the cosmic rays (see also Appendix B). The resulting Kp,disk is plotted asthe “solar target” curve in Fig. 2.

Another “indirect” contribution would be solar wind implantation of 10 Be produced onthe proto-Sun (Bricker & Caffee 2010). Although Bricker & Caffee (2010) originally envisioned implantation on bare solids, one can equally envision implantation on the gas disk surface (followed by vertical mixing); this would amount to adding to the right-hand side of Eq. (14), with the young Sun’s 10Be (number) production rate and fSW the counterpart of f. This would amount to an effective increase of Kp,disk of . If I set equal to the present-day ratio of long-term average 1 AU 10Be (Nishiizumi & Caffee 2001)3 and >10 MeV proton (Reedy 1996) energy omnidirectional fluxes, this evaluates to 6 × 10−20 s2 m−2. From Fig. 2, at face value, solar wind implantation thus appears negligible compared to local (in-disk) spallation (except for steep p ≳ 4), but since we do not really know how toextrapolate 10Be production on the Sun to its very active early times, the nominal one order-of-magnitude deficit may not warrant definitive conclusions yet.

Returning to Eq. (14), the requirement that 10Be/9Be vanishes at infinity leads to the first integration (15)

Provided the radial Schmidt number ScRνDR has a finite lower bound, the requirement that 10Be/9Be does not diverge at the inner edge of the disk (since uRDR would not be integrable there from Eq. (3)) leads to the unique solution (16)

where the second equality assumes uniform and ScR, with (17)

which is plotted in Fig. 3.

Interestingly, the result is weakly sensitive to the details of the protoplanetary disk (e.g., turbulence level) with the reference to heliocentric distance being only implicit in feff. This is because the approximate ∝ R−2 decrease of the flux density is essentially compensated by the R2 factor in the radial transport timescale tvis. The above expression can be interpreted as a sum of contributions advected from outer regions (see next subsection, in particular Eq. (19)) and contributions diffused back outward. 10Be/9Be decreases monotonically outward, falling off roughly as R−3ScR∕2 outside the irradiated region, conforming to the equilibrium gradient of a passive scalar (Clarke & Pringle 1988; Jacquet & Robert 2013). In the limit ScR → 0, it converges pointwise toward the inner edge value f(R0, +) and a flat profile, but ScR is probably of order unity in magnetohydrodynamical turbulence (Johansen et al. 2006), not to mention the possibility of laminar wind-driven accretion (and thus even higher ScR) further out (e.g., Bai 2016). So contrary to earlier statements (Gounelle et al. 2013; Kööp et al. 2018a), 10 Be production in the gas would entail no spatial uniformity of the 10Be/9Be ratio.

thumbnail Fig. 3

feff (Eq. (17)) as a function of heliocentric distances for two values of Rmax for my toy parameterization of the particle flux distribution (see end of Sect. 2.2) and two values of the radial Schmidt number, with the lower one (higher relative diffusivity) in dotted lines. The parameter feff is proportional to the 10Be/9Be ratio (whose absolute magnitude is discussed at Eq. (22)) in the bulk (gas+solids) disk for given mass accretion rate and proton output and energy distributions.

Open with DEXTER

3.2 In situ 10Be production in free-floating solids

I now consider a (sub-) millimeter rocky solid, say a CAI, formed (or more precisely, whose beryllium last equilibrated with the gas) at an heliocentric distance R1. It will certainly inherit the 10Be/9Be ratio of the local reservoir calculated above, but it should also acquire additional 10 Be so long as itwanders through the irradiated region of the disk. The purpose of this subsection is to estimate this contribution.

Since the inner disk (where irradiation may occur) is dense, I assume the solid to be tightly coupled to the gas (that is, a gas-grain decoupling parameter S ≪ 1; Jacquet et al. 2012; this is consistent with the lack of bulk chemical fractionation assumed previously). It is also small enough for the thin target approximation (e.g., Sossi et al. 2017) to apply, that is, for Eq. (5) to apply at the scale of the whole particle as a function of the local (outside) monoenergetic intensities. Since the radial transport timescale tvis is much longer than the vertical transport timescale H2ν with H = cs∕ΩK the pressure scale height, I can use theaverage rate of 10Be production given by Eq. (9). That is, I view the random vertical motion of the solids due to turbulence as an ergodic process (see, e.g., Fig. 8 in Ciesla 2010) and identify the time average of this rate for a given solid to the instantaneous average over the long-term probability distribution of same-size solids, which here follows the distribution of the gas (but see Appendix C for the case of finite settling)4. For this refractory target, the sum on the right-hand side of Eq. (10) is essentially restricted to k = O and i = 1H, 4 He (with no “indirect reaction” contribution). This corresponds to the “rocky target” curve in Fig. 2.

As a result of the turbulent diffusion superimposed on advection by the gas mean velocity, the radial motion of an individual solid has a stochastic character. I will be content to set the “typical timescale” spent crossing a radial bin of width ΔRR to |ΔRuR|. Indeed, even when the flow is opposite to the transport envisioned (as in this subsection), the fraction of a population of solids at R at t = 0, which has diffused upstream beyond that distance at time t (, with erfc the complementary error function), is maximum for t = − ΔRuR. The total 10Be produced in the solid during transport between an heliocentric distance R1 and its exit from the irradiated region is then given by (19)

where the last equality assumes a uniform (R). Figure 4 plots f(R, +). Unlike the inherited disk value (which benefited from outward diffusion), nonzero values do not extend to source radii outside the irradiated region.

thumbnail Fig. 4

f(R, +) as a functionof heliocentric distances for two values of Rmax for my toy parameterization of the particle flux distribution (see end of Sect. 2.2). This is proportional to the 10Be/9Be ratio produced by in situ irradiation of a refractory solid as a function of its heliocentric distance of origin.

Open with DEXTER

4 Discussion

4.1 Magnitude of 10Be production

Provided a given CAI underwent isotopic equilibration after the last production of 10 Be, the 10Be/9Be indicated by the slope of an internal isochron is the sum of that inherited from the disk gas and that produced in situ in the solid (Eqs. (16) and (19), respectively), which may be written as (20)

with (21)

The contribution of in situ production (the term with the “CAI” subscripts in the above equation) is subdominant with respect to spallation in the solar gas (that is, fsumfeff(R1)) since, not to mention the ignored incipient settling (see Appendix C), (i) Kp,disk outweighs Kp,CAI by a factor ~3 (Fig. 2) as it includes additional production channels (e.g., target 12 C, and indirect reactions) and (ii) Be, as a refractory element (half-condensation temperature of 1452 K according to Lodders 2003), should be concentrated relative to O in a CAI with respect to the overall disk. Indeed Be concentrations typically are of order 0.1–1 ppm in the data of McKeegan et al. (2000) and Gounelle et al. (2013), compared to the CI chondrite value of 25.2 ppb (the latter amounting to O∕Be = 107 ≈ half the solar value; Palme et al. 2014). Nevertheless, the fine-grained CAIs of Sossi et al. (2017) exhibit concentrations down to 6 ppb5. So, at least for some Be-poor CAIs having spent a random walk in the irradiated region of the disk longer than the “typical timescale” used in Sect. 3.2, in situ production may not be entirely negligible. The deviation of the initial boron isotopic composition from solar in some CAIs, too large for an irradiated solar gas, may be a collateral effect of such a contribution (Liu et al. 2010; Gounelle et al. 2013); the same applies for cosmogenic 3 He and 21 Ne excesses found by Kööp et al. (2018a) in the same type of inclusions. This indicates that, at least for some objects, heliocentric distance R1 was inside the proton-irradiated region of the disk.

In order to numerically evaluate Eq. (20), I will, as previous authors (Lee et al. 1998; Gounelle et al. 2001; Sossi et al. 2017), use the X-ray luminosity LX as a proxy for L10, since only the former can be measured from distant young stellar objects (although Ceccarelli et al. (2014) estimated a proton output equivalent to L10 ≳ 1027 W, for LX < 7 × 1024 W, from the ionization level seen toward the single source OMC-2 FIR 4). Since their luminosities are higher than the most powerful flares of the contemporaneous Sun (e.g., Feigelson et al. 2002), solar flares are the least improper analogs in that calibration. Lee et al. (1998) derived a scaling L10 = 0.09LX by ratioing the total proton and X-ray fluences of the impulsive flares of solar cycle 21 (peaking around 1980). However, while impulsive flares correlate best with X-ray outputs (Lee et al. 1998), gradual flares actually dominate the proton fluences at 1 AU (Desai & Giacalone 2016), whatever the ratio at the X point considered by Lee et al. (1998) may be. This warrants re-examination of the scaling between X-rays and solar energetic particles (SEP) without restricting to impulsive flares. Emslie et al. (2012) evaluated the energy outputs of 38 solar eruptive events between 2002 and 2006. Those 21 with reported SEP outputs totaled a SEP output of 1.5 × 1025 J (dominated by >10 MeV protons since the energy spectra is shallower below this; Mewaldt 2006) and 2.6 × 1024 J in X-rays, hence a ratio of average L10LX ≈ 6 which I adopt as a normalizing value6. Another relevant – if more indirectly for our purpose – dataset is the catalog of 314 SEP events between 1984 and 2013 of Papaioannou et al. (2016). They represent a total >10 MeV proton 1 AU energy fluence of 5.1 × 103 J m−2 and total X-ray fluence of 54 J m−2 hence (L10LX)∂f ln R∕sin ϕ ≈ 23 (which averages different geometries of the flare sources vis-à-vis the Earth) assuming isotropic distribution of SEP at 1 AU, suggestive of the same order of magnitude for L10LX. Equation (20) then becomes (22)

The nominal value obligingly coincides with the average one measured for regular CAIs in CV chondrites (Davis & McKeegan 2014). This, however, should not hide the considerable uncertainties in several factors, in addition to the L10LX discussed above. In particular, Kp,disk, normalized to the value for p = 2.5 typical of gradual flares favored by Sossi et al. (2017), may lose one order of magnitude or so (depending on the possible implantation contributions) for a steeper energy distribution; on the other hand, the normalization value for fsum, inspired from the case of ballistic emission, could be an underestimate if the cosmic rays are focused toward the disk (see Fig. 3), but it depends on the essentially unknown magnetic field configuration and turbulent diffusion efficiency. While the predictive value of the model should not thus be overrated, this result does show that, in the current state of our knowledge, the evidence of extinct 10 Be in refractory inclusions does not mandate an origin at the very inner edge of the disk, and formation over a wider range of heliocentric distance, say 0.1–1 AU, can be envisioned. This is the main point of this work.

4.2 Spatio-temporal variations of the 10Be/ 9Be ratio

In this model, 10Be/9Be is proportional to the LX that is an observable. Since LX has only a shallow dependence on time (e.g., ∝t−0.36 for one solar mass according to Telleschi et al. 2007), compared to the decrease of the accretion rate (∝ t−1.4 according to Hartmann et al. 1998), LX should increase over time (as t1.04 if I combine these examples although derived from different data sources; see Fig. 5). It may indeed be verified in Fig. 6 that LX anticorrelates with , if with a fair amount of scatter, which reminds us of the elusive determinants of LX. Since with decreasing isotherms should recede toward the Sun (see Eq. (4)), fsum may also be expected to increase for a given T(R1) (e.g., a condensation front ) for a fixed magnetic field configuration. So, barring a strong decrease of L10LX, 10Be/9Be should increase with time for a given formation and/or equilibration temperature. That is, mass loss from the disk would increase the relative importance of the surficial layers prone to spallation, overcoming the decline in proton luminosity.

This expectation would provide an explanation for the lower 10Be/9Be of FUN CAIs in CV chondrites (MacPherson et al. 2003; Wielandt et al. 2012) and PLACs in CM chondrites (Liu et al. 2010) compared to “regular” CAIs. Indeed these CAIs are characterized by nucleosynthetic stable isotopic anomalies (e.g., in Ti, Ca; MacPherson 2014; Kööp et al. 2018b) much in excess of those seen in bulk chondrules (Gerber et al. 2017) or meteorites (Trinquier et al. 2009). This suggests they formed early, perhaps during infall of an isotopically heterogeneous protosolar cloud, before turbulent mixing essentially suppressed the heterogeneities (e.g., Boss 2012). Perhaps the LX was so low that their 10Be/9Be was in fact dominated by the background inherited from the protosolar cloud (with perhaps some irradiation from the protostar directly on the envelope before arrival on the disk; Ceccarelli et al. 2014) rather than later dominating local spallation (Liu et al. 2010; Wielandt et al. 2012).

At the other end of the spectrum, a relatively late formation might explain the high 10Be/9Be (7 × 10−3) of fine-grained CV chondrite CAIs analyzed by Sossi et al. (2017) or that of Isheyevo CAI 411 (10−2; Gounelle et al. 2013). Indeed the former CAIs have group II rare earth element patterns that are commonly ascribed to condensation in a reservoir previously depleted in an ultrarefractory condensate (e.g., Boynton 1989). One may then speculate that this preliminary fractionation took some time to complete (given the tight coupling of millimeter-sized solids and gas) so group II CAIs had to be of relatively late formation. However, while Isheyevo CAI 411 shows 26 Al/27Al < 10−6 – an 26Al depletion common to many CH chondrite CAIs –, the fine-grained CV chondrite CAIs tend to exhibit 26 Al/27Al ratios no lower than their melted counterparts (MacPherson et al. 2012; Kawasaki et al. 2019), unlike what a later formation would suggest in view of the 26Al half-life of 0.72 Ma (e.g., MacPherson et al. 2012). This conclusion is, however, predicated on an homogeneous 26 Al distribution, which is far from a foregone assumption in early times. In fact, the nucleosynthetic anomaly-bearing CAIs discussed above show reduced 26 Al/27Al (~ 10−5; MacPherson et al. 2014; Park et al. 2017) indicating an increase of that ratio in the CAI-forming region between their formation epoch and that of melted CAIs, and perhaps somewhat beyond. Interestingly also, some 26 Al should be produced by spallation along with 10Be; for their fine-grained CAIs, Sossi et al. (2017) estimated an increase of 26 Al/27Al of (0.3 − 1.2) × 10−5, comparable to the spread resolved by MacPherson et al. (2012) and Kawasaki et al. (2019). At any rate, this calls into question the usefulness of 26 Al as a chronometer during CAI formation, although it may have become more homogenized afterward. Alternatively, the coarse-grained CAIs may have acquired 10Be/9Be lower than their fine-grained counterparts because of equilibration with their surroundings at higher heliocentric distances, during the localized heating events that melted them.

This very dependence of 10Be/9Be on heliocentric distance also implies that the comparatively younger chondrules (or other meteoritic material), which likely formed several AUs away from the Sun (so probably outside the range for in situ irradiation), need not have high 10Be/9Be, not to mention effects of finite time of transport or radioactive decay, as timescales for the chondrule formation epoch (0–3 Ma after CAIs; Nagashima et al. 2018; Connelly & Bizzarro 2018) or their transport are comparable to the half-life of 10 Be. For a Schmidt number of order unity, one might expect a 2 ± 1 orders of magnitude depletion with respect to the CAI-forming region (one order of magnitude closer to the Sun), but it is difficult to narrow down. Still, no chondrule 10Be/9Be values are known yet since silicate crystallization there hardly fractionates Be from B (Davis & McKeegan 2014).

thumbnail Fig. 5

LX for Taurusas a function of age (data from Güdel et al. 2007). Although a least-square fit (dashed line) shows a systematic increase with age, scatter is evident, in part owing to the diversity of stellar masses. Overplotted is the curve (red, continuous) corresponding to one solar mass (Telleschi et al. 2007; see beginning of Sect. 4.2). For reference, the normalizations used in Eq. (22) amount to LX = 107.7 J kg−1.

Open with DEXTER
thumbnail Fig. 6

LX for young stellar objects as a function of . Taurus and Orion data are taken from Güdel et al. (2007) and Bustamante et al. (2016), respectively (where LX also shows little systematic trend with stellar mass).

Open with DEXTER

4.3 Implications for vanadium isotopic systematics

While it is beyond the scope of the current work to calculate the effects of proton irradiation on other isotopic systems, some comments on the reinterpretation of the vanadium-50 excesses (up to 4.4o) reported by Sossi et al. (2017) are in order. In this framework, spallogenic 50V would also be dominated by inheritance from the nebula, rather than in situ production, but to a lesser extent as the target nuclei/element ratio (essentially Ti/V playing the role of O/Be; Sossi et al. 2017), between two refractory elements, would not markedly differ between a solar gas and a condensate (see second paragraph of Sect. 4.1). So the contention by Kööp et al. (2018a) that lack of evidence of in situ production of noble gases by solar energetic production in fine-grained CAIs analyzed by Vogel et al. (2004)7 prevents a spallogenic origin for their 50V excesses andfavors mass-dependent fractionation upon condensation is not necessarily valid. Since differential attenuation of the energetic particles modifies their energy distribution at finite penetration columns, the production ratio (independent of the fluence, and only dependent on p) of two nuclides would be affected. In particular, for a given p, since the reactions producing 50V are efficient at lower energies than for 10Be (Lee et al. 1998), where attenuation is strongest, 50V production would be comparatively lower. Thus, for given (measured) 10Be/9Be and δ51V, a steeper slope than in the formalism of Sossi et al. (2017) would be indicated (for a given target composition, which however would have to be revised to solar). Since, if we ratio Eq. (8) with its 50 V counterpart, the only difference with the Sossi et al. (2017) formalism is the presence of Σp (E) in the integrals in the denominator and the numerator, one may surmise that a shift of p comparable to the power-law exponent of Σp(E) in the Continuous Slowing Down Approximation (about 1.82; see Appendix B) would largely cancel out its effect and restore the observed ratio. Although a steeper p would diminish 10 Be production for a given L10, this would hardly affect our conclusions given the other uncertainties discussed at the end of Sect. 4.1. This also does not take into account the neutron contribution alluded to in Sect. 2.2, which may alter inferences on the energy distribution. Clearly, a dedicated study on the simultaneous effects expected for different isotopic systems in this scenario, whose proportions are independent of the uncertainties on the fluences, would be worthwhile.

4.4 X-rays and D/H fractionation

As mentioned above, X-ray emission should accompany cosmic ray flares from the Sun. These could also leave isotopic fingerprints (though unrelated to spallation) in meteorites. Indeed Gavilan et al. (2017) linked deuterium enrichment of chondritic organic matter (previously ascribed by Remusat et al. 2006, 2009 to ionizing radiation) to X-ray irradiation near the surface of the disk. The energy fluence F expected for organic matter wandering in the outer disk can be calculated with the same formalism as above if I replace by E, so I end this discussion with this short aside. For an attenuation column with n = 2.485 and an energy distribution with TX the emission temperature (Igea & Glassgold 1999), the energy-weighted average of the former (the equivalent of Kp ) is ΣX (kBTX)Γ(n + 1). For distant part of the (flared) disk, I can use f(R0, R) = HXR (with HX the X-ray absorption height, likely a few times H by analogy with Appendix C) so that the net fluence has the fairly well-constrained expression (modified after Eq. (19)) (23)

with the difference of HXR between the locus of formation of the organic matter (or the location inward of which settling ceases to be important, i.e. the S = 1 line of Jacquet et al. 2012; see Appendix C) and that of accretion in the chondrite parent body. This is five orders of magnitude short of the critical fluence of 5 × 1027 eV cm−2 = 8 × 1012 J m−2 indicated by the experiments of Gavilan et al. (2017) for 0.5–1.3 keV photons (and is also much shorter than the astrophysical estimate of Gavilan et al. 2017, which erroneously used X-ray attenuation at surface layers as representative of the bulk). While the steady-state approximation may overestimate the surface density in the outer regions of interest (which nevertheless should not allow efficient settling of the grains), alleviating it would hardly bridge the gap. Nevertheless, the quantitative assessment of isotopic effects of irradiation, whether X-rays (in particular at higher energies) or other parts of the spectrum, is still in its infancy and this additional calculation is essentially intended for future applications in similar scenarios.

5 Conclusion

I have analytically investigated the production of short-lived radionuclide beryllium-10 in surface layers of the disk irradiated by protosolar flares. I found that 10Be production in the gas outweighs 10Be production in solids after condensation because the gas contains a greater breadth of suitable target nuclides (e.g., 12 C, 1 H) and incurs less dilution by stable refractory Be. Taking into account incipient settling and possible implantation of solar wind-borne 10 Be would further widen the difference. Although many uncertainties remain on the magnetic field configuration, the scaling of cosmic rays with X-ray luminosities and other factors, it does appear that this model can reproduce 10Be/9Be ratios measured in CAIs. Therefore, the past presence of 10 Be does not require (at least at present) that CAIs formed at the inner edge of the disk and allows formation at a fraction of an AU, as thermal models would suggest, and more in line with evidence for a genetic link with their host carbonaceous chondrites (abundance,fraction of the refractory budget; Jacquet et al. 2012). If this model holds true, an interesting corollary (barring strong variations of the energetic protons/X-ray ratio) is that the oldest CAIs should have the lowest 10Be/9Be ratios, which would explain those of nucleosynthetic anomalies-bearing CAIs ((F)UN, PLAC). This would also suggest that the fine-grained group II CAIs in CV chondrites measured by Sossi et al. (2017) were a relatively late generation of refractory inclusions, a possibility that remains to be explored. This does not mean that chondrules, which formed at comparatively much larger heliocentric distances, should have high 10Be/9Be, since this ratio decreases outward, following passive diffusion outside the irradiated inner disk. I finally note that the same formalism allows an estimate of fluences of X-rays (produced in the same protosolar flares) or other types of radiations, on aggregates freely floating in the disk which can be, for example, compared to experimental evidence of D/H fractionation in meteoritic organic matter by irradiation.


I am grateful to the organizers and participants of the workshop “Core to Disk”, a program of the ψ2 initiative that took place at the Institut d’Astrophysique Spatiale in Orsay from May 14 to June 22, 2018, in particular Patrick Hennebelle and Matthieu Gounelle, who elicited my interest in the subject. I also thank Dr Ming-Chang Liu, Manuel Güdel, and Thomas Preibisch who kindly answered my requests for information. Comments by the anonymous reviewer greatly improved the clarity of the result sections as well as the discussion of other isotopic systems. This work is supported by ANR-15-CE31-004-1 (ANR CRADLE).

Appendix A Steady-state disk with infall

In this appendix, I study the effects of infall on the steady-state solution for the disk. I adopt the Cassen & Moosman (1981) expression for the infall source term per unit area for a cloud in solid-body rotation, assuming bipolar jets do not suppress it in the inner disk. The mass accretion rate is no longer constant but obeys (A.1)

with RC the instantaneous centrifugal radius, in the total infall rate on the star+disk system, and θ the Heaviside function. Treating all matter arriving inside R = R0 as directly accreted by the star, this can be integrated as (A.2)

with the mass accretion rate of the star. If the whole disk can be treated as in steady state, we should have = in but I refrain from making this identification immediately as the steady-state approximation may only hold locally in the inner regions and/or in the absence of infall (in = 0) altogether (as assumed in the main text).

Angular momentum conservation reads (Cassen & Moosman 1981) (A.3)

with M the stellar mass. If the latter is much bigger than the mass of the disk, the middle term on the right-hand side can be neglected upon integration (after division by R1∕2), which yields (A.4)

For RRC this approximates to (A.5)

which isthus a modest correction to Eq. (2).

For RR0 and = in, I have (A.6)

The factor between brackets decreases from 1 to 2/3 when R increases from 0 to RC. Beyond the centrifugal radius, = 0 and, from Eq. (A.3), R1∕2Σν becomes constant so that (A.7)

For bounded αT, 2πRΣ is not integrable, but this steady-state solution, which merely sets an upper bound to the disk viscous expansion, can only hold for distances of tvis (R) shorter than the timescale of evolution of the system8.

Appendix B Calculation of Kp

I take the attenuation column for protons as the minimum between the continuous slowing-down approximation (CSDA) and diffusion regimes of Padovani et al. (2018), after integrating their Eqs. (E.3) and (21)9, respectively: (B.2)

where A = 1.77 × 10−6 eV1+sm2, s = 0.82, αp = 1.28 (noted α in Padovani et al. 2018), σMT = 10−30 m2, Lp (E) = 10−21 m2 eV are defined in Padovani et al. (2018) and Γ is Euler’s gamma function. For the other species, I adopt , with Zi, Ai its atomic and mass numbers of species i (e.g., Gounelle et al. 2001). The transition between the CSDA and the diffusion regime thus occurs at (B.3)

Table B.1

Cross-section parameters.

As to the cross sections, similar to Desch et al. (2004), I adopt the Yanasak et al. (2001) formulation (B.4)

with σ0, Eth, E0, x constants (for a given (k, i)) and θ is again the Heaviside function. I fit those for i = 4He using the Lange et al. (1995) data (unknown to Yanasak et al. 2001), with power-law fitting below 30 MeV nucleon−1 and the flat regime being assigned the maximum cross section measured (rather than a fit beyond the aforementioned threshold, as this regime is only incipient in the energy range studied). The resulting parameters are presented in Table B.1.

Then I have for Eth ≤ E0Ecrit (B.5)

and for EthEcritE0 (B.6)

whose contributions can be summed to yield Kp (Eq. (10)).

Appendix C Height ofspallation layer and effect of settling on 10Be production in solids

The spallation layer column Σspall may be defined as (C.1)

Assuming a vertically isothermal disk (and thus a gaussian gas density ), the height Hspall of this layer obeys (C.2)

For Σ ~ 105±1 kg m−2 and Σspall = 101±1 kg m−2, I have χHspallH ≈ 3 ± 1, so that (C.3)

Since, for x ≳ 2, I have10 (C.4)

the density at z = Hspall may be approximated as (C.5)

The nominal value (corresponding to a number density of 2 × 1012 cm−3) is halfway in the range studied by Umebayashi & Nakano (1981) so the neutron contribution (ignored in this paper) can be seen in their Figs. 2 and 5 to be limited. Then, the settling parameter (Jacquet et al. 2012) for a solid of density ρs and radius a there is (C.6)

with δz the vertical diffusion coefficient normalized to (which should be of order α).

Simulations of ideal MHD turbulence indicate δz ~ αρ−1 over a large extent of the disk thickness, although it should stall in the corona (e.g., Jacquet 2013, and references therein). Since Sz may be expected to be ≲1 there from the above evaluation, it may not be a dramatic underestimate to treat it as constant throughout the vertical thickness of the disk. Under that assumption, a population of identical solids has a density (Jacquet 2013) (C.7)

with Σc the surface density of the population. The 10Be production rate averaged over this distribution is then (C.8)

with and Σ′ the columns of the solids and gas, respectively, integrated vertically from the surface. If I let (C.9)

and since (C.10)

this becomes (C.11)

where I have used approximation (C.4), with the corollary for z ≪ 1. Roughly speaking, the integrand only contributes (~1) for a g argument ≲ 1, that is so that (C.12)

So 10Be production should quickly dropas Sz approaches unity and become negligible for Sz ≳ 1.



Assuming symmetry of irradiation about the midplane, although this is not a prerequisite to the final result (Eq. (9)).


If we follow a field line from (r, θ) = (R0, θ0) in polar coordinates, we have (11)

so that upon crossing the midplane r = R = R0∕sin2 θ0 with f(R0, R) = cos θ0.


Taking into account the factor of four between the mean flux of 10Be on a randomly oriented surface ((2.9 ± 1.2) × 10−2 m−2s−1 on the Moon) and its omnidirectional flux in vacuo.


It may also be noted that since the portion of tvis spent within a column Σspall of the surface (18)

is much longer than the weekly periodicity of protostellar flares (Wolk et al. 2005), the relevant proton luminosity L10 in Eq. (9) must be the long-term average luminosity rather than the “characteristic” baseline or the typical flare peak value.


It may nonetheless be wondered whether this average of SIMS points does not underestimate the bulk concentrations given the incompatible behavior of Be in melilite (Paque et al. 2014).


Of course, this may shorten the timescales required by Sossi et al. (2017) at the disk inner edge as well, but they would still represent dozens of orbits.


These incidentally do show evidence for spallogenic He and Ne, even though it is unclear there is any excess over that due to the recent exposure of the meteoroid to galactic cosmic rays.


Mathematically, inclusion of the hitherto neglected middle term of the right-hand side of Eq. (A.3) would suffice to make the disk mass finite but on a too large radial scale for it to be relevant.


For the latter, after a change of variable E0 = E∕sinθ, I use the identity

(for derivation, multiply the left-hand side by and rewrite the resulting double integral in Cartesian coordinates).


This is because erfc.

All Tables

Table B.1

Cross-section parameters.

All Figures

thumbnail Fig. 1

Sketch of the scenario of 10Be production developed in this paper. Solar cosmic rays cause spallation reactions in the upper layers of the disk, in the gas and/or condensates, within a region (outlined by a dashed separatrix) where magnetic field lines originating from the proto-Sun and its outskirts hit the disk. The field and particle emission geometries are purely illustrative and do not detract from the generality of the formalism.

Open with DEXTER
In the text
thumbnail Fig. 2

Kp as a function of the power-law exponent of the incoming cosmic ray energy distribution, for a solar (without possible implantation contributions) and a rocky target (see Sects. 3.1 and 3.2, respectively).

Open with DEXTER
In the text
thumbnail Fig. 3

feff (Eq. (17)) as a function of heliocentric distances for two values of Rmax for my toy parameterization of the particle flux distribution (see end of Sect. 2.2) and two values of the radial Schmidt number, with the lower one (higher relative diffusivity) in dotted lines. The parameter feff is proportional to the 10Be/9Be ratio (whose absolute magnitude is discussed at Eq. (22)) in the bulk (gas+solids) disk for given mass accretion rate and proton output and energy distributions.

Open with DEXTER
In the text
thumbnail Fig. 4

f(R, +) as a functionof heliocentric distances for two values of Rmax for my toy parameterization of the particle flux distribution (see end of Sect. 2.2). This is proportional to the 10Be/9Be ratio produced by in situ irradiation of a refractory solid as a function of its heliocentric distance of origin.

Open with DEXTER
In the text
thumbnail Fig. 5

LX for Taurusas a function of age (data from Güdel et al. 2007). Although a least-square fit (dashed line) shows a systematic increase with age, scatter is evident, in part owing to the diversity of stellar masses. Overplotted is the curve (red, continuous) corresponding to one solar mass (Telleschi et al. 2007; see beginning of Sect. 4.2). For reference, the normalizations used in Eq. (22) amount to LX = 107.7 J kg−1.

Open with DEXTER
In the text
thumbnail Fig. 6

LX for young stellar objects as a function of . Taurus and Orion data are taken from Güdel et al. (2007) and Bustamante et al. (2016), respectively (where LX also shows little systematic trend with stellar mass).

Open with DEXTER
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.