通过第一性原理计算研究内地核的结构与成分

The determination of the composition and structure of the Earth’s inner core has long been the major subject in the study of the Earth’s deep interior. It’s widely believed that the Earth’s core is formed by iron with a fraction of nickel. However, light elements must exist in the inner core because the earth core is less dense than pure iron-nickel alloy (~2-3% in the solid inner core and ~6-7% in the liquid outer core). The questions are what and how much light element is there in the iron-nickel alloy. Besides the composition, the crystal structure of the iron with or without light element is also not well known. According to the seismological observations, the sound waves propagate 3-4% faster along the spin axis than in the equatorial plane. That means the inner core is anisotropic. The densest structure of iron-nickel alloy should be h.c.p structure under the very high pressures. However, the h,c,p structure does not propagate waves anisotropic ally. Then what is the structure of the iron-nickel alloy or the iron-nickle-light element alloy. In this study, we tried to predict the composition and the structure of the inner core through ab initio calculation of the Gibbs free energy, which is a function of internal energy, density and entropy. We conclude that the h.c.p structure is more stable than the b.c.c structure under high pressure and 0 K, but with the increase of temperature, the free energy of the b.c.c structure is decreasing much faster than the h.c.p structure caused by the vibration of the atomics, so the b.c.c structure is more stable at high temperatures. With the addition of light elements (S or Si or both), the free energy of b.c.c. decreases even faster, about 3at% of Si not only explains why the inner core is about 2-3 % lighter than the iron-nickle alloy, but also reasons why the inner core is anisotropic, since the b.c.c. structure becomes more stable than the h.c.p structure at 5500-6000K and b.c.c. is anisotropic in propagating seismic waves. Therefore, we infer that the inner core of the earth is formed by b.c.c iron and a fraction of nickel plus ~3at.% Si, with a temperature higher than 5500K, which is consistent with the studies from other approaches.