26Al as a heat source for early melting of planetesimals: Results from isotopic studies of meteorites

Ion microprobe studies of magnesium isotopic composition in igneous components from several chondritic meteorites have been carried out to look for²⁶Mg excess that may be attributed to the presence of the now-extinct radionuclide²⁶Al(τ ∼ 1 Ma) at the time of formation of these objects. A positive evidence for the presence of²⁶Al in the analysed objects will strengthen its case as the primary heat source for the early thermal metamorphism/melting of meteorite parent bodies. Based on calculated temperature profiles inside chondritic objects of different sizes and initial²⁶Al/²⁷Al ratios, we have estimated the initial abundances of²⁶Al needed to provide the heat necessary for the wide range of thermal processing seen in various types of meteorites. The magnesium isotopic data obtained by us do not provide definitive evidence for the presence of²⁶Al at the time of formation of the analysed igneous phases in different chondritic meteorites. Experimental evidence for a planetary scale distribution of²⁶Al in the early solar system to serve as a significant heat source for the thermal metamorphism and melting of meteorite parent bodies (planetesimals) remains elusive.     

Thermal evolution of the early Moon

The early thermal evolution of Moon has been numerically simulated to understand the magnitude of the impact-induced heating and the initially stored thermal energy of the accreting moonlets. The main objective of the present study was to understand the nature of processes leading to core–mantle differentiation and the production and cooling of the initial convective magma ocean. The accretion of Moon was commenced over a time scale of 100 yr after the giant impact event around 30–100 million years in the early solar system. We studied the dependence of the planetary processes on the impact scenarios, the initial average temperature of the accreting moonlets, and the size of the protomoon that accreted rapidly beyond the Roche limit within the initial 1 yr after the giant impact. The simulations indicate that the accreting moonlets should have a minimum initial averaged temperature around 1600 K. The impacts would provide additional thermal energy. The initial thermal state of the moonlets depends upon the environment prevailing within the Roche limit that experienced episodes of extensive vaporization and recondensation of silicates. The initial convective magma ocean of depth more than 1000 km is produced in the majority of simulations along with the global core–mantle differentiation in case the melt percolation of the molten metal through porous flow from bulk silicates was not the major mode of core–mantle differentiation. The possibility of shallow magma oceans cannot be ruled out in the presence of the porous flow. Our simulations indicate the core–mantle differentiation within the initial 10² to 10³ yr of the Moon accretion. The majority of the convective magma ocean cooled down for crystallization within the initial 10³ to 10⁴ yr.     

Inhomogeneous Chemical Evolution of the Galaxy in the Solar Neighbourhood

N-body numerical simulations of an inhomogeneous Galactic Chemical Evolution (GCE) of the solar neighbourhood with a high temporal resolution are presented. The solar annular ring is divided into distinct spatial grids of area ~1–2 kpc². Each grid evolves distinctly in terms of star formation and nucleosynthetic yields from numerous generations of stars. The evolution of the galaxy is simulated by considering discrete episodes of star formation. Subsequent to the evolution of the simulated stars within each grid the stellar nucleosynthetic yields are homogenized within the grid rather than the traditionally adopted criteria of homogenizing over the entire solar annular ring. This provides a natural mechanism of generating heterogeneities in the elemental abundance distribution of stars. A complex chemical evolutionary history is inferred that registers episodes of time-dependent contributions from SN II+Ib/c with respect to SN Ia. It was observed that heterogeneities can remerge even after episodes of large scale homogenizations on scales larger than the grid size. However, a comparison of the deduced heterogeneities with the observed scatter in the elemental abundances of the dwarf stars suggest only a partial match, specifically, for [Fe/H] > ?0.5. The deduced heterogeneities in the case of carbon, oxygen, magnesium, silicon, sulphur, calcium and titanium can explain the observed heterogeneities for [Fe/H] < ?0.5. It may not be possible to explain the entire observed spread exclusively on the basis of the inhomogeneous GCE.     

Determination of rare earth and refractory trace element abundances in early solar system objects by ion microprobe

Experimental and analytical procedures devised for measurement of rare earth element (REE) abundances using a secondary ion mass spectrometer (ion microprobe) are described. This approach is more versatile than the conventional techniques such as neutron activation analysis and isotope dilution mass spectrometry by virtue of its high spatial resolution that allows determination of REE abundances in small domains (10-20 micron) within individual mineral phases. The ion microprobe measurements are performed at a low mass-resolving power adopting the energy-filtering technique (Zinner and Crozaz 1986) for removal and suppression of unresolved complex molecular interferences in the REE masses of interest. Synthetic standards are used for determining various instrument specific parameters needed in the data deconvolution procedure adopted for obtaining REE abundances. Results obtained from analysis of standards show that our ion microprobe may be used for determining REE abundances down to ppm range with uncertainties of ∼ 10 to 15%. Abundances of rare earth and several other refractory trace elements in a set of early solar system objects isolated from two primitive carbonaceous chondrites were determined using the procedures devised by us. The results suggest that some of these objects could be high temperature nebular condensates, while others are products of melting and recrystallization of precursor nebular solids in a high temperature environment.     

Did 26Al and impact‐induced heating differentiate Mercury?

Numerical models dealing with the planetary scale differentiation of Mercury are presented with the short‐lived nuclide, ²⁶Al, as the major heat source along with the impact‐induced heating during the accretion of planets. These two heat sources are considered to have caused differentiation of Mars, a planet with size comparable to Mercury. The chronological records and the thermal modeling of Mars indicate an early differentiation during the initial ~1 million years (Ma) of the formation of the solar system. We theorize that in case Mercury also accreted over an identical time scale, the two heat sources could have differentiated the planets. Although unlike Mars there is no chronological record of Mercury's differentiation, the proposed mechanism is worth investigation. We demonstrate distinct viable scenarios for a wide range of planetary compositions that could have produced the internal structure of Mercury as deduced by the MESSENGER mission, with a metallic iron (Fe‐Ni‐FeS) core of radius ~2000 km and a silicate mantle thickness of ~400 km. The initial compositions were derived from the enstatite and CB (Bencubbin) chondrites that were formed in the reducing environments of the early solar system. We have also considered distinct planetary accretion scenarios to understand their influence on thermal processing. The majority of our models would require impact‐induced mantle stripping of Mercury by hit and run mechanism with a protoplanet subsequent to its differentiation in order to produce the right size of mantle. However, this can be avoided if we increase the Fe‐Ni‐FeS contents to ~71% by weight. Finally, the models presented here can be used to understand the differentiation of Mercury‐like exoplanets and the planetary embryos of Venus and Earth.     

Stellar nucleosynthetic contribution of extinct short‐lived nuclei in the early solar system and the associated isotopic effects

Abstract— A wide range of stellar nucleosynthetic sources has been analyzed to derive their contributions of short‐lived and stable nuclei to the presolar cloud. This detailed study is required to infer the most plausible source(s) of short‐lived nuclei through a critical comparison among the various stellar sources that include AGB stars, novae, supernovae II, Ia, and Wolf‐Rayet stars that evolved to supernovae Ib/c. In order to produce the canonical value of ²⁶Al/²⁷Al in the early solar system, almost all stellar sources except low‐mass AGB stars imply large isotopic anomalies in Ca‐Al‐rich inclusions (CAIs). This is contrary to the observed isotopic compositions of CAIs. The discrepancy could impose stringent constraints on the formation and thermal evolution of CAIs from different chondrites. Among the various stellar scenarios, the injection of short‐lived nuclei into the previously formed solar protoplanetary disc by a massive star of an ad hoc chosen high‐injection mass cut is a possible scenario. There is a possibility of the contribution of short‐lived nuclides by a 1.5–3 M_⊙ AGB star as it implies the smallest shift in stable isotopes. A low‐mass AGB star of relatively low metallicity would be even a better source of short‐lived nuclei. However, this scenario would require independent gravitational collapse of the presolar cloud coupled with ambipolar diffusion of magnetic flux. Alternatively, numerous scenarios can be postulated that involve distant (≥10 pc) massive stars can contribute ⁶⁰Fe to the presolar cloud and can trigger its gravitational collapse. These scenarios would require production of ²⁶Al and ⁴¹Ca by irradiation in the early solar system. Significant production of ²⁶Al and ⁶⁰Fe can be explained if massive, rotating Wolf‐Rayet stars that evolved to supernovae Ib/c were involved.     

Numerical simulations of the differentiation of accreting planetesimals with 26Al and 60Fe as the heat sources

Abstract— Numerical simulations have been performed for the differentiation of planetesimals undergoing linear accretion growth with ²⁶Al and ⁶⁰Fe as the heat sources. Planetesimal accretion was started at chosen times up to 3 Ma after Ca‐Al‐rich inclusions (CAIs) were formed, and was continued for periods of 0.001–1 Ma. The planetesimals were initially porous, unconsolidated bodies at 250 K, but became sintered at around 700 K, ending up as compact bodies whose final radii were 20, 50, 100, or 270 km. With further heating, the planetesimals underwent melting and igneous differentiation. Two approaches to core segregation were tried. In the first, labelled A, the core grew gradually before silicate began to melt, and in the second, labelled B, the core segregated once the silicate had become 40% molten. In A, when the silicate had become 20% molten, the basaltic melt fraction began migrating upward to the surface, carrying ²⁶Al with it. The ⁶⁰Fe partitioned between core and mantle. The results show that the rate and timing of core and crust formation depend mainly on the time after CAIs when planetesimal accretion started. They imply significant melting where accretion was complete before 2 Ma, and a little melting in the deep interiors of planetesimals that accreted as late as 3 Ma. The latest melting would have occurred at <10 Ma. The effect on core and crust formation of the planetesimal's final size, the duration of accretion, and the choice of (⁶⁰Fe/⁵⁶Fe)_initial were also found to be important, particularly where accretion was late. The results are consistent with the isotopic ages of differentiated meteorites, and they suggest that the accretion of chondritic parent bodies began more than 2 or 3 Ma after CAIs.     

The role of impact and radiogenic heating in the early thermal evolution of Mars

Journal of Earth System Science - The planetary differentiation models of Mars are proposed that take into account core–mantle and core–mantle–crust differentiation. The numerical...     

Thermal evolution of trans‐Neptunian objects,icy satellites,and minor icy planets in the early solar system

Numerical simulations are performed to understand the early thermal evolution and planetary scale differentiation of icy bodies with the radii in the range of 100–2500 km. These icy bodies include trans‐Neptunian objects, minor icy planets (e.g., Ceres, Pluto); the icy satellites of Jupiter, Saturn, Uranus, and Neptune; and probably the icy‐rocky cores of these planets. The decay energy of the radionuclides, ²⁶Al, ⁶⁰Fe, ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th, along with the impact‐induced heating during the accretion of icy bodies were taken into account to thermally evolve these planetary bodies. The simulations were performed for a wide range of initial ice and rock (dust) mass fractions of the icy bodies. Three distinct accretion scenarios were used. The sinking of the rock mass fraction in primitive water oceans produced by the substantial melting of ice could lead to planetary scale differentiation with the formation of a rocky core that is surrounded by a water ocean and an icy crust within the initial tens of millions of years of the solar system in case the planetary bodies accreted prior to the substantial decay of ²⁶Al. However, over the course of billions of years, the heat produced due to ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th could have raised the temperature of the interiors of the icy bodies to the melting point of iron and silicates, thereby leading to the formation of an iron core. Our simulations indicate the presence of an iron core even at the center of icy bodies with radii ≥500 km for different ice mass fractions.