HETEROGENEOUS DISTRIBUTION OF 26Al AT THE BIRTH OF THE SOLAR SYSTEM

, , , , , , , and

Published 2011 May 10 © 2011. The American Astronomical Society. All rights reserved.
, , Citation Kentaro Makide et al 2011 ApJL 733 L31 DOI 10.1088/2041-8205/733/2/L31

2041-8205/733/2/L31

ABSTRACT

It is believed that 26Al, a short-lived (t1/2 = 0.73 Ma) and now extinct radionuclide, was uniformly distributed in the nascent solar system (SS) with the initial 26Al/27Al ratio of ∼5.2 × 10−5, suggesting an external, stellar origin rather than local, solar source. However, the stellar source of 26Al and the manner in which it was injected into the SS remain controversial: the 26Al could have been produced by an asymptotic giant branch star, a supernova, or a Wolf–Rayet star and injected either into the protosolar molecular cloud, protosolar cloud core, or protoplanetary disk. Corundum (Al2O3) is predicted to be the first condensate from a cooling gas of solar composition. Here we show that micron-sized corundum condensates from 16O-rich (Δ17O ∼ −25‰) gas of solar composition recorded heterogeneous distribution of 26Al at the birth of the SS: the inferred initial 26Al/27Al ratio ranges from ∼6.5×10−5 to <2×10−6; 52% of corundum grains measured are 26Al-poor. Abundant 26Al-poor, 16O-rich refractory objects include grossite- and hibonite-rich calcium–aluminum-rich inclusions (CAIs) in CH (high metal abundance and high iron concentration) chondrites, platy hibonite crystals in CM (Mighei-like) chondrites, and CAIs with fractionation and unidentified nuclear effects CAIs chondrites. Considering the apparently early and short duration (<0.3 Ma) of condensation of refractory 16O-rich solids in the SS, we infer that 26Al was injected into the collapsing protosolar molecular cloud and later homogenized in the protoplanetary disk. The apparent lack of correlation between 26Al abundance and O-isotope composition of corundum grains constrains the stellar source of 26Al in the SS.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

26Al (which decays to 26Mg with a half-life, t1/2 = 0.73 Ma), 10Be (t1/2 = 1.5 Ma), 36Cl (t1/2 = 0.3 Ma), 41Ca (t1/2 = 0.1 Ma), 53Mn (t1/2 = 3.7 Ma), and 60Fe (t1/2 = 2.3 Ma) are short-lived radionuclides that were present in the early solar system (SS; McKeegan & Davis 2003). The recently inferred uniform distribution of 26Al in the early SS (Jacobsen et al. 2008; Villeneuve et al. 2009; Davis et al. 2010) makes the 26Al–26Mg system one of the most important relative chronometers of the early SS processes, such as the evaporation, condensation, and melting of solids in the protoplanetary disk (McKeegan & Davis 2003; Kita et al. 2005; Krot et al. 2009). The maximum amount of 26Al that could have been produced by solar energetic particle irradiation is insufficient to explain the canonical 26Al/27Al ratio (∼5.2×10−5) and suggests a stellar, nucleosynthetic origin of 26Al (Duprat & Tatischeff 2007). The stellar source of 26Al and the mechanism by which it was introduced into the SS, however, remain controversial. Possible sources of freshly synthesized 26Al are a neighboring asymptotic giant branch (AGB) star (Wasserburg et al. 2006), a supernova (SN; Boss et al. 2008), or a wind from a Wolf–Rayet (WR) star (Arnould et al. 2006; Gaidos et al. 2009; Tatischeff et al. 2010), any of which could have injected 26Al into the parent molecular cloud (Gaidos et al. 2009), protosolar cloud core (Boss et al. 2008), or protoplanetary disk (Ouellette et al. 2007).

AGB stars are not associated with star-forming regions and the probability of an encounter between an AGB star and a star-forming molecular cloud is very low, ∼2.5×10−6 (Huss et al. 2009). Although massive SN progenitors are commonly associated with star-forming regions (Lada & Lada 2003), the probability of injection of 26Al into the protosolar cloud core or protoplanetary disk from a nearby SN is also low, ⩽3×10−3 (Williams & Gaidos 2007; Gounelle & Meibom 2008). To overcome this problem, it has been recently suggested that the Sun was born in a stellar cluster of second generation that was self-enriched in short-lived radionuclides (60Fe and possibly 26Al) by a single or multiple 10–15 solar mass (M) SN several million years before our SS formed (Gounelle et al. 2009). Finally, it was proposed that 26Al could have been added to the progenitor molecular cloud by WR wind (Arnould et al. 2006; Gaidos et al. 2009; Tatischeff et al. 2010).

CAIs are the oldest SS solids dated (Amelin et al. 2002; Bouvier & Wadhwa 2010) and are believed to have formed when the Sun was a class 0 or class I protostar (e.g., Krot et al. 2009; Tscharnuter et al. 2009), although Ciesla (2010) argues that most CAIs found in meteorites formed at the time that the solar nebula stopped accreting significant amounts of material from the protosolar parent cloud, and thus corresponds roughly to the class I–class II transition. The majority of CAIs in unmetamorphosed chondrites have the initial 26Al/27Al ratio ((26Al/27Al)0) of ∼(4–5)×10−5 (MacPherson et al. 1995; Jacobsen et al. 2008; Makide et al. 2009b). CAIs with low (26Al/27Al)0, ≪5×10−6, are very rare (Makide et al. 2009b; Sahijpal & Goswami 1998; Liu et al. 2009), except in CH chondrites where they are dominant (Krot et al. 2008). CH chondrites escaped thermal metamorphism and aqueous alteration, and are considered to be one of the most primitive meteorites in our collections. These 26Al-poor CAIs are interpreted either to be the result of early formation, prior to injection and homogenization of 26Al in the SS (Makide et al. 2009b; Sahijpal & Goswami 1998; Liu et al. 2009; Krot et al. 2008) or as a result of preferential loss of 26Al carrier from CAI precursors by sublimation (Thrane et al. 2008). The late-stage formation of 26Al-poor CAIs, after decay of 26Al, can also not be excluded, because CAIs are known to have experienced multistage thermal processing in the solar nebula and on chondrite parent asteroids (Krot et al. 2005) that could have erased possible evidence for radiogenic 26Mg (MacPherson et al. 1995).

Corundum (Al2O3) is predicted to be the first condensate from a cooling gas of solar composition (Ebel & Grossman 2000). Therefore, corundum condensates can potentially provide additional constraints on the origin of 26Al and degree of its heterogeneity in the nascent SS. The pioneering work on solar corundum grains was done by Virag et al. (1991), who studied 26 individual corundum grains 3–20 μm grains from the unmetamorphosed CM carbonaceous chondrite (CC) Murchison. On the basis of the oxygen and magnesium-isotopic compositions, the grains were divided into three groups. Group 1 (n = 17) and group 2 (n = 5) grains show 26Mg excesses corresponding to (26Al/27Al)0 of 5×10−5 and 5×10−6, respectively; group 3 grains (n = 4) show no resolvable 26Mg excesses ((26Al/27Al)0 <3×10−7).

To understand the origin and distribution of 26Al in the early SS, we measured 26Al–26Mg isotope systematics and oxygen-isotope compositions of micron-sized corundum grains from acid-resistant residues of unequilibrated ordinary chondrites (UOCs; H (high total iron contents) and LL (low metallic iron and low total iron contents) of low petrologic type, ⩽3.15) and unmetamorphosed CCs (CI1 (Ivuna-like), CM2, CR2 (Renazzo-like), and CO3.0 (Ornans-like)) using the University of Hawai'i (UH) Cameca ims-1280 ion microprobe. Oxygen isotopes were used to distinguish solar from presolar corundum grains; the latter are characterized by isotopically anomalous oxygen compositions (Nittler 2003) compared to the solar and terrestrial values, Δ17O ∼ −25‰ and 0‰, respectively (Makide et al. 2009b; McKeegan et al. 2010), where Δ17O = δ17O−0.52×δ18O, δ17,18O = ((17,18O/16O)sample)/(17,18O/16O)SMOW−1)×1000, and SMOW is the Standard Mean Ocean Water.

2. ANALYTICAL PROCEDURES

2.1. Sample Preparation

Acid-resistant residues of ordinary and CCs used in this study were prepared by Huss & Lewis (1995) and Huss et al. (2003). The residues consist of SiC, spinel (MgAl2O4), hibonite (CaAl12O19), and corundum and range in size from 0.5 to 5 μm. The residues were diluted by a mixture of 90% isopropanol and 10% water and mixed in an ultrasonic bath. About 0.2 μl of the mixture was siphoned using a micro-aspirator and dispersed onto a gold substrate under a stereomicroscope. More than 100 grains were dispersed on each gold substrate. Prior to the dispersion, the gold substrate was examined for possible contamination using the UH JEOL JXA-8500F field emission electron microprobe (FE-EPMA) equipped with a cathodoluminescence detector (CL) and an energy dispersive spectrometer (EDS). Several relatively large (>5 μm) corundum grains embedded in the substrates were identified and avoided during subsequent oxygen- and magnesium-isotopic measurements. To locate the chondritic corundum grains comprising less than 1% of the dispersed grains dominated by spinel, we used CL and EDS of the FE-EPMA. Among the dispersed grains, only corundum and SiC show intense CL; the latter is very rare and can easily be distinguished by EDS. The corundum grains identified were imaged in secondary and backscattered electrons.

2.2. Oxygen- and Magnesium-isotope Measurements

Oxygen-isotope compositions of individual corundum grains were measured with the UH Cameca ims-1280 (for details see Makide et al. 2009a).

Aluminum- and magnesium-isotope measurements with UH Cameca ims-1280 were conducted on the corundum grains which have been previously analyzed for oxygen isotopes. An 16O primary beam, either defocused (∼30 μm spotsize) or focused (∼3 μm), was used to generate secondary Mg+ and Al+ ions. A field aperture of 1000×1000 μm2 corresponding to ∼7 μm on the sample was used to minimize the contribution of Mg signals from the substrate and other dispersed grains surrounding the grain of interest. The contribution of the Mg signal from the substrate was estimated to be <1% of Mg from the sample. The mass resolving power was set to ∼3500–3800, sufficient to separate interfering hydrides and 48Ca++. Secondary ions of 24Mg+, 25Mg+, and 26Mg+ were measured simultaneously using the monocollector electron measure (EM) and two multicollector EMs; 27Al+ was subsequently measured with the monocollector EM by peak-jumping. The primary beam current was adjusted so that the 27Al+ count rate did not exceed 4.5×105 cps. Acquisition times for Mg isotopes and 27Al were 15 s and 1 s, respectively. Each measurement consisted of 20 cycles. Excess of radiogenic 26Mg (δ26Mg*) was calculated using the following equation: δ26Mg* = Δ26Mg − 2×Δ25Mg, where Δ25,26Mg = ((25,26Mg/24Mg)sample)/(25,26Mg/24Mg)ref−1) × 1000 and reference Mg-isotope ratios are from Catanzaro et al. (1966). We note that Δ25,26Mg values wre calculated using the mean of ratios for all cycles. Since Mg count rates are low, it results in bias in 25Mg/24Mg and 26Mg/24Mg ratios (Ogliore et al. 2011). The typical reproducibility of standard measurements was ∼10‰ (2σ) in Δ25Mg and Δ26Mg. To correct for instrumental mass fractionation effects and estimate the relative sensitivity factor for the 27Al/24Mg ratio, we used a Montana sapphire (α-Al2O3) standard (NMNH 126321 from the Smithsonian Institution). The instrumental fractionation correction was made based on the measurements of micron-sized Montana sapphire grains dispersed on the gold substrate. To correct the 27Al+/24Mg+ ratios measured by SIMS in micron-sized corundum grains, we measured the Mg concentration in the Montana sapphire grain by the UH FE-EPMA; a piece of the same grain was also measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS) by E. D. Young (UCLA). To measure trace amount of Mg in the sapphire by FE-EPMA, we used defocused (10 μm) electron beam at ∼300 nA current and 15 keV accelerating voltage. Madagascar hibonite was used as a standard. The 27Al/24Mg ratios in the sapphire obtained by FE-EPMA and LA MC-ICP-MS are 6084±498 (2σ, 45 spots) and 5830 ± 510 (2σ, 4 spots), respectively, and are the same within the uncertainties and show no evidence for heterogeneity.

3. RESULTS AND DISCUSSION

3.1. Oxygen- and Magnesium-isotope Compositions of Micron-sized Corundum Grains

Oxygen-isotope compositions of micron-sized corundum grains from acid-resistant residues of UOCs (Semarkona (LL3.0), Bishunpur (LL3.1), Roosevelt County (RC) 075 (H3.1)) and unmetamorphosed CCs (Orgueil (CI1), Murray (CM2), Renazzo (CR2), and Allan Hills (ALH) A77307 (CO3.0)) together with compositions of the mineralogically pristine CAIs from CR2 CCs (Makide et al. 2009b) and the solar wind returned by the Genesis spacecraft (McKeegan et al. 2010) are listed in online Table 1 and plotted in Figure 1. All but three corundum grains have 16O-rich compositions similar to those of CR CAIs and solar wind, consistent with a condensation origin from an 16O-rich gas of solar composition. Three corundum grains, not shown in Figure 1, have highly anomalous O-isotope compositions, and are probably presolar in origin.

Figure 1.

Figure 1. Δ17O values of micron-sized corundum grains from unequilibrated ordinary chondrites (UOCs) and unmetamorphosed carbonaceous chondrites (CCs), CAIs from CR chondrites (Makide et al. 2009b), and solar wind (McKeegan et al. 2010). Errors are 2σ. The 16O-rich composition of corundum grains similar to those of CR CAIs and solar wind is consistent with a condensation origin from a gas of solar composition. Petrologic types of ordinary chondrites are from Grossman & Brearley (2005).

Standard image High-resolution image

Table 1. Oxygen-isotope Compositions (in ‰) and the Inferred (26Al/27Al)0 Ratios from Individual Micron-sized Corundum Grain in UOCs and Unmetamorphosed CCs

Grain δ18O δ17O Δ17O (26Al/27Al)0 × 10−5
Carbonaceous            
chondrites            
ALHA 03-03 −35.9 2.3 −43.4 4.5 −24.7 4.7 N.D.  
...                

Notes. Sem, Semarkona; Bis, Bishunpur; RC075, Roosevelt County 075; ALHA, Alan Hills A77307; Murr, Murray; Org, Orgueil; Ren, Renazzo; *P, presolar grains; *U, isotopically unusual grains, possibly an analytical artifact. The *P and *U grains are not plotted in Figure 1. N.D. = not determined. ^grains with O-isotope compositions reported by Makide et al. (2009a).

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

The solar corundum grains were subsequently measured for magnesium-isotope compositions and 27Al/24Mg ratio; the results are listed in Table 1 and plotted in Figure 2. The grains show large variations of the inferred initial 26Al/27Al ratio ((26Al/27Al)0): 52% of grains have no resolvable 26Mg*: an upper limit on (26Al/27Al)0 is 2×10−6; 40% of grains have high (26Al/27Al)0, (3.0–6.5)×10−5; 8% of grains have intermediate values of (26Al/27Al)0, (1–2)×10−5.

Figure 2.

Figure 2. (a) Aluminum–magnesium evolutionary isotope diagram for the micron-sized corundum grains #02–01 and #01–10 from Orgueil. The grain #02–01 has excesses of 26Mg corresponding to the initial 26Al/27Al ratio ((26Al/27Al)0) of (4.7±1.0)×10−5. No resolvable excess of 26Mg is found in the grain #01–10; an upper limit on (26Al/27Al)0 is 4.8×10−7. Errors are 2σ. (b) Δ17O values vs. (26Al/27Al)0 in corundum grains from UOCs and CCs. Errors are 2σ. There is a bi-modal distribution of (26Al/27Al)0: 52% of grains have no detectable excess of 26Mg (an upper limit on (26Al/27Al)0 is 2×10−6); 40% of grains have large excesses of 26Mg corresponding to (26Al/27Al)0 from ∼3.0×10−5 to ∼6.5×10−5; 8% of grains have intermediate values of (26Al/27Al)0. There is no correlation between O- and 26Al–26Mg systematics of corundum grains. In contrast to the corundum grains, the (26Al/27Al)0 in whole-rock CV CAIs shows small variations, (5.23 ± 0.13)×10−5. (c) The distribution of (26Al/27Al)0 in solar corundum grains.

Standard image High-resolution image

The coexistence of the 26Al-rich and 26Al-poor corundum grains in the same primitive meteorite precludes late-stage (after decay of 26Al) resetting of the 26Al–26Mg systematics of the 26Al-poor corundum grains during thermal metamorphism in the host chondrite parent body. Late-stage resetting in the solar nebula during thermal processing associated with chondrule formation is also excluded. Chondrule formation occurred in an 16O-poor (Δ17O > −10‰) nebular reservoir (e.g., Krot et al. 2006) when 26Al was still alive (Kita et al. 2005; Kurahashi et al. 2008; Krot et al. 2009). We infer that 16O-rich corundum grains with low initial 26Al/27Al ratios never contained the canonical abundance of 26Al, i.e., the lack of 26Mg excess in 52% of 16O-rich corundum grains is a primary characteristic.

3.2. Multiple Generations of CAIs and Refractory Grains and Their Distribution in the Protoplanetary Disk

A high abundance (>50%) of 26Al-poor, 16O-rich CAIs has been previously reported in CH CCs (Krot et al. 2006). Similar to corundum grains, the CH CAIs are among the smallest (<20 μm) and the most refractory (rich in hibonite and grossite, CaAl4O7) objects known (Krot et al. 2002). The 26Al heterogeneity and high abundance of 26Al-poor, 16O-rich refractory objects in a fine fraction of solar dust can be explained in the context of the formation setting of refractory objects (MacPherson et al. 1995; Krot et al. 2009) and theoretical modeling of their radial transport in the protoplanetary disk (Boss 2008; Ciesla 2009, 2010; Yang & Ciesla 2010). Refractory objects are thought to have formed in region(s) with high ambient temperature (at or above condensation temperature of forsterite (Mg2SiO4), ∼ 1350 K at total pressure of <10−4 bar), most likely within 1–2 astronomical units from the Sun (Krot et al. 2009, and references therein). They were subsequently transported outward to the accretion regions of chondrite parent asteroids and cometary nuclei (Gounelle & Meibom 2007; Boss 2008; Ciesla 2009, 2010). The oldest 207Pb–206Pb absolute ages of CAIs from CV (Vigarano-like) chondrites (Amelin et al. 2002; Bouvier & Wadhwa 2010), a narrow range of their initial 26Al/27Al ratios (Jacobsen et al. 2008), corresponding to a range of crystallization ages of <0.1 Ma, and the high ambient temperature in the CAI-forming region(s), are all consistent with an early and brief duration of CAI formation, during a period of high mass-accretion rate of dust and gas to the proto-Sun (∼10−6 M yr−1; Krot et al. 2009; Ciesla 2009). We propose that the 26Al-poor and 26Al-rich 16O-rich corundum grains represent different generations of refractory objects formed during a brief epoch of CAI formation, and, therefore, reflect heterogeneous or variable 26Al at the birth of the SS.

To understand the distribution of refractory objects with the different (26Al/27Al)0 in the early SS, we modeled the evolution of a viscous disk and redistribution of refractory objects formed in the high-temperature (above 1400 K) region within it during the first million years of its evolution (Figure 3). We consider both the mass and angular momentum transport in the disk as well as the addition of mass to both the Sun and the disk due to infall from the parent molecular cloud (Yang & Ciesla 2010; Hueso & Guillot 2005). The parent cloud is assumed to be 1 M, spherical, at a temperature of 15 K, rotating with a uniform angular velocity of Ωc = 10−14 s−1, and we adopt a mass accretion rate of ∼3×10−6 M yr−1; the disk viscosity is parameterized with a turbulent viscosity coefficient of α = 10−3. The accreted mass falls onto the disk at locations where the specific angular momentum in the disk equals that in the parent cloud (Yang & Ciesla 2010; Hueso & Guillot 2005). We track populations of fine refractory objects as they are redistributed throughout the disk by viscous expansion, dividing them into batches depending on the timing of when they are last exposed to the high-temperature region (Figures 3(a) and (b)).

Figure 3

Figure 3 (a) Masses of the solar mass star, disk, and star+disk as a function of time. (b) Fraction of dust processed at high temperature (> 1400 K) and survived in a 1 Myr old disk, plotted in different batches (numbered from 1 to 12) according to when they were last exposed to those high temperatures.

Standard image High-resolution image

We find that the mass of the disk reaches its maximum by 300,000 years of its evolution (Figure 3(a)), and that the most dominant population of refractory objects in a 1 million year old disk consists of those that formed immediately before and after cessation of infall from the parent molecular cloud (Figure 3(b)). This prediction, combined with the observed heterogeneity of (26Al/27Al)0, suggests that 26Al was introduced into the SS contemporaneously with the collapse of the protosolar cloud core; injection into the disk (Ouellette et al. 2007) is much less likely. Alternatively, if 26Al was heterogeneously distributed in the parent molecular cloud, mixing in the disk did not homogenize it until the infall and formation of refractory solids had ceased. Thus, refractory objects formed prior to the end of collapse would not record the average solar nebula 26Al/27Al ratio, but rather the ratio that existed in the hot nebular region at the time of their formation. This would allow refractory objects with different initial 26Al/27Al ratios (26Al-poor and 26Al-rich) to form and to be preserved in the disk. It also implies that the initial 26Al/27Al ratio of (5.23 ± 0.13)×10−5 recorded by CV CAIs (Jacobsen et al. 2008; Davis et al. 2010) may not necessarily represent the initial 26Al/27Al ratio in the SS.

Our modeling indicates that a significant population of refractory objects that formed prior to the cessation of infall survives in the disk (Figure 3(b)). The survival time of dust particles (up to ∼1 m in size) in the disk is inversely proportional to their sizes: fine (5 μm) dust is well-coupled to the nebular gas and less affected by head wind than the coarse (∼1 mm) dust (Ciesla 2010). These observations may explain the high proportion of 26Al-poor refractory objects in a fine fraction of solar dust, e.g., as we have observed in corundum grains and CH CAIs, as well as the presence of an 26Al-poor, 16O-rich CAI in comet Wild2 (Matzel et al. 2010).

3.3. Stellar Source of  26Al

The heterogeneous distribution of 26Al at the birth of the SS may preclude inheritance of 26Al from a molecular cloud enriched by a previous generation of stars (Gounelle et al. 2009), which is expected to have a uniform distribution of 26Al. Instead, it indicates injection of 26Al nearly contemporaneously with the collapse molecular cloud core that formed the SS. The similar oxygen-isotope compositions of the 26Al-rich and 26Al-poor corundum grains and the Sun suggest injection of 26Al by a massive (>20 M) Type II SN, by a wind from a WR star or by a low-mass (1.5 M) AGB star (Sahijpal & Soni 2006; Krot et al. 2008; Ellinger et al. 2010). Injection of 26Al by a less massive SN or by a 3–6 M AGB star is unlikely, because it is expected to produce significant change in oxygen-isotope composition of the SS that has not been observed (Sahijpal & Soni 2006; Gounelle & Meibom 2007; Ellinger et al. 2010).

To distinguish between an explosive SN and a WR wind of 26Al, measurements of other isotopes (e.g., 60Fe–60Ni isotope systematics) of 26Al-rich and 26Al-poor 16O-rich CAIs or refractory grains are required: in contrast to SNe, the WR stars do not produce 60Fe (Sahijpal & Soni 2006). Therefore, if 26Al was injected by a WR star (Gaidos et al. 2009; Tatischeff et al. 2010) and 60Fe was inherited from a molecular cloud of an earlier generation (Gounelle et al. 2009), the two isotopes would be decoupled, i.e., the 26Al-rich and 26Al-poor CAIs will have similar initial abundance of 60Fe. In contrast, if 26Al and 60Fe were injected by the same SN, the 26Al-poor CAIs would have lower initial abundance of 60Fe than the 26Al-rich CAIs.

We note that the initial abundance of 60Fe in the SS remains controversial. The inferred initial 60Fe/56Fe ratios range from (2–5)×10−7 (Telus et al. 2011) to <1×10−8 (Tang & Dauphas 2011). Due to the presence of nucleosynthetic isotopic anomalies in Ni, the detection of radiogenic 60Ni in bulk CAIs is difficult (Quitté et al. 2007) and may require high-precision Fe–Ni isotope measurements of CAI mineral separates to construct an internal Fe–Ni isochron.

Constructive review by R. Gallino is highly appreciated. We thank E. Young (UCLA) for LA-ICPMS measurements of sapphire standard kindly provided by G. MacPherson (Smithsonian Institution). This research has received support from NASA grants NNX08AH91G and NNX07AI81G (A.N.K.).

Please wait… references are loading.
10.1088/2041-8205/733/2/L31