Arxiv:1908.10927V1 [Cond-Mat.Mtrl-Sci] 28 Aug 2019 Calculations, and a Variety of Alloying Elements and Crys- 2 Tal Structures Has Been Considered
Total Page:16
File Type:pdf, Size:1020Kb
Impact of anharmonicity on sound wave velocities at extreme conditions Olle Hellman1 and S. I. Simak1 1Department of Physics, Chemistry and Biology (IFM), Link¨opingUniversity, SE-581 83, Link¨oping,Sweden. Theoretical calculations of sound-wave velocities of materials at extreme conditions are of great importance to various fields, in particular geophysics. For example, the seismic data on sound-wave propagation through the solid iron-rich Earth's inner core have been the main source for elucidating its properties and building models. As the laboratory experiments at very high temperatures and pressures are non-trivial, ab initio predictions are invaluable. The latter, however, tend to disagree with experiment. We notice that many attempts to calculate sound-wave velocities of matter at extreme conditions in the framework of quantum-mechanics based methods have not been taking into account the effect of anharmonic atomic vibrations. We show how anharmonic effects can be incorporated into ab initio calculations and demonstrate that in particular they might be non- negligible for iron in Earth's core. Therefore, we open an avenue to reconcile experiment and ab initio theory. Sound-wave velocities of materials at extreme condi- mological values are well constrained, the difference be- tions are important parameters for building models in tween the observations and theory suggests that a simple different physical areas, in particular in geophysics. For model for the inner core based on the commonly assumed example, modeling of Earth's inner core has to rely on phases is wrong. Different and rather complex explana- sound-velocities from seismic data. They, together with tions of this puzzling result have been suggested, includ- other cosmochemical and geochemical data, indicate that ing the partial melting of the inner core.[13] It has been, the inner core should be made out mostly of iron with however, noticed by Voˇcadlo et al.[13] that the disagree- approximately 10 at. % of Ni and some light alloying ment between theory and experiment could in principle elements.[1] The choice and amount of the light elements be accounted for by anelasticity. The simple estimates in is a rather controversial issue.[2, 3] Ref. 13 suggest that the effect of anelasticity should not The structural state of iron at Earth's core conditions overcome 1.5 % and therefore has to be ruled out with (temperature around 6000 K and pressure around 360 the reservation that anelasticity is a very complex issue GPa) is also a matter of debate. For a long time it has that requires material data at the conditions of Earth's widely been accepted that the hexagonal close-packed inner core in order to draw irrefutable conclusions. (hcp) phase of iron is the most likely stable phase in In this letter we show that in contrast to previous as- the inner core (see Tateno et al.[4] and Tateno et al.[5] sumptions and simple models the effect of anelasticity for corresponding experimental claims). Although some due to anharmonic lattice vibrations may have dramatic studies suggest that the face-centered cubic (fcc) [6, 7] impact on the sound-wave anisotropy of crystalline mat- or body-centered cubic (bcc) [3, 8, 9] structures may also ter at high temperatures, in particular in iron at Earth's be stabilized at inner core conditions. The bcc structure core conditions. This is proved via calculations involving has been strongly supported by recent theoretical works the temperature-dependent effective potential (TDEP) [10, 11] and has also been claimed in the latest high- method [14{16], which is based on ab initio molecular dy- pressure experiment employing a laser-heated diamond namics simulations and provides accurate temperature- anvil cell and microstructure analysis combined with in- dependent properties of solids. situ x-ray diffraction.[12] We start from the notion that the standard procedure A possible route to build a model of the solid inner core for calculating the sound-wave velocity in a system like is to match sound velocities in laboratory experiments or high-pressure high-temperature Fe is to obtain second- computation to the observed velocities from seismic data. order elastic constants from molecular dynamics simu- There have been many attempts employing such an ap- lations at the proper temperature and density, and use proach in the framework of quantum-mechanics based these elastic constants in the Christoffel equations arXiv:1908.10927v1 [cond-mat.mtrl-sci] 28 Aug 2019 calculations, and a variety of alloying elements and crys- 2 tal structures has been considered. Remarkably, none of Cqq = ρωqq; (1) these attempts has led to results matching all the exper- where imentally known peculiarities. In particular, the seismi- X cally observed velocities of the S-waves are lower than the Cαγ = c q q : (2) calculated ones (by more than 10 %). The theoretical ve- q αβγδ β δ βδ locities of the P-waves, however, are in good agreement with experiment. Here q is the wavevector of the acoustic wave, cαβγδ elas- As discussed in Ref. 13, if the uncertainties in the seis- tic constants (any index α, β, γ or δ runs over Cartesian 2 a b 0 0 2000 2000 T (K) T (K) 4000 4000 0 0 5 5 V (km/s) 10 6000 V (km/s) 10 6000 FIG. 1. Distribution of spherically averaged sound velocities in hcp Fe at Earth's inner core pressure (360 GPa) as a function of temperature. In panel a the anharmonic contributions are not taken into account, and in panel b they are included. We identify the two regions as primary (Vp) and secondary (Vs) waves. One may see that the anharmonicity has little effect on Vp, but heavily influences the distribution of Vs. x, y and z), and ρ is the density. The sound velocities where ri locates atom i. Inserting Eq. 4 into Eq. 3 and are obtained from the frequencies ! as v = d!=dq, and collecting terms in powers of η gives (to first order) are classified via the polarization vectors . Eq. 1 is in- 2 deed valid, but formally only for a perfectly harmonic ^ X pi X αβ α β X αβ,γδ α β H = + Φij ui uj + Φij ηγδui uj 2mi crystal in the long wavelength limit. In the general case, i ijαβ ijαβγδ the sound waves may be perturbed, diffused and eventu- | {z } | {z } ^ ^ 0 ally changed by anharmonic lattice vibrations. In other H0 H words, the conventional approach assumes a perfect prin- (5) ciple of superposition: any sound wave is unaffected by where the intrinsic vibrations already present in the crystal { X Φαβ,γδ = Φαβγ rδ (6) something that is generally not true. At ambient con- ij ijk k k ditions such effects are usually ignored in all but some pathological cases.[17] However, at Earth's core condi- αβγ Here Φijk are the third order interatomic force con- tions ignoring the effects of anharmonicity can become a stants. H^ 0 results in a correction to the ideal harmonic dangerously crude approximation. waves from a strain, such as an acoustic wave. This cor- Let us consider how anharmonicity may be accurately rection expresses the change in second order force con- taken into account. One starts with the general form of stants upon a strain. The correction is nontrivial and the Hamiltonian describing lattice dynamics[18] is best treated using perturbation theory. First, one ex- 2 presses the reciprocal-space matrix elements correspond- ^ X pi X αβ α β H = + Φij ui uj + ing to a strain perturbation:[19] 2mi i ijαβ (3) iγ jδ X αβγ α β γ X qsq0s0 0 Φ u u u + :::; αβ ~ αβ,γδ iqri+iq rj ijk i j k Vqsq0s0 = p Φij e (7) 4 !qs!q0s0m m ijkαβγ ij,γδ i j where p is the momentum, m the mass, and ui the dis- Here !qs and qs stand for phonon frequencies and polar- placement of atom i. Φ designate interatomic force con- ization vectors at reciprocal momentum q for the phonon stants of increasing order. When this expansion is trun- branch s, respectively. As can be shown (see Ref. 19 and cated at the second order terms, the harmonic approxi- Ref. 20), the Cq tensor in Eq. (2) should now contain 3 mation is obtained, where in the long wavelength limit additional terms due to anharmonicity the Christoffel equations (Eq. 1) can be derived. Let us now introduce strain parameters η and write 1 8 X X αβ γδ !qs(2nq¯s0 + 1) ∆ cαβ,γδ = Vqsq¯s0 Vqsq¯s0 2 2 (8) V ! 0 − ! the displacement ui of an atom i from its equilibrium ~ q ss0 q¯s qs position as 2 X X 2nqs + 1 ∆2c = − V αβ V γδ (9) αβ,γδ V qsq¯s qsq¯s ! ~ q s qs ui = ηri (4) 3 14 the full derivation of general results). However, its prac- a tical usage has been hindered by the inability to calculate 12 accurate interatomic force constants up to the third or- der. The temperature dependent effective potential tech- 10 nique (TDEP)[14{16] allows one to achieve this goal in a 8 rigorous and efficient way. Actual calculations have been done for hcp Fe, which 6 is expected to be the least anharmonic phase out of pos- sible candidates for Earth's core conditions. That is if Sound velocity (km/s) velocity Sound 4 any effect on sound velocities is observed for hcp Fe, it should be expected to be even larger for other iron V 2 p phases. We also notice that the hcp case is not compli- V s cated by collective self-diffusion events predicted for bcc 0 1000 2000 3000 4000 5000 6000 Fe.[10] The calculations have been performed at a pres- Temperature (K) sure of 360 GPa corresponding to the one in Earth's inner core[21] in the temperature range from 0 to 6000 K.