Next Article in Journal / Special Issue
Elastic, Mechanical and Phonon Behavior of Wurtzite Cadmium Sulfide under Pressure
Previous Article in Journal
Synthesis, Crystal Structure and Hirshfeld Topology Analysis of Polymeric Silver(I) Complex with s-Triazine-Type Ligand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Analysis of Elastic, Mechanical and Phonon Properties of Wurtzite Zinc Sulfide under Pressure

Department of Physics, Hitit University, 19030 Corum, Turkey
*
Author to whom correspondence should be addressed.
Crystals 2017, 7(6), 161; https://doi.org/10.3390/cryst7060161
Submission received: 18 April 2017 / Revised: 23 May 2017 / Accepted: 31 May 2017 / Published: 4 June 2017
(This article belongs to the Special Issue Crystal Chemistry of Zinc, Cadmium and Mercury)

Abstract

:
We report for the first time the application of a mixed-type interatomic potential to determine the high-pressure elastic, mechanical, and phonon properties of wurtzite zinc sulfide (WZ-ZnS) with geometry optimization calculations under pressures up to 12 GPa. Pressure dependency of typical elastic constants, bulk, shear, and Young moduli, both longitudinal and shear wave elastic wave velocities, stability, as well as phonon dispersions and corresponding phonon density of states of WZ-ZnS were surveyed. Our results for the ground state elastic and mechanical quantities of WZ-ZnS are about experiments and better than those of some published theoretical data. Obtained phonon-related results are also satisfactory when compared with experiments and other theoretical findings.

1. Introduction

The computational modeling of materials has been a successful and rapid tool to address the unclear points of physical interest. In addition, predicting good elastic and thermodynamic properties of materials is a necessary demand for present-day solid state science and industry. In particular, these properties at high pressure and temperature are important for the development of modern technologies [1,2,3,4,5].
As emphasized in the authors’ recent work, wide-band-gap II–VI semiconductors such as ZnX (X: S, O, Se, Te) and CdX (X: S, Se, Te) are remarkable materials for the design of high-performance opto-electronic devices including light-emitting diodes and laser diodes in the blue or ultraviolet region [6].
According to its crystallographic structure, ZnS compounds can crystallize in either zinc blende (ZB) crystal structure with space group F-43m or wurtzite (WZ) crystal structure with the P63mc space group at ambient conditions [7,8,9,10].
To date, there have been a number of studies performed especially for the ground state (T = 0 K and P = 0 GPa) structural, elastic, thermodynamic, and optical properties of ZB and WZ crystal structures of ZnS compounds—not only with experiments, but also with theoretical works because of their technological importance. In 2004, Wright and Gale [11] reported new interatomic potentials to model the structures and stabilities of the ZB and WZ ZnS phases in which their theoretical results are within experiments. Later, in 2008, Bilge et al. [8], performed ab initio calculations based on projector augmented wave pseudo potential (PAW). They employed generalized gradient approximation (GGA) of density functional theory (DFT) to investigate the ground state mechanical and elastic properties of ZB, rock salt (NaCl), and WZ phases of ZnS. They concluded that the mechanical properties of ZnS under high pressure are quite different from those at ambient conditions. At the same time, Rong et al. [10] reported the pressure dependence of the elastic properties of ZB and WZ crystals of ZnS by the GGA within the plane-wave pseudopotential (PWP) of DFT and found reasonable results which are consistent with former experimental and theoretical data. Further, in 2009, Cheng et al. [12] conducted an experimental and theoretical study on first and second-order Raman scattering of ZB and WZ ZnS and found satisfactory experimental and theoretical results. Afterwards, in 2012, Grünwald et al. [9] established transferable pair potentials for ZnS crystals to accurately describe the ground state properties of ZB and WZ phases of ZnS compounds. Recently, Yu et al. [13] performed DFT calculations by using both the local density approximation (LDA) and GGA for the exchange-correlation potential and calculated the phonon dispersion curves and the phonon density of states of WZ ZnS in which the calculated values display good agreement with earlier experimental and theoretical data. In another recent study, Ferahtia et al. [14] published the first-principles plane-wave-based pseudopotential method calculations of the structural, elastic, and piezoelectric parameters of WZ ZnS and concluded a reasonable degree of agreement between their results and data available in the literature.
The above recent and increasing theoretical efforts on ZnS motivated us to perform this work by addressing the lesser-known high-pressure elastic, mechanical, and phonon properties of WZ ZnS with a different method. In contrast to the above applied methods and interatomic potentials of literature, this is the first report regarding the application of a mixed-type potential to determine the mentioned high-pressure quantities of WZ ZnS with geometry optimization calculations. During calculations, we concentrated on the pressure behavior of five independent elastic constants, bulk, shear, and Young moduli, elastic wave velocities, mechanical stability conditions, and phonon properties of WZ ZnS under pressures between 0 GPa and 12 GPa at T = 0 K.
The following part of the article provides details of our theoretical calculations with the employed interatomic potential in Section 2. In Section 3, we compare our results with the available experimental results and other theoretical data. Finally, Section 4 presents the main findings of this work in the conclusions.

2. Details of Theoretical Calculations

The most significant feature of materials modeling is the choice of the proper interatomic potentials that sufficiently and accurately describe the physical properties of the concerned problem. Simple empirical potentials are such modelling tools for materials, since they can yield reasonable results. These potentials can also provide good explanations of the defect energies, surface energies, elastic properties, and mechanical aspects of oxides [15], fluorides [16] and many other materials [17,18]. Most of these potentials include Coulomb interactions, short-range pair interactions, and ionic polarization treated by the shell model of Dick and Overhauser [19]. The sum of the Coulomb terms, short-range interactions, and ionic polarization expresses the total energy for these potentials. If we presume that the electron cloud of an ion is simulated by a massless shell charge qs and the nucleus by a core of charge qc, then the total charge becomes q = qs + qc. In the shell model interatomic potential approach, a harmonic force with spring constant K couples the core and the shell of the ion. So, for modelling the short-range pair interactions acting between the shells, we can use a typical Buckingham potential, as presented in Equation (1):
V i j B u c k i n g h a m = A exp ( r i j ρ ) C r i j 6
The first part of Equation (1) indicates the Born–Mayer term, whereas the second part represents the Van der Waals energies. Further, in our work, we applied the original form of a mixed-type interatomic potential of Hamad et al. [20] for short-range interactions, which incorporates the Buckingham and Lennard–Jones 9–6 potentials form as in Equation (2):
V i j s h o r t = A exp ( r i j ρ ) + B r i j 9 C r i j 6
We also considered the semi-covalent nature of ZnS by using a three-body potential for S-Zn-S angel as in its original form [20] and expressed by Equation (3):
V i j k = ( 1 2 ) K T B ( θ i j k θ 0 ) 2
In Equation (3), θ 0 and K T B show the equilibrium constant angle between S-Zn-S and fitting constant of Hamad et al. [20], respectively. Lastly, sulphur anion polarizibility treated by the shell model of Reference [19] can be written as in Equation (4):
V i j c o r e s h e l l = ( 1 2 ) K r i j 2
where r i j describes the core–shell separation, and K indicates the spring constant. Although extra details of presently employed potential and its parameterization procedure can be found in Reference [20], Table 1 lists the original potential parameters of Reference [20] used in our calculations.
All theoretical calculations in this work were carried out with the General Utility Lattice Program (GULP) 4.2 molecular dynamics code [21,22]. This versatile code allows the concerned structures to be optimized at constant pressure (all internal and cell variables are included) or at constant volume (unit cell remains frozen). To avoid the constraints, constant pressure optimization was applied to the geometry of WZ ZnS cell with the Newton–Raphson method based on the Hessian matrix calculated from the second derivatives. The cell geometry of WZ ZnS was assigned as a = b = 3.811   Å , c = 6.234   Å , u = 0.375, and α = β = 90 ° and γ = 120 ° with space group P 6 3 m c . During the present geometry optimization calculations, the Hessian matrix was recursively updated with the BFGS [23,24,25,26] algorithm. After setting the necessities for the geometry optimization of WZ ZnS, we devised multiple runs at zero Kelvin temperature and checked the pressure ranges between 0 GPa and 12 GPa in steps of 3 GPa. Further, phonon and associated properties of WZ ZnS were also addressed after constant pressure geometry optimization calculations as a function of pressure within the quasiharmonic approximation under zero Kelvin temperature, as implemented in GULP code. It is possible to capture the phonon density of states (PDOS) and dispersions for a material after specifying a shrinking factor with GULP phonon computations. Further, phonons are described by calculating their values at points in reciprocal space within the first Brillouin zone of the given crystal. To achieve the Brillouin zone integration and obtain the PDOS, we have used a standard and reliable scheme developed by Monkhorts and Pack [27] with 8 × 8 × 8 k-point mesh.

3. Results and Discussion

Figure 1 indicates the density behavior of WZ ZnS under pressure. As is well known, the density of many materials exhibits clear increments under pressure [6,15,17,18,28]. This is also the case for WZ ZnS under pressure, as seen in Figure 1. The minimum and maximum values of the density of WZ ZnS was found to be 4.08 g/cm3 and 4.59 g/cm3 for 0 GPa and 12 GPa, respectively, at zero temperature.
The elastic constants of materials not only provide precise and essential information about the materials, but also explain many mechanical and physical properties. Once the elastic constants are determined, one may get a deeper insight into the stability of the concerned material [6,15,17,18,28]. These constants can be also helpful to predict the properties of materials (i.e., interatomic bonding, equation of state, and phonon spectra). They also link to several thermodynamic parameters, such as the specific heat, thermal expansion, Debye temperature, Grüniesen parameter, etc. However, in general, elastic constants derived from the total energy calculations correspond to single-crystal elastic properties. So, the Voigt–Reuss–Hill approximation is a confident scheme for polycrystalline materials [6,15,17,18,28]. To obtain the accurate values of elastic constants and other analyzed parameters of WZ ZnS, the Voigt–Reuss–Hill values were taken into account for this research. For wurtzite crystal structures, five independent elastic constants exist as: C 11 , C 12 , C 13 , C 33 ,   and   C 44 [28].
The plot in Figure 2 shows our results for C 11 , C 12 , C 13 , C 33   and   C 44 elastic constants of WZ ZnS under pressures between 0 GPa and 12 GPa. All obtained elastic constants of WZ ZnS increase with the applied pressure except C 33   and   C 44 . Beyond these increments, a closer examination of Figure 2 reveals that the magnitudes of the elastic constants are in the range of: C 33 > C 11 > C 12 > C 13 > C 44 . Both the range of elastic constants and slight decrement of C 44 constant under pressure mimic the DFT findings of Tan et al. [7]. However, our results for the ground state parameters of WZ ZnS are obviously better than those of Tan et al. [7] as well as Wright and Gale [11], and are much closer to the experimental results of Neumann et al. [29]. Table 2 gives a numerical comparison of present results for the elastic constants and other calculated quantities of WZ ZnS with some previously published experimental and theoretical data.
From a stability outlook, the proverbial Born mechanical stability condition for a hexagonal structure also holds for wurtzite crystal, and must satisfy [28]:
C 44 > 0 ,   C 11 > C 12   and   ( C 11 + 2 C 12 ) C 33 2 C 13 2 > 0
The present results for the obtained elastic constants of WZ ZnS satisfy the mechanical stability condition (Figure 2 and Table 2), and this result points out that WZ ZnS is mechanically stable at 0 K temperature and 0 GPa pressure.
The bulk modulus (B) is the only elastic constant of a material that conveys much information about the bonding strength. Moreover, it is a measure of the material’s resistance to external deformation, and occurs in many formulae describing diverse mechanical–physical characteristics. The shear modulus (G), however, portrays the resistance to shape change caused by a shearing force. In addition to B and G, Young’s modulus (E) is the resistance to uniaxial tensions. These three distinct moduli (B, G, and E) are other valuable parameters for identifying the mechanical properties of materials [6,15,17,18,28].
Figure 3 represents the pressure behavior of B , G ,   and   E of WZ ZnS for the entire pressure range. From the prevalent physical definition of bulk modulus ( B = Δ P / Δ V ) is expected an increment for B because of its direct proportion to applied pressure. Hence, bulk modulus of WZ ZnS exemplifies a straight increment in Figure 3. Contrary to sharp increment in B, E and G moduli have slight decrements under pressure again similar to the DFT findings of Tan et al. [7]. Apart from the pressure behavior of these three moduli, Table 2 also summarizes another numerical comparison for B with former experimental and theoretical data. The present result for the bulk modulus of WZ ZnS with 79.5 GPa approximately agrees with experimental values and is better than other DFT results, as seen in Table 2.
The adjectives “brittle” and “ductile” signify the two distinct mechanical characters of solids when they are exposed to stress. Ductile and brittle features of materials play a key role during the manufacturing of materials [6,15,17,18,28]. So, we also evaluated the ductile (brittle) nature of WZ ZnS under pressure. In general, brittle materials are not deformable or are less deformable before fracture, whereas ductile materials are very deformable before fracture. At this point, the Pugh ratio is a determinative limit for the ductile (brittle) behavior of materials, and has popular use in literature. If the B/G ratio is about 1.75 and higher, the material is accepted to be ductile; otherwise, the material becomes brittle [6,15,17,18,28]. After a careful evaluation, we determined B/G values changing from 2.38 to 4.95, with a monotonous increment between P = 0 GPa and P = 12 GPa at zero temperature, respectively, for WZ ZnS (as seen in Figure 4). So, this result manifests the ductile character of WZ ZnS in both ground state and under pressure. As another result, all values of the B/G are higher than 1.75 and increase with pressure, which suggests that pressure can improve ductility (Figure 4).
In solids, low temperature (T = 0 K in our case) acoustic modes can trigger the vibrational excitations. Depending on this fact, two typical (longitudinal and shear) elastic waves exist. The velocity V L represents the longitudinal elastic wave velocity and V S denotes shear wave velocity [6,15,17,18,28]. Figure 5 shows the behavior of V L and V S elastic wave velocities for WZ ZnS as a function of pressure. V L has a substantial increment compared to V S , which is the most observed circumstance for many materials.
Figure 6 shows the phonon dispersion of WZ ZnS along the chosen Γ–A path. In addition, Table 3 compares the numerical values of the zone–center (Γ points) phonon frequencies of this work with earlier experimental and theoretical data of WZ ZnS. As seen in Table 3, the agreement with experiment is quite good, especially for A1 (TO), E1 (TO) and A2 (TO), E2 (TO) phonon modes and competing with recent GGA and LDA DFT data of Reference [14].
Further, we have focused on the ground state phonon density of states (PDOS) of WZ ZnS to clarify the contribution of the both elements (Zn and S) to the phonon properties of the compound. Figure 7 demonstrates the partial PDOS of WZ ZnS under zero pressure and temperature. As is apparent in Figure 7, the contribution of the Zn element to acoustic phonon modes is higher than S, whereas the opposite case is valid for the S element due to its dominant contribution to the optical modes of the compound.
Figure 8 shows the phonon dispersion curves of WZ ZnS under different pressures. Each pressure value above 0 GPa shifts the dispersion curves of WZ ZnS up slightly to the higher frequency values, as is clear in Figure 8. This behavior strictly originates from the atoms of WZ ZnS which are getting closer to each other under pressure, since they sit in steeper potential wells. The effect of pressure also causes the same behavior in the total density of states curves of WZ ZnS as shown in Figure 9.
Overall, the obtained results display a fair agreement with the experiments—in particular in elastic constants, bulk modulus, and phonon properties of WZ ZnS. Finally, the presented results for all considered quantities of WZ ZnS through this research are not only consistent with experiments, but also better than those of some published theoretical data.

4. Conclusions

In summary, we applied a mixed-type interatomic potential for the first time in conjunction with geometry optimization calculations to predict the high-pressure elastic, mechanical, and phonon properties of WZ ZnS. As our results demonstrate, the application of a mixed-type interatomic potential that is originally used for predicting the surface energies of zinc blende-type ZnS crystal structure well capture the elastic, mechanical and phonon features of WZ ZnS under pressure. From a quantitative evaluation, the obtained results of this work approximate the former experimental values—in particular for the elastic constants, bulk modulus, and phonon modes, and better than those of some other published data thanks to the sensitivity of the applied potential. The results of the present work may especially be helpful to future studies regarding the high-pressure elastic, mechanical, and other related properties of WZ ZnS.

Author Contributions

Melek Güler and Emre Güler conceived and designed the theoretical calculations; Melek Güler performed the all calculations; Melek Güler analyzed the all obtained data of work; both Melek Güler and Emre Güler wrote the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bouhemadoua, A.; Allali, D.; Bin-Orman, S.; Muhammad Abud Al Safi, E.; Khenata, R.; Al-Douri, Y. Elastic and thermodynamic properties of the SiB2O4 (B=Mg, Zn and Cd) cubic spinels: An ab initio FP-LAPW study. Mater. Sci. Semicond. Process. 2015, 38, 192–202. [Google Scholar] [CrossRef]
  2. Boudrifa, O.; Bouhemadou, A.; Guechi, N.; Bin-Omran, S.; Al-Douri, Y.; Khenata, R. First-principles prediction of the structural, elastic, thermodynamic, electronic and optical properties of Li4Sr3Ge2N6 quaternary nitride. J. Alloy. Compd. 2015, 618, 84–94. [Google Scholar] [CrossRef]
  3. Boulechfar, R.; Meradji, H.; Chouahda, Z.; Ghemid, S.; Drablia, S.; Khenata, R. FP-LAPW investigation of the structural, electronic and thermodynamic properties of Al3Ta compound. Int. J. Mod. Phys. B 2015, 29, 1450244. [Google Scholar] [CrossRef]
  4. Murtaza, G.; Gupta, S.K.; Seddik, T.; Khenata, R.; Alahmed, Z.A.; Ahmed, R.; Khachai, H.; Jha, P.K.; Bin Omran, S. Structural, electronic, optical and thermodynamic properties of cubic REGa3 (RE = Sc or Lu) compounds: Ab initio study. J. Alloy. Compd. 2014, 597, 36–44. [Google Scholar] [CrossRef]
  5. Zeng, Z.; Christos, S.G.; Baskoutas, S.; Bester, G. Excitonic optical properties of wurtzite ZnS quantum dots under pressure. J. Chem. Phys. 2015, 142, 114305. [Google Scholar] [CrossRef] [PubMed]
  6. Güler, E.; Güler, M. A Theoretical Investigation of the Effect of Pressure on the Structural, Elastic and Mechanical Properties of ZnS Crystals. Braz. J. Phys. 2015, 45, 296–301. [Google Scholar] [CrossRef]
  7. Tan, J.J.; Li, Y.; Ji, G.F. High-Pressure Phase Transitions and Thermodynamic Behaviors of Cadmium Sulfide. Acta Phys. Pol. A 2011, 120, 501–506. [Google Scholar] [CrossRef]
  8. Bilge, M.; Kart, S.Ö.; Kart, H.H.; Çağın, T. Mechanical and electronical properties of ZnS under pressure. J. Ach. Mat. Man. Eng. 2008, 31, 29–34. [Google Scholar]
  9. Grünwald, M.; Zayak, A.; Neaton, J.B.; Geissler, P.L. Transferable pair potentials for CdS and ZnS crystals. J. Chem. Phys. 2012, 136, 234111. [Google Scholar] [CrossRef] [PubMed]
  10. Rong, C.X.; Hu, C.E.; Yeng, Z.Y.; Cai, L.C. First-Principles Calculations for Elastic Properties of ZnS under Pressure. Chin. Phys. Lett. 2008, 25, 1064. [Google Scholar] [CrossRef]
  11. Wright, K.; Gale, J.D. Interatomic potentials for the simulation of the zinc-blende and wurtzite forms of ZnS and CdS: Bulk structure, properties, and phase stability. Phys. Rev. B 2004, 70, 035211. [Google Scholar] [CrossRef]
  12. Cheng, Y.C.; Jin, C.Q.; Gao, F.; Wu, X.L.; Zhong, W.; Li, S.H.; Chu, P.K. Raman scattering study of zinc blende and wurtzite ZnS. J. App. Phys. 2009, 106, 123505. [Google Scholar] [CrossRef]
  13. Yu, Y.; Fang, D.; Zhao, G.D.; Zheng, X.L. Ab initio Calculation of the thermodynamic properties of of wurtzite ZnS: Performance of the LDA and GGA. Chalcogenide Lett. 2014, 11, 619–628. [Google Scholar]
  14. Ferahtia, S.; Saib, S.; Bouarissa, N.; Benyettou, S. Structural parameters, elastic properties and piezoelectric constants of wurtzite ZnS and ZnSe under pressure. Superlattices Microstruct. 2014, 67, 88–96. [Google Scholar] [CrossRef]
  15. Güler, M.; Güler, E. High pressure phase transition and elastic behavior of europium oxide. J. Optoelectron. Adv. Mater. 2014, 16, 1222–1227. [Google Scholar]
  16. Ayala, A.P. Atomistic simulations of the pressure-induced phase transitions in BaF2 crystals. J. Phys. Condens. Matter 2001, 13, 11741. [Google Scholar] [CrossRef]
  17. Güler, E.; Güler, M. Phase transition and elasticity of gallium arsenide under pressure. Mater. Res. 2014, 17, 1268–1272. [Google Scholar] [CrossRef]
  18. Güler, E.; Güler, M. Elastic and mechanical properties of hexagonal diamond under pressure. Appl. Phys. A 2015, 119, 721–726. [Google Scholar] [CrossRef]
  19. Dick, B.G.; Overhauser, A.W. Theory of the Dielectric Constants of Alkali Halide Crystals. Phys. Rev. 1958, 112, 90. [Google Scholar] [CrossRef]
  20. Hamad, S.; Cristol, S.; Catlow, C.R.A. Surface Structures and Crystal Morphology of ZnS: Computational Study. J. Phys. Chem. B 2012, 106, 11002–11008. [Google Scholar] [CrossRef]
  21. Gale, J.D. GULP: A computer program for the symmetry-adapted simulation of solids. J. Chem. Soc. Faraday Trans. 1997, 93, 629. [Google Scholar] [CrossRef]
  22. Gale, J.D.; Rohl, A.L. The General Utility Lattice Program (GULP). Mol. Simul. 2003, 29, 291. [Google Scholar] [CrossRef]
  23. Broyden, G.C. The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations. J. Inst. Math. Appl. 1970, 6, 76. [Google Scholar] [CrossRef]
  24. Fletcher, R. A new approach to variable metric algorithms. Comput. J. 1970, 13, 317. [Google Scholar] [CrossRef]
  25. Goldfarb, D. A family of variable-metric methods derived by variational means. Math. Comput. 1970, 24, 23. [Google Scholar] [CrossRef]
  26. Shanno, D.F. Conditioning of quasi-Newton methods for function minimization. Math. Comput. 1970, 24, 647. [Google Scholar] [CrossRef]
  27. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188. [Google Scholar] [CrossRef]
  28. Güler, E.; Güler, M. High pressure elastic properties of wurtzite aluminum nitrate. Chin. J. Phys. 2014, 52, 1625–1635. [Google Scholar] [CrossRef]
  29. Neumann, H. Lattice dynamics and related properties of AIBIII and AIIBIV compounds, I. Elastic constants. Cryst. Res. Technol. 2004, 39, 939–958. [Google Scholar] [CrossRef]
  30. Hu, E.C.; Sun, L.L.; Zeng, Z.Y.; Chen, X.R. Pressure and Temperature Induced Phase Transition of ZnS from First-Principles Calculations. Chin. Phys. Lett. 2008, 25, 675–678. [Google Scholar]
  31. Brafman, O.; Mitra, S.S. Raman Effect in Wurtzite- and Zinc-Blende-Type ZnS Single Crystals. Phys. Rev. 1968, 171, 931. [Google Scholar] [CrossRef]
Figure 1. Density behavior of wurtzite (WZ) ZnS under pressure.
Figure 1. Density behavior of wurtzite (WZ) ZnS under pressure.
Crystals 07 00161 g001
Figure 2. Elastic constants of WZ ZnS under pressure.
Figure 2. Elastic constants of WZ ZnS under pressure.
Crystals 07 00161 g002
Figure 3. Elastic moduli (B, G, and E) behavior of WZ ZnS versus pressure.
Figure 3. Elastic moduli (B, G, and E) behavior of WZ ZnS versus pressure.
Crystals 07 00161 g003
Figure 4. B/G ratio of WZ ZnS against pressure.
Figure 4. B/G ratio of WZ ZnS against pressure.
Crystals 07 00161 g004
Figure 5. Behavior of longitudinal (VL) and shear wave (VS) velocities of WZ ZnS under pressure.
Figure 5. Behavior of longitudinal (VL) and shear wave (VS) velocities of WZ ZnS under pressure.
Crystals 07 00161 g005
Figure 6. Phonon dispersion curve of WZ ZnS at zero temperature along the chosen Γ–Α path.
Figure 6. Phonon dispersion curve of WZ ZnS at zero temperature along the chosen Γ–Α path.
Crystals 07 00161 g006
Figure 7. Partial and total phonon density of states of (PDOS) curve of WZ ZnS under zero temperature and pressure.
Figure 7. Partial and total phonon density of states of (PDOS) curve of WZ ZnS under zero temperature and pressure.
Crystals 07 00161 g007
Figure 8. Phonon dispersion curves of WZ ZnS at zero temperature along the chosen Γ–Α path under different pressures.
Figure 8. Phonon dispersion curves of WZ ZnS at zero temperature along the chosen Γ–Α path under different pressures.
Crystals 07 00161 g008
Figure 9. Total phonon density of states of (PDOS) curves of WZ ZnS under zero temperature with pressures 3 GPa, 6 GPa, 9 GPa, and 12 GPa.
Figure 9. Total phonon density of states of (PDOS) curves of WZ ZnS under zero temperature with pressures 3 GPa, 6 GPa, 9 GPa, and 12 GPa.
Crystals 07 00161 g009
Table 1. Mixed type theoretical interatomic potential used in this work. Parameters of the applied potential were taken from its original reference [20].
Table 1. Mixed type theoretical interatomic potential used in this work. Parameters of the applied potential were taken from its original reference [20].
Mixed Potential Parameters from Ref. [20]
General potentialA (eV)ρ (Å)B (eV. Å9)C (eV. Å6)
Zn-S213.200.475664.3510.54
S-S11,4130.1530.0129.18
Three-body potentialθ0 (degree)KTB (ev.rad−2)
Zn core-S Shell-S Shell109. 470.778
Spring potentialK (eV. Å−2)
S core-S Shell27.690
Ion chargesCharge (e)
Zn core2.000
S core1.357
S Shell−3.357
Table 2. Comparing our results with previous experimental and theoretical data for the considered parameters of WZ ZnS under zero pressure and temperature. * Bulk modulus experimental value was taken from Reference [30].
Table 2. Comparing our results with previous experimental and theoretical data for the considered parameters of WZ ZnS under zero pressure and temperature. * Bulk modulus experimental value was taken from Reference [30].
ParameterExp [29]PresentOther Theoretical
Ref [8]Ref [14]
C11 (GPa)124.2124.4115.6135.4
C12 (GPa)60.159.848.965.8
C13 (GPa)45.558.937.151.6
C33 (GPa)140.0113.0132.5160.7
C44 (GPa)28.637.327.232.4
B (GPa)75.8 *79.568.584.9
E (GPa) 75.3
G (GPa) 33.4
VS (km/s) 2.86
VL (km/s) 5.51
Table 3. Comparing our results with previous experimental and theoretical data for the phonon frequencies (cm−1) of WZ ZnS under zero pressure and temperature.
Table 3. Comparing our results with previous experimental and theoretical data for the phonon frequencies (cm−1) of WZ ZnS under zero pressure and temperature.
SymmetryExp [31] (cm−1)Present (cm−1)Ref [13]
GGA (cm−1)LDA (cm−1)
E2 (low)72986869
B1 (low) 230186199
A1 (TO)273279257290
E1 (TO)273253257293
E2 (high)286260261297
B1 (high) 279318347
A1 (LO)351345327350
E1 (LO)351352324349

Share and Cite

MDPI and ACS Style

Güler, M.; Güler, E. Theoretical Analysis of Elastic, Mechanical and Phonon Properties of Wurtzite Zinc Sulfide under Pressure. Crystals 2017, 7, 161. https://doi.org/10.3390/cryst7060161

AMA Style

Güler M, Güler E. Theoretical Analysis of Elastic, Mechanical and Phonon Properties of Wurtzite Zinc Sulfide under Pressure. Crystals. 2017; 7(6):161. https://doi.org/10.3390/cryst7060161

Chicago/Turabian Style

Güler, Melek, and Emre Güler. 2017. "Theoretical Analysis of Elastic, Mechanical and Phonon Properties of Wurtzite Zinc Sulfide under Pressure" Crystals 7, no. 6: 161. https://doi.org/10.3390/cryst7060161

APA Style

Güler, M., & Güler, E. (2017). Theoretical Analysis of Elastic, Mechanical and Phonon Properties of Wurtzite Zinc Sulfide under Pressure. Crystals, 7(6), 161. https://doi.org/10.3390/cryst7060161

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop