Research Article  Open Access
N. K. Omboga, C. O. Otieno, P. W. O. Nyawere, "First Principle Study of the Mechanical Properties and Phonon Dispersion of the Iron Pnictide Compound EuFe_{2}As_{2}", The Scientific World Journal, vol. 2020, Article ID 5986073, 6 pages, 2020. https://doi.org/10.1155/2020/5986073
First Principle Study of the Mechanical Properties and Phonon Dispersion of the Iron Pnictide Compound EuFe_{2}As_{2}
Abstract
We present results on the first principle study of the elastic constants and the phonon dispersion of EuFe_{2}As_{2} at zero pressure. The groundstate energy calculations were performed within Density Functional Theory (DFT) and the generalized gradient approximation using the pseudopotential method with planewave basis sets. The projector augmentedwave (PAW) pseudopotentials were used in our calculation. The open source code QUANTUM ESPRESSSO was used with its pseudopotential database. The study on the elastic constants at zero pressure was a clear indication that the compound is mechanically stable, and the phonon dispersion study also indicated that the compound is dynamically stable. The elastic constants and mechanical properties also led to the conclusion that the compound is ductile and anisotropic.
1. Introduction
Iron pnictides belong to the “122” family of compounds and have in studies of the recent past shown to exhibit high temperature superconductivity [1–4]. These are basically type II superconductors and they consist of iron arsenide layers. By substitution of europium with sodium or potassium, the material becomes superconducting. Also, in the application of external pressure or doping, superconductivity is achieved [4]. There is a reduction in phase transition temperature when the compound is doped or adequate pressure is applied on it, making it superconducting [5].
The iron pnictide compounds are also unique because their phase transition to superconductivity and magnetic properties are connected to structural properties and conceal the key to comprehending the basic properties of these materials [6]. According to [7], EuFe_{2}As_{2} experiences a phase transition from the tetragonal phase to orthorhombic phase at a pressure of 4.3 GPa, which is maintained until a pressure of 11 GPa where it transits to a collapsed tetragonal phase. This pressure is maintained up to around 35 GPa. The conductivity of iron pnictides shows a highly anisotropic characteristic [8].
Phase transition at low temperature is common in the iron pnictide family of compounds. Experiments performed at high pressure reveal characteristics which link the electronic and structural properties to the superconducting properties of the compound [6].
Superconductivity and ferromagnetism coexist in europium diiron diarsenide doped with phosphorous [5]. The iron spin density and the europium () antiferromagnetically make europium diiron diarsenide outstanding among the 122 iron pnictide superconducting compounds. The FeAs layers in the compound europium diiron diarsenide are seen to be the source of the superconductivity [4]. This compound is also unique because of the additional magnetic moment of the localized given that other 122 compounds undergo only the SDW transition.
Application of hydrostatic pressure on indicates that the SDW transition is continuously suppressed on pressure application. Superconductivity is achieved on suppression of the magnetism. Superconductivity transition temperature occurs at 30 K at a pressure of 2.8 GPa [9].
Other studies in similar iron pnictides have been done and documented, although mostly experimental. Barium diiron diarsenide exhibits superconductivity on distorting the structure through application of pressure [10]. Superconductivity is achieved at a pressure of 5.5 GPa and temperature of 30.5 K. Studies on CaFe_{2}As_{2} have also revealed that it exhibits superconductivity on application of pressure and up to temperatures of 40 K [11]. Experimental results on europium diiron diarsenide have also indicated that superconductivity occurs below temperatures of 33 K [12].
In the ground state, pnictides show spin density wave. The main differentiating factor between Cuprates (which were initially the major materials under the study of superconductivity) is the behavior of the electrons that conduct [13]. In iron pnictides, the Fe 3d prevails over the electrons near the Fermi level. The character is largely contained in the states near the Fermi level [13]. EuFe_{2}As_{2} has unique characteristics, since the Eu^{2+} ions have 4f electrons and the combined spin of the electrons is 7/2 [14]. In this article, we present the elastic constants and phonon dispersion of the iron pnictide EuFe_{2}As_{2} as studied from first principles. Seismic wave velocities can be explained using elastic constants and also the comprehending of chemical bonds is facilitated by the elastic constants [15]. This paper is arranged such that Section 2 presents computational methodology, Section 3 reports on results and discussions, and Section 4 is conclusion.
2. Computational Details
The open source code QUANTUM ESPRESSO was used in this study. Projector augmentedwave (PAW) pseudopotentials of the Perdew, Erzenhoff, and Burke (PBE) were employed to cater for the core electrons [16]. The generalized gradient approximation (GGA) [17] functional was employed for the exchange correlation. The Kohn–Sham equations were solved iteratively by running a selfconsistent cycle. The unit cell used to generate the super cells comprised 5 atoms: iron and arsenic had two atoms each and europium had one atom. Nine super cells of 2 × 2× 2 dimensions were created for the phonon dispersion study, and the number of atoms in each generated super cell was 40. The phonon calculation of finite differences was done using Phonopy code. A smearing of 0.2 eV of the MethfesselPaxton scheme was used. We used an energy cutoff of 30 Ry, which dictated the number of plane waves to be employed. The elastic constants were obtained using the postprocessing code THERMO_PW within Quantum ESPRESSO. Variable cell relaxation was done before undertaking all the above calculations to ensure the material is stressfree.
Bulk modulus is a measure of the ability of a material to withstand variations in volume when compressed on all sides [18]. The bulk modulus in this study was computed usingwhere .
And, .
Then, .
The shear modulus (G) is the measure of resistance to reversible deformations upon shear stress [18].
For the shear modulus,
Therefore, .
The Young’s modulus E [18] and Poisson’s ratio [18] were obtained, respectively, using the equations below:
The study of phonons assists in understanding properties of materials such as thermal and electrical conductivity. It enables quantization of the energy of vibration. A phonon refers to a kind of lattice vibration in a crystal, where the particles vibrate at the same single frequency [19]. From the phonon dispersion calculation, one can be able to calculate the critical temperature of the material in question. A dispersion relation refers to the relationship between the frequency of vibration and the wave vector and this relationship is given as follows:where k is wave vector, is frequency of vibration, and is velocity of sound [19, 20].
3. Results and Discussion
Structural optimization was done by doing a variable cell relaxation calculation and also optimizing Kpoints and cutoff energy. The lattice parameters of the tetragonal crystal structure were also calculated being guided by parameters from previous experimental studies on the compound and optimized to obtain a relaxed structure, [21], and the parameters are as shown in Table 1.
From Table 1, it is realized that optimized and experimental lattice parameters are in good agreement [22].
The crystal structure of this iron pnictide compound EuFe_{2}As_{2} is as shown in Figure 1:
Mechanical properties such as elastic properties guide in understanding more about the nature of the force in the crystal. The elastic constants c_{ij} of the compound were computed at T = 0 K and = 0 GPa. There exist six elastic constants in the body centered tetragonal structure [22], and they are computed and reported for the first time in Table 2.

Given the BornHuang criteria [24] shown below, a stable tetragonal structure should satisfy the following criteria:(i)C_{ii} > 0 (i = 1, 3, 4, 6)(ii)C_{11} + C_{33} − 2C_{13} > 0(iii)2 (C_{11} + C_{12}) + C_{33} + 4C_{13} > 0(iv)C_{11} − C_{12} > 0
The values reported in Table 2 are all positive and satisfy all the above four criteria, hence proving that this tetragonal compound is stable. C_{11} > C_{33} as seen from Table 2, implying that the crystal is more stiff along 100 and 010 axes, results which are in agreement with the study on a similar iron pnictide BaPd_{2}As_{2} [25]. C_{12}–C_{44} gives the Cauchy pressure, which, if positive, makes the compound ductile. From Table 2, the Cauchy pressure is positive; hence the material is ductile, a result similar to that of a related iron pnictide [25].
In our earlier studies, we reported the bulk, shear, and Young’s modulus, and Poisson’s ratio values of EuFe_{2}As_{2} at zero pressure; we tabulate these earlier results here for completeness [21] in Table 3. Bulk, shear, and Young’s modulus are obtained using the Voigt–Reuss–Hill approximation method.

Pugh’s ratio, [26] indicates the following:
for ductile materials and for brittle materials.
The value from Table 3 gives a value of 2.53, implying that the compound is ductile. The value of Pugh’s ratio for this compound is similar to that of BaPd_{2}As_{2}, an iron pnictide of similar structure [25].
Using the Frantserich method, ductile and brittle materials are classified on the basis of their Poisson’s ratio. The value acts as the border of ductility and brittleness.
If , then the material is ductile and if , the material is brittle. Since , as seen from Table 3, this is further proof that the compound is indeed ductile. Also, Poisson’s ratio is in close range with that obtained in the study of the iron pnictide BaPd_{2}As_{2} [25]. Poisson’s ratio also suggests that the compound is metallic, which is in good agreement with a similar pnictide in [25].
When Poisson’s ratio is in the range 0.25–0.5, it implies that the forces present in the compound are central [25]. The anisotropic character was also calculated usingwhere if = 0, the compound is totally isotropic.
In this case, the calculations indicate that the material is anisotropic since the value is 0.408. This is in good agreement with past studies of compounds in the pnictide family that have also been proved to be anisotropic as in [25]. The bulk, shear, and Young’s moduli values increased upon application of increasing pressure on the compound. For instance, Poisson’s ratio, bulk modulus, shear modulus, and Young’s modulus at a pressure of 1 GPa were 0.36100, 1196.9 GPa, 326.6 GPa, and 889.0 GPa, respectively, showing a progressive increase up to 1 GPa.
The Fermi energy of increasing pressure up to 0.8 GPa was obtained and plotted, and the result is as seen in Figure 2.
Figure 2 indicates that the population of charge carries with respect to density of states increases; hence, more electrons are made available for electrical conductivity, enhancing electrical conductivity in the compound. The Fermi energy also is important in describing the electronic structure properties of a material.
The change in volume of the crystal upon application of pressure was also obtained, and a graph of volume against pressure was plotted as seen in Figure 3.
4. Phonon Dispersion
There are two branches in phonon dispersion relation: the acoustic mode which is the lower mode and the optical mode which is the upper mode [27]. The acoustic mode refers to the inphase vibration mode while the optical mode refers to the outofphase vibration. The optical phonons get their name from the fact that they get excited by the radiation of infrared in crystals that are ionic [28]. There are two types of phonon calculation methods: the frozen phonon method and the Density Functional Perturbation Theory [29, 30] as implemented in PWSCF. Figure 4 shows the phonon dispersion for the iron pnictide from the current study.
In the study of phonons, there exist three modes that are associated with each mode number n [19], taking into account the x, y, and z axes [31]. The number of acoustic modes is usually three for crystals whose number of atoms is equal to or greater than two, and the optical modes are given by 3N3 [28]. Therefore, given that europium diiron diarsenide has got five atoms, there are 15 modes of vibration expected. As seen in Figure 4, there are 15 vibration modes with twelve being optical modes and three being acoustic modes. The acoustic modes converge at the gamma high symmetry point. Acoustic modes vibrate at a slower frequency and are in the same phase with the unit cell. Optical modes of vibration have a higher frequency compared to acoustic modes and two neighboring atoms vibrate in a direction opposite to each other. In the acoustic mode, the two adjacent atoms will vibrate together in the same direction. Phonon dispersions are computed along a given line of high symmetry points. The above information therefore confirms that the compound is dynamically stable.
The phonon density of states was also plotted and it is as shown in Figure 5.
A system is considered to be dynamically stable at equilibrium if the potential energy is always increasing for any combination of displacement of atoms; therefore, phonons should have nonnegative and real frequencies for stability [32]. Negative frequencies imply that the potential energy reduces; hence, the system is unstable. Phonon frequencies arise as a result of the displacement of atoms in a given crystal from the rest position, which in turn makes the forces rise [32]. It is important to establish the number of normal modes that are neighboring a certain phonon energy; these details are necessary when studying thermal and electrical conductivity and also establishing the critical temperature of superconducting materials [33]. The Debye temperature is a constant that is associated with the highest allowed mode of vibration [34]. The Debye temperature in this study was 436.454 K and the average Debye sound velocity was 3330.336 m/s. A Debye temperature of above 400 K implies that the crystal’s thermal conductivity is high [35], and since the compound we were studying had a Debye temperature of above 400 K, we concluded that the thermal conductivity is high. Note that, for temperatures that are below the Debye temperature, the heat capacity of the compound rises with the temperature cube and for temperatures above the Debye temperature, the heat capacity of the crystal remains constant; that is to say, it no longer depends on temperature. The Debye temperature and the heat capacity are directly proportional [36].
5. Conclusion
We presented results on the elastic constants and phonon dispersion of the iron pnictide EuFe_{2}As_{2}, studied from first principles. For the first time, we have reported elastic constants and Debye temperature of this iron pnictide. From the study, it can be concluded that the compound is mechanically and dynamically stable. The iron pnictide is ductile and anisotropic. The compound is also a good thermal conductor deduced through its Debye temperature.
Data Availability
All the data in this study are available and can be made available on request with the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to thank the Center for High Performance Computing (CHPC), South Africa, for providing the computational resources which contributed to the success of this study. The authors also recognize the efforts of Dr. George Manyali of Kaimosi Friends University College for the insightful discussion offered.
References
 K. Igawa, H. Okada, H. Takahashi et al., “Pressureinduced superconductivity in iron pnictide compound SrFe_{2}As_{2},” Journal of the Physical Society of Japan, vol. 78, no. 2, Article ID 025001, 2009. View at: Publisher Site  Google Scholar
 W. Li, J.X. Zhu, Y. Chen, and C. S. Ting, “Firstprinciples calculations of the electronic structure of ironpnictide EuFe_{2}AsP_{2} superconductors: evidence for antiferromagnetic spin order,” Physical Review B, vol. 86, no. 15, Article ID 155119, 2012. View at: Publisher Site  Google Scholar
 I. I. Mazin, “Superconductivity gets an iron boost,” Nature, vol. 464, no. 7286, pp. 183–186, 2010. View at: Publisher Site  Google Scholar
 C. MicleaM. Nicklas, H. S. Jeevan et al., “Evidence for a reentrant superconducting state in EuFe_{2}As_{2} under pressure,” Physical Review B, vol. 79, no. 21, Article ID 212509, 2009. View at: Publisher Site  Google Scholar
 T. Terashima, M. Kimata, H. Satsukawa et al., “EuFe_{2}As_{2} under high pressure: an antiferromagnetic bulk superconductor,” Journal of the Physical Society of Japan, vol. 78, no. 8, Article ID 083701, 2009. View at: Publisher Site  Google Scholar
 N. Colonna, G. Profeta, A. Continenza, and S. Massidda, “Structural and magnetic properties of CaFe_{2}As_{2} and BaFe_{2}As_{2} from firstprinciples density functional theory,” Physical Review B, vol. 83, no. 9, Article ID 094529, 2011. View at: Publisher Site  Google Scholar
 W. O. Uhoya, G. M. Tsoi, Y. K. Vohra, M. A. McGuire, and A. S. Sefat, “Structural phase transitions in EuFe_{2}As_{2} superconductor at low temperatures and high pressures,” Journal of Physics: Condensed Matter, vol. 23, no. 36, Article ID 365703, 2011. View at: Publisher Site  Google Scholar
 A. L. Ivanovskii, “New superconductors based on (Ca, Sr, Ba) Fe_{2}As_{2} ternary arsenides: synthesis, properties, and simulation,” Journal of Structural Chemistry, vol. 50, no. 3, pp. 539–551, 2009. View at: Publisher Site  Google Scholar
 N. NiS. Nandi, A. Kreyssig et al., “Firstorder structural phase transition in CaFe_{2}As_{2},” Physical Review B, vol. 78, no. 1, Article ID 014523, 2008. View at: Publisher Site  Google Scholar
 S. A. J. Kimber, A. Kreyssig, Y.Z. Zhang et al., “Similarities between structural distortions under pressure and chemical doping in superconducting BaFe_{2}As_{2},” Nature Materials, vol. 8, no. 6, pp. 471–475, 2009. View at: Publisher Site  Google Scholar
 W. Uhoya, G. Tsoi, Y. K. Vohra et al., “Anomalous compressibility effects and superconductivity of EuFe_{2}As_{2} under high pressures,” Journal of Physics: Condensed Matter, vol. 22, no. 29, Article ID 292202, 2010. View at: Publisher Site  Google Scholar
 P. L. Paulose, H. S. Jeevan, C. Geibel, Z. Hossain, and Anupam, “Superconductivity and magnetism in Kdoped EuFe_{2}As_{2},” Journal of Physics: Condensed Matter, vol. 21, no. 26, Article ID 265701, 2009. View at: Publisher Site  Google Scholar
 K. Maiti, “Electronic structure of Febased superconductors,” Pramana, vol. 84, no. 6, pp. 947–956, 2015. View at: Publisher Site  Google Scholar
 A. Oleaga, A. Salazar, A. Thamizhavel, and S. K. Dhar, “Thermal properties and Ising critical behavior in EuFe_{2}As_{2},” Journal of Alloys and Compounds, vol. 617, pp. 534–537, 2014. View at: Publisher Site  Google Scholar
 P. W. O. Nyawere, N. W. Makau, and G. O. Amolo, “Firstprinciples calculations of the elastic constants of the cubic, orthorhombic and hexagonal phases of BaF_{2},” Physica B: Condensed Matter, vol. 434, pp. 122–128, 2014. View at: Publisher Site  Google Scholar
 P. Jomo, C. Otieno, and P. Nyawere, “Ab initio highpressure study of semiconductormetal phase transition of the chalcogenide compound KPSe_{6},” Advances in Condensed Matter Physics, vol. 2020, Article ID 3407141, p. 9, 2020. View at: Publisher Site  Google Scholar
 J. O. Agora, C. Otieno, P. W. O. Nyawere, and G. S. Manyali, “Ab initio study of pressure induced phase transition, structural and electronic structure properties of superconducting perovskite compound GdBa_{2}Cu_{3}O_{7x},” Computational Condensed Matter, vol. 23, Article ID e00461, 2020. View at: Publisher Site  Google Scholar
 M. H. Sadd, Elasticity: Theory, Applications, and Numerics, Academic Press, Cambridge, MA, USA, 2009.
 S. H. Simon, The Oxford Solid State Basics, Oxford University, England, UK, 2013.
 P. Y. Yu and M. Cardona, Fundamentals of Semiconductors: Physics and Materials Properties, Springer, New York, NY, USA, 1996.
 N. K. Omboga and C. Otieno, “Structural and electronic properties of the iron pnictide EuFe_{2}As_{2} from first principles,” Journal of Physics Communications, vol. 4, no. 2, Article ID 025007, 2019. View at: Publisher Site  Google Scholar
 M. S. Ali, M. A. Rahman, and M. Z. Rahaman, “A theoretical investigation of ThCr_{2}Si_{2}type Pdbased superconductors XPd_{2}Ge_{2} (X = Ca, Sr, La, Nd),” Physica C: Superconductivity and its Applications, vol. 561, pp. 35–44, 2019. View at: Publisher Site  Google Scholar
 P. L. Alireza, Y. T. C. Ko, J. Gillett et al., “Superconductivity up to 29 K in SrFe_{2}As_{2} and BaFe_{2}As_{2} at high pressures,” Journal of Physics: Condensed Matter, vol. 21, no. 1, Article ID 012208, 2009. View at: Publisher Site  Google Scholar
 S. Piskunov, E. Heifets, R. I. Eglitis, and G. Borstel, “Bulk properties and electronic structure of SrTiO_{3}, BaTiO_{3}, PbTiO_{3} perovskites: an ab initio HF/DFT study,” Computational Materials Science, vol. 29, no. 2, pp. 165–178, 2004. View at: Publisher Site  Google Scholar
 F. Parvin and S. H. Naqib, “Structural, elastic, electronic, thermodynamic, and optical properties of layered BaPd_{2}As_{2} pnictide superconductor: a first principles investigation,” Journal of Alloys and Compounds, vol. 780, pp. 452–460, 2019. View at: Publisher Site  Google Scholar
 S. F. Pugh, “XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals,” The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, vol. 45, no. 367, pp. 823–843, 1954. View at: Publisher Site  Google Scholar
 P. Misra, Physics of Condensed Matter, Academic Press, Cambridge, MA, USA, 2011.
 G. D. Mahan, Condensed Matter in a Nutshell, Princeton University Press, New Jersey, NY, USA, 2011.
 A. Dal Corso, S. Baroni, R. Resta, and S. de Gironcoli, “Ab initiocalculation of phonon dispersions in II–VI semiconductors,” Physical Review B, vol. 47, no. 7, pp. 3588–3592, 1993. View at: Publisher Site  Google Scholar
 K. Karch, P. Pavone, W. Windl, O. Schütt, and D. Strauch, “Ab initiocalculation of structural and latticedynamical properties of silicon carbide,” Physical Review B, vol. 50, no. 23, pp. 17054–17063, 1994. View at: Publisher Site  Google Scholar
 J. Birman, Theory of Crystal Space Groups and Lattice Dynamics, Springer, Berlin, Germany, 1984.
 A. Togo and I. Tanaka, “First principles phonon calculations in materials science,” Scripta Materialia, vol. 108, pp. 1–5, 2015. View at: Publisher Site  Google Scholar
 F. Giustino, Materials Modelling Using Density Functional Theory: Properties and Predictions, Oxford University Press, England, UK, 2014.
 T. L. Hill, An Introduction to Statistical Thermodynamics, Courier Corporation, Chicago, IL, USA, 1986.
 I.M. Low, Advances in Science and Technology of Mn+ 1AXn Phases, Elsevier, Cambridge, MA, USA, 2012.
 P. Debye, “Zur Theorie der spezifischen Wärmen,” Annalen der Physik, vol. 344, no. 14, pp. 789–839, 1912. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 N. K. Omboga et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.