Iron silicate perovskite and postperovskite in the deep lower mantle.

Ferromagnesian silicates are the dominant constituents of the Earth's mantle, which comprise more than 80% of our planet by volume. To interpret the low shear-velocity anomalies in the lower mantle, we need to construct a reliable transformation diagram of ferromagnesian silicates over a wide pressure-temperature (P-T) range. While MgSiO3 in the perovskite structure has been extensively studied due to its dominance on Earth, phase transformations of iron silicates under the lower mantle conditions remain unresolved. In this study, we have obtained an iron silicate phase in the perovskite (Pv) structure using synthetic fayalite (Fe2SiO4) as the starting material under P-T conditions of the lower mantle. Chemical analyses revealed an unexpectedly high Fe/Si ratio of 1.72(3) for the Pv phase in coexistence with metallic iron particles, indicating incorporation of about 25 mol% Fe2O3 in the Pv phase with an approximate chemical formula (Fe2+0.75Fe3+0.25)(Fe3+0.25Si0.75)O3. We further obtained an iron silicate phase in the postperovskite (PPv) structure above 95 GPa. The calculated curves of compressional (VP) and shear velocity (VS) of iron silicate Pv and PPv as a function of pressure are nearly parallel to those of MgSiO3, respectively. To the best of our knowledge, the iron silicate Pv and PPv are the densest phases among all the reported silicates stable at P-T conditions of the lower mantle. The high ferric iron content in the silicate phase and the spin-crossover of ferric iron at the Si-site above ~55 GPa should be taken into account in order to interpret the seismic observations. Our results would provide crucial information for constraining the geophysical and geochemical models of the lower mantle.

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity.

The compositional and thermal state of Earth's mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earth's early history.

Hydrogen-bearing iron peroxide and the origin of ultralow-velocity zones.

Ultralow-velocity zones (ULVZs) at Earth's core-mantle boundary region have important implications for the chemical composition and thermal structure of our planet, but their origin has long been debated. Hydrogen-bearing iron peroxide (FeO2Hx) in the pyrite-type crystal structure was recently found to be stable under the conditions of the lowermost mantle. Using high-pressure experiments and theoretical calculations, we find that iron peroxide with a varying amount of hydrogen has a high density and high Poisson ratio as well as extremely low sound velocities consistent with ULVZs. Here we also report a reaction between iron and water at 86 gigapascals and 2,200 kelvin that produces FeO2Hx. This would provide a mechanism for generating the observed volume occupied by ULVZs through the reaction of about one-tenth the mass of Earth's ocean water in subducted hydrous minerals with the effectively unlimited reservoir of iron in Earth's core. Unlike other candidates for the composition of ULVZs, FeO2Hx synthesized from the superoxidation of iron by water would not require an extra transportation mechanism to migrate to the core-mantle boundary. These dense FeO2Hx-rich domains would be expected to form directly in the core-mantle boundary region and their properties would provide an explanation for the many enigmatic seismic features that are observed in ULVZs.

Genkwanin suppresses mitochondrial dysfunction to alleviate IL-1ß-elicited inflammation, apoptosis, and degradation of extracellular matrix in chondrocytes through upregulating DUSP1.

Osteoarthritis (OA) is a form of chronic degenerative disease contributing to elevated disability rate among the elderly. Genkwanin is an active component extracted from Daphne genkwa possessing pharmacologic effects. Here, this study is designed to expound the specific role of genkwanin in OA and elaborate the probable downstream mechanism. First, the viability of chondrocytes in the presence or absence of interleukin-1 beta (IL-1ß) treatment was detected by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay. Terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling assay was used to assess cell apoptosis. Inflammatory response was estimated through enzyme-linked immunosorbent assay and Western blot. In addition, immunofluorescence staining and Western blot were utilized to measure the expression of extracellular matrix (ECM)-associated proteins. Dual-specificity protein phosphatase-1 (DUSP1) expression was tested by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and Western blot. Following DUSP1 elevation in genkwanin-treated chondrocytes exposed to IL-1ß, inflammatory response and ECM-associated factors were evaluated as forementioned. In addition, 5,5',6,6'-tetrachloro-1,1',3,3'-tetraethylbenzimidazolocarbocyanine iodide staining was to assess the mitochondrial membrane potential. Adenosine triphosphate (ATP) level was examined with ATP assay kit, and RT-qPCR was used to test mitochondrial DNA expression. Results indicated that genkwanin administration enhanced the viability while ameliorated the apoptosis, inflammatory response, and ECM degradation in IL-1ß-induced chondrocytes. Besides, genkwanin treatment fortified DUSP1 expression in IL-1ß-exposed chondrocytes. DUSP1 interference further offsets the impacts of genkwanin on the inflammation, ECM degradation, and mitochondrial dysfunction in IL-1ß-challenged chondrocytes. In short, genkwanin enhanced DUSP1 expression to mitigate mitochondrial dysfunction, thus ameliorating IL-1ß-elicited inflammation, apoptosis, and degradation of ECM in chondrocytes.

Spin crossover in ferropericlase and velocity heterogeneities in the lower mantle.

Deciphering the origin of seismic velocity heterogeneities in the mantle is crucial to understanding internal structures and processes at work in the Earth. The spin crossover in iron in ferropericlase (Fp), the second most abundant phase in the lower mantle, introduces unfamiliar effects on seismic velocities. First-principles calculations indicate that anticorrelation between shear velocity (VS) and bulk sound velocity (Vφ) in the mantle, usually interpreted as compositional heterogeneity, can also be produced in homogeneous aggregates containing Fp. The spin crossover also suppresses thermally induced heterogeneity in longitudinal velocity (VP) at certain depths but not in VS. This effect is observed in tomography models at conditions where the spin crossover in Fp is expected in the lower mantle. In addition, the one-of-a-kind signature of this spin crossover in the RS/P (∂ ln VS/∂ ln VP) heterogeneity ratio might be a useful fingerprint to detect the presence of Fp in the lower mantle.

Elastic anomalies in a spin-crossover system: ferropericlase at lower mantle conditions.

The discovery of a pressure induced iron-related spin crossover in Mg((1-x))Fe(x)O ferropericlase (Fp) and Mg-silicate perovskite, the major phases of Earth's lower mantle, has raised new questions about mantle properties which are of central importance to seismology. Despite extensive experimental work on the anomalous elasticity of Fp throughout the crossover, inconsistencies reported in the literature are still unexplained. Here we introduce a formulation for thermoelasticity of spin crossover systems, apply it to Fp by combining it with predictive first principles density-functional theory with on-site repulsion parameter U calculations, and contrast results with available data on samples with various iron concentrations. We explain why the shear modulus of Fp should not soften along the crossover, as observed in some experiments, predict its velocities at lower mantle conditions, and show the importance of constraining the elastic properties of minerals without extrapolations for analyses of the thermochemical state of this region.

CKLF1 interference alleviates IL­1ß­induced inflammation, apoptosis and degradation of the extracellular matrix in chondrocytes via CCR5.

Osteoarthritis (OA) is a type of joint disease with a rising prevalence and incidence among the elderly across the global population. Chemokine-like factor 1 (CKLF1) is a human cytokine, which has been demonstrated to be involved in the progression of multiple human diseases. However, little attention has been paid to the impact of CKLF1 on OA. The present study was designed to identify the role of CKLF1 in OA and to clarify the regulatory mechanism. The expression levels of CKLF1 and its receptor CC chemokine receptor 5 (CCR5) were examined by reverse transcription-quantitative PCR (RT-qPCR) and western blotting. A Cell Counting Kit-8 assay was used to estimate cell viability. The levels and expression of inflammatory factors were determined by ELISA and RT-qPCR, respectively. Apoptosis was investigated by TUNEL assays and the protein levels of apoptosis-related factors were analyzed by western blotting. RT-qPCR and western blotting were used to examine the expression of extracellular matrix (ECM) degradation-associated proteins and ECM components. Dimethylmethylene blue analysis was used to analyze the production of soluble glycosamine sulfate additive. A co-immunoprecipitation assay was used to confirm the protein interaction between CKLF1 and CCR5. The results revealed that CKLF1 expression was increased in IL-1ß-exposed murine chondrogenic ATDC5 cells. Furthermore, CKLF1 silencing enhanced the viability of IL-1ß-induced ATDC5 cells, while inflammation, apoptosis and degradation of the ECM were reduced. Additionally, CKLF1 knockdown led to decreased CCR5 expression in IL-1ß-challenged ATDC5 cells, and CKLF1 bound with CCR5. The enhanced viability, as well as the suppressed inflammation, apoptosis and degradation of the ECM, following CKLF1 knockdown in the IL-1ß-induced ATDC5 cells were all restored after CCR5 was overexpressed. In conclusion, CKLF1 might serve a detrimental role in the development of OA by targeting its receptor CCR5.

Small-scale layered structures at the inner core boundary.

The fine-scale seismic features near the inner core boundary (ICB) provide critical insights into the thermal, chemical, and geodynamical interactions between liquid and solid cores, and may shed light on the evolution mechanism of the Earth's core. Here, we utilize a dataset of pre-critical PKiKP waveforms to constrain the fine structure at the ICB, considering the influence of various factors such as source complexity, structural anomalies in the mantle, and properties at the ICB. Our modeling suggests a sharp ICB beneath Mongolia and most of Northeast Asia, but a locally laminated ICB structure beneath Central Asia, Siberia, and part of Northeast Asia. The complex ICB structure might be explained by either the existence of a kilometer-scale thickness of mushy zone, or the localized coexistence of bcc and hcp iron phase at the ICB. We infer that there may be considerable lateral variations in the dendrites growing process at ICB, probably due to the complicated thermochemical and geodynamical interaction between the outer and inner core.

Inelastic Neutron Scattering Study of Phonon Density of States of Iodine Oxides and First-Principles Calculations.

Iodine oxides I2Oy (y = 4, 5, 6) crystallize into atypical structures that fall between molecular- and framework-base types and exhibit high reactivity in an ambient environment, a property highly desired in the so-called "agent defeat materials". Inelastic neutron scattering experiments were performed to determine the phonon density of states of the newly synthesized I2O5 and I2O6 samples. First-principles calculations were carried out for I2O4, I2O5, and I2O6 to predict their thermodynamic properties and phonon density of states. Comparison of the INS data with the Raman and infrared measurements as well as the first-principles calculations sheds light on their distinctive, anisotropic thermomechanical properties.

[Treatment of severe medial tibial bone defect in primary total knee arthroplasty with autogenous bone graft and plate fixation].

OBJECTIVE: To explore the technique of autogenous bone graft combined with plate fixation in total knee arthroplasty(TKA) with severe proximal medial tibial bone defect. METHODS: From March 2012 to October 2018, 21 patients (9 males and 12 females) with severe bone defects in the proximal medial tibia during primary total knee arthroplasty were treated with autogenous structural bone grafting and steel plate fixation, with an age of 61 to 77 years old with an average of (69.6±9.1) years and a course of 64 to 257 months with an average of (73.6±170.7) months. According to Rand classification, there were 13 cases of type Ⅲb and 8 cases of type Ⅳb. Postoperative complications were observed, and knee joint function was evaluated by the Hospital for Special Surgery (HSS) score and SF-36 quality of life score. RESULTS: All 21 patients were followed up for 37 to 64 months with an average of (49.5±13.7) months. The incisions of all patients healed smoothly, and 2 patients developed lower limb intermuscular venous plexus thrombosis after operation. There were no periprosthetic infection, loosening of prosthesis and other complications. The autogenous bone grafts of all patients achieved bony healing during postoperative X-ray follow-up, and the healing time was 8 to 13 months with an average of (10.1±2.3) months. The HSS score of patients increased significantly from 30 to 48 with an average of (53.4±4.2) before operation to 75 to 92 with an average of (81.2±8.4) at the final follow-up (P<0.05). The SF-36 quality of life score of patients after operation was significantly different from that before operation (P<0.05). CONCLUSION: The technique of autogenous bone graft combined with steel plate fixation can achieve satisfactory osseointegration effect in the treatment of severe proximal tibial bone defects in primary knee arthroplasty, with less complications and obvious improvement in knee function.

Elastic properties of Fe-bearing Akimotoite at mantle conditions: Implications for composition and temperature in lower mantle transition zone.

The pyrolite model, which can reproduce the upper-mantle seismic velocity and density profiles, was suggested to have significantly lower velocities and density than seismic models in the lower mantle transition zone (MTZ). This argument has been taken as mineral-physics evidence for a compositionally distinct lower MTZ. However, previous studies only estimated the pyrolite velocities and density along a one-dimension (1D) geotherm and never considered the effect of lateral temperature heterogeneity. Because the majorite-perovskite-akimotoite triple point is close to the normal mantle geotherm in the lower MTZ, the lateral low-temperature anomaly can result in the presence of a significant fraction of akimotoite in pyrolitic lower MTZ. In this study, we reported the elastic properties of Fe-bearing akimotoite based on first-principles calculations. Combining with literature data, we found that the seismic velocities and density of the pyrolite model can match well those in the lower MTZ when the lateral temperature heterogeneity is modeled by a Gaussian distribution with a standard deviation of ∼100 K and an average temperature of dozens of K higher than the triple point of MgSiO3. We suggest that a harzburgite-rich lower MTZ is not required and the whole mantle convection is expected to be more favorable globally.

[Case-control study on carpal canal endoscopy and arthroscopy for the treatment of plantar fasciopathy].

OBJECTIVE: To explore clinical effects of carpal canal endoscopy in treating patients with plantar fasciopathy who failed by conservative treatment. METHODS: From August 2018 to August 2019, 50 patients with plantar fascia were divided into two groups and 25 patients in each group. In carpal canal endoscopy group, included 11 males and 14 females, aged from 39 to 67 years old with an average of(57.7±6.4) years old;carpal canal endoscopy was used to plantar fascia release. In arthroscopy group, included 9 males and 16 females, aged from 41 to 73 years old with an average of (58.1±7.2) years old;conventional 4.0 mm arthroscopy Instruments was used to plantar fascia release. Operation time, hospitalization expense and postoperative complications between two groups were observed and compared. Postoperative visual analogue scale(VAS) and American Orthopedic Foot Ankle Society (AOFAS) score were used to evaluate clinical function. RESULTS: All patients were followed up from 12 to 18 months with an average of (14.3±2.1) months. There were significant differentces in operation time and hospitalization expense between two groups (P<0.05). Surgical incision healed well in carpal canal endoscopy group, and 2 patients delayed union in arthroscopy group, and no difference between two groups (P>0.05). There were no statistical differences in VAS, AOFAS and grading between two groups at 12 months after operation(P>0.05). CONCLUSION: The outcome of carpal canal endoscopy and arthroscopy has similar effects in treating plantar fascia. While carpal canal endoscopy has advantages of need not perfusion during opertaion, protect soft tissue well, less operation time, and lower cost.

First-principles calculations of the structural and dynamic properties, and the equation of state of crystalline iodine oxides I2O4, I2O5, and I2O6.

The structural and dynamical correlations, and the equation of state of crystalline I(2)O(4), I(2)O(5), and I(2)O(6) are investigated by first-principles calculations based on the density functional theory (DFT). The lattice dynamics results reveal distinctive features in the phonon density of states among the three crystals. The frequencies of the stretch modes in I(2)O(4) and I(2)O(5) are clearly separated from those of the other (e.g., bending) modes by a gap, with all stretch modes above the gap. In contrast, the gap in I(2)O(6) separates the highest-frequency stretch modes with other stretch modes, and there is no gap between the stretch and the other modes in I(2)O(6). The motion of iodine atoms is involved in all vibrational modes in I(2)O(5), but only in low-frequency lattice modes in I(2)O(6). In I(2)O(4), iodine atoms are involved in modes with frequency below 700 cm(-1). Van der Waals correction within our DFT calculations is found to reduce the overestimation of the equilibrium volume, with its effect on structure similar to the pressure effect. Namely, both effects significantly decrease the inter-molecular distances, while slightly increasing the bond lengths within the molecules. This causes the frequencies of some vibrational modes to decrease with pressure, resulting in negative "modes Grüneisen parameters" for those modes. Thermodynamic properties, derived from the equation of state, of crystalline I(2)O(4), I(2)O(5), and I(2)O(6) are discussed within the quasi-harmonic approximation.

Vibrational and thermodynamic properties of ß-HMX: a first-principles investigation.

Thermodynamic properties of ß-HMX crystal are investigated using the quasi-harmonic approximation and density functional theory within the local density approximation (LDA), generalized gradient approximation (GGA), and GGA + empirical van der Waals (vdW) correction. It is found that GGA well describes the thermal expansion coefficient and heat capacity but fails to produce correct bulk modulus and equilibrium volume. The vdW correction improves the bulk modulus and volume, but worsens the thermal expansion coefficient and heat capacity. In contrast, LDA describes all thermodynamic properties with reasonable accuracy, and overall is a good exchange-correlation functional for ß-HMX molecular crystal. The results also demonstrate significant contributions of phonons to the equation of state. The static calculation of equilibrium volume for ß-HMX differs from the room-temperature value incorporating lattice vibrations by over 5%. Therefore, for molecular crystals, it is essential to include phonon contributions when calculated equation of state is compared with experimental data at ambient condition.

Coupled deep-mantle carbon-water cycle: Evidence from lower-mantle diamonds.

Diamonds form in a variety of environments between subducted crust, lithospheric and deep mantle. Recently, deep source diamonds with inclusions of the high-pressure H2O-phase ice-VII were discovered. By correlating the pressures of ice-VII inclusions with those of other high-pressure inclusions, we assess quantitatively the pressures and temperatures of their entrapment. We show that the ice-VII-bearing diamonds formed at depths down to 800 ± 60 km but at temperatures 200-500 K below average mantle temperature that match the pressure-temperature conditions of decomposing dense hydrous mantle silicates. Our work presents strong evidence for coupled recycling of water and carbon in the deep mantle based on natural samples.

Biomimetic Polydopamine-Modified Silk Fibroin/Curcumin Nanofibrous Scaffolds for Chemo-photothermal Therapy of Bone Tumor.

The simultaneous therapy of tumor recurrence and bone defects resulting from surgical resection of osteosarcoma is still a challenge in the clinic. Combination therapy based on a localized drug-delivery system shows great promise in the treatment of osteosarcoma. Herein, bifunctional polydopamine (PDA)-modified curcumin (CM)-loaded silk fibroin (SF) composite (SF/CM-PDA) nanofibrous scaffolds, which combined photothermal therapy with chemotherapy to synergistically enhance osteosarcoma therapy, were prepared by PDA coating of the SF/CM nanofibrous scaffolds fabricated by supercritical carbon dioxide (SC-CO2) technology. The PDA coating improved hydrophilicity and mechanical strength of the SF/CM scaffolds. The SF/CM-PDA scaffolds present good photothermal conversion capacity and excellent photostability. The low pH and near-infrared (NIR) irradiation could effectively accelerate release of CM in the SF/CM-PDA scaffolds. The in vitro anticancer results indicated that the biocompatible SF/CM-PDA scaffolds had a long-term, stable, and superior anticancer effect compared to pure CM. Furthermore, the SF/CM-PDA scaffolds significantly increased the growth inhibition of osteosarcoma MG-63 cells under NIR irradiation (808 nm and 1.3 W/cm2). Besides, the SF/CM-PDA scaffolds could enhance osteoblast MC3T3-E1 cell proliferation in vitro when the mass ratio of CM was 0.05-0.5%. This work has therefore demonstrated that the bifunctional SF/CM-PDA scaffolds provide a competitive strategy for local osteosarcoma therapy and bone regeneration.

Formation of large low shear velocity provinces through the decomposition of oxidized mantle.

Large Low Shear Velocity Provinces (LLSVPs) in the lowermost mantle are key to understanding the chemical composition and thermal structure of the deep Earth, but their origins have long been debated. Bridgmanite, the most abundant lower-mantle mineral, can incorporate extensive amounts of iron (Fe) with effects on various geophysical properties. Here our high-pressure experiments and ab initio calculations reveal that a ferric-iron-rich bridgmanite coexists with an Fe-poor bridgmanite in the 90 mol% MgSiO3-10 mol% Fe2O3 system, rather than forming a homogeneous single phase. The Fe3+-rich bridgmanite has substantially lower velocities and a higher VP/VS ratio than MgSiO3 bridgmanite under lowermost-mantle conditions. Our modeling shows that the enrichment of Fe3+-rich bridgmanite in a pyrolitic composition can explain the observed features of the LLSVPs. The presence of Fe3+-rich materials within LLSVPs may have profound effects on the deep reservoirs of redox-sensitive elements and their isotopes.

High-temperature inter-mineral potassium isotope fractionation: implications for K-Ca-Ar chronology.

Recent advances in high-precision potassium (K) isotopic analysis have found considerable isotopic variation in rock samples of the Earth's continental and oceanic crusts; however, it is still uncertain whether there is any resolvable inter-mineral and mineral-melt K isotopic fractionation during igneous and metamorphic processes. Here, we report K isotope compositions of mineral separates from three extremely well preserved igneous rocks (intrusive/extrusive and mafic/intermediate/felsic) in order to investigate possible inter-mineral and mineral-melt K isotopic fractionation at magmatic temperatures. For the first time, we found large inter-mineral fractionation of K isotopes in natural samples (up to 1.072‰), where plagioclase displays a significant enrichment of heavier K isotopes compared to potassium feldspar and biotite in a granite. In addition, we also observed smaller but measurable K isotope fractionation (0.280‰±0.030‰) between ternary feldspar phenocrysts and matrix in a trachyandesite, as well as a comparable isotope fractionation (0.331‰±0.010‰) between plagioclase and the bulk in a gabbroic intrusive rock. We also evaluated such results by comparing the theoretically calculated equilibrium K isotope fractionation factors between relevant igneous minerals in literature and this study. In general, the measured inter-mineral fractionations are consistent with the theoretical calculations (i.e., plagioclase is enriched in heavier isotopes compared to potassium feldspar). Specifically, the measured K isotope fractionation between phenocryst rim and matrix in the trachyandesite agrees well with the calculated equilibrium isotope fractionation. However, the measured K isotope fractionations between phenocryst core and matrix as well as between plagioclase and K-feldspar are significantly larger (by a factor of ~2-3) than the calculated isotope fractionations, which suggest isotopic disequilibrium due to kinetic processes. Using a range of plagioclase-melt isotope fractionation factors inferred from the theoretical calculations in this study, we modeled the K isotopic fractionation during the formation of lunar anorthositic crust, and the result shows a negligible effect on the K isotopic compositions in both lunar crust and mantle. The K isotopic difference between Earth and Moon, therefore, cannot be the result of Lunar Magma Ocean differentiation. Finally, we evaluate the effect of observed inter-mineral fractionations on K-Ar and 40Ar-39Ar dating. This study indicates the variation of 40K/K ratio would contribute a maximum 0.08% error to the K-Ar and 40Ar-39Ar age uncertainties. We propose a refined 40K/total K ratio as 0.00011664±0.00000011 (116.64±0.11ppm) instead of the conventional value, 0.0001167(2) for the present Earth. Because some minerals fractionate K isotopes, ultrahigh precision age dating with the K-Ca-Ar dating systems must measure the K isotope fractionation in the same mineral fractions used for age dating.

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.

Density functional study of 1,3,5-trinitro-1,3,5-triazine molecular crystal with van der Waals interactions.

Volume dependence of the total energy and vibrational properties of crystalline l,3,5-trinitro-l,3,5-triazine (RDX) are calculated using the density functional theory (DFT). For this molecular crystal, properties calculated with a generalized gradient approximation to the exchange-correlation energy differ drastically from experimental values. This discrepancy arises from the inadequacy in treating weak van der Waals (vdW) interactions between molecules in the crystal, and an empirical vdW correction to DFT (DFT-D approach by Grimme) is shown to account for the dispersion effects accurately for the RDX crystal, while incurring little computational overhead. The nonempirical van der Waals density-functional (vdW-DF) method also provides an accurate description of the vdW corrections but with orders-of-magnitude more computation. We find that the vibrational properties of RDX are affected in a nontrivial manner by the vdW correction due to its dual role--reduction of the equilibrium volume and additional atomic forces.