The Scientific World Journal

The Scientific World Journal / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 5986073 |

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 EuFe2As2", The Scientific World Journal, vol. 2020, Article ID 5986073, 6 pages, 2020.

First Principle Study of the Mechanical Properties and Phonon Dispersion of the Iron Pnictide Compound EuFe2As2

Academic Editor: Aggelos Avramopoulos
Received19 Jun 2020
Revised14 Sep 2020
Accepted16 Sep 2020
Published26 Sep 2020


We present results on the first principle study of the elastic constants and the phonon dispersion of EuFe2As2 at zero pressure. The ground-state energy calculations were performed within Density Functional Theory (DFT) and the generalized gradient approximation using the pseudopotential method with plane-wave basis sets. The projector augmented-wave (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 [14]. 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], EuFe2As2 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 CaFe2As2 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]. EuFe2As2 has unique characteristics, since the Eu2+ 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 EuFe2As2 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 augmented-wave (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 self-consistent 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 Methfessel-Paxton scheme was used. We used an energy cut-off of 30 Ry, which dictated the number of plane waves to be employed. The elastic constants were obtained using the post-processing code THERMO_PW within Quantum ESPRESSO. Variable cell relaxation was done before undertaking all the above calculations to ensure the material is stress-free.

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 K-points and cut-off 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.

ParameterThis workExperimentalReference

a0 = b0 (ang)3.9453.916[11]
c0 (ang)12.84612.052[11]

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 EuFe2As2 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 cij 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.

CijValue (GPa)


Given the Born-Huang criteria [24] shown below, a stable tetragonal structure should satisfy the following criteria:(i)Cii > 0 (i = 1, 3, 4, 6)(ii)C11 + C33 − 2C13 > 0(iii)2 (C11 + C12) + C33 + 4C13 > 0(iv)C11 − C12 > 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. C11 > C33 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 BaPd2As2 [25]. C12C44 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 EuFe2As2 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.

PropertyVoigt approximationReuss approximationVoigt–Reuss–Hill average

Bulk modulus (B)346.939 GPa346.352 GPa346.646 GPa
Young’s modulus (E)375.378 GPa350.224 GPa362.801 GPa
Shear modulus (G)142.224 GPa131.518 GPa136.871 GPa
Poisson ratio (n)0.319670.331470.3254

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 BaPd2As2, 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 BaPd2As2 [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 in-phase vibration mode while the optical mode refers to the out-of-phase 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 3N-3 [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 non-negative 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 EuFe2As2, 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.


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.


  1. K. Igawa, H. Okada, H. Takahashi et al., “Pressure-induced superconductivity in iron pnictide compound SrFe2As2,” Journal of the Physical Society of Japan, vol. 78, no. 2, Article ID 025001, 2009. View at: Publisher Site | Google Scholar
  2. W. Li, J.-X. Zhu, Y. Chen, and C. S. Ting, “First-principles calculations of the electronic structure of iron-pnictide EuFe2AsP2 superconductors: evidence for antiferromagnetic spin order,” Physical Review B, vol. 86, no. 15, Article ID 155119, 2012. View at: Publisher Site | Google Scholar
  3. I. I. Mazin, “Superconductivity gets an iron boost,” Nature, vol. 464, no. 7286, pp. 183–186, 2010. View at: Publisher Site | Google Scholar
  4. C. MicleaM. Nicklas, H. S. Jeevan et al., “Evidence for a reentrant superconducting state in EuFe2As2 under pressure,” Physical Review B, vol. 79, no. 21, Article ID 212509, 2009. View at: Publisher Site | Google Scholar
  5. T. Terashima, M. Kimata, H. Satsukawa et al., “EuFe2As2 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
  6. N. Colonna, G. Profeta, A. Continenza, and S. Massidda, “Structural and magnetic properties of CaFe2As2 and BaFe2As2 from first-principles density functional theory,” Physical Review B, vol. 83, no. 9, Article ID 094529, 2011. View at: Publisher Site | Google Scholar
  7. W. O. Uhoya, G. M. Tsoi, Y. K. Vohra, M. A. McGuire, and A. S. Sefat, “Structural phase transitions in EuFe2As2 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
  8. A. L. Ivanovskii, “New superconductors based on (Ca, Sr, Ba) Fe2As2 ternary arsenides: synthesis, properties, and simulation,” Journal of Structural Chemistry, vol. 50, no. 3, pp. 539–551, 2009. View at: Publisher Site | Google Scholar
  9. N. NiS. Nandi, A. Kreyssig et al., “First-order structural phase transition in CaFe2As2,” Physical Review B, vol. 78, no. 1, Article ID 014523, 2008. View at: Publisher Site | Google Scholar
  10. S. A. J. Kimber, A. Kreyssig, Y.-Z. Zhang et al., “Similarities between structural distortions under pressure and chemical doping in superconducting BaFe2As2,” Nature Materials, vol. 8, no. 6, pp. 471–475, 2009. View at: Publisher Site | Google Scholar
  11. W. Uhoya, G. Tsoi, Y. K. Vohra et al., “Anomalous compressibility effects and superconductivity of EuFe2As2 under high pressures,” Journal of Physics: Condensed Matter, vol. 22, no. 29, Article ID 292202, 2010. View at: Publisher Site | Google Scholar
  12. P. L. Paulose, H. S. Jeevan, C. Geibel, Z. Hossain, and Anupam, “Superconductivity and magnetism in K-doped EuFe2As2,” Journal of Physics: Condensed Matter, vol. 21, no. 26, Article ID 265701, 2009. View at: Publisher Site | Google Scholar
  13. K. Maiti, “Electronic structure of Fe-based superconductors,” Pramana, vol. 84, no. 6, pp. 947–956, 2015. View at: Publisher Site | Google Scholar
  14. A. Oleaga, A. Salazar, A. Thamizhavel, and S. K. Dhar, “Thermal properties and Ising critical behavior in EuFe2As2,” Journal of Alloys and Compounds, vol. 617, pp. 534–537, 2014. View at: Publisher Site | Google Scholar
  15. P. W. O. Nyawere, N. W. Makau, and G. O. Amolo, “First-principles calculations of the elastic constants of the cubic, orthorhombic and hexagonal phases of BaF2,” Physica B: Condensed Matter, vol. 434, pp. 122–128, 2014. View at: Publisher Site | Google Scholar
  16. P. Jomo, C. Otieno, and P. Nyawere, “Ab initio high-pressure study of semiconductor-metal phase transition of the chalcogenide compound KPSe6,” Advances in Condensed Matter Physics, vol. 2020, Article ID 3407141, p. 9, 2020. View at: Publisher Site | Google Scholar
  17. 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 GdBa2Cu3O7-x,” Computational Condensed Matter, vol. 23, Article ID e00461, 2020. View at: Publisher Site | Google Scholar
  18. M. H. Sadd, Elasticity: Theory, Applications, and Numerics, Academic Press, Cambridge, MA, USA, 2009.
  19. S. H. Simon, The Oxford Solid State Basics, Oxford University, England, UK, 2013.
  20. P. Y. Yu and M. Cardona, Fundamentals of Semiconductors: Physics and Materials Properties, Springer, New York, NY, USA, 1996.
  21. N. K. Omboga and C. Otieno, “Structural and electronic properties of the iron pnictide EuFe2As2 from first principles,” Journal of Physics Communications, vol. 4, no. 2, Article ID 025007, 2019. View at: Publisher Site | Google Scholar
  22. M. S. Ali, M. A. Rahman, and M. Z. Rahaman, “A theoretical investigation of ThCr2Si2-type Pd-based superconductors XPd2Ge2 (X = Ca, Sr, La, Nd),” Physica C: Superconductivity and its Applications, vol. 561, pp. 35–44, 2019. View at: Publisher Site | Google Scholar
  23. P. L. Alireza, Y. T. C. Ko, J. Gillett et al., “Superconductivity up to 29 K in SrFe2As2 and BaFe2As2 at high pressures,” Journal of Physics: Condensed Matter, vol. 21, no. 1, Article ID 012208, 2009. View at: Publisher Site | Google Scholar
  24. S. Piskunov, E. Heifets, R. I. Eglitis, and G. Borstel, “Bulk properties and electronic structure of SrTiO3, BaTiO3, PbTiO3 perovskites: an ab initio HF/DFT study,” Computational Materials Science, vol. 29, no. 2, pp. 165–178, 2004. View at: Publisher Site | Google Scholar
  25. F. Parvin and S. H. Naqib, “Structural, elastic, electronic, thermodynamic, and optical properties of layered BaPd2As2 pnictide superconductor: a first principles investigation,” Journal of Alloys and Compounds, vol. 780, pp. 452–460, 2019. View at: Publisher Site | Google Scholar
  26. 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
  27. P. Misra, Physics of Condensed Matter, Academic Press, Cambridge, MA, USA, 2011.
  28. G. D. Mahan, Condensed Matter in a Nutshell, Princeton University Press, New Jersey, NY, USA, 2011.
  29. 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
  30. K. Karch, P. Pavone, W. Windl, O. Schütt, and D. Strauch, “Ab initiocalculation of structural and lattice-dynamical properties of silicon carbide,” Physical Review B, vol. 50, no. 23, pp. 17054–17063, 1994. View at: Publisher Site | Google Scholar
  31. J. Birman, Theory of Crystal Space Groups and Lattice Dynamics, Springer, Berlin, Germany, 1984.
  32. 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
  33. F. Giustino, Materials Modelling Using Density Functional Theory: Properties and Predictions, Oxford University Press, England, UK, 2014.
  34. T. L. Hill, An Introduction to Statistical Thermodynamics, Courier Corporation, Chicago, IL, USA, 1986.
  35. I.-M. Low, Advances in Science and Technology of Mn+ 1AXn Phases, Elsevier, Cambridge, MA, USA, 2012.
  36. 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 © 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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.