3.1. Elastic Constants under Pressure
Elastic constants are important elasticity characteristics of a crystal’s unit cell that provide knowledge about its elasticity, its physical properties, and structural anisotropy. The hexagonal crystal symmetry of sH gas hydrate requires five independent elastic constants to determine its elastic properties which are C
11, C
12, C
13, C
33, and C
44. Its elastic constants matrix also has the dependent elastic constant C
66 which equals 0.5(C
11–C
12). The elastic constants of different sH gas hydrate structures were determined using first-principles computations at 0 K and stress-free conditions in our previous work [
3]. In this work, pressure-dependent elastic constants of empty, xenon-neohexane (Xe-NH), and carbon dioxide-neohexane (CO
2-NH) sH gas hydrates are investigated using first-principles at 0 K. The elastic constants of CO
2-NH sH gas hydrate at 0 GPa are also presented in this work in
Figure 1.
Comparing CO
2-NH sH hydrate to other sH gas hydrate structures at 0 GPa available in our previous work [
3] shows that CO
2-NH has C
11 value close to those reported for multiple filled sH gas hydrate structures. However, C
33 of CO
2-NH sH gas hydrate is lower than those reported for different sH gas hydrates including the empty metastable structure. This reflects the impact of the shape and chemistry of carbon dioxide in lowering the sH gas hydrate’s resistance to principal strain in the z-axis direction. Similarly the resistance to basal shear deformation is the lowest for CO
2-NH sH gas hydrate that has the smallest C
44 value. The elastic constants of CO
2-NH sH gas hydrate at 0 GPa are close to those reported for methane-propane sH gas hydrate from first-principles [
61], but with higher C
11. Jendi et al. [
9] reported the elastic constants of CO
2 sI gas hydrate from first principles. Comparing their values of C
11, C
12, and C
44 to those reported for CO
2-NH sH gas hydrate of this work shows that CO
2-NH sH hydrate has slightly higher resistance to principal strain in the x-axis direction (larger C
11) and a higher resistance to basal shear deformation (larger C
44). Using the structure stability criteria presented in
Section 2.2, all studied sH gas hydrates’ structures are found stable in the selected pressure range of investigation.
Figure 1 shows the elastic constants of CO
2-NH sH gas hydrate versus pressure and those of empty and Xe-NH sH hydrates are presented in
Table 2. What all systems share is the increase of C
11, C
33, C
12, and C
13 with increased applied pressure. C
44 of the empty sH gas hydrate shows very small dependence on pressure, while that of filled systems increases with applied pressure. This indicates the reduced resistance to shear deformation with increased pressure for empty sH gas hydrate. The rate of change of elastic constants of filled sH gas hydrates with pressure (ΔC
αβ/ΔP) is the highest for C
11 and C
33 and it is smaller for the other elastic constants. C
12, C
13 and C
44 are less affected by applied pressure as they reflect the shear elasticity [
23] and not the resistance to compression.
C
11 of empty sH hydrate is less than those of filled systems under all pressures considered in this work. This agrees with our previous findings [
3] of different sH gas hydrates at 0 K and 0 GPa. Xe-NH gas sH hydrate has a higher C
11 compared to CO
2-NH at 1 GPa, but at higher pressures CO
2-NH sH hydrate experiences higher C
11 values. C
33 of empty sH hydrate comes in the middle between those of Xe-NH and CO
2-NH at all pressure points. Xe-NH sH gas hydrate has the largest C
33 values at all pressures while CO
2-NH sH gas hydrate has the smallest C
33 values. Changing the help gas that occupies the small and medium cages of sH gas hydrate affects the bonding characteristics in the z-axis direction and can strengthen or weakens them as compared to those of empty sH structure. As clear from
Table 2, C
11 of empty sH gas hydrate is lower than its C
33 under all pressures. The same is observed for Xe-NH sH gas hydrate (except at 1 GPa) which reflects that atomic bonding is stronger in the z-axis (c-axis) direction as compared to the x-axis for those structures. However, C
11 of CO
2-NH sH gas hydrate is higher than its C
33 at all pressures. This indicates the stronger atomic bonding in the x-axis direction as compared to the z-axis for CO
2-NH sH gas hydrate [
18,
21]. This signifies the role of guest type and shape in characterizing sH gas hydrate’s properties in the different directions.
C
44 of the Xe-NH sH gas hydrate is the highest at all pressure points and it increases with compression, indicating the increased hydrate’s resistivity to basal shear deformation with pressure. The effect of pressure on sH gas hydrate’s resistivity to shear is obvious as CO
2-NH is less resistive to shear compared to empty structure at stress-free conditions (Ref. [
3]), however, this changes for higher pressures where C
44 values of CO
2-NH sH gas hydrate become larger than those of the empty structure for P ≥ 1 GPa. C
12 and C
13 elastic constants of empty sH gas hydrate are generally higher than those of Xe-NH and CO
2-NH sH gas hydrates (except for C
12 of Xe-NH at 1 GPa).
The rate of change of elastic constants with pressure for sH gas hydrates investigated in this work depends on the nature of guests inside the hydrate cages. The filled sH gas hydrates experience larger rates of change with pressure for C11 and C33 compared to the empty structure. However, empty sH gas hydrate experiences larger rates of change with pressure for C12 and C13 compared to the Xe-NH and the CO2-NH sH gas hydrates. CO2-NH has the largest rate of change with pressure for C11, while Xe-NH has the largest change of C33 with pressure. Those differences highlight the impact of guest chemistry on sH gas hydrate’s structure bonding in different directions.
The empty sH gas hydrate has close rates of change for C
11, C
33, C
12, and C
13. This agrees with the finding of Gagnon et al. [
36] who studied the effect of pressure on adiabatic elastic constants of ice (Ih). The percentage of change of C
12 and C
13 of ice (Ih) were found to be larger than those of C
11 and C
33. Gagnon et al. [
36] also found that C
44 of ice (Ih) undergoes a small negative change with applied pressure. This agrees with the small negative rate of change of C
44 with pressure of empty sH gas hydrate investigated in this work. This softening of C
44 could indicate the increased instability of the empty sH gas hydrate structure with applied pressure. It also reflects the increased tendency to structural phase transitions as pressure increases [
18]. This observation can be explained by the absence of guest molecules that provide the needed support for the hydrate’s structure via cage occupancy and guest-host interactions.
3.2. Polycrystalline Properties under Pressure
The polycrystalline properties of different sH gas hydrates—including empty and Xe-NH structures—at 0 GPa are available in our previous work [
3]. In this work, the CO
2-NH sH gas hydrate’s mechanical properties have been determined at 0 K and 0 GPa and are presented in
Table 3. The formulas used to compute those properties are available in
Appendix A. CO
2-NH sH gas hydrate has the lowest resistance to volume change (bulk modulus) among all structures previously investigated (Ref. [
3]). The bulk modulus of this hydrate is close to the isothermal bulk modulus for ice (Ih) at −35.5 °C [
36]. Its shear and Young’s moduli are close to those reported for empty and nitrogen-neohexane sH gas hydrates (Ref. [
3]). Poisson’s and shear to bulk moduli ratios of CO
2-NH sH gas hydrate are very close to those reported for Xe-NH sH hydrate (Ref. [
3]). The universal elastic anisotropy index of CO
2-NH sH gas hydrate is high, compared to that of empty, methane-neohexane, xenon-neohexane, argon-neohexane, and hydrogen-neohexane previously reported, but remains lower than that of nitrogen-neohexane sH gas hydrate (Ref. [
3]).
The linear shape and chemistry of carbon dioxide imposes extra anisotropy on the hydrate’s structure as shown by the three anisotropy measures (A
U, f
E, f
G) that deviate from isotropic values. The wave velocities of CO
2-NH sH gas hydrate are within the same range of values reported for other gas hydrate structures [
3,
9,
15]. The properties of CO
2 sI gas hydrate from first-principles were reported by Jendi et al. [
9]. Compared to CO
2 sI gas hydrate [
9], CO
2-NH sH gas hydrate has a close density value, a lower bulk modulus, higher Young’s and shear moduli, a smaller Poisson’s ratio, a close longitudinal wave velocity, and a lager transverse wave velocity. Jendi et al. [
9] showed that CO
2 sI gas hydrate has larger bulk modulus than methane sI gas hydrate. However, for sH gas hydrates of neohexane it is the opposite. Using CO
2 as a help gas in neohexane sH gas hydrate slightly lowers the structure’s bulk modulus (this work) compared to using methane (Ref. [
3]) as a help gas. Carbon dioxide’s quadrupole moment contributes to its intermolecular interactions with hydrate’s water molecules and other guests which impacts the overall elasticity of the hydrate’s structure.
To understand the effect of applied pressure on the mechanical properties of sH gas hydrates, the empty, Xe-NH, and CO
2-NH sH gas hydrate structures were investigated under multiple pressures. The pressure-dependent bulk (B), shear (G), and Young’s (E) moduli of CO
2-NH sH gas hydrate are presented in
Figure 2 and those of empty and Xe-NH sH gas hydrates are available in
Table 4. In general, the moduli of filled sH gas hydrates increase with increased pressure, however, the chemistry of guests play a role in the rate of change of moduli with pressure. The moduli of sH gas hydrates are affected by the hydrate’s hydrogen bonds and its strength that affects the structure’s compressibility and elasticity. The increase in bulk modulus with compression of empty, Xe-NH, and CO
2-NH sH gas hydrates can be linked to the strengthened hydrogen bonds under pressure which was discussed in our previous work [
7]. The rate of change of bulk modulus of sH gas hydrates with pressure is comparable to that of the pressure-dependent isothermal bulk modulus of ice (Ih) reported by Gagnon et al. [
36] at −35.5 °C. Manakov et al. [
62] computed the pressure-dependent bulk modulus from experimentally obtained pressure-volume data of xenon sI, ethane sI, and THF-xenon sII gas hydrates. Their results show the increase in bulk modulus of the three gas hydrates with increased pressure which agrees with this work’s findings of sH gas hydrate pressure-dependent bulk modulus at 0 K.
The pressure-dependent Young’s modulus computed from elastic constants in this work can be compared with that presented in our previous work [
7] which was obtained using IR-based technique. Both methods agree on the increase of sH gas hydrate’s Young’s modulus with compression. Young’s modulus of empty sH gas hydrate under compression obtained using elastic constants is smaller than that approximated using the hydrogen bonds’ IR stretching frequencies (Ref. [
7]). The difference between the two tend to increase with increased compression for empty sH gas hydrate. The opposite is observed for the Xe-NH and CO
2-NH sH gas hydrates with higher Young’s moduli obtained using elastic constants compared to the IR-based ones (Ref. [
7]). In most cases investigated here, the increased pressure of Xe-NH and CO
2-NH sH gas hydrates is accompanied by an improved agreement between Young’s modulus values obtained using elastic constants and those from the hydrate’s IR spectra relations (Ref. [
7]). This adds valuable knowledge on the limits and potentials of approximating sH gas hydrate’s mechanical properties using the vibrational characteristics of its hydrogen bonds. However, expanding the investigation to a wider pressure range is essential to better understand and evaluate the contribution of the hydrogen bond to the sH gas hydrate’s elasticity.
Young’s modulus of filled sH gas hydrate structures is higher than their bulk moduli under pressure. This agrees with our previous investigation of different sH gas hydrates at 0 K and 0 GPa (Ref. [
3]). However, the empty sH gas hydrate that has a larger Young’s modulus compared to bulk modulus at 0 GPa (Ref. [
3]) experiences the opposite under compression. The rates of Young’s modulus change of Xe-NH and CO
2-NH are much higher than those of empty sH gas hydrate. The empty structure’s bulk modulus is more pressure-dependent than its Young’s modulus. This is analogous to the findings of Gagnon et al. [
36] for ice (Ih). The bulk modulus is most sensitive to pressure for the CO
2-NH sH gas hydrate while the Xe-NH sH gas hydrate comes last in bulk modulus pressure-dependency. Moreover, Young’s moduli of Xe-NH and CO
2-NH sH gas hydrates are slightly more pressure-dependent compared to their bulk moduli. The shear modulus of sH gas hydrates is less sensitive to pressure change compared to other moduli, however, the presence of guest molecules tends to increase the shear modulus dependency on applied pressure. The resistance to shear increases under compression for both Xe-NH and CO
2-NH sH gas hydrates. However, the empty sH gas hydrate experiences a slight decrease in shear modulus with compression which is related to the softening of C
44 value as pressure increases. A similar result was obtained for the pressure-dependent shear modulus of hexagonal ice (Ih) by Gagnon et al. [
36].
To measure the degree of brittleness and ductility and their pressure-dependency, the Poisson’s ratio (ν), and the shear to bulk moduli (G/B) ratio are studied.
Figure 3 shows the pressure-dependent Poisson’s ratio (ν) and (G/B) for CO
2-NH sH gas hydrate and
Table 5 presents those of empty and Xe-NH sH gas hydrates. The ductile-brittle nature of sH gas hydrates is found to change slightly with pressure. Poisson’s ratio of empty and CO
2-NH sH gas hydrates slightly increase with compression, however, their G/B ratios slightly decrease with pressure. On the other hand, the Xe-NH sH gas hydrate experience slight fluctuations in ν and G/B with increased pressure. Empty sH gas hydrate is found to be more ductile under pressure compared to filled hydrate structures investigated here with higher Poisson’s ratio and lower G/B ratio. This agrees with our previous findings [
3] for different sH gas hydrates of neohexane.
The detection of gas hydrates’ location requires proper measurements of wave velocities or elastic speeds of sound. The importance of obtaining the longitudinal (V
L) and transverse (V
T) wave velocities is clear in the relationship between those velocities and other critical elastic properties of a material. The longitudinal (V
L) and transverse (V
T) wave velocities of different sH gas hydrate structures were computed from first principles at 0 GPa and are presented in our previous work [
3]. The pressure-dependency of those velocities for empty, Xe-NH, and CO
2-NH sH gas hydrates is discussed here.
Table 6 has V
T and V
L of empty and Xe-NH sH gas hydrates at three pressure points and 0 K and
Figure 4 shows the pressure-dependent velocities of CO
2-NH sH gas hydrate. For filled sH gas hydrates (Xe-NH, CO
2-NH) there is a clear increase of both velocities with increased pressure. However, the longitudinal wave velocity (V
L) experiences a larger rate of change with pressure compared to the transverse wave velocity (V
T). The transverse wave velocity of empty sH gas hydrate is weakly dependent on pressure and experience a slight decrease unlike that of filled sH gas hydrates. This is due to the slight decrease of empty sH gas hydrate’s shear modulus and its density increase with compression. Empty sH gas hydrate has the highest longitudinal wave velocity among the other systems presented here at all pressures. This can be explained by the higher density of filled sH gas hydrates compared to empty sH. Transverse wave velocity depends on shear modulus and density, while longitudinal wave velocity depends on Young’s modulus (and hence bulk and shear moduli), Poisson’s ratio, and density. The pressure-dependence of those velocities is basically understood through the pressure-dependence of Poisson’s ratio, density, and shear and bulk moduli.
Helgerud et al. [
32] investigated using experimental data the pressure-dependent wave speeds for ice (Ih), methane sI, and methane-ethane sII gas hydrates. They found that the compressional wave speed increases with pressure for the three structures which agrees with the findings presented here. However, they determined that the transverse wave speed of the three structures experience an uncommon decrease with applied pressure. This decrease in transverse wave speed with pressure was previously reported based on experimental findings for ice (Ih) by Gagnon et al. [
36]. Suwa et al. [
34] presented acoustic velocities of argon sII and methane sI gas hydrates at different pressures based on experimental data. The trend they presented shows the slight increase of longitudinal wave velocity with pressure with almost pressure-independent transverse wave velocity. The computed pressure-dependent longitudinal and transverse wave velocities from experimental data of cubic and hexagonal phases of methane hydrates of Beam et al. [
33] show an increase with increased pressure for both velocities. This agrees with the findings in this work for filled sH gas hydrates.
3.3. Direction-Dependent Elastic Properties under Pressure
Anisotropy of a material describes the direction dependence of its properties. It is essential to determine the anisotropy of gas hydrates to properly understand their direction-dependent mechanical properties which is important for gas hydrates’ detection and recovery. There are different anisotropy measurements available in the literature. In this work, the universal elastic anisotropy index (A
U) [
56] is computed for the empty, Xe-NH, and CO
2-NH sH gas hydrates under different pressures. This factor is general for any crystal symmetry and it is a generalization of the Zener anisotropy factor. For isotropic materials A
U equals zero and a deviation from this reflects the material’s structural anisotropy. The Young’s and shear moduli anisotropy factors can be obtained by taking the ratio of the moduli in the z-axis direction and that in the x-axis direction. A value of one for both moduli anisotropy factors is expected for an isotropic structure and a deviation from one indicates the material’s moduli anisotropy. This measurement was previously utilized by Tromans [
57] for computing different HCP metals’ anisotropy. The formulas of those factors are available in
Appendix A.
As observed by
Figure 5 and
Table 7, the universal elastic anisotropy index (A
U) increases with applied pressure for empty and CO
2-NH sH gas hydrates. This reflects the slight increase in difference between Voigt and Reuss limits of each moduli with compression. The filling of hydrate cages provides support to the structure via the guest-host intermolecular interactions; however, the shape and chemistry of guests might introduce additional anisotropy under pressure. As clear by the Xe-NH sH gas hydrate’s data, the universal elastic anisotropy index decreases and stabilizes with pressure up to 2.9 GPa, which might be due to the spherical shape of xenon. However, the CO
2 molecule is linear in shape which can impact the cages’ symmetry and hence overall structure’s anisotropy. A
U can reflect the structure’s bonding characteristics [
31] and its variation with applied pressure indicates the structure’s pressure-dependent anisotropy.
Murayama et al. [
63] investigated sH gas hydrate of neohexane and different help gases experimentally and using molecular dynamics. Their work discussed the impact of guest molecules’ type on sH gas hydrate’s anisotropic compressibility. The experimental findings of Takeya et al. [
64] revealed the anisotropy of sH gas hydrate unit cell parameters’ response to change in large guest molecule. In addition, Arai et al. [
65] used powder X-ray diffraction and molecular dynamics simulations to investigate the impact of guest molecules on sH gas hydrate structure. Their findings highlighted the anisotropy of the sH gas hydrate structure due to the inclusion of nitrogen molecule as a small guest. The anisotropy of the lattice constants temperature-dependence is explained by the rotational motion of the linear nitrogen molecule in addition to its anisotropic distribution inside the nonspherical medium cages. The findings of Arai et al. [
65] coincide with the outcomes of this work on the imposed anisotropy on neohexane sH gas hydrate due to the encapsulation of the linear carbon dioxide molecule in the small and medium cages of sH gas hydrate at 0 K and different pressures.
The Young’s and shear moduli anisotropy factors (f
E and f
G, respectively) vary with pressure and with guest type. For empty sH gas hydrate, Young’s modulus anisotropy factor (f
E) gets closer to a value of one with increased pressure (
Table 7). This means that as compression increases, Young’s modulus on the prismatic and basal planes become closer to each other. The opposite is observed for empty sH gas hydrate’s shear modulus that experiences an increase in its shear modulus anisotropy factor (f
G) which reflects the larger difference between shear modulus on the basal and prismatic planes as pressure increases. The values of f
E and f
G of Xe-NH sH gas hydrate (
Table 7) show the integration of guest chemistry and pressure in affecting the moduli anisotropy factors. Young’s modulus anisotropy factor (f
E) of Xe-NH sH gas hydrate is found to be most deviated from a value of one at 2 GPa, followed by 2.9 GPa, and 1 GPa. However, its shear modulus anisotropy factor (f
G) deviation from one at 2.9 GPa is bigger than those at 1 and 2 GPa.
Figure 5 shows that the pressure-dependent f
E and f
G of CO
2-NH sH gas hydrate decrease below one with increased pressure and hence reflect the increased structure’s moduli anisotropy. As evident by
Figure 5, the Young’s modulus anisotropy factor (f
E) of CO
2-NH sH gas hydrate is more pressure-dependent compared to its shear modulus anisotropy factor (f
G).
Elastic constants of sH gas hydrates can also reflect the structure’s anisotropy. C
11 and C
33 both define the material’s resistance to strain in the x-axis and the z-axis directions, respectively. The larger the difference between these two constants is, the higher is the elastic anisotropy of the structure. As pressure increases, the difference between C
11 and C
33 of CO
2-NH sH gas hydrate increases (
Figure 1) which explains its increased elastic anisotropy. In the same manner, the difference between C
44 and C
66 contributes to the anisotropy in shear modulus of sH gas hydrate. C
44 is related to the structure’s resistance to basal shear deformation, while C
66 is related the structure’s resistance to prismatic shear deformation in a hexagonal crystal [
66]. Computing the difference between C
44 and C
66 of CO
2-NH sH gas hydrate investigated in this work confirms the increased shear anisotropy of the structure with compression.
Different anisotropy factors indicate the level of anisotropy of a certain structure; however, the angular dependence of properties shows their variation with orientation and indicates the directions along which a property is maximum/minimum. Equations (4) and (5) compute the reciprocal of the Young’s and shear moduli, respectively, using the crystal’s compliances (S
αβ).
Figure 6 presents the relationship between the hexagonal unit cell axes (a, b, c) and the x, y, and z-axes. The vector (z
1) is the normal to a plane that intercepts the x, y, and z-axes and θ is its direction with respect to the z-axis (c-axis) of the unit cell. Further details on axes rotation are provided in reference [
57].
To characterize the variation of Young’s and shear moduli of sH gas hydrates with direction, the moduli versus the angle θ of empty, Xe-NH, and CO
2-NH sH gas hydrates are presented in
Figure 7,
Figure 8, and
Figure 9, respectively. The trends of the direction dependent moduli are the same for the same hydrate structure at different pressures.
Figure 7,
Figure 8 and
Figure 9 clearly show the cylindrical symmetry of sH gas hydrates’ Young’s and shear moduli with respect to the z-axis (c-axis, [0001]) of the hexagonal unit cell of sH gas hydrate.
As
Figure 7 clearly shows, the Young’s modulus of empty sH gas hydrate is maximum at θ = 0° or in the z-axis direction (on the basal plane) and another maximum is observed at θ = 90° (on the prismatic planes). Young’s modulus is minimum at 30° < θ < 60° for empty sH gas hydrate at different pressures. On the other hand, empty sH gas hydrate’s shear modulus has a maximum at 30° < θ < 70°. It is minimum in the z-axis direction (θ = 0°) and another minimum is observed at θ = 90°. The angular moduli curves also reflect the anisotropy of the structure at certain pressure. Angular Young’s modulus of empty sH gas hydrate on the prismatic and basal planes become closer to each other with increased pressure which agrees with the findings from the pressure dependency of f
E for this structure. However, the angular shear modulus on the prismatic and basal planes tend to deviate from each other with increased compression and that agrees with the behavior of f
G of empty sH gas hydrate under compression.
Figure 7 shows that Young’s modulus has little dependency on pressure, while shear modulus of empty sH gas hydrate slightly decreases with increased pressure. This is clear by the slight downshift of the shear modulus curve with increased compression (
Figure 7b). This observation agrees with the pressure-dependent polycrystalline moduli presented in
Table 4.
The Xe-NH sH gas hydrate’s direction-dependent moduli are presented in
Figure 8 which shows a similar behavior of both moduli with direction to that of empty sH gas hydrate. The Young’s modulus of Xe-NH sH gas hydrate is maximum on the basal (θ = 0°) and prismatic (θ = 90°) planes and is minimum for 30° < θ < 60°. However, the shear modulus is minimum on the basal and prismatic planes and has a maximum value in the direction corresponding to 20° < θ < 70°. The pressure also affects the location of the smallest value of shear modulus for Xe-NH sH gas hydrate. The smallest shear modulus for this structure is observed at θ = 90° when the structure is compressed to 1 GPa. However, the shear modulus is smallest in the z-axis direction (θ = 0°) for higher pressures. Angular Young’s modulus of Xe-NH sH gas hydrate at 2 GPa has the highest difference between its values on the basal and the prismatic planes, while that at 1 GPa has the smallest difference which agrees with f
E values from
Table 7. The outcomes from the angular shear modulus curves (
Figure 8b) agree with the f
G values change with pressure for Xe-NH sH gas hydrate (
Table 7). The effect of pressure on the direction-dependent Young’s and shear moduli of Xe-NH is observed by the upward shift in their curves as pressure increases. This coincides with the pressure dependence of the polycrystalline moduli presented in
Table 4 for Xe-NH sH gas hydrate.
The impact of guest type and chemistry on the angular properties of sH gas hydrates is apparent through the direction-dependent Young’s modulus of CO
2-NH sH gas hydrate (
Figure 9a). The general behavior of the direction-dependent Young’s modulus of CO
2-NH in terms of the location of moduli maxima and minima is comparable with those of empty and Xe-NH sH gas hydrates. However, Young’s modulus presented in
Figure 9a has the biggest value on the prismatic planes (θ = 90°) and has another (smaller) maximum on the basal plane (θ = 0°) which is opposite to what is observed in
Figure 7a and
Figure 8a where the biggest Young’s modulus is in the z-axis direction. The shear modulus of CO
2-NH sH gas hydrate is minimum on the basal and the prismatic planes and it is maximum for 30° < θ < 70° which resembles the shear moduli curves of empty (
Figure 7b) and Xe-NH (
Figure 8b) sH gas hydrates. The increased anisotropy of both Young’s and shear moduli of CO
2-NH sH gas hydrate with compression is noticeable in the rising difference between the moduli at θ = 0° and θ = 90° with compression which acknowledges the pressure-dependent anisotropy factors presented in
Figure 5. In agreement with the polycrystalline pressure dependence of CO
2-NH sH gas hydrate’s moduli (
Figure 2), the angular Young’s and shear moduli (
Figure 9) increase with compression.
To better characterize the anisotropy in Young’s and shear moduli of sH gas hydrates under pressure, the projections of the moduli on polar plots are used where the y-axis of the plot is the vector parallel to the z-axis (c-axis) and the x-axis of the plot is the vector perpendicular to the z-axis (c-axis) of the sH gas hydrate’s unit cell. Isotropic materials have a plot of circular shape and the deviation from this constant curvature shape reflects the material’s moduli anisotropy [
42].
Figure 10 shows the Young’s (a) and shear (b) moduli projections or polar plots of empty sH gas hydrate. The empty structure experiences a distorted circular symmetry of both moduli at all pressures. However, the deviation from the circular shape is more pronounced for the Young’s modulus (
Figure 10a). The shear modulus of empty sH gas hydrate in polar coordinates slightly deviates from the circular shape reflecting the small shear moduli anisotropy that marginally increases with pressure (
Figure 10b). The higher anisotropic characteristics of Young’s modulus compared to shear modulus anisotropy of empty sH gas hydrate agrees with the findings for hexagonal ice (Ih) [
38]. Franco Pinheiro Moreira et al. [
38] evaluated the Young’s and shear moduli of ice (Ih) using molecular dynamics and different water models at zero temperature and zero pressure. They also discussed the impact of chosen water model on elastic anisotropy. Their results show that angular Young’s modulus of ice (Ih) deviates from the circular shape more than its shear modulus does. Notably this particular behavior of ice (Ih) moduli [
38] is equivalent to that of the hexagonal empty sH gas hydrate of this work, highlighting the impact of symmetry.
As previously discussed, the polar-coordinate plots of moduli projections (
Figure 10) of empty sH gas hydrate demonstrate that its Young’s modulus is more anisotropic than its shear modulus. This agrees with the moduli anisotropy factors (
Table 7) of this structure except at 3 GPa. It is critical to highlight that using the moduli anisotropy factors does not always accurately reflect the overall moduli anisotropic behavior. For example, a system can have equal Young’s modulus in the z-axis direction and that in the direction perpendicular to it (E(θ = 90°) = E(θ = 0°) → f
E = 1) which differs from its Young’s modulus at θ = 45°. In that case using the moduli anisotropy factor alone can give inaccurate indication of a structure’s elastic anisotropy.
Figure 11 presents the moduli polar-coordinate plots of Xe-NH sH gas hydrates and shows the impact of pressure on moduli anisotropy. The higher the pressure, the more anisotropic the moduli are for Xe-NH sH gas hydrate. Young’s modulus of Xe-NH sH gas hydrate is more anisotropic than its shear modulus with higher distortion of the circular symmetry. At 1 GPa, the shear modulus of Xe-NH is in fact isotropic with an almost circular projection of its shear modulus. Despite that f
G values of the Xe-NH sH gas hydrate at 1 and 2 GPa are almost equally different from one, the 2 GPa projection of shear modulus appears more anisotropic compared to that at 1 GPa. This once again illustrates that the moduli projection or polar plot is a better representation of the overall moduli direction dependency. The polar-coordinate plots also reflect the effect of compression on moduli values and proves the higher dependency of Xe-NH sH gas hydrate’s moduli on applied pressure compared to that of the empty sH gas hydrate.
As discussed earlier; the linear shape and chemistry of the CO
2 molecule play a role in defining the properties of gas hydrates encapsulating it. Understanding the factors affecting carbon dioxide gas hydrates is essential for a better interpretation of their mechanical properties and how they behave under different conditions. The angular dependence of this hydrate’s properties is crucial for engineering applications related to carbon dioxide sequestration.
Figure 12 presents the polar-coordinate Young’s and shear moduli of CO
2-NH sH gas hydrate at equilibrium (0 GPa) and at different pressures. As evident by
Figure 12, Young’s modulus of CO
2-NH sH gas hydrate is more anisotropic than its shear modulus which agrees with the earlier findings of empty and Xe-NH sH gas hydrates.
The shear and Young’s moduli anisotropy of the CO
2-NH sH gas hydrate increases with pressure and is more pressure-dependent compared to those of empty and Xe-NH sH gas hydrates. In general, the filling of sH gas hydrate changes its moduli anisotropy which is tunable by the guest shape and chemistry. This was previously concluded from our first-principles study of different sH gas hydrates at 0 GPa (Ref. [
3]). This work demonstrates that the compression of neohexane structure-H gas hydrates containing xenon and carbon dioxide as guests is accompanied by an increased anisotropy of the structures’ Young’s and shear moduli. This calls for expanding the investigation of sH gas hydrate’s mechanical properties in different chemical environments using different guest molecules.
Panda and Ravi Chandran [
43] explained the moduli anisotropy of hexagonal TiB
2 by interpreting its chemical bonding. They attributed the elastic anisotropy on the YZ plane to the difference in the structure’s bond strengths. In a similar manner, the Young’s modulus anisotropy of sH gas hydrate in the z-axis direction can be explained by the variation of bond strength among the OH covalent bonds, the hydrogen bonds, and the guest-host interactions. As previously discussed for sH gas hydrates [
3,
7], the structures’ hydrogen bonds contribute to their moduli. The difference in hydrogen bond density and strength across the different planes of sH gas hydrate can influence the structure’s bulk and Young’s moduli. This could explain the more pronounced Young’s modulus anisotropy compared to that of shear modulus for sH gas hydrates. The guest-host interactions add to the chemical bond environment variations in sH gas hydrate’s structure and should be carefully studied.
Investigating the sH gas hydrate’s unit cell elastic constants, polycrystalline properties, and direction-dependent moduli under compression highlights noticeable changes in its stiffness and resistance to changes in shape and volume. It is of interest to examine sH gas hydrate’s elastic and direction dependent behavior under tensile stress as well. For this purpose, the empty sH gas hydrate structure was studied under a pressure of −1 GPa. The different elastic properties of this structure under a pressure of −1 GPa are listed in
Table 8. The elastic constants are lower than those reported in
Table 2 for empty sH gas hydrate under pressure of 1 to 3 GPa. They are also lower than those previously reported for empty sH gas hydrate under zero temperature and pressure (Ref. [
3]). The bulk modulus is noticeably reduced under tensile, while shear and Young’s modulus are less affected. Brittle nature is higher compared to the system under zero stress (Ref. [
3]) and under compression (
Table 5) as observed by the lower Poisson’s ratio (ν) and the higher G/B ratio. The longitudinal wave velocity (V
L) is more affected by tensile stress compared to the transverse wave velocity (V
T), which was previously confirmed for sH hydrates under compression.
The values of the elastic constants reflect the increased isotropic features in empty sH gas hydrate under tensile conditions. C
11 and C
33 are very close to each other indicating isotropy in elastic moduli. In addition, C
66 (= 0.5(C
11-C
12) = 4.80 GPa) is very close to C
44 reflecting the isotropy in shear. This is further confirmed by the moduli anisotropy factors f
E and f
G with values almost equal to one, and the universal elastic anisotropy index A
U that has a nearly zero value.
Figure 13 also confirms the isotropic nature of Young’s and shear moduli of empty sH gas hydrate under a pressure of −1 GPa.