# Ab Initio Study of the Crystalline Structure of HgS under Low and High Pressure

Ahmed Amine Aidouni, Abdelkader Aissat, Mounir Ould-Mohamed,

Mohamed El Amine Benamar, Samuel Dupont, Jean Pierre Vilcot

► To cite this version:

Ahmed Amine Aidouni, Abdelkader Aissat, Mounir Ould-Mohamed, Mohamed El Amine Benamar, Samuel Dupont, et al. Ab Initio Study of the Crystalline Structure of HgS under Low and High Pressure. Crystals, 2024, 14 (9), pp.780. 10.3390/cryst14090780. hal-04722216

HAL Id: hal-04722216

https://hal.science/hal-04722216v1

Submitted on 4 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Ab Initio Study of the Crystalline Structure of HgS under Low and High Pressure

Ahmed Amine Aidouni <fnref n="1" />,<fnref n="2" />, Abdelkader Aissat <fnref n="1" />,<fnref n="3" />,* ID, Mounir Ould-Mohamed <fnref n="2" />, Mohamed El Amine Benamar <fnref n="1" />, Samuel Dupont <fnref n="3" /> ID and Jean Pierre Vilcot <fnref n="3" /> ID

**Abstract:** This study analyzes the lattice dynamics of HgS under various pressures using ab initio self-consistent calculations based on the plane-wave method (PW) and generalized gradient approximation (GGA). The static study, performed by enthalpy calculations, predicts that the transition from the cinnabar phase (α-HgS) to the zinc-blende $B_3$ (β-HgS) or wurtzite (2H) structures occurs at very low pressures, at 0.65 or 0.70 GPa, respectively. Furthermore, the transition from β-HgS to the rocksalt ($B_1$) phase occurs at 7 GPa, and at high pressure, specifically at 110 GPa, HgS can adopt the CsCl ($B_2$) phase. The mechanical study confirms the stability of the β and 2H phases at 0 GPa. Phonon calculations corroborate the results of the static and mechanical studies regarding stability ($α \xrightarrow{0.7\text{GPa}} 2H \xrightarrow{0.9\text{GPa}} β$), and the results indicate that the instabilities of the transverse acoustic (TA) modes, induced by the application of pressures of 10.5 GPa, 21 GPa, and 190 GPa, are responsible for the observed phase transitions in part of the Brillouin.

**Keywords:** cinnabar; phonon dispersion curves; high pressure; structural properties

## 1. Introduction

Binary Hg-X systems (where X = S, Se, Te) have garnered significant interest in electronic and optical technologies [1–5]. While HgS typically crystallizes in the trigonal cinnabar structure (α-HgS), it may also be found in the β phase (metacinnabar), either naturally or in the laboratory [6,7]. However, the existence of a hexagonal close-packed (HCP) modification of HgS is still debated, and there is insufficient theoretical work on this topic in the existing literature. Experimental investigations have discussed the presence or absence of the hexagonal phase at high or low pressures, but with mixed results [8,9]. In this article, we present new theoretical findings on the structural properties of HgS, specifically focusing on the potential existence of an HCP modification. Our study provides insights into the nature of HgS and could have significant implications for its electronic and optical applications.

At high pressures up to 20 GPa, no sign of phase transition is mentioned in studies carried out using X-rays and Raman spectroscopy [10]. In parallel, Bridgeman failed in the compression of cubic HgS (β-HgS), which transformed into the cinnabar form at only 5 GPa [11]. However, using energy dispersive X-ray diffraction, Huang et al. were able to determine that in a pressure range between 13 and 30 GPa, the cinnabar structure of HgS changes to the rocksalt structure ($B_1$) [12]. Nelmes et al. found 20.57 GPa to be the transition pressure [13].

<figure><img src="image_2.png" /></figure>

### check for updates

**Citation:** Aidouni, A.A.; Aissat, A.; Ould-Mohamed, M.; Benamar, M.E.A.; Dupont, S.; Vilcot, J.P. Ab Initio Study of the Crystalline Structure of HgS under Low and High Pressure. *Crystals* **2024**, *14*, 780. https://doi.org/10.3390/cryst14090780

Academic Editor: Thomas M. Klapötke

Received: 12 July 2024

Revised: 28 July 2024

Accepted: 29 August 2024

Published: 31 August 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Department of Material Science, Faculty of Material Science, Mathematics and Computer Science, University of Ahmed Draya, Adrar 01000, Algeria; aa_aidouni@hotmail.com (A.A.A.); benamardz64@gmail.com (M.E.A.B.)</footnote>

<footnote n="2">² LPTHIRM-Physics Department, Faculty of Sciences, University of Blida 1, P.O. Box 270, Blida 09000, Algeria; ouldmohamedmounir79@gmail.com</footnote>

<footnote n="3">³ Institute of Electronics, Microelectronics and Nanotechnology (IEMN)—UMR CNRS 8520, University of Sciences and Technologies of Lille 1, Avenue Poincare, BP 60069, 59652 Villeneuve of Ascq, France; jean-pierre.vilcot@univ-lille.fr (J.P.V.)</footnote>

<footnote>* Correspondence: sakre23@yahoo.fr or aissat_abdelkader@uni-blida.dz; Tel.: +00213-559604539</footnote>
Theoretically, all calculations based on density functional theory (DFT) (the majority of which are based on static calculations of Gibbs free enthalpy, which do not take into account the vibrations of the atoms in the crystal (dynamic study of phonons)) show the appearance of the rocksalt phase at different pressures: 20 GPa [14], 19.9 GPa [15], 16.1 GPa [16], and 26.75 GPa [17]. Using first-principles calculations, Yang et al. [18] determined that at very low pressure (3.4 GPa), the transition of the β-HgS ($B_3$) phase (space group $F\bar{4}3m$) to cinnabar α-HgS (space group $P3_121$) occurs, followed by a transition to the NaCl ($B_1$) structure (space group $Fm\bar{3}m$) at a pressure of approximately 14.3 GPa. Recently, Cihan et al. [14] showed that a phase transition can occur from the NaCl phase to the CsCl phase ($B_2$) (space group $Pm\bar{3}m$) at a pressure of 28 GPa, using calculations based on DFT (density functional theory), with the prediction of hydrostatic pressure conducted using enthalpy calculations.

In summary, the experimental study of the structural properties of HgS is not really rich, particularly in terms of the presence or absence of the hexagonal phase at low pressure, as well as the phase transitions at high pressure. First-principle studies on the structural properties of HgS are primarily based on static calculations (the Murnaghan equation of state EOS and Gibbs free energy), not accounting for phonon vibrations. Recent work has addressed HgS's dynamic properties, providing vibrational frequency data [19]. We sought the most stable phase among candidate structures (cinnabar α-HgS, wurtzite (2H), and β-HgS zinc blende) and determined transition pressures using enthalpy calculations under high pressure. To ensure structural stability, we conducted a dynamic study, confirming static results. This method, used in Zunger's work [20], helps identify new phases at various pressures. This work focuses on exploring the structural and dynamic properties of HgS, with an emphasis on its hexagonal phase at low pressure. Understanding these characteristics is crucial for optoelectronic applications, due to the unique electronic and optical properties of HgS. This material is promising for high-sensitivity sensors and detectors, owing to its phase transitions. The knowledge gained could also foster technological innovations in various fields, such as energy and communication technologies.

## 2. Materials and Methods

In this study, we used the linear response of the pseudopotentials method (PWscf) implemented in the Quantum Espresso code [21]. The exchange and correlation effects were treated using the generalized gradient approximation (GGA) of the PBE (Perdew-Burke-Ernzerhof) type [22]. The dynamical matrix is obtained by modifying the electron density through the inverse of the dielectric matrix, which describes the response of the valence electron density to a periodic lattice perturbation. The dielectric matrix is calculated from the eigenfunctions and energy levels of the unperturbed system [21].

For the dynamical properties, we used ultra-soft Rappe-Rabe-Kaxiras-Joannopoulos (rrkjus) pseudopotentials [23], developed for the PWscf code [21]. The electron wave functions were expressed using a plane-wave basis set with a cutoff for kinetic energy that varies between 70 and 100 Ry and a cutoff energy of charge density worth 10 times the ecutwfc. The standard method of Monkhorst and Pack (MP) [24] special k-points was employed to perform accurate Brillouin-zone integrations. We used 8 × 8 × 8 MP meshes for the zinc-blende ($B_3$) phase, rocksalt phase ($B_1$) and the CsCl phase ($B_2$); 9 × 9 × 4 MP meshes for the cinnabar phase; and 8 × 8 × 6 MP meshes for the wurtzite phase.

The corresponding integrating points covering the irreducible Brillouin zone were 29 k points for the $B_3$ and $B_1$ phase, 35 k-points for the $B_2$ phase, 163 k-points for the cinnabar phase, and 60 k-points for the $B_1$ phase. The self-consistent electronic iterations were kept until the calculated total energy of the crystal converged with a tolerance less than $10^{-7}$ Ryd. The lattice dynamics properties were calculated using the density functional perturbation theory (DFTP) [25] by means of a 4 × 4 × 4 q-point MP mesh for cubic systems and 4 × 4 × 2 q-point MP mesh for noncubic systems. A Fourier interpolation of these matrices allowed for obtaining the phonon dispersion curves. The convergence of the
total energy in relation to the k-points and basis was tested and found to be adequate to ensure convergence.

**3. Results and Discussion**

**3.1. *Structural Properties***

The initial calculations focused on the structural parameters of HgS, namely the equilibrium lattice parameter, bulk modulus, and its pressure derivative. These structural properties are studied in the following different structures: the cinnabar α-HgS structure, which corresponds to space group *P*3₁2₁; the zinc-blende *B*₃ structure corresponding to space group *F*43*m* (No.: 216); the wurtzite (2H) phase with space group *P*6₃*m*c; the rocksalt (*B*₁) phase corresponding to space group *F*m3*m*; and finally the CsCl (*B*₂) phase in space group *P*m3*m*. Figure 1 illustrates the total energy (Ryd) as a function of the volume (Bohr) E(V) ofHgS, using GGA for the cinnabar, zinc-blende, wurtzite, rocksalt, and CsCl phases. The calculated total energies are adjusted using Murnaghan's equation of state to obtain the structural parameters.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Total energy (Ryd) as a function of volume (au³) for HgS.</figcaption></figure>

In Table 1 are listed the calculated structural parameters for the ground-state phase, namely the equilibrium lattice parameter (a), the bulk modulus (B) and its pressure derivative (B₀), along with available theoretical and experimental values [14,15,19,26–28]. According to Table 1, our calculated values align well with the available experimental and theoretical values. The E(V) curves in Figure 1 clearly show that, under ambient conditions, the cinnabar structure is the fundamental phase of HgS. It is also evident that three structures (cinnabar, B₃, and 2H) exhibit almost equal minimal energies, with a difference of a few meV. Although the hcp structure of HgS is not observed in nature, its study is justified by several aspects, including the structural analogy between sphalerite and metacinnabar, as well as the reference to a natural form of HgS (hypercinnabar) with hexagonal symmetry [7]. In high-pressure conditions, the NaCl structure would be favored, which is consistent with experimental findings [12,13]. With greater compression, a transition to the CsCl structure can occur, which aligns well with theoretical calculations [14,15].

The structural phase transition occurs at the pressure where the Gibbs free energies of the two structures are equal. Given that all our calculations are conducted at *T* = 0 K, the Gibbs free energy *G* = *E*₀ + *PV* − *TS* simplifies to the enthalpy *H* = *E*₀ + *PV*. Therefore, the pressure at which the phase transition takes place is identified when the enthalpies of the two phases are the same. Figure 2 illustrates the enthalpy *H* as a function of pressure *p*, indicating that the static phase transition from the cinnabar phase to the zinc-blende (or wurtzite) phase of HgS occurs at 0.65 (0.70) GPa.
<table><caption>**Table 1.** The calculated equilibrium lattice parameter (a), bulk modulus ($B_0$), and the pressure derivative of the bulk modulus ($B'$) for HgS are compared with various experimental and theoretical values.</caption><thead><tr><th>Structure</th><th>Reference</th><th>a (Å)</th><th>C (Å)</th><th>B (GPa)</th><th>B<sub>0</sub> (GPa)</th></tr></thead><tbody><tr><td rowspan="4">α-HgS</td><td>This work</td><td>4.28</td><td>9.71</td><td>63.96</td><td>3.73</td></tr><tr><td>Exp [28]</td><td>4.15</td><td>9.50</td><td>-</td><td>4.00</td></tr><tr><td>Cal [19]</td><td>4.36</td><td>9.70</td><td>16.00</td><td>5.40</td></tr><tr><td>Cal [19]</td><td>4.29</td><td>9.69</td><td>22.50</td><td>4.80</td></tr><tr><td rowspan="4">β-HgS</td><td>This work</td><td>6.00</td><td></td><td>70.01</td><td>4.91</td></tr><tr><td>Exp [27]</td><td>5.85</td><td></td><td>68.60</td><td></td></tr><tr><td>Cal [15]</td><td>5.97</td><td></td><td>65.40</td><td></td></tr><tr><td>Cal [14]</td><td>6.04</td><td></td><td>59.08</td><td>4.69</td></tr><tr><td>Wurtzite</td><td>This work</td><td>4.22</td><td>6.91</td><td>50.50</td><td>4.93</td></tr><tr><td rowspan="3">B<sub>1</sub></td><td>This work</td><td>5.52</td><td></td><td>92.30</td><td>3.58</td></tr><tr><td>Cal [14]</td><td>5.59</td><td></td><td>98.90</td><td>4.65</td></tr><tr><td>Cal [15]</td><td>5.49</td><td></td><td>86.30</td><td></td></tr><tr><td rowspan="3">B<sub>2</sub></td><td>This work</td><td>3.45</td><td></td><td>53.80</td><td>5.00</td></tr><tr><td>Cal [14]</td><td>3.43</td><td></td><td>94.31</td><td>5.44</td></tr><tr><td>Cal [15]</td><td>3.52</td><td></td><td>67.40</td><td>4.80</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Enthalpy as a function of pressure (P) for HgS.</figcaption></figure>

It turns out that the α-HgS cinnabar phase has the same enthalpy as the β-HgS zinc-blende phase and the hexagonal phase (within a few meV), which explains the possibility of obtaining both phases in the laboratory (the cinnabar phase and one of the other two phases), although this requires the addition of a very small amount of iron [29] to promote the growth of the zinc-blende phase. However, the zinc-blende (wurtzite) phase is calculated to be stable at zero pressure with a very slight enthalpy difference compared to the cinnabar phase. Therefore, this cubic phase must be stable from 0.7 GPa until the appearance of another cubic phase—the rocksalt phase. Although there are many experimental measurements of these transitions (α-HgS-β-HgS) at very low pressures (see Introduction), no theoretical calculation seems to have been reported, except for that of Cordona et al. [15], who found the same transition pressure value of 0.7 GPa, although they did not study the hexagonal structure.

The values obtained by the Gibbs static method using the first-principles pseudopotential method are given in Table 2. At a high pressure, the transition goes to the rocksalt phase. The value obtained seems far from the values obtained dynamically in our phonon calcula-
tions (discussed in the phonon section) and other theoretical and experimental values. This difference is due to the fact that the $B_3$ to $B_1$ transition starts at 6 GPa but can only occur beyond 20 GPa. In mercury compounds, the series of phase transitions ($B_3$-cinnabar-to-$B_1$) suggests that the sp3 bond first transitions to the spbond, then to the NaCl structure, unlike other IB-VI compounds that transform directly from sp3 to the NaCl structure [12]. This transition might involve a change in the coordination or bonding environment from a more complex geometry (such as tetrahedral, found in sp3 hybridization) to a simpler, linear one before finally transitioning to the NaCl structure.

<table><caption>**Table 2.** Comparison of transition pressures (Pt) between Gibbs statistical method and dynamic method with experimental measurements and theoretical calculations.</caption><thead><tr><th>Structure</th><th>Reference</th><th>Pt (GPa)</th></tr></thead><tbody><tr><td rowspan="4">α-HgS</td><td>Stat</td><td>0</td></tr><tr><td>Dyn</td><td>0</td></tr><tr><td>Exp [28]</td><td>0</td></tr><tr><td>Cal [19]</td><td>0</td></tr><tr><td rowspan="4">β-HgS</td><td>Stat</td><td>0.65</td></tr><tr><td>Dyn</td><td>0</td></tr><tr><td>Exp [28]</td><td>0</td></tr><tr><td>Cal-LDA [19]</td><td>0.70</td></tr><tr><td rowspan="2">Wurtzite</td><td>Stat</td><td>0.70</td></tr><tr><td>Dyn</td><td>0</td></tr><tr><td rowspan="5">B<sub>1</sub></td><td>Stat</td><td>6.07</td></tr><tr><td>Dyn</td><td>21</td></tr><tr><td>Exp [12,13]</td><td>13, 20.57</td></tr><tr><td>Cal-GGA [14–17]</td><td>20, 19.9, 16.1, 26.75</td></tr><tr><td>Cal-LDA [19]</td><td>22.2</td></tr><tr><td rowspan="3">B<sub>2</sub></td><td>Dyn</td><td>190</td></tr><tr><td>Exp</td><td>-</td></tr><tr><td>Cal-GGA [14,30]</td><td>28, 210</td></tr></tbody></table>

HgS can adopt the CsCl ($B_2$) phase from 112 GPa, which does not agree with the theoretical value also obtained by a statistical enthalpy calculation [14]; however, it confirms the result obtained by [30], which predicts that HgS will transition to a simple cubic structure beyond 100 GPa. The Cmcm and Pmmn structures were subjected to calculations. This confirms the stability of the Cmcm phase compared to Pmmn. These results are in agreement with the conclusions of Nelmes and McMahon (52 GPa), who suggested “a distortion of the RS structure, possibly to Cmcm” [31]. Lastly, Biering predicted a transition from the Pmmn structure to the CsCl structure at a pressure of 210 GPa [30].

## 3.2. Elastic Properties

The mechanical properties of the binary HgS compound have been calculated for the cubic (zinc blende) and hexagonal (wurtzite) phases. For these structures, the elastic constants $C_{11}$, $C_{12}$, $C_{13}$, $C_{33}$ and $C_{44}$ were used to determine the shear modulus (G), Young's modulus (E), Poisson's ratio (ν) and Debye temperature. The mechanical stability of a cubic system necessitates that three independent elastic constants—$C_{11}$, $C_{12}$, and $C_{44}$—satisfy the Born conditions [32–36] as follows:

$$C_{11} > 0, C_{44} > 0, C_{11} - C_{12} > 0, C_{11} + 2C_{12} > 0, C_{12} < B < C_{11} \quad (1)$$

For a hexagonal structure, the stability criteria are

$$C_{11} > |c_{12}|, C_{44} > 0, (C_{11} + C_{12})C_{33} > 2C_{13}^2 \quad (2)$$
In the case of the B₃ structure, the elastic constants satisfy these conditions, which reveal stability at zero pressure. For the hexagonal structure, the elastic constants also meet these criteria, which show that it is a metastable phase at ambient pressure. The mechanical properties calculated for these two structures are reported in Table 3.

<table><caption>**Table 3.** Elastic constant in GPa, Shear modulus, G (GPa), Young's modulus E (GPa), Poisson's ratio v and Debye temperature θ<sub>D</sub> (K).</caption><thead><tr><th>Ref.</th><th>C<sub>11</sub></th><th>C<sub>12</sub></th><th>C<sub>13</sub></th><th>C<sub>33</sub></th><th>C<sub>44</sub></th><th>E</th><th>G</th><th>v</th><th>θ<sub>D</sub></th></tr></thead><tbody><tr><td>This work</td><td>82.50</td><td>71.10</td><td>-</td><td>-</td><td>18.20</td><td>32.70</td><td>11.40</td><td>0.42</td><td>138</td></tr><tr><td>[16]</td><td>70.40</td><td>46.50</td><td>-</td><td>-</td><td>43.30</td><td>76.90</td><td>30.30</td><td>0.26</td><td>-</td></tr><tr><td>[32]</td><td>61.80</td><td>44.90</td><td>-</td><td>-</td><td>25.20</td><td>44.10</td><td>-</td><td>0.33</td><td>213</td></tr><tr><td>This work</td><td>86.70</td><td>37.80</td><td>27.60</td><td>93.80</td><td>16.80</td><td>59.50</td><td>23.40</td><td>0.27</td><td>196</td></tr></tbody></table>

### 3.3. Dynamical Stability

In order to statically and mechanically study the stable phase discussed earlier, we conducted a dynamic study based on phonon calculations at different pressures. The phonon dispersions along the high-symmetry lines in different structural phases illustrate the various symmetries of the phonon branches. The point groups can have up to 12 irreducible representations, which are differentiated in the phonon dispersion plots by distinct colors. Figure 3 illustrates the phonon dispersion along the high-symmetry lines for the cinnabar phase at 0 GPa (a), the wurtzite phase at 0.7 GPa (b), and the zinc-blende phase at 0.9 GPa (c).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Dispersion of vibrational modes in the cinnabar, wurtzite and zinc-blende phases of HgS.</figcaption></figure>

It is evident that HgS is dynamically stable in the α phase, with a first transition to the hexagonal phase at 0.7 GPa. At 0.9 GPa, HgS stabilizes in the β phase. These structures differ only slightly in volume and energy, as the zinc blende can be considered as an ABCABC stacking and wurtzite as an ABAB stacking. These results are consistent with the experimental observations of hexagonal stacking in the works of KovácsKis and Pandey [8,9].

Figure 4 illustrates the phonon dispersion in the zinc-blende phase of HgS for various pressure values. As the volume decreases, the TO and LO phonon modes of the zinc-blende phase increase in frequency, while the TA mode decreases. At a volume of 0.909 V₀ (where V₀ indicates the calculated equilibrium volume of β - HgS), corresponding to a pressure of
10.2 GPa, the phonon modes along the [0e0] direction become imaginary, indicating the dynamic instability of the zinc-blende phase. To properly identify these instabilities and obtain accurate pressure values, the frequency dependence of the phonons TA for pressures between 4 and 10.50 GPa and for several points ($q_i$) of the direction [e00] were studied.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Phonons dispersion curve of the zinc-blende phase of HgS at different pressures.</figcaption></figure>

Figure 5 shows that all TA ($q_i$) modes decrease sharply with increasing pressure, with an imaginary TA phonon mode identified at a pressure of 5.7 GPa. This TA phonon mode begins to become imaginary only at point X at a pressure of around 10.2 GPa, indicating the start of the $B_3$ phase transition to phase $B_1$ of HgS. The instabilities of the TA modes at points $q_i$, particularly at point X, are accountable for the phase transition. The instabilities of the phonons predicted at a pressure of 10.2 GPa are consistent with the transition pressure of 13 GPa measured by Huang [12].

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Calculated TA phonons frequencies at the $q_i(e, 0, 0)$ points for the $B_3$ phase of HgS as a function of pressure.</figcaption></figure>

Additionally, Figure 6 demonstrates that HgS is dynamically stable in phase $B_1$ only from a pressure of 21 GPa, which aligns with the value obtained from experimental mea-
surements [13] and ab initio calculations [14]. At very high pressure, HgS is stable in the CsCl phase ($B_2$). By studying the TA modes of phase $B_1$ at various pressures, the instabilities of these modes at points qi, specifically at point X, are the cause of the phase transition at a pressure of 96 GPa. This value marks the beginning of the phase transition to another structure. Although theoretical simulations (static method) suggest that the transition from the NaCl structure to the CsCl structure could occur at a pressure of 28 GPa [14] and from $B_1$ $57.9\text{GPa}$ $\rightarrow$ $Pmmn$ $210\text{GPa}$ $\rightarrow$ $B_3$ [30], experimental studies have not yet examined HgS at very high temperatures and pressures above 60 GPa.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Phonon dispersion curve of the NaCl phase of HgS at different pressures.</figcaption></figure>

Figures 7 and 8 show that our dynamic calculations predicted a value around 109 GPa, which aligns well with [30]. However, the static method founded on enthalpy calculations [14] cannot confirm this stabilization since it does not account for atomic vibrations. Additionally, it is possible that another structure with lower total energy could exist but has not yet been investigated.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Calculated TA phonons frequencies at the $q_i(e, 0, 0)$ points for the $B_1$ phase of HgS as a function of pressure.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Phonon dispersion curve of the CsCl phase of HgS at different pressures.</figcaption></figure>

**4. Conclusions**

In this study, the static, mechanical, and dynamic properties of HgS were investigated using first-principles calculations. The results of the static study showed that HgS can crystallize in the cinnabar, wurtzite, or zinc-blende structures with a difference of a few meV. At a high pressure, the NaCl and CsCl structures are favored. The study of mechanical properties confirmed the stability of the cubic and hexagonal phases. The dynamic analysis, carried out through phonon calculations, confirmed the previous results on the existence of the hexagonal phase between the alpha and beta phases and indicated that a phase transition occurs towards the NaCl structure and that at very high pressure (>190 GPa), HgS will adopt the CsCl phase. Our results are in agreement with experimental measurements and theoretical data. Our study represents the first theoretical proposal of the existence of the hexagonal phase and the first dynamic analysis at high pressure.

**Author Contributions:** Conceptualization, A.A.A. and A.A.; methodology, A.A.A.; software, A.A.A.; validation, A.A.; formal analysis, A.A.A. and A.A.; investigation, A.A.A. and A.A.; resources, J.P.V.; data curation, A.A.A. and A.A.; writing—original draft preparation, A.A.; writing—review and editing, A.A. and J.P.V.; visualization, A.A., M.O.-M., M.E.A.B. and S.D.; supervision, A.A., M.O.-M., J.P.V. and S.D.; project administration, M.E.A.B.; funding acquisition, A.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Data Availability Statement:** The original contributions presented in the study are included in the article.

**Acknowledgments:** The authors would like to thank the optoelectronics team, Institute of Electronics, Microelectronics and Nanotechnology (IEMN), UMR CNRS 8520. University of Sciences and Technologies of Lille 1, Avenue Poincare, 60069, 59652 Villeneuve of Ascq, France, for the helpful support during the accomplishment of this work.

**Conflicts of Interest:** The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

**References**

1. Wissmann, H.; Anh, T.T.; Rogaschewski, S.; Von Ortenberg, M. Self-organized MBE growth of II-VI epilayer son patterned GaSb substrates. *J. Cryst. Growth* **1999**, *201*, 619–622. [CrossRef]

2. Hu, T.-J.; Cui, X.-Y.; Li, X.-F.; Wang, J.-S.; Yang, J.-H.; Gao, C.-X. In Situ Electrical Resistivity and Hall Effect Measurement of beta-Hg Sunder High Pressure. *Chin. Phys. Lett.* **2015**, *32*, 3. [CrossRef]
3. Ren, T.; Xu, S.; Zhao, W.-B.; Zhu, J.-J. A surfactant-assisted photochemical route to single crystalline HgS nanotubes. *J. Photochem. Photobiol. A Chem.* **2005**, *173*, 93–98. [CrossRef]

4. Kershaw, S.V.; Yiu, W.K.; Sergeev, A.; Rogach, A.L. Development of synthetic methods to grow long-wavelength infrared-emitting HgTe quantum dots in dimethyl form amide. *Chem. Mater.* **2020**, *32*, 3930–3943. [CrossRef]

5. Mews, A.; Kadavanich, A.V.; Banin, U.; Alivisatos, A.P. tructural and spectroscopic investigations of CdS/HgS/CdS quantum-dot quantum wells. *Phys. Rev. B* **1996**, *53*, R13242. [CrossRef]

6. Zallen, R.; Michael, E.S. Plasma edge and band structure of cubic HgS. *Solid State Commun.* **1970**, *8*, 1291–1294. [CrossRef]

7. Szuszkiewicz, W.; Witkowska, B.; Jouanne, M.; Balkanski, M. Raman spectroscopy of cubic Hg<sub>1-x</sub>F<sub>x</sub>S. *Mater. Sci. Forum.* **1995**, *182–184*, 711–714. [CrossRef]

8. Kovácskis, V.; Dódony, I. Structural disorder in natural cubic HgS. *Acta Mineral. Petrogr.* **1999**, *40*, 3–10.

9. Chatterjee, P.M.; Pandey, D.; Raval, R.; Dubey, A.K. Coherent Aspects of Multifaceted Eco-friendly Biopolymer-Polyglutamic Acid from the Microbes. *J. Pure Appl. Microbiol.* **2019**, *13*, 741–756. [CrossRef]

10. Werner, A.; Hochheimer, H.D.; Strössner, K.; Jayaraman, A. High-pressure X-ray diffraction studies on HgTe and HgS to 20 GPa. *Phys. Rev. B* **1983**, *28*, 3330. [CrossRef]

11. Bridgman, P.W. *The Compression of 46 Substances to 50,000 kg/cm²*; American Academy of Arts& Sciences: Boston, MA, USA, 1940; pp. 21–51.

12. Huang, T.; Ruoff, A.L. Pressure induced phase transition of HgS. *J. Appl. Phys.* **1983**, *54*, 5459–5461. [CrossRef]

13. Nelmes, R.J.; McMahon, M.I. Structural transitions in the group IV, III-V, and II-VI semiconductors under pressure. *Semicond. Semimet.* **1998**, *54*, 145–246.

14. Kürkçü, C.; Selgin, A.L.; Merdan, Z.; Yamçiçier, Ç.; Öztürk, H. Investigation of structural and electronic properties of β-HgS: Molecular dynamics simulations. *Chin. J. Phys.* **2018**, *56*, 783–792. [CrossRef]

15. Ullah, N.; Murtaza, G.; Khenata, R.; Wong, K.M. Phase transition, electronic and optical properties of mercury chalcogenides under pressure. *Phase Transit.* **2014**, *87*, 571–581. [CrossRef]

16. Saini, P.K.; Ahlawat, D.S.; Singh, D. Structural, electronic and optical properties of HgS under pressure using FP-LAPW method. *Int. Res. J. Eng. Technol. (IRJET)* **2017**, *4*, 1080–1086.

17. Vodolazov, D.Y.; Peeters, F.M. Superconducting rectifier based on the asymmetric surface barrier effect. *Phys. Rev. B* **2005**, *72*, 172508. [CrossRef]

18. Yang, X.C.; Yang, J.; Zhang, E.-J.; Gao, C.-X. Study of electronic and elastic properties of beta-HgS under high pressure viafirst-principles calculations. *Phys. Status Solidi C Curr. Top. Solid State Phys.* **2011**, *8*, 1703–1707.

19. Cardona, M.; Kremer, R.K.; Lauck, R.; Siegle, G.; Muñoz, A.; Romero, A.H. Electronic, vibrational, and thermodynamic properties of metacinnabar β-HgS, HgSe, and HgTe. *Phys. Rev.* **2009**, *B 80*, 195204. [CrossRef]

20. Ozolinş, V.; Zunger, A. The oryofsy stematicabsence of NaCl-Type(β-Sn-Type)high pressure phases in covalent(ionic) semiconductors. *Phys. Rev. Lett.* **1999**, *82*, 767. [CrossRef]

21. Baroni, S.; de Gironcoli, S.; Stefano, Corso, D. Phonons and related crystal properties from density functional perturbationtheory. *Rev. Mod. Phys.* **2001**, *73*, 515–562. [CrossRef]

22. Perdew, J.P.; Burke, K.; Ernzerhof, M. Matthias, Generalized gradient approximation made simple. *Phys. Rev. Lett.* **1996**, *77*, 3865. [CrossRef] [PubMed]

23. Andrew, M.R.; Karin, M.R.; Kaxiras, E.; Joannopoulos, J.D. Optimised pseudopotential. *Phys. Rev. B* **1990**, *41*, 1227, Erratum in *Phys. Rev.* **1991**, *44*, 13175–13176.

24. Pack, J.D.; Hendrik, M.J. Special points for Brillouin-zone integrations, a reply. *Phys. Rev. B* **1977**, *16*, 1748. [CrossRef]

25. Dominic, F.M. The compressibility of media under extreme pressures. *Proc. Natl. Acad. Sci. USA* **1944**, *30*, 244–247.

26. Madelung, O. *Landolt Bornste: Numerical Data and Functional Relationships in Science and Technology*; New Series, Group III; Springer: Berlin/Heidelberg, Germany, 1982; Volume 22, pp. 244–247.

27. Strobel, S.; Schleid, T. Quaternäre Caesium-Kupfer(I)-Lanthanoid (III)-Selenide Vom Typ CsCu<sub>3</sub>M<sub>2</sub>Se<sub>5</sub>(M=Sm, GdLu). *Z. Für Anorg. Und Allg. Chem.* **2014**, *630*, 706–711. [CrossRef]

28. Földvári, M. *Hand book of Thermogravimetric System of Minerals and Its Use in Geological Practice*; Geological Institute of Hungary: Budapest, Hungary, 2011.

29. Biering, S.; Schwerdtfeger, P.A. comparative density functional study of the high-pressure phases of solid ZnX, CdX, and HgX (X = S, Se, and Te): Trends and relativistic effects. *J. Chem. Phys.* **2012**, *137*, 034705. [CrossRef] [PubMed]

30. Nelmes, R.J.; McMahon, M.I. *High Pressure Semiconductor Physics I*; Suski, T., Paul, W., Eds.; Academic Press: Cambridge, MA, USA, 1998; pp. 145–247.

31. McMahon, M.I.; Nelmes, R.J. Observation of a cinnabar phase in GaAs at high pressure. *Phys. Rev. Lett.* **1997**, *19*, 3697. [CrossRef]

32. Born, M.; Huang, K. *Dynamical Theory of Crystal Lattices*; Oxford University Press: Oxford, UK, 1954.

33. Eshelby, J.D. The determination of the elastic field o fan ellipsoidal inclusion, and related problems. *Proc. R. Soc. A* **1957**, *241*, 376–396.

34. Reuss, A. Berechnung Der Fließgrenzen Mischkristallen Aufgrund Der Plastizität Bedingung Für Einkristalle. *ZAMM-J. Appl. Math. Mech./Zeitschrift für Angewandte Mathematik und Mechanik* **1929**, *9*, 49–58. [CrossRef]
35. Morse, P.M. The elastic constants of isotropic materials. *Phys. Rev.* **1929**, *34*, 57–64. [CrossRef]

36. Nayfeh, A.H. *Introduction to Perturbation Techniques*; John Wiley&Sons: Hoboken, NJ, USA, 1981.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.