short communications
Pressure effects on the
Debye–Waller factor of ironaInstitute of Research and Development, Duy Tan University, 03 Quang Trung, Hai Chau, Da Nang 550000, Vietnam, bFaculty of Natural Sciences, Duy Tan University, 03 Quang Trung, Hai Chau, Da Nang 550000, Vietnam, cTan Trao University, Km 6, Yen Son, Tuyen Quang 301910, Vietnam, dVNU University of Science, 334 Nguyen Trai, Thanh Xuan, Ha Noi 100000, Vietnam, eHanoi Pedagogical University No 2, Vinh Phuc 280000, Vietnam, and fHong Duc University, 565 Quang Trung, Dong Ve, Thanh Hoa 441430, Vietnam
*Correspondence e-mail: hieuhk@duytan.edu.vn
The pressure effects on atomic mean-square relative displacement characterizing the extended X-ray absorption fine structure (EXAFS) Debye–Waller factor of iron metal have been investigated based on the Debye model. The analytical expressions of the Debye frequency and
Debye–Waller factor have been derived as functions of crystal volume compressibility. Based on the well established equation-of-state including the contributions of the anharmonic and electronic thermal pressures, numerical calculations have been performed for iron up to a pressure of 220 GPa and compared with experimental data when possible. These results show that the Debye frequency increases rapidly with compression, and beyond 150 GPa it behaves as a linear function of pressure. Meanwhile the mean-square relative displacement curve drops robustly with pressure, especially at pressures smaller than 100 GPa. This phenomenon causes the enhancement of signals at high pressure. Reversely, the increasing of temperature will reduce the amplitude of spectra.Keywords: Debye model; EXAFS; Debye–Waller factor; Debye frequency.
1. Introduction
Iron is a main component in the Earth's core. The high-pressure thermodynamic properties of iron are primarily important information for explaining geochemical observations and seismic data, and for modeling the dynamic properties, evolution process, structure and dynamo of the Earth. A number of theoretical and experimental works have been performed to study the structure and thermodynamic properties of iron at extreme conditions, such as first-principles calculations (Stixrude, 2012), simulations (Alfè, 2010), laser-heated diamond-anvil cells (DACs) (Errandonea, 2013), shock-wave experiments (Nguyen & Holmes, 2004) and X-ray absorption spectroscopy (Wang & Ingalls, 1998; Ishimatsu et al., 2016; Torchio et al., 2016).
Extended X-ray absorption fine-structure spectroscopy (EXAFS) (Crozier et al., 1988) is a powerful method for investigation of local structure and thermodynamic properties of materials. It is a result of the interference between outgoing photoelectron waves and back-scattering due to the presence of neighboring atoms. The Debye–Waller factor (DWF) W(k) = (where k is the wavenumber and σ2 is the second cumulant) is one important factor of the oscillation function. The second cumulant σ2 or DWF accounts for the effects of atomic thermal vibrations and influences sensitively the amplitudes. Previously, compression effects on the DWF have been considered in some theoretical approaches (Hung et al., 2011; Hieu & Hung, 2013). Using dynamic compression with multiple shocks, measurements on the compressed iron were performed up to 560 GPa by taking into account anharmonic effects (Ping et al., 2013). In spite of this, there is still no complete theory that allows us to theoretically predict the parameters at high pressure.
In this work, we investigate the pressure effects on the Debye frequency and
DWF of iron metal based on the Debye model. Numerical calculations are performed up to a pressure of 220 GPa and compared with those of previous works when possible.2. Theoretical approach
The Debye–Waller factor W(k) = accounts for the effects of atomic thermal vibrations in The second cumulant σ2 characterizes the variance of the distance distribution corresponding to the parallel mean-square relative displacement (MSRD) of the bond between absorber and backscatterer atoms. Hung et al. (2010), using the anharmonic correlated Debye model, derived an expression for the DWF as a function of temperature T as
where a is the lattice constant, q is the phonon wavenumber, M is the mass of composite atoms, z(q) = , kB is the and ω(q) is the dispersion relation which achieves the maximum value of the Debye frequency ωD at the boundaries q = ±π/a of the first of the linear chain. The phonon frequency ω(q) is expressed as (Hung et al., 2010)
From equation (1), the pressure effects on the MSRD can be investigated indirectly through the evaluation of lattice parameter a and the Debye frequency ωD under compression. The volume-dependent lattice parameter can be determined as a = a0ζ1/3, where ζ is the volume compression ratio ζ = V/V0 (V is the volume of the crystal; a0 and V0 are values of a and V at ambient conditions). Meanwhile the volume dependence of the Debye frequency can be considered through the Grüneisen parameter γG as (Hieu, 2015)
Previous works (Boehler, 1983; Hieu, 2014) showed that at low pressure the Grüneisen parameter reduces gradually when pressure increases. A simple power-law form of the volume-dependent Grüneisen parameter suggested by Graf et al. (2004) is γG = γ0ζp, where p > 0 is a material constant, and γ0 is the Grüneisen parameter at ambient conditions (Errandonea, 2013; Hieu & Ha, 2013). However, this power law does not describe well the Grüneisen parameter of materials at high compression. Another well described form of the Grüneisen parameter following the Al'tshuler et al. (1987) proposition has been chosen as (Dewaele et al., 2006)
where γ∞ is the value of the Grüneisen parameter γG at infinite compression.
By substituting equation (4) into equation (3) and taking the integral, we find the Debye frequency as a function of volume compression ζ = V/V0 as
where ω0 is the value of the Debye frequency of iron at ambient conditions. This quantity can be measured from experiments or calculated from the correlated Debye model (Sevillano et al., 1979).
By substituting equations (5) and (2) into equation (1), we finally derive the volume-dependent MSRD expression as
In order to study the pressure effects on the Debye frequency and atomic MSRD of iron, we need an accurate equation-of-state (EOS) relating the pressure–volume–temperature of the system. In the work of Dewaele et al. (2006), based on the DAC experiments, combined with X-ray diffraction, Hugoniot data and ab initio modeling, a high-pressure high-temperature quasi-hydrostatic EOS of h.c.p.-iron was constructed as
where PV(V, 300 K) is the isothermal (cold) pressure at temperature 300 K and PTH(V, T) is the temperature-dependent (thermal) pressure. These two pressures, PV(V, 300 K) and PTH(V, T), are expressed, using the fitted Vinet EOS (Anderson, 1995) and the formalism derived by Dorogokupets & Oganov (2007), by equations (8) and (9), respectively, as
where B0 and are the isothermal bulk modulus and its first-pressure derivative, respectively; R denotes the The first term of the right-hand side of equation (9) represents the main part of PTH, which is the quasi-harmonic Debye thermal pressure (Anderson, 1995). The second term includes the contributions of the anharmonic and electronic thermal pressure (Alfè et al., 2001).
3. Results and discussion
In this section, numerical calculations are performed for iron up to 220 GPa. Here it should be noted that experimental evidence shows an h.c.p.–b.c.c. transition line in pressure–temperature space from about 95 ± 2 GPa and 2986 ± 79 K to at least 222 ± 6 GPa and 4192 ± 104 K (Hrubiak et al., 2018). The parameters b0, e0, m and n of equation (7) have been derived for iron by fitting first-principles anharmonic and electronic thermal pressures as (Dewaele et al., 2006) b0 = 3.7 × 10−5 K−1, e0 = 1.95 × 10−4 K−1, m = 1.87 and n = 1.339. The fitted values of γ0 and γ∞ are 1.710 and 1.305, respectively. The Debye temperature θ0 at ambient conditions is fixed as 417 K (the Debye frequency ω0 = kBθ0/ℏ = 5.46 × 1013 Hz), the isothermal bulk modulus B0 and its first-pressure derivative are 163.4 GPa and 5.38, respectively (Dewaele et al., 2006).
In Fig. 1, we show the pressure dependence of the Debye frequency of iron up to 220 GPa. Experimental results of Wang & Ingalls (1998) and Anderson et al. (2001) have also been presented for comparison. Anderson et al. (2001) have extracted this thermodynamic quantity from the atomic mean-square displacement which was obtained by X-ray structural using the on powder diffraction measurements. As can be observed from Fig. 1, our theoretical calculation follows the experimental data (Wang & Ingalls, 1998; Anderson et al., 2001) up to a pressure of about 200 GPa. Beyond 150 GPa, the slope of the Debye frequency curve becomes reducible and it is almost linear proportional to pressure. The initial slope of this curve is about 4.4 K GPa−1, but, from a pressure of 150 GPa and beyond, the average slope is about 1.2 K GPa−1.
Using the derived Debye frequency ωD, we can investigate the pressure effects on the parallel MSRD of iron. The pressure dependence of the atomic MSRD characterizing the DWF of iron metal at different temperatures (T = 300 K, 500 K, 700 K, 900 K) is shown in Fig. 2. From this figure we can see that the MSRD curves of iron drop rapidly (especially at pressures smaller than 100 GPa). This phenomenon forecasts the enhancement of signals at high pressure. Above 200 GPa, atomic MSRD curves are readily seen to be in approximately linear proportion to pressure. The initial slopes of MSRD curves dσ2/dP at ambient pressure are −1.12 × 10−4, −1.88 × 10−4, −2.63 × 10−4 and −3.39 × 10−4 Å2 GPa−1 corresponding to temperatures of 300, 500, 700 and 900 K. Meanwhile, the slopes of MSRD curves at temperature 300 K, 500 K, 700 K and 900 K and at pressure 200 GPa are −4.87 × 10−6, −8.23 × 10−6, −1.17 × 10−5 and −1.52 × 10−5 Å2 GPa−1, respectively. Due to the lack of available experimental data and theoretical calculations, the comparison for pressure-dependent MSRD of iron has been neglected.
Our calculations predict that the Debye frequency increases robustly with pressure, while the et al., 2019). As can be realized from Fig. 2, the DWF enhances robustly when increasing the temperature. This behavior of atomic MSRD is palpable evidence for significant contribution of thermal disorder in the high-temperature region.
DWF of iron drops rapidly. The rise of the phonon frequency and the reduction of atomic MSRD can be explained by the constraint of atomic vibrations due to the increasing of compression. It is worth mentioning that the DWF is susceptible to both static disorder and thermal disorder. The first one represents the structural disorder caused by alloying or strain while the second one represents the disorder due to the thermal vibrations of atoms. The latter gives us the information of the bond-stretching force constant between absorber and backscatterer atoms coming close to the dynamic properties of them (HongUsing the derived . Here it should be noted that, for simplicity, we have neglected the pressure-induced and temperature-induced phase shift of spectroscopy caused by the variation of temperature and pressure. As can be clearly observed from this figure, at a fixed pressure, increasing the temperature will reduce the amplitude of spectra. Reversely, at a specified temperature, an increasing of pressure will intensify the signals.
DWF, we calculate the theoretical spectra of iron at temperatures 300 K and 900 K, and at pressures 50 GPa and 200 GPa. These four calculated curves are shown in Fig. 34. Conclusions
In conclusion, the Debye model has been developed to investigate the pressure dependence of the
DWF of iron based on the volume-dependent Grüneisen parameter. We derived the analytical expressions of the Debye frequency and parallel MSRD as functions of volume compression. Numerical calculations for iron up to a pressure of 220 GPa show that the Debye frequency increases rapidly with compression, and beyond 150 GPa it behaves like a linear function of pressure. Meanwhile, the MSRD curve drops robustly with pressure, especially at pressures smaller than 100 GPa. The decreasing of MSRD of iron causes the enhancement of signals at high pressure. In contrast, the increasing of temperature will fall off the spectra. This work can be seen as a useful reference for future experimental measurements of high-pressure MSRD factors of materials.Acknowledgements
The authors would like to acknowledge Professor Anatoly Belonoshko and an anonymous reviewer for their useful comments and suggestions.
Funding information
Funding for this research was provided by: the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.01-2017.343 (to Hieu Ho).
References
Alfè, D. (2010). Rev. Mineral. Geochem. 71, 337–354. Google Scholar
Alfè, D., Price, G. D. & Gillan, M. J. (2001). Phys. Rev. B, 64, 045123. Google Scholar
Al'tshuler, L. V., Brusnikin, S. E. & Kuz'menkov, E. A. (1987). J Appl. Mech. Tech. Phys. 28, 129–141. Google Scholar
Anderson, L., Dubrovinsky, L., Saxena, S. K. & LeBihan, T. (2001). Geophys. Res. Lett. 28, 399–402. Web of Science CrossRef CAS Google Scholar
Anderson, O. L. (1995). Equations of State for Solids in Geophysics and Ceramic Science, Oxford Monographs on Geology and Geogphysics. Oxford University Press. Google Scholar
Boehler, R. (1983). Phys. Rev. B, 27, 6754–6762. CrossRef CAS Web of Science Google Scholar
Crozier, E. D., Rehr, J. J. & Ingalls, R. (1988). X-ray Absorption: Principles, Applications, Techniques of EXAFS, SEXAFS and XANES, 1st ed., edited by D. C. Koningsberger and R. Prins. Wiley-Interscience. Google Scholar
Dewaele, P., Loubeyre, F., Occelli, M., Mezouar, M., Dorogokupets, P. I. & Torrent, M. (2006). Phys. Rev. Lett. 97, 215504. Web of Science CrossRef PubMed Google Scholar
Dorogokupets, P. I. & Oganov, A. R. (2007). Phys. Rev. B, 75, 024115. Web of Science CrossRef Google Scholar
Errandonea, D. (2013). Phys. Rev. B, 87, 054108. Web of Science CrossRef Google Scholar
Graf, M. J., Gree, C. W. & Boettger, J. C. (2004). AIP Conf. Proc. 706, 65–68. CrossRef CAS Google Scholar
Hieu, H. K. (2014). Vacuum, 109, 184–186. Web of Science CrossRef CAS Google Scholar
Hieu, H. K. (2015). Vacuum, 120, 13–16. Web of Science CrossRef CAS Google Scholar
Hieu, H. K. & Ha, N. N. (2013). AIP Adv. 3, 112125. Google Scholar
Hieu, H. K. & Hung, V. V. (2013). High. Press. Res. 33, 768–776. Web of Science CrossRef CAS Google Scholar
Hong, N. T., Hieu, H. K., Duc, N. B., Phuong, D. D., Tuyen, N. V. & Khoa, D. Q. (2019). Vacuum, 163, 210–215. Web of Science CrossRef CAS Google Scholar
Hrubiak, R., Meng, Y. & Shen, G. (2018). arXiv:1804.05109 [physics. geo-ph]. Google Scholar
Hung, N. B., Bao Trung, N. & Kirchner, B. (2010). Physica B, 405, 2519–2525. Web of Science CrossRef Google Scholar
Hung, N. V., Hung, H. K., Hieu, H. K. & Frahm, R. R. (2011). Physica B, 406, 456–460. Google Scholar
Ishimatsu, N., Kawamura, M., Mizumaki, H., Maruyama, H., Sumiya, H. & Irifune, T. (2016). High. Press. Res. 36, 381–390. Web of Science CrossRef CAS Google Scholar
Nguyen, J. H. & Holmes, N. C. (2004). Nature, 427, 339–342. Web of Science CrossRef PubMed CAS Google Scholar
Ping, F., Coppari, F., Hicks, B., Yaakobi, B., Fratanduono, S., Hamel, S., Eggert, J. R., Rygg, R. F., Smith, D. C., Swift, D. G., Braun, T. R., Boehly, T. R. & Collins, G. W. (2013). Phys. Rev. Lett. 111, 065501. Web of Science CrossRef PubMed Google Scholar
Sevillano, H., Meuth, H. & Rehr, J. J. (1979). Phys. Rev. B, 20, 4908–4911. CrossRef CAS Web of Science Google Scholar
Stixrude, L. (2012). Phys. Rev. Lett. 108, 055505. Web of Science CrossRef PubMed Google Scholar
Torchio, F., Occelli, O., Mathon, A., Sollier, E., Lescoute, L., Videau, T., Vinci, A., Benuzzi-Mounaix, J., Headspith, W., Helsby, S., Bland, D., Eakins, D., Chapman, S., Pascarelli, S. & Loubeyre, P. (2016). Sci. Rep. 6, 26402. Web of Science CrossRef PubMed Google Scholar
Wang, F. M. & Ingalls, R. (1998). Phys. Rev. B, 57, 5647–5654. Web of Science CrossRef CAS Google Scholar
© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.