Evidence for an oxygen-depleted liquid outer core of the Earth.

On the basis of geophysical observations, cosmochemical constraints, and high-pressure experimental data, the Earth's liquid outer core consists of mainly liquid iron alloyed with about ten per cent (by weight) of light elements. Although the concentrations of the light elements are small, they nevertheless affect the Earth's core: its rate of cooling, the growth of the inner core, the dynamics of core convection, and the evolution of the geodynamo. Several light elements-including sulphur, oxygen, silicon, carbon and hydrogen-have been suggested, but the precise identity of the light elements in the Earth's core is still unclear. Oxygen has been proposed as a major light element in the core on the basis of cosmochemical arguments and chemical reactions during accretion. Its presence in the core has direct implications for Earth accretion conditions of oxidation state, pressure and temperature. Here we report new shockwave data in the Fe-S-O system that are directly applicable to the outer core. The data include both density and sound velocity measurements, which we compare with the observed density and velocity profiles of the liquid outer core. The results show that we can rule out oxygen as a major light element in the liquid outer core because adding oxygen into liquid iron would not reproduce simultaneously the observed density and sound velocity profiles of the outer core. An oxygen-depleted core would imply a more reduced environment during early Earth accretion.

Density functional theory investigation of the phonon instability, thermal equation of state and melting curve of Mo.

The phonon instability and thermal equation of state of Mo are extensively investigated using density functional theory. The calculated phonon dispersion curves agree well with experiments. Under compression, we captured a large softening in the transverse acoustic (TA) branches of body-centred cubic Mo. When the pressure is raised to 716 GPa, the frequencies along Γ-N in the TA branches soften to imaginary frequencies, indicating structural instability. For face-centred cubic Mo, the phonon calculations predicted the stability by promoting the frequencies from imaginary to real. Within quasi-harmonic approximation, we predicted the thermal equation of state and some other properties including the thermal expansion coefficient α, product αK(T), heat capacity C(V), entropy S, Grüneisen parameter γ and Debye temperature Θ(D). The melting curves of Mo were also obtained successfully.

Lattice dynamics and thermodynamics of molybdenum from first-principles calculations.

We calculated the phase transition, elastic constants, full phonon dispersion curves, and thermal properties of molybdenum (Mo) for a wide range of pressures using density functional theory. Mo is stable in the body-centered-cubic (bcc) structure up to 703 +/- 19 GPa and then transforms to the face-centered close-packed (fcc) structure at zero temperature. Under high temperature and pressure, the fcc phase of Mo is more stable than the bcc phase. The calculated phonon dispersion curves accord excellently with experiments. Under pressure, we captured a large softening along H-P in the TA branches. When the volume is compressed to 7.69 A(3), the frequencies along H-P in the TA branches soften to imaginary frequencies, indicating a structural instability. When the pressure increases, the phonon calculations on the fcc Mo predict the stability by promoting the frequencies along Gamma to X and Gamma to L symmetry lines from imaginary to real. The thermal equation of state was also investigated. From the thermal expansion coefficient and the heat capacity, we found that the quasiharmonic approximation was valid only up to about melting point at zero pressure. However, under pressure, the validity can be extended to a much higher temperature.

A flyer-impact technique for measuring viscosity of metal under shock compression.

A flyer-impact technique, different from the explosive method of [Sakharov et al., Sov. Phys. Dokl. 9, 1091 (1965)], is developed to investigate the viscosity of shocked metals. The shock wave with a front of sinusoidal perturbation is induced by the sinusoidal profile of the impact surface of the sample by use of two-stage light-gas gun. The oscillatory damping process of the perturbation amplitude is monitored by electric pins. Two damping curves (perturbation amplitude relative to its initial value versus propagated distance relative to the wavelength of sinusoidal perturbation) of aluminum are determined at 78 and 101 GPa. The effective shear viscosity coefficients are deduced to be about 1300 and 800 Pa s based on the Miller and Ahrens analytic solution for viscous fluid.

Ab initio refinement of the thermal equation of state for bcc tantalum: the effect of bonding on anharmonicity.

We report a detailed ab initio study for body-centered-cubic (bcc) Ta within the framework of the quasiharmonic approximation (QHA) to refine its thermal equation of state and thermodynamic properties. Based on the excellent agreement of our calculated phonon dispersion curve with experiment, the accurate thermal equations of state and thermodynamic properties are well reproduced. The thermal equation of state (EOS) and EOS parameters are considerably improved in our work compared with previous results by others. Furthermore, at high temperatures, the excellent agreement of our obtained thermal expansion and Hugoniot curves with experiments greatly verifies the validity of the quasiharmonic approximation at higher temperatures. It is known that pressure suppresses the vibrations of atoms from their equilibrium positions, i.e. the bondings among atoms are strengthened by pressure; for the same temperature, anharmonicity becomes less important at high pressure. Thus the highest valid temperature of the QHA can be reasonably extended to the larger range.

Brillouin scattering studies of liquid argon at high temperatures and high pressures.

Brillouin scattering measurements were performed on liquid argon in a diamond anvil cell at various solidification points up to 503 K. With the measured results from the 60 degree platelet- and 180 degree back-scattering geometries, the sound velocity, refractive index, experimental equation of state, and adiabatic bulk modulus of liquid argon as a function of pressure were determined. The discrepancy between experimental and previous calculated equation of state indicates that the many-body contribution to the density of liquid argon increases with increasing pressure and decreases with increasing temperature. By analyzing the Brillouin spectra in the coexistence of liquid and solid phase regions, the volume change and latent heat of solid-liquid transformation along the equilibrium curve have been also obtained for the first time.

Yield strength of molybdenum at high pressures.

In the diamond anvil cell technology, the pressure gradient approach is one of the three major methods in determining the yield strength for various materials at high pressures. In the present work, by in situ measuring the thickness of the sample foil, we have improved the traditional technique in this method. Based on this modification, the yield strength of molybdenum at pressures has been measured. Our main experimental conclusions are as follows: (1) The measured yield strength data for three samples with different initial thickness (100, 250, and 500 microm) are in good agreement above a peak pressure of 10 GPa. (2) The measured yield strength can be fitted into a linear formula Y=0.48(+/-0.19)+0.14(+/-0.01)P (Y and P denote the yield strength and local pressure, respectively, both of them are in gigapascals) in the local pressure range of 8-21 GPa. This result is in good agreement with both Y=0.46+0.13P determined in the pressure range of 5-24 GPa measured by the radial x-ray diffraction technique and the previous shock wave data below 10 GPa. (3) The zero-pressure yield strength of Mo is 0.5 GPa when we extrapolate our experimental data into the ambient pressure. It is close to the tensile strength of 0.7 GPa determined by Bridgman [Phys. Rev. 48, 825 (1934)] previously. The modified method described in this article therefore provides the confidence in determination of the yield strength at high pressures.

Self-consistent fluid variational theory for the dissociation of dense nitrogen.

The self-consistent fluid variational theory is used to calculate the pressure dissociation of dense nitrogen at high temperatures. The accurate high-pressure and high-temperature effective pair potentials are adopted to describe the intermolecular interactions, which are made to consider molecular dissociation. This paper focuses on a mixture of nitrogen atoms and molecules and is devoted to the study of the phenomenon of pressure dissociation at finite temperature. The equation of state and dissociation degree are calculated from the free-energy functions in the range of temperature of 2000-15 000 K and density of 0.2-3.0 gcm(3), which can be compared with other approaches and experiments.

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer.

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.