Pressure-Driven Responses in Cd₂SiO₄ and Hg₂GeO₄ Minerals: A Comparative Study

1. Introduction

Thenardite is a mineral consisting of anhydrous sodium sulfate and is commonly found in dry evaporite environments. Gmelin [1] identifies eight distinct anhydrous phases of sodium sulfate. The phase commonly referred to as thenardite, named after the mineral, is designated as Na₂SO₄(V). It is documented to exhibit stability within the temperature range of 32°C to approximately 180°C. The thenardite mineral crystallizes in an orthorhombic structure with Fddd space group. The atomic arrangement in ternary compounds of the AB₂O₄ type is known to be influenced by the relative sizes of the A and B cations. Muller and Roy [2] demonstrated that the stability of the structure type relies on the comparative sizes of these cations. They found that within silicate and germanate compounds, the olivine structure type can accommodate a wider spectrum of octahedral cation sizes in contrast to other structure types. Despite earlier expectations from cation radius assessments indicating that chromous orthosilicate would adopt the olivine structure, Cr₂SiO₄ was observed to crystallize in the thenardite-type structure instead [3,4,5]. Typically, the latter is associated with cations that are larger than those present in silicates and germanates exhibiting the olivine structure type. Mercury orthogermanate, Hg₂GeO₄ [6], and cadmium orthosilicate, Cd₂SiO₄ [7], represent the sole other known compounds featuring this particular structure type.

Mehrotra et al. [5] conducted a comprehensive review of compounds exhibiting isomorphism with thenardite in 1978. Apart from Na₂SO₄ and Na₂SeO₄, the structural configuration resembling thenardite is also observed in Ag₂SO₄, various mixed phosphates and arsenates (such as AgHgPO₄, NaCdAsO₄, AgCdAsO₄). In such materials, non-tetrahedral elements (Na, Ag, Cd, Hg) are situated close to the centres of distorted six-cornered polyhedra, maintaining significant separation. Cr₂SiO₄ exhibits the Fddd space group, akin cell dimensions, and comparable oxygen atomic coordinates. The primary distinction lies in the z coordinate of Cr in contrast to that of Na, Ag, Hg, or Cd atoms. In Na₂SO₄, Na atoms within highly distorted octahedra are spaced 3.60 Å apart. Contrastingly, in the comparable observation of Cr₂SiO₄, the movement of O atoms, notably the displacement of the Cr atom along the c axis, results in (1) flattening one side of the octahedron to form a (distorted) equatorial plane, and (2) relocating the Cr atom from its position near the midpoint of the octahedron to the face of the equatorial plane. This motion simultaneously extends the distance to the oxygens in the tetrahedral edge while reducing the distance to the Cr atom in the neighboring equatorial plane.

The exploration of cadmium silicate, encompassing its three distinct stoichiometries: CdSiO₃ [8], Cd₂SiO₄, and Cd₃SiO₅ [7], has been ongoing since the early 20^th century, primarily due to its luminescent properties [9]. This luminescence, often referred to as persistent luminescence, has found applications in various fields, including bioimaging, optical data storage, as well as decorative and emergency signalling [10,11,12]. Notably, Dy³⁺ doped cadmium silicate materials have demonstrated persistent luminescence exceeding five hours [12]. Cadmium silicates adopt specific structures at room temperature which are: monoclinic CdSiO₃, orthorhombic Cd₂SiO₄, and tetragonal Cd₃SiO₅. The initial structure (CdSiO₃) is classified within the P121/c1 space group, characterized by three asymmetrically distorted (CdO₆)^4- octahedra arranged in slabs parallel to the (0 1 0) plane, with unbranched silicate chains serving as separators. These chains are composed of regular corner-sharing (SiO4)4- tetrahedra [8]. The second structure (Cd₂SiO₄) is classified within the Fddd space group and comprises solely one asymmetrically distorted (CdO₆)^4- octahedron amidst highly regular (SiO₄)^4- tetrahedra [7]. The final structure (Cd₃SiO₅) is categorized under the P4/nmm space group, featuring two distinctly symmetric (CdO₆)^4- octahedra positioned between regular (SiO₄)^4- tetrahedra [7,13]. Furthermore, the thenardite-type compound Hg₂GeO₄ was synthesized by Rfpke and Eysel (1978) through hydrothermal methods at 673 K and 2 x 10⁸ Pa in a welded gold capsule. Hg₂GeO₄ is the sole known compound within the HgO-SiO₂-GeO₂ system. Mehrotra et al. [5] proposed a tentative structure based on geometrical considerations, a proposal later generally confirmed by Hesse in 1981 [6]. The inclination of Hg²⁺ toward twofold coordination is apparent from the varying short and long Hg-O distances, leading to a markedly distorted structure within the thenardite family [14].

Cr₂SiO₄ and Cd₂SiO₄ possessing the orthorhombic thenardite structure have undergone examination under elevated pressures [15,16]. The Cd₂SiO₄ structure underwent thorough refinement under different pressures, reaching up to 9.5 GPa. This analysis unveiled a bulk modulus of 119.5±0.5 GPa, accompanied by a pressure derivative B₀’ valued at 6.17(4). On the contrary, the structure of Cr₂SiO₄, refined up to 9.2 GPa, demonstrated a bulk modulus of 94.7±0.5 GPa, alongside a B₀’ value of 8.32(14). Miletech et al. [16] attributed the relatively lower bulk modulus and higher B₀’ of the chromous structure to the compression of the unusually long Cr-O bond and the comparatively small size of the Cr²⁺ ion relative to the coordination polyhedron’s size. Meanwhile, there have been no reports on the structural properties of a similar compound, Hg₂GeO₄ of the thenardite-type, under high pressure to date.

In this study, our objective is to examine the distinct physical characteristics of Cd₂SiO₄ and Hg₂GeO₄ and investigate the effects of pressure variations on their behaviour. The subsequent section of our research focuses on the computational methodology employed. Transitioning to the third segment, we present a thorough overview of our major discoveries and engage in discussions regarding structural properties, elasticity, mechanical responses, dynamic features, and electronic properties. We will analyse these aspects under both ambient and pressured conditions. Lastly, we will synthesize our findings and present concluding remarks in the final section.

2. Materials and Methods

The molecular geometry optimization was carried out using the Projector-Augmented Wave (PAW) method within the density-functional theory (DFT) framework, as implemented in VASP [17]. Specifically, the Projector-Augmented Wave formalism-based pseudopotentials were employed. The Perdew-Burke-Ernzerhof for solids (PBEsol) [18] functional, within the Generalized-Gradient Approximation (GGA), was utilized for generating exchange-correlation functionals. Additionally, PBE [19] calculations were performed, revealing that PBEsol provides a better description of the crystal structure of Be₂SiO₄. Consequently, this study primarily focuses on the results obtained from PBEsol calculations. For all computations, a plane wave energy cutoff of 600 eV was set, with a selected energy convergence criterion of 10^-8 eV. Geometry optimization employed a dense k-mesh based on the Monkhorst-Pack technique [20]. Phonon dispersion calculations were conducted using density-functional perturbation theory with VASP and Phonopy [21]. The phonon calculations utilized a 2×2×2 supercell (112 atoms), employing a 6×6×6 k-point mesh and an energy convergence criterion of 10^-8 eV. To address the underestimation of the band gap for semiconductors/insulators in GGA approximation, the hybrid functional HSE06 was employed to compute electronic characteristics. The Hartree-Fock screening value was set at 0.2 Å [22]. For analyzing the band gap response to pressure variations, calculations were carried out using VASP with HSE06 functional. The bulk modulus and its pressure dependency were determined by fitting the pressure-volume (P-V) data to a third-order Birch-Murnaghan equation of state (EOS), elucidating the unit-cell volume’s response to compression [23].

3. Results

3.1. Structural Properties

Both mineral compounds feature an orthorhombic crystalline structure with the space group Fddd which is illustrated in Figure 1. In our case, calculations were conducted using GGA-PBE [19] and GGA-PBEsol [18] functionals. Our computational analysis indicates that the GGA-PBEsol functional accurately describes the ground state of the investigated compound, as summarized in Table 1. The predicted cell parameters closely match experimental values [6,15]. Table 1 reveals that PBE tends to slightly overestimate the cell parameters. The obtained atomic positions are listed in Table 2.

The thenardite-type minerals’ structure, as described in Figure 1, features tetrahedrally coordinated Si (Ge) and octahedrally coordinated Cd (Hg) atoms in Cd₂SiO₄ (Hg₂GeO₄). In Cd₂SiO₄, the O atom is coordinated by three Cd atoms and one Si atom. Similarly, in Hg₂GeO₄, the O atom forms bonds in a tetrahedral coordination with three identical Hg atoms and one Ge atom. The CdO₆ (HgO₆) octahedra link together via shared edges, creating zigzag chains aligned with the crystallographic [110] and [-110] directions. These chains establish a three-dimensional network, further connected through shared edges among CdO₆ (HgO₆) octahedra. Additionally, the SiO₄ (GeO₄) tetrahedra play a role in this connectivity, sharing edges with CdO₆ (HgO₆) octahedra. Both octahedra and tetrahedra exhibit significant distortion due to the unique connectivity of the SiO₄ (GeO₄) groups. In Cd₂SiO₄, the tetrahedron undergoes uniaxial elongation along its two-fold axis parallel to the c axis due to the sharing of two edges of the silicate tetrahedron with edges of the CdO₆ polyhedra. Conversely, in Hg₂GeO₄, the elongation of the tetrahedron along its two-fold axis, aligned with the b axis, is induced by the sharing of two edges between the germanate tetrahedron and the HgO₆ polyhedra.

Furthermore, we have computed the structural parameters as variables dependent on pressure, illustrated in Figure 2(a) and Figure 2(b). For Cd₂SiO₄, we observed a close correspondence between the volume (V/V₀) and lattice constants (a/a₀, b/b₀, c/c₀) with experimental data, as indicated in Figure 2(a). The response of lattice parameters to pressure exhibits significant anisotropy in both compounds. As pressure increases, all lattice parameters decrease in Cd₂SiO₄, whereas in Hg₂GeO₄, parameters a and b decrease with increasing pressure, while parameter c shows an increase under pressure. From Figure 2(a) and Figure 2(b), it is evident that both compounds exhibit less compressibility along the b axis. Additionally, the volume (V₀), bulk modulus (B₀), and its pressure derivative (B₀’) at zero pressure were determined through least-squares analysis of pressure-volume data. To calculate the bulk modulus for both compounds, we utilized the third-order Birch-Murnaghan equation of state (EOS). The resulting values for Cd₂SiO₄ are 702.03 ų, 120.53 GPa, and 4.43 for V₀, B₀, and B₀’ respectively, and for Hg₂GeO₄, they are 812.95 ų, 54.95 GPa, and 7.50. In Figure 2(c) and Figure 2(d), the unit-cell volume data plotted against pressure is shown, along with the pressure-volume curve determined using these fitted parameters. R Miletich obtained B₀ = 119.2(5) GPa and B₀’ = 6.17(4) for Cd₂SiO₄, indicating agreement with our theoretical results [15]. The compressibility of Cd₂SiO₄, with B₀ = 120.53 GPa, is lower than Hg₂GeO₄ and Cr₂SiO₄ (B₀ = 94.7(4) GPa) [16]. Other known cadmium oxides with compression data include CdO and CdWO₄, for which B₀ = 108 GPa [24] and 123 GPa [25] were reported, respectively. This is a consequence of the fact that in the three compounds compressibility is dominated by changes induced by pressure in the coordination polyhedral of Cd. Conversely, the compressibility of Cd₂SiO₄ is quite similar that of olivine-type M₂SiO₄ compounds containing large M cations, such as Fe₂SiO₄ (B₀ = 123.9 GPa) [26] or CaMgSiO₄ (B₀ = 113 GPa) [27].

The compression of both the CdO₆ and SiO₄ coordination polyhedra in Cd₂SiO₄ is affected by their respective polyhedral geometries. The CdO₆ polyhedra exhibit an increase in angular and bond-length distortion as pressure increases, as depicted in Figure 3(a). Figure 3(c) illustrates the Cd-O bond lengths as pressure varies, demonstrating anisotropic polyhedral compression. The longest bonds, Cd-O(5,6) share edge with the SiO₄ tetrahedra, resulting in the shortest inter-cation distance (Cd-Si distance of 3.098 Å). Conversely, the shorter Cd-O(1,2) and Cd-O(3,4) bonds engage in weaker cation-cation interactions across shared O-O edges. The atypical angular distortion of the CdO₆ octahedron, arising from polyhedral connections and shared edges, elucidates the alterations in polyhedral geometry induced by pressure. Significantly, the axes of the O-Cd-O octahedron exhibit substantial deviations from the ideal 180˚, with O(3)-Cd-O(4) measuring 141.5˚ and O(1)-Cd-O(6) at 113.39˚ as illustrated in Supplementary Figure S1. The displacement between Cd and O atoms along the c axis influences the O(3)-Cd-O(4) angle, whereas displacement along the b direction affects bond angles like O(5)-Cd-O(6), O(1)-Cd-O(4), and O(3)-Cd-O(5). Furthermore, the displacement of the Cd atom concerning the surrounding O atoms elucidates the distinction in compression between Cd-O(1) and Cd-O(2) bonds under pressure, leading to variations in compression along the crystallographic a and b axes, as depicted in Figure 2(a). In the SiO₄ tetrahedron, polyhedral volumes and Si-O distances remain relatively unchanged with pressure, as depicted in Figure 3(b) and Figure 3(c) respectively, while angular distortion notably decreases. The decrease in distortion signifies the diminishing uniaxial elongation caused by repulsion between Cd and Si atoms along shared edges. Likewise, O-Si-O angles tend toward the ideal tetrahedral angle with increasing pressure. The compression mechanism of Cd₂SiO₄ structure is mainly governed by cation-cation repulsions across shared O-O polyhedral edges. Symmetrically distinct metal-metal distances exhibit minimal compression up to 10 GPa, indicating stiffness relative to the overall structure. These repulsions induce distortions in CdO₆ octahedra and displace Cd atoms, leading to rapid compression of Cd-Cd(3) distance between opposing octahedra, even surpassing compression along the parallel c axis.

Likewise, we investigated the HgO₆ and GeO₄ polyhedra within the Hg₂GeO₄ compound. Both the HgO₆ and GeO₄ polyhedra demonstrate heightened angular distortion and bond-length variation with increasing pressure, as shown in Figure 4(a) and Figure 4(b). Figure 4(c) depicts the fluctuation in Hg-O bond lengths under varying pressures, highlighting anisotropic compression of the polyhedra. The longest bonds, Hg-O(3), share edges with the GeO₄ tetrahedra, resulting in the shortest inter-cation distance (Hg-Ge distance of 3.403 Å). Conversely, the shorter Hg-O(1) and Hg-O(2) bonds engage in weaker cation-cation interactions across shared O-O edges. The atypical angular deformation of the HgO₆ octahedron, stemming from polyhedral connections and shared edges, elucidates pressure-induced alterations in polyhedral geometry. Notably, the O-Hg-O octahedron axes exhibit significant deviations from the ideal 180˚, with O(3)-Hg-O(4) measuring 145.78˚ and O(1)-Hg-O(6) at 112.26˚ as shown in Supplementary Figure S2. The displacement of Hg atoms relative to O atoms along the b axis impacts the O(3)-Hg-O(4) angle, whereas displacement along the c direction alters bond angles such as O(5)-Hg-O(6), O(1)-Hg-O(4), and O(3)-Hg-O(5). Moreover, the displacement of the Hg atom relative to surrounding O atoms accounts for the discrepancy between Hg-O(1) and Hg-O(2) bond compressions with pressure, contributing to variations in compression along crystallographic a and b axes, as depicted in Figure 2(b). In the GeO₄ tetrahedron, polyhedral volumes and Ge-O distances remain relatively stable under pressure, as shown in Figure 4(b) and Figure 4(c), respectively. Unlike SiO₄, angular and bond-length distortion in GeO₄ increases with pressure.

3.2. Elastic and Mechanical Properties

A material’s elastic properties govern its reaction to stress, encompassing both deformation and the subsequent restoration to its original form when stress is relieved. These properties play a crucial role in revealing the bonding dynamics between neighboring atomic layers, the directional characteristics of binding, and the overall structural integrity. Elastic constants of solids serve as a bridge between their mechanical and dynamical behaviours, offering crucial insights into the forces operating within them. Furthermore, these constants serve as predictive tools for determining the structural stability of materials. In the case of orthorhombic symmetry, there are nine distinct elastic constants: C₁₁, C₂₂, C₃₃, C₄₄, C₅₅, C₆₆, C₁₂, C₁₃, and C₂₃ [28]. These elastic constants adhere to the generalized lattice stability criteria [29] across various pressure ranges, signifying the mechanical robustness of both compounds up to 10 GPa. With increasing pressure, almost all elastic constants experience growth, reflecting strong interactions between atoms. Consequently, the compounds exhibit enhanced strength, as depicted in the Figure 5(a) and Figure 5(b).

Employing the elastic constants acquired, we calculated the bulk modulus (B) and shear modulus (G) utilizing the Voigt-Reuss-Hill (VRH) approximation [30,31]. For Cd₂SiO₄, the calculated B and G values are 124.41 GPa and 42.05 GPa, respectively. Conversely, for Hg₂GeO₄, the B and G values are 73.06 GPa and 21.54 GPa, respectively. Notably, Cd₂SiO₄ exhibits a comparatively larger bulk modulus and shear modulus, indicating greater stiffness and resistance to deformation under stress. The determined B values at ambient pressure for both compounds align closely with the B₀ value derived from the Birch-Murnaghan equation of state (EOS). Additionally, Young’s modulus (E) can be derived from the bulk and shear moduli. For Cd₂SiO₄, the obtained value of E is 113.37 GPa, while for Hg₂GeO₄, it is 58.78 GPa. We have calculated the elastic moduli upto 10 GPa. As pressure increases, the values of B, G, and E also increase, as depicted in Figure 5(c) and Figure 5(d).

The values of ν and the B/G ratio serve to characterize the brittle or ductile nature of a structure. If ν and B/G are both less than 0.26 and 1.75, respectively, the structure is considered brittle; otherwise, it is regarded as ductile [32,33]. Our findings suggest that both proposed structures exhibit ductile behavior. Poisson’s ratio serves as an indicator of volume alteration during uniaxial deformation, with a value of ν = 0.5 indicating no volume change during elastic deformation. The low values observed for both compounds imply significant volume changes during their deformation. Additionally, ν provides insights into the bonding forces’ characteristics more effectively than other elastic constants [34]. It has been established that ν = 0.25 represents the lower limit for central-force solids, while 0.5 indicates infinite elastic anisotropy [35]. The low ν values observed for both structures suggest central interatomic forces within the compounds. In addition, we have determined the ν and B/G ratio under various pressures for both compounds. Both the ν and B/G ratio exhibit an increase as pressure increases, as illustrated in Figure 5(e) and Figure 5(f).

To assess the elastic anisotropy of the compounds under investigation, we acquired shear anisotropic factors, which measure the degree of anisotropy in atomic bonding across various planes. These factors play a critical role in evaluating the durability of materials. We calculated shear anisotropic factors for the {100} (A₁), {010} (A₂), and {100} (A₃) crystallographic planes, as well as percentages of anisotropy in compression (A_B) and shear (A_G) [36,37]. In a crystal displaying isotropy, the values of A₁, A₂, and A₃ should be one; any deviation from this indicates the presence of elastic anisotropy. A percentage anisotropy of 0% signifies perfect isotropy. For both compounds, the calculated shear anisotropy values under different external pressures are detailed in Table 3. For Cd₂SiO₄, anisotropies increase in the {010}, {001}, and {100} planes (A₂ < A₃ < A₁) at zero pressure, with the {100} plane exhibiting the highest anisotropy. Conversely, Hg₂GeO₄ displays an anisotropy sequence of A₃ < A₂ < A₁. Percentage anisotropies in compression and shear were approximately 5% and 8%, respectively for Cd compound whereas for Hg compound they are 32% and 22%, respectively. In addition, the universal anisotropy index (A^U) also provides the information of anisotropy in crystals [38]. The departure of A^U from zero denotes the level of single crystal anisotropy, incorporating both shear and bulk contributions, which sets it apart from other established metrics. Thus, A^U serves as a universal metric for quantifying single crystal elastic anisotropy. Cd₂SiO₄ exhibits an anisotropy index of 0.96, while Hg₂GeO₄ displays a value of 3.78.

Elastic wave velocities describe the rate at which waves travel through a substance when it undergoes elastic deformation. These waves consist of compression (longitudinal) waves (v_l), and shear (transverse) waves (v_t). For both compounds, the transverse (v_t) and longitudinal (v_l) mode velocities can be derived from the elastic constants [39]. The findings for both compounds are delineated in Table 4, indicating that the v_l increases with pressure, whereas the v_t initially rises before declining. Moreover, the Debye temperature (θ_D), a fundamental parameter, is linked to various solid-state properties like specific heat, elastic constants, and melting temperature. In this investigation, θ_D was determined from the mean elastic wave velocity (v_m) [40] and shown in Table 4. For both compounds, the calculated value of θ_D found to increase with pressure.

3.3. Lattice Dynamics

We conducted phonon dispersion analysis along the high-symmetry path (Γ-Y-T-Z-Γ-X) under zero pressure to assess the dynamical stability of the material. Figure 6(a) and Figure 6(b) display the calculated phonon dispersions for both compounds. Our findings indicate that the acoustic branches exhibit positive frequencies, confirming the dynamical stability of the system under analysis. Notably, there is a significant interaction between the optical and acoustic modes. Additionally, we computed the total and partial phonon density of states (PhDOS) for both compounds, as depicted in Figure 6(c) and Figure 6(d) for the Cd and Hg compounds, respectively. The PhDOS profiles of both compounds are primarily influenced by oxygen atoms. Furthermore, significant contributions of Cd and Hg atoms are evident in the lower frequency region (0-200 cm⁻¹), indicating that Cd and Hg atoms have lower frequencies with increasing atomic weight: Cd < Hg.

Moreover, we have computed the vibrational modes at zero pressure. Both compounds exhibit 42 vibrational modes. The mechanical representation of 39 optical modes is: Γ_M = 4A_g + 6B_1g + 5B_2g + 6B_3g + 4A_u + 5B_1u + 4B_2u + 5B_3u. The calculated optical modes contain the 21 Raman active (4A_g + 6B_1g + 5B_2g + 6B_3g), and 18 IR active (4A_u + 5B_1u + 4B_2u + 5B_3u) modes. The obtained frequencies of these modes are given in Table 5.

3.4. Electronic Properties

To thoroughly characterize the physical properties of the compounds under investigation, we employed the GGA-PBE functional to calculate their electronic band structures. The resulting band gap is 1.54 eV for Cd₂SiO₄ and 0.82 eV for Hg₂GeO₄, as shown in Figure 7(a) and Figure 7(b) respectively. Considering the known tendency of GGA to underestimate band gaps in insulators and semiconductors [41], we applied the HSE06 functional to improve the accuracy of the band gap determination. Figure 7(c) and Figure 7(d) illustrate the revised band gaps, which are 3.34 eV for the Cd compound and 2.09 eV for the Hg compound. Both compounds exhibit indirect gap characteristics, with the valence band’s highest point and the conduction band’s lowest point located at different positions. Furthermore, the analysis of the band structures indicates minimal dispersion in the valence band, while significant dispersion is observed in the conduction band. This implies that electrons will have a significant smaller effective mass than holes.

To gain a deeper insight into the electronic properties, we analyzed the total and partial density of states (DOS), as depicted in Figure 7(e) and Figure 7(f) using the HSE06 functional. From these figures, it is evident that the primary contributors to the highest peak of the valence bands are the O-3d states. Additionally, minor contributions from Cd-4d and Si-3s states are observed in the valence bands of the Cd compound, and Hg-5d and Ge-4s states in the Hg compound. The conduction bands of the Cd and Hg compounds are primarily dominated by Cd-5s and Hg-6s states, respectively. The Hg-6s orbitals are the responsible of the smaller band gap of Hg₂GeO₄. The same phenomenon is observed when PbWO₄ [42] is compared to other tungstates [43] due to the role of Pb-6s states.

Moreover, we carried out high-pressure calculations to investigate the electronic properties under different pressure conditions. The calculated electronic band gaps under pressure for both structures are depicted in Figure 8. As pressure rises, there is a notable expansion in the band gap of Cd₂SiO₄ (see Figure 8(a)), whereas in the case of the Hg compound depicted in Figure 8(b), the band gap decreases with increasing pressure. The widening of the band gap in Cd₂SiO₄ results from the intensified crystal field splitting between bonding and antibonding states under compression [44]. Conversely, the distinctive behavior of the band gap in Hg₂GeO₄ is attributed to the involvement of Hg-6s states in the lower portion of the conduction band, causing it to decrease under compression. Similar phenomena occur with Pb-6s states in PbWO₄ [42], in PbMoO₄ [45], and in PbCrO₄ [46].

4. Conclusions

In summary, we examined the structural, elastic, and electronic properties of Cd₂SiO₄ and Hg₂GeO₄ across different pressures using first-principles calculations based on density functional theory. Our analysis revealed that the computed lattice parameters at ambient pressure closely matched experimental data. The investigated compounds were determined to be mechanically and dynamically stable. We explored the pressure-dependent behavior of structural parameters. The compression of lattice parameters displayed an anisotropic behavior in both compounds. Notably, the compression mechanism in these structures is mainly governed by cation-cation repulsions across shared O-O polyhedral edges. The calculated elastic and mechanical properties reveal that Cd₂SiO₄ exhibits a comparatively larger bulk modulus and shear modulus, indicating greater stiffness and resistance to deformation under stress. We have computed the frequencies of vibrational modes, which include 21 Raman active modes, 3 acoustic modes, and 18 IR active modes. Both minerals were identified as indirect band gap insulators. Notably, Cd₂SiO₄ exhibited a significant increase in band gap with increasing pressure, while the band gap of Hg₂GeO₄ decreased under pressure. These findings revealed distinct structural and electronic responses despite the similar structures of the two compounds.