Nickel isotopic evidence for late-stage accretion of Mercury-like differentiated planetary embryos

Earth’s habitability is closely tied to its late-stage accretion, during which impactors delivered the majority of life-essential volatiles. However, the nature of these final building blocks remains poorly constrained. Nickel (Ni) can be a useful tracer in characterizing this accretion as most Ni in the bulk silicate Earth (BSE) comes from the late-stage impactors. Here, we apply Ni stable isotope analysis to a large number of meteorites and terrestrial rocks, and find that the BSE has a lighter Ni isotopic composition compared to chondrites. Using first-principles calculations based on density functional theory, we show that core-mantle differentiation cannot produce the observed light Ni isotopic composition of the BSE. Rather, the sub-chondritic Ni isotopic signature was established during Earth’s late-stage accretion, probably through the Moon-forming giant impact. We propose that a highly reduced sulfide-rich, Mercury-like body, whose mantle is characterized by light Ni isotopic composition, collided with and merged into the proto-Earth during the Moon-forming giant impact, producing the sub-chondritic Ni isotopic signature of the BSE, while delivering sulfur and probably other volatiles to the Earth. Based on Nickel isotope analysis of meteorites and terrestrial rocks, the authors suggest that the Bulk Silicate Earth has a sub-chondritic Nickel isotope composition. This signature is thought to result from the impact and accretion of a Mercury-like impactor which originated from the innermost Solar System.

T he Earth experienced a protracted accretion history over several tens up to 100 million years, which proceeded by the collision of numerous planetesimals and planetary embryos 1,2 . A fundamental assumption was that the Earth's building blocks as a whole were compositionally similar to undifferentiated chondritic meteorites. Researchers have looked among different classes of chondrites for the closest representative of the accreting materials that formed Earth [3][4][5][6] . However, emerging evidence points to a mismatch in many crucial elemental and isotopic ratios between chondritic meteorites and the accessible Earth, arguing for the possible accretion of additional materials that are chemically and isotopically different from extant meteorite collections [7][8][9][10][11] . Constraining the nature of these building blocks of Earth is important, because they not only provide fundamental information on terrestrial planet formation, but also help understand how the Earth evolved into its current habitable status.
The late accretion stages, including the Moon-forming giant impact and the late veneer event, likely account for only <10% of Earth's total mass 12 , but they represent a critical step for Earth to build its life-essential volatile budgets [13][14][15][16][17][18][19][20] . Dynamical models of Earth's growth suggest that the late accretion stages were highly heterogeneous, consisting of a mixture of materials from two genetically distinct reservoirs in the Solar nebula 20,21 . One endmember may originate from the inner Solar system and contain a reduced, non-carbonaceous component that is probably 'missing' in known meteorites 11,[20][21][22][23][24] . The other may be oxidized, carbonaceous chondrite-like material from the outer Solar system 17,[19][20][21]25 . When these materials were added to Earth is still debated [19][20][21][22][23][24][25][26][27] . The carbonaceous chondrite-like materials are commonly thought to be the source of major volatiles in Earth 28,29 . Recent high pressure-temperature experiments on metal alloy-silicate partitioning of volatiles (e.g., carbon, sulfur, and nitrogen), however, suggest that Earth's volatile abundance patterns could have been largely established by impact of a sulfurrich, differentiated planetary body with minimal contributions from carbonaceous chondrite-like materials 15,16 . Due to the lack of proper meteorite proxies, the nature of late-stage impactors remains poorly known.
Nickel isotopic compositions of meteorites and terrestrial rocks may hold important clues. Nickel in the bulk silicate Earth (BSE) was mostly derived from late-stage impactors, as that from earlier stages was largely segregated into the core due to its moderately siderophile nature 23 . Models predict that~95% of Ni in the BSE was derived from the last~35% of mass that accreted to Earth 23,30 . Nickel is non-volatile and partitions compatibly into the mantle dominant phaseolivinefollowing accretion, so that the BSE can potentially capture the Ni isotopic signature of late-stage accreting materials.
Nickel isotopic variations in meteorites have been well documented as shown in Fig. 1. Mass-independent nucleosynthetic Ni isotope anomalies arise from the heterogeneous distribution of presolar matters in the Solar protoplanetary disk, and thus trace the provenance of Earth's building blocks. The nucleosynthetic anomalies are present in carbonaceous and ordinary chondrites, with enstatite chondrites largely within error of the BSE 23,31-35 , supporting the general idea that the late-stage accreting materials mainly originated from an enstatite-like source region in the inner Solar system 23 . Iron meteorites display similar anomalies, together with different groups of chondrites, forming a dichotomy between carbonaceous and non-carbonaceous meteorites as found in many other isotope systems (e.g., Mo, Cr, Ru, Ti) 36 . Mass-dependent isotopic variations stem from physico-chemical processes in the Solar nebula and on the planetary parent bodies. Nickel isotopic compositions of enstatite, ordinary, and most carbonaceous chondrites exhibit a common value, expressed as δ 60/58 Ni (the 60 Ni/ 58 Ni ratio in parts per thousand, relative to the SRM986 standard; δ 60/58 Ni = ( 60/58 Ni sample / 60/58 Ni SRM986 − 1) × 1000), with an average of +0.23 ± 0.11‰ (2 SD, n = 34) 32,37-40 . The small isotopic variation observed in the carbonaceous chondrites most likely reflects the heterogeneous distribution of an isotopically light sulfide component 41,42 , supported by the roughly negative correlation between δ 60/58 Ni and sulfur content (Fig. 1). It is not surprising that iron meteorites have δ 60/58 Ni values within the 'chondritic' range ( Fig. 1), because they represent fragments of the disrupted cores of planetary bodies, and dominate the Ni budget.
The Ni isotopic composition of present BSE is poorly constrained. An earlier report of a few ultramafic rocks yielded δ 60/58 Ni values indistinguishable from the chondritic average, and they concluded that the BSE has a chondritic Ni isotopic composition 39 . This conclusion is questioned in a recent study combining new and reported peridotite samples 40 , which suggested that the BSE has a δ 60/58 Ni~0.1‰ lower than the chondrite average, a difference that was attributed to Earth's core formation 40,43 . Central to this debate is the limited Ni isotope data for terrestrial silicate rocks and scant information on Ni isotope fractionation during igneous and core-mantle differentiations.
Here, we show that the BSE has a sub-chondritic Ni isotopic composition by applying Ni isotope analysis to meteorites and terrestrial rocks. Our first-principles calculations further suggest that the light Ni isotopic signature of the BSE is not a result of  61 Ni/ 58 Ni), are from literature [31][32][33]35,36 , and the δ 60 Ni of chondrites and iron meteorites are from this study and literature [37][38][39][40][41]70 . The carbonaceous chondrites (CC), enstatite chondrites (EC), and ordinary chondrites (OC) have average δ 60 Ni values of 0.23 ± 0.14‰ (2 SD), 0.22 ± 0.02‰ (2 SD) and 0.24 ± 0.02‰ (2 SD), respectively. A roughly negative correlation between δ 60 Ni and sulfur abundance is observed in carbonaceous chondrites (upper panel), which may be caused by the presence of various abundances of sulfides with δ 60 Ni values as low as −1‰ 41,42 . The gray areas represent the dichotomy between carbonaceous and non-carbonaceous meteorites 36 . The sulfur abundances are from ref. 71  core-mantle differentiation. Rather, the signature was established during Earth's late-stage accretion, via impact and accretion of a highly reduced, Mercury-like impactor that likely originated from the innermost Solar system.
First-principles calculations on Ni isotope fractionation during core formation. The sub-chondritic Ni isotopic composition of the BSE could have resulted from two possible processes: isotope fractionation associated with Earth's differentiation, or the accretion of non-chondritic materials.
The former hypothesis is examined using first-principles calculations on Ni isotope fractionation factors (10 3 lnα of 60 Ni/ 58 Ni) among Earth's major Ni-bearing phases: olivine, wadsleyite, ringwoodite, bridgmanite, and Fe-Ni alloy. Limited differences in 10 3 lnα are found between olivine and wadsleyite/ringwoodite in the mantle transition zone and bridgmanite in the lower mantle (e.g., 10 3 lnα < 0.05‰ at 1500 K and <0.03‰ at 2000 K; Fig. 3a), which excludes the possibility of a hidden reservoir enriched in heavy Ni isotopes in the mantle. This lends credence to the use of accessible mantle and mantle-derived samples as representative of the present BSE Ni isotopic signature. Nickel isotope fractionation between Fe-Ni alloys and silicates (e.g., bridgmanite) under core-formation conditions is also negligible (P = 25-130 GPa; Fig. 3b and Supplementary Note 2). Notably, incorporation of sulfur into the Fe-Ni alloy slightly reduces the force constant of Ni, leading to the enrichment of light Ni isotopes in Fe-Ni alloys relative to the silicates (Fig. 3b).
To directly assess the equilibrium Ni isotope fractionation between silicate and metallic melts during core-mantle differentiation, we performed first-principles molecular dynamic simulations on melt phases of Fe 92 Ni 5 S 3 and Mg 30 NiSi 32 O 96 based on the density functional theory. The 10 3 lnα between Fe 92 Ni 5 S 3 and Mg 30 NiSi 32 O 96 melts is −0.011‰ at~38 GPa and 3500 K ( Fig. 3b; Supplementary Note 2), further confirming the conclusion based on crystals that core-mantle differentiation does not significantly fractionate Ni isotopes.
Two experimental studies investigated equilibrium Ni isotope fractionation between metal and silicate mineral or melt at low pressures (≤1.3 GPa) and temperatures (≤1623 K) 45,46 . Both studies predict limited Ni isotope fractionation under coreformation temperatures (<0.01‰ in terms of δ 60/58 Ni at T > 3000 K), a result that is consistent with our first-principles calculations at higher pressures, implying negligible pressure effect on silicate-metal Ni isotope fractionation. If the bulk Earth (BE) has a chondritic Ni isotopic composition, mass balance calculations using high-pressure metal-silicate Ni elemental partition coefficients 47 and isotope fractionation factors obtained from our first-principles calculation demonstrate that core-mantle differentiation cannot explain the sub-chondritic Ni isotopic composition of the BSE (Fig. 3c).
Other possible events including evaporative loss, collision erosion, and core-mantle chemical diffusion can also be discounted as causes for the sub-chondritic Ni isotopic signature of BSE 40 . Evaporative loss of Ni is unlikely given the relatively refractory nature of Ni. In addition, kinetic isotope fractionation associated with evaporation would lead to a heavy BSE Ni isotopic composition, opposite to observations (Fig. 2). Collisional erosion during Earth's formation preferentially removed early formed basaltic crust 48 . The terrestrial oceanic basalts have an average δ 60/58 Ni (0.03 ± 0.16‰; 2 SD) slightly lower than the BSE value, and thus collisional erosion cannot explain the The nature of late-stage accreting materials and its implications. Therefore, the Earth's mantle must have accreted subchondritic materials during its growth. In the early stages, metal and silicate melts equilibrate completely in the magma ocean 49 , so that the proto-BSE likely has a low Ni concentration and a chondritic Ni isotopic composition. Because of the moderately siderophile nature of Ni 23 and a possible disequilibrium scenario for the late-stage accretion 49,50 , the BSE's sub-chondritic Ni isotopic signature was likely established in the late stages. The late veneer following the main growth stage added the last <0.5% of mass to Earth and contributed <5% of the Ni budget of the BSE 18 . Hence, it is unlikely to be the event that produced the subchondritic Ni isotopic composition of BSE. To account for the observed Ni isotopic value of the present BSE, the late-veneer material accreted to the BSE would have to have had extremely low δ 60/58 Ni of around −2.5‰ ( Supplementary Fig. 8), a value that has not been found in any natural rocks, and would be inconsistent with the Ni isotopic composition of an average carbonaceous chondrite-like material for the late veneer 17,19,24,25 . The last significant stage of Earth's accretion was the Moonforming giant impact, contributing >20% of Ni budget of the BSE 49 . The Ni isotopic composition of the BSE could have been strongly influenced by the Moon-forming impactor 18,23 . Assuming a chondritic Ni isotopic composition for the proto-BSE as discussed above, mass balance calculation suggests that materials accreted to the proto-BSE have δ 60/58 Ni values as low as −0.35‰ (Supplementary Fig. 8). Accordingly, the Moon-forming impactor is unlikely to have a composition represented by chondrites.
Instead, we hypothesize that the sub-chondritic Ni isotopic composition of the BSE resulted from the impact and accretion of the sulfide-rich mantle of a highly reduced, differentiated planetary body. It has long been recognized that accretion of planetary embryos that were already differentiated into cores and mantles contributed significantly to the growth of Earth 51-53 . The Moon-forming impactor has been suggested to be a sulfur-rich, differentiated planetary body 15,16,27 ; but uncertainties remain as to whether it is a highly reduced, Mercury-like impactor 23,24,27 , or a relatively more oxidized body 20,21,26 . Mercury is the most reduced planet in the inner Solar system and has an abnormally high abundance of sulfides in its mantle [54][55][56] , whereas oxidized planetary embryos have sulfur segregated into their cores 57 . This is because sulfur is highly siderophile at high oxygen fugacity (fO 2 ) and partitions into the metallic core, but becomes lithophile and enters into the silicate melt as sulfide species under low fO 2 (e.g., five units below the iron-wüstite buffer; IW-5) 58,59 . Magmatic sulfides are the only major Ni-bearing phases that are isotopically much lighter than silicates (δ 60/58 Ni sulfide down to −1‰) 41 . Rocks with high sulfide/silicate ratios have light Ni isotopic compositions, which is most evident in magmatic Nisulfide deposits where the bulk δ 60/58 Ni values are negatively correlated with the sulfur content 41,42 . Therefore, when small planetary embryos (the proto-impactor) were formed in a sulfurrich early Solar nebula and differentiated into core and mantle under highly reduced environment similar to the Mercury (mean IW-5.4) 58 , the mantle would be sulfur-rich and have a light Ni isotopic composition (Fig. 4). By contrast, in the large proto-Earth, core-mantle differentiation proceeded under much higher pressure and likely more oxidizing condition (>IW-3; Fig. 4) 60 , in which sulfur behaves as a siderophile element 61,62 , leading to a sulfur-poor mantle. In this case, limited silicate-metal Ni isotope fractionation is expected (see discussion above), and thus the proto-Earth mantle likely has a chondritic Ni isotopic  composition. During the Moon-forming impact, the impactor's core merged directly into the proto-Earth's core due to its limited emulsification, while the remaining parts of the impactor were incorporated into the Earth's mantle 49,50 . The Mercury-like impactor's sulfur-rich mantle would have been completely dissolved in the planet-wide, more oxidizing terrestrial magma ocean 63 , and produced the sub-chondritic Ni isotopic signature of Earth's mantle (Fig. 4).
Our findings imply that, instead of an outer Solar system origin 20,21 , the Moon-forming impactor may represent a 'missing' Earth's building component originated from a highly reduced reservoir in the inner Solar system. While the absence of nucleosynthetic Ni isotopic anomalies in enstatite chondrites is consistent with an inner Solar system provenance for the Moonforming impactor 23 , the sub-chondritic mass-dependent Ni isotopic composition of the BSE further points towards a sulfide-rich, Mercury-like impactor, likely from closer to the Sun. This is in accordance with variations in nucleosynthetic isotope anomalies of other elements in meteorites, suggesting that the Earth possesses the most s-process enriched materials from the inner Solar system 11,22,24,64 . A most recent study found sprocess enriched ruthenium (Ru) isotopic signatures in Eoarchean rocks, supporting the idea that the pre-late veneer Earth incorporated building materials from the innermost region of the Solar system, most likely through the Moon-forming giant impact 24 . The impactor might be sulfur-rich, such that the highly siderophile Ru was partially retained in the mantle without being completely extracted to the core. Later addition of a carbonaceous chondrite-like late veneer with s-process Ru deficits from the outer Solar system ultimately built up the modern mantle Ru isotopic composition 24 . The proposed late-stage accretion of the highly reduced, Mercury-like planetary body may not only explain the broad geochemical similarity between Earth and Moon 27 , but also account for the volatile abundance patterns in the BSE 16 . Our study highlights the importance of inner planets, e.g., Mercury and Venus, in searching for the Earth's 'missing' building blocks that are not present in extant meteorite collections. Future studies on achondrites from the inner Solar system and samples from Venus and Mercury as well as experimental work will shed more light on these issues.

Methods
Nickel isotope analyses. We undertook an inter-laboratory comparison of geological reference materials using different analytical protocols in two labs: Indiana University (IU) and China University of Geosciences, Beijing (CUGB). Despite different double spike solutions and different column chemistry methods, the results of standards from two labs agree with each other within analytical uncertainty (Supplementary Table 1). A comparison of our data with published values is shown in Supplementary Fig. 1.
Ni isotope analyses at IU:. Sample powders were digested in a mixture of distilled HF + HNO 3 + HCl. After complete dissolution, aliquots of sample solutions containing 1.5 μg Ni were spiked with a 61 Ni -62 Ni double spike to reach an optimal spike-sample ratio of 64:36. The mixtures were refluxed on a hotplate to ensure sample-spike equilibration before column chemistry. Separation of Ni from the matrices was achieved using a three-stage, cation exchange chromatography procedure using Bio-Rad 200-400 mesh AG 50W-X8 resin. Briefly, the first column applies mixture of 20% 10 M HCl and 80% acetone to separate Ni from Fe, Mn, and Cr. The second column uses 15% 10 M HCl and 85% acetic acid to separate Ni from elements such as Mg, Al, Ca, and Ti, and the last column further purifies Ni using 0.9 M HNO 3 to remove Na and K. The Ni isotopic ratios were measured using Nu Plasma II MC-ICPMS at Indiana University.
Ni isotope analyses at CUGB:. Sample powders were dissolved and spiked follow the same protocol at IU, but with a different 61 Ni- 62 Ni spike solution to obtain optimal ratio of 62 Ni spike / 58 Ni sample = 1.3. The spike-sample solution was then passed through four-stage column chemistry. The step 1 column uses AG 50W-X8 and AG 1-X8 resins to separate Ni from Fe and Ca.
Step 2 column uses AG 50W-X8 resin to separate Ni from Mg, Ti, and Al in a media of 0.15 M HNO 3  First-principles calculations. Ab initio calculations were performed using the software "Quantum Espresso" 65 , which is based on the density functional theory (DFT), plane wave, and pseudopotentials. The generalized gradient approximation (GGA) was adopted to describe the exchange-correlation functional. The pseudopotential for magnesium was generated using von Barth and Car's methods with a cutoff radius 2.5 Bohr. The electron configurations are 3s 2 3p 0 , 3s 1 3p 1 , 3s 1 3p 0.5 3d 0.5 , 3s 1 3p 0.5 , and 3s 1 3d 1 with decreasing weights of 1.5, 0.6, 0.3, 0.3, and 0.2, respectively. The pseudopotentials for nickel, silicon, and oxygen were generated by the method in Troullier and Martins 66 . The cutoff radius are 1.45 Bohr with the electron configuration of 2s 2 2p 4 for oxygen and 1.47 Bohr with the electron configuration of 3s 2 3p 4 3d 0 for silicon. The cutoff radius for nickel is 2.1 Bohr with the electron configuration of 4s 2 3d 8 4p 0 . The pseudopotential for Fe was generated using the Vanderbilt method 67 with a valence configuration of 3s 2 3p 6 3d 6.5 4s 1 4p 0 and a cutoff radius of 2.0 Bohr for Fe.
We first optimized all crystal structures of Ni-bearing minerals using the variable cell shape molecular dynamics method 68 with different k-point grids dependent on the sizes of unit cells (Supplementary Table 7). The energy cutoff for plane wave and charge density are set to 70 Ry and 700 Ry, respectively. The residual forces converge within 10 −4 Ry/Bohr. After the relaxed structures were obtained, we then calculated vibrational frequencies using the finite displacement method as implemented in the open-source code PHONOPY 69 . Consequently, the reduced partition function ratios β of 60 Ni/ 58 Ni for all phases can be calculated from the equation: 1Àe Àu il e À 1 2 u il . In order to directly estimate the equilibrium Ni isotope fractionation between silicate and metallic melts, we conducted first-principles molecular dynamics (FPMD) simulations on Mg 30 Ni 2 Si 32 O 96 and Fe 92 Ni 5 S 3 melts based on the DFT using the Vienna ab initio simulation package (VASP). The GGA was adopted for the exchange-correlation functional and the projector-augmented-wave (PAW) pseudopotentials were used. The energy cutoff for the plane wave was 600 eV. The Brillouin zone summations over the electronic states were performed at gamma Fig. 4 Cartoon showing the merger of a highly reduced, Mercury-like planetary body with the relatively more oxidizing proto-Earth, and schematic evolution of Ni isotopic composition of the BSE. a The small, proto-impactor from the highly reduced innermost region of the Solar system differentiated into a core, a sulfur-rich mantle, and likely a sulfide layer at the mantle-core boundary, because sulfur behaves more lithophile at low fO 2 and pressure 58,59,61,62 . The impactor's mantle likely had a subchondritic Ni isotopic composition due to the enrichment of sulfides. Differentiation on the large, proto-Earth partitioned sulfur dominantly into the core because sulfur is more siderophile at relatively high fO 2 and pressure 58,59,61,62 . Therefore, the mantle of the proto-Earth had a Ni isotopic composition close to the chondritic value. b The Moon-forming giant impact would have completely melted the Earth, forming a planetwide, more oxidizing magma ocean with high solubility of sulfur 60 . The sulfides in the impactor's mantle were dissolved in the terrestrial magma ocean, producing a sub-chondritic Ni isotopic composition for the Earth's mantle. c The late veneer following the main growth stage added carbonaceous chondrite-like materials to the Earth with limited effects on the Ni isotopic systematics of the BSE as discussed in the text.
point. The FPMD simulations were performed in the NVT thermodynamic ensemble with a fixed temperature of 3000 K and the Nosé thermostat was used. The time step was set to be 1 fs, and the total running time is up to 60 ps. The initial liquid configurations were prepared by conducting simulations on the structures at 6000 K. The cell parameter of the cubic box is 11.15 Å for Mg 30 Supplementary Fig. 5), respectively. After equilibration, we extracted 66 snapshots from the FPMD trajectory every 200 steps and only optimized the atomic positions of Ni with fixed cubic boxes. This strategy makes the Ni atoms in each snapshot are at the local equilibrium positions. Then we estimated the force constant matrix of Ni atoms in all snapshots using the small displacement method based on the harmonic approximation. The force constants <F > of the Ni atom in Mg 30 Supplementary Fig. 6), respectively. Therefore, the 10 3 lnα between Mg 30 Ni 2 Si 32 O 96 and Fe 92 Ni 5 S 3 melts is 0.016 ± 0.003 ‰ at 3000 K and 0.011 ± 0.002 ‰ at 3500 K.

Data availability
All data in this study are included in the supplementary information files and are available from the corresponding author.