1. Introduction
Ice has a diverse presure temperature (P–T) diagram, with currently more than 17 crystal phases known under a variety of pressures and temperatures, including clathrate ice under negative pressure [1,2,3,4,5,6,7,8,9,10,11]. Natural gas hydrates are crystalline compounds composed of clathrate ice and natural gas. Compared with petroleum and coal, natural gas produces less carbon dioxide and pollution. The amount of gas hydrate stored in the earth is twice that of the existing oil and gas reserves and will be an important energy source in the future [12,13,14,15,16]. In the past, it was thought that empty hydrated lattices could not exist because guest molecules were needed to stabilize the host framework [17,18,19]. Later, Wooldridge and Jacobson et al. suggested that these guests could be removed, leaving a metastable, empty inclusion compound [20,21]. In 2014, Falenty et al. experimentally prepared a clathrate ice crystal by evacuating the S-II type cage hydrate [10]. The resulting crystal, named ice XVI, has the lowest density of any member of the ice family found to date (0.81 g/cm2) and its discovery opened up research into negative-pressure clathrate hydrates. Ice XVI was later reported to have good mechanical stability and negative thermal expansion [22]. However, information on its vibrational spectrum is lacking in the literature. This work focuses on the phonon statistics of ice XVI and prediction of its IR absorption, Raman scattering, and inelastic neutron scattering (INS) spectra. The typical vibration modes are discussed to clarify the dynamic process of lattice vibrations.
2. Computational Methodology
Using the CASTEP code [23], a first-principles density functional theory (DFT) method, we calculated the lattice vibrations of the ice XVI phase. On the basis of our preliminary calculations, the exchange-correlation (XC) function of revised Perdew-Burke-Ernzerhof (RPBE) [24], one of the generalized gradient approximation (GGA) method, was selected for geometry optimization. The energy convergence threshold of the self-consistent field (SCF) for geometric optimization was set as 1.0 × 10−9 hartree. The energy cut-off was set as 750 eV and the K-point separation was 0.07/Å. Norm-conserving pseudopotentials were used to calculate the phonon density of states (PDOS). We took the periodic structure of ice XVI from Ref. [10]. The original crystal data only included the positions of oxygen. To obtain a complete structure, a program was designed to position hydrogen atoms randomly within a 34-molecule supercell. The number of possible permutations and combinations of hydrogens is very large. In practice, we reduced the number of possible states by restricting the search to those with a uniformly distributed dipole moment. Thus, the hydrogen atoms were disordered, and the water molecules were linked by hydrogen bonds (H-bonds) to form a clathrate structure as shown in Figure 1. The supercell contained a cage composed of 28 molecules, forming 12 5-membered H-bonded water rings and 4 6-membered rings, denoted 51264. These rings were capable of enclosing large guest molecules. The geometry optimization was performed under zero pressure. Because RPBE slightly underestimates H-bond interactions, the density of the calculated crystal was 0.77 g/cm2, 5% lower than the experimental value. The H-bond lengths ranged from 1.797 Å to 1.845 Å, the length of O–H covalent bonds ranged between 0.984 and 0.990 Å, and the distance between adjacent oxygen ranged from 2.768 to 2.83 Å.
In addition, we counted the numbers of four-H-bond and two-H-bond normal modes based on the calculated results using a purpose-designed, self-compiled program. To clearly characterize the complex vibration modes, which are a key property of ice hydrates, we distinguished the mode types by one molecule in the supercell possessing the largest amplitude. The standard classification of stretching modes in terms of the vibrational direction was applied. Finally, we plotted a distribution diagram to mimic the PDOS of these two types of H-bonds.
3. Results and Discussion
Based on the harmonic approximation, the PDOS, infrared (IR) absorption, and Raman scattering spectra of ice XVI were calculated as shown in Figure 2. In theory, the PDOS spectrum integrates all phonons throughout the first Brillouin zone (BZ). The INS signals are proportional to the PDOS and can be compared with the PDOS spectrum. However, the interactions between photons and phonons occur only near the center of the BZ. The peaks in the Raman and IR spectra correspond to the frequencies of the calculated normal modes at the gamma point. Thus, we used these three simulated spectra as reference data for the INS, Raman, and IR experiments.
Because of the wide range of spectral intensities of Raman and IR, Figure 2 presents four sections of each of the spectra separately for comparison, on appropriately adjusted intensity scales. For a unit cell containing 34 molecules, there are 34 × 3 × 3 − 3 = 303 optical normal vibration modes. To illustrate the dynamic process of the normal modes, some typical modes of each section are presented for discussion.
The most interesting results are in the translational band. For a water molecule, there are three spatial vibration modes: stretching, wagging, and rocking. The mechanism of H-bond vibrations was the central topic in our previous studies [25,26,27,28,29]. We found two basic H-bond vibration modes in ice Ic, which gave rise to two triangular peaks in INS [27]. In all ice phases, each water molecule is linked with four neighbors via H-bonds to form a tetrahedral structure. According to our investigations, in such a local structure, when the central molecule vibrates along its angle bisector, the four connected H-bonds oscillate simultaneously to yield a four-bond mode. In contrast, there also exist vibrational motions along two H-bonds while the other two remain unchanged, corresponding to a two-bond mode. The vibration modes of H-bonds can thus be classified into two groups, four-bond modes and two-bond modes. We integrated these two groups in the region from 160 cm−1 to 320 cm−1 to mimic the PDOS spectrum of the H-bonds, as shown in Figure 3. The resulting spectrum matches the H-bond peaks in Figure 2 very well. The two-bond modes are distributed in a lower energy range while the four-bond modes are much stronger. In ideal ice Ic, the strength ratio of these two modes is 2 . Considering the distance between nearest neighbors in the hydrogen-disordered structure of ice XVI, the lengths of the H-bonds are calculated to range from 1.797 Å to 1.845 Å. Therefore, the strengths of the two kinds of modes have a distribution that gives rise to two main peaks, in which some values overlap. The simulated data confirmed that the two types of H-bonds also occur in clathrate ice, which may be of interest to researchers into natural gas hydrates.
For example, in the normal mode at 314 cm−1 shown in Figure 4, the molecules mainly stretch along the bisector of the angle in-plane. The molecule oscillates relative to its four neighbors in this mode. Under this situation, the four H-bonds act in synchrony. However, for the mode at 220 cm−1, the molecules undergo a combination of rocking and wagging along the H-bonded linkages. The vibration of each molecule only aligns with two H-bonds. Herein, the strength ratio of these two modes is approximately equal to 2 . We provide a video that depicts the normal modes at 220 cm−1 and 314 cm−1 to illustrate the dynamic process (Supplementary Materials). As shown in Figure 2, in the translational band, the Raman spectrum has a distinct peak at 212 cm−1. There is also a distinct peak at 175 cm−1 in the IR spectrum. Both modes correspond to the weak H-bond vibration and together make up the weak peak at around 200 cm−1. The main peak of the PDOS at around 300 cm−1 is inactive for IR absorption and very weak for Raman scattering.
For the intermolecular librational band, the frequency range is from 592 cm−1 to 1049 cm−1 with 102 normal mode patterns. This number is three times the number of molecules within a unit cell because there are three vibration modes for each molecule: rocking, wagging, and twisting. Note that the molecules move through space in the translational band, while the oxygen atoms remain static in the librational band. Very interestingly, one uniform vibration was found in this region for this ice phase. As shown in Figure 4, the absorption at the frequency of 592 cm−1 is a rocking mode for all molecules, and that at 1049 cm−1 is a twisting mode for all molecules. All of the peaks between these frequencies are mixtures of the three vibration types. This indicates that the rocking mode possesses the lowest energy while the twisting vibration has the highest. Three main peaks of Raman scattering are seen in Figure 2 at 648 cm−1, 737 cm−1, and 968 cm−1. The analysis of vibration modes for these three peaks shows that they are all mixtures of rocking, wagging, and twisting vibrations.
There are 34 normal modes in the bending zone, with a range from 1643 cm−1 to 1709 cm−1. In our previous studies, it was found that the vibrational energy increases from in-phase bending to out-of-phase bending. Due to the disordered structure of hydrogen, the synchronous opening and closing of all molecules is impossible. Figure 5 illustrates the minimum frequency of bending at 1643 cm−1. In this vibration mode the numbers of in-phase and out-of-phase motions are almost identical. The greatest frequency is at 1709 cm−1 and is an equal mixture of the two bending mode types. The Raman spectrum shows only one peak at 1682 cm−1, and the intensities of both Raman and IR are almost undetectably weak. However, peaks in this region can be detected in the INS experiment, as shown in the PDOS spectrum.
In the intramolecular O–H stretching region, each molecule has two kinds of vibration mode: symmetric stretching and asymmetric stretching. There are 68 normal modes in the stretching zone, with a range from 3105 cm−1 to 3368 cm−1. In Figure 5, four representative modes are selected, with the typical vibrations of the molecule shown in gold. To maintain a static center of mass, the vibrations of molecules are correlated in a crystal lattice. The minimum-energy mode at 3105 cm−1 is predominantly symmetric stretching with some coupled asymmetric stretching. Four peaks were detected by Raman scattering as shown in Figure 2. The minimum-energy peak is at 3180 cm−1 and is still mainly symmetric stretching. The highest-energy peak in the Raman spectrum in this region is at 3236 cm−1, in which most molecules undergo asymmetric stretching. The highest-energy peak of the entire Raman spectrum is at 3363 cm−1, in which almost all of the vibrations are asymmetric. The increasing trend of vibrational energies in this region from symmetric to asymmetric stretching is consistent with the literature [28]. Owing to selection rules, the IR absorptions are mostly symmetric stretching. In Figure 5, one can see that the highest-intensity peak of IR is at around 3180 cm−1.
4. Conclusions
Ice XVI, the first negative-pressure clathrate ice produced in the lab, was simulated using the first-principles DFT method. The vibrational spectra of INS, IR absorption, and Raman scattering were calculated theoretically. Four regions of the spectra are discussed individually together with some typical vibration modes. Compared with calculated results of ice Ih, its density is 0.77 kg/cm3, much lower than 0.88 kg/cm3 of Ih. Accordingly, the lengths of the H-bonds in ice XVI are bigger than Ih, manifesting an obvious red-shift of the H-bond (for the upmost mode, 314 cm−1 vs. 327 cm−1), and a slight blue-shift in the stretching region [25]. In the translational region, two H-bond peaks were found similar to the spectrum of ice Ic [27]. These H-bond vibrations could be classified into two groups based on their physical motions. The strong peak originated from four H-bonds of one molecule vibrating together while the weak peak corresponded to two H-bonds vibrating. We confirmed that this is a general rule that also applies to clathrate hydrates.
Based on this study, we also found a rule for the number of vibration modes in the four vibration regions in the ice. Assuming that the number of molecules in one unit cell of ice is N, the total number of optical normal modes should be 3N × 3 − 3. The O–H stretching region contains 2N modes corresponding to the two kinds of molecular stretching: symmetric and asymmetric. The bending region contains N modes as there is only one such mode for each molecule. The librational band has 3N modes representing the twisting, wagging, and rocking vibrations. Finally, the number of optical normal modes in the translational band is 3N − 3. The coefficient is 3 because there are three spatial vibration modes for each molecule, while 3 acoustic modes must be subtracted. This is a general formula governing any ice phase.
Figure 1. Diagram of ice XVI structure. (a) Unit cell of 34 molecules. The red balls indicate the oxygen atoms and the hydrogen atoms are in gray. The hydrogen bonds are shown with blue dashed lines; (b) Side view of a 4 × 4 × 4 supercell.
Figure 2. Vibration spectra of ice XVI simulated with the CASTEP code. From top to bottom: Raman scattering, IR absorption, and PDOS, which is proportional to inelastic neutron scattering. (a) Intermolecular translation band; (b) Intermolecular libration band; (c) Intramolecular bending band; (d) Intramolecular stretching band.
Figure 3. The distribution diagrams of four-bond modes (blue) and two-bond modes (red) of ice XVI in the range from 160 cm-1 to 320 cm-1. The standard classification of stretching modes in terms of the vibrational direction is applied.
Figure 4. Examples of normal modes of intermolecular vibrations. Representative vibrating molecules are shown in gold. For H-bond vibrations, the molecules present rocking or wagging vibration modes at 220 cm-1, while the molecules vibrate along the bisector of the angle for the mode at 314 cm-1. The mode at 592 cm-1 is the minimum-energy mode in the librational band and involves rocking of all of the molecules. The mode at 1049 cm-1 is the maximum-energy mode in this region and involves twisting vibrations for all molecules.
Figure 5. Four examples of normal vibration modes in the intramolecular O-H stretching region (3105 cm-1, 3180 cm-1, 3236 cm-1, 3363 cm-1) and two examples for bending (1643 cm-1, 1709 cm-1). The green arrows indicate the direction of vibration and their thickness is proportional to the amplitude. Typical vibration modes are shown in gold.
Supplementary Materials
The following are available online https://www.mdpi.com/1996-1944/12/2/246/s1, Video S1: Dynamic process of normal modes at 314 cm-1. Video S2: Dynamic process of three normal modes at 220 cm-1.
Author Contributions
Z.-R.W. performed simulations and drafted the manuscript. X.-L.Z. designed a special program. L.J. and K.Z. assisted data processing. H.-W.L. and Y.G. participated discussions. P.Z. conducted simulations, data analysis and edited the manuscript. All the authors gave final approval for publication.
Funding
This research received no external funding.
Acknowledgments
The numerical calculations were done on the supercomputing system in the Supercomputing Center, Shandong University, Weihai.
Conflicts of Interest
The authors declare no conflict of interest.
© 2019 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 (http://creativecommons.org/licenses/by/4.0/).
1. Bertie, J.E.; Calvert, L.D.; Whalley, E. Transformations of ice II, ice III, and ice V at atmospheric pressure. J. Chem. Phys. 1963, 38, 840.
2. Bertie, J.E.; Whalley, E. Optical spectra of orientationally disordered crystals. 2. Infrared spectrum of ice Ih and ice Ic from 360 to 50 cm-1. J. Chem. Phys. 1967, 46, 1271.
3. Whalley, E.; Heath, J.B.R.; Davidson, D.W. Ice 9-An antiferroelectric phase related to ice 3. J. Chem. Phys. 1968, 48, 2362.
4. Tajima, Y.; Matsuo, T.; Suga, H. Phase-transition in KOH-doped hexagonal ice. Nature 1982, 299, 810-812.
5. Kuhs, W.F.; Finney, J.L.; Vettier, C.; Bliss, D.V. Structure and hydrogen ordering in ices VI, VII, and VIII by neutron powder diffraction. J. Chem. Phys. 1984, 81, 3612-3623.
6. Hirsch, K.R.; Holzapfel, W.B. Effect of high-pressure on the Raman-spectra of ice-VIII and evidence for ice-X. J. Chem. Phys. 1986, 84, 2771-2775.
7. Lobban, C.; Finney, J.L.; Kuhs, W.F. The structure of a new phase of ice. Nature 1998, 391, 268-270.
8. Salzmann, C.G.; Radaelli, P.G.; Hallbrucker, A.; Mayer, E.; Finney, J.L. The preparation and structures of hydrogen ordered phases of ice. Science 2006, 311, 1758-1761.
9. Salzmann, C.G.; Radaelli, P.G.; Mayer, E.; Finney, J.L. Ice XV: A new thermodynamically stable phase of ice. Phys. Rev. Lett. 2009, 103, 105701.
10. Falenty, A.; Hansen, T.C.; Kuhs, W.F. Formation and properties of ice XVI obtained by emptying a type sII clathrate hydrate. Nature 2014, 516, 213-233.
11. Del Rosso, L.; Celli, M.; Ulivi, L. New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen-filled ice. Nat. Commun. 2016, 7, 13394. [Green Version]
12. Sloan, E.D. Fundamental principles and applications of natural gas hydrates. Nature 2003, 426, 353-359.
13. Yin, Z.Y.; Chong, Z.R.; Tan, H.K.; Linga, P. Review of gas hydrate dissociation kinetic models for energy recovery. J. Nat. Gas Sci. Eng. 2016, 35, 1362-1387.
14. Chong, Z.R.; Yang, S.H.B.; Babu, P.; Linga, P.; Li, X.S. Review of natural gas hydrates as an energy resource: Prospects and challenges. Appl. Energy 2016, 162, 1633-1652.
15. Koh, C.A. Towards a fundamental understanding of natural gas hydrates. Chem. Soc. Rev. 2002, 31, 157-167.
16. Davidson, D.W.; Handa, Y.P.; Ratcliffe, C.I.; Tse, J.S.; Powell, B.M. The ability of small molecules to form clathrate hydrates of structure-II. Nature 1984, 311, 142-143.
17. VanderWaals, J.H.; Platteeuw, J.C. Clathrate solutions. Adv. Chem. Phys. 1959, 2, 1-57.
18. Belosludov, V.R.; Inerbaev, T.M.; Subbotin, O.S.; Belosludov, R.V.; Kudoh, J.; Kawazoe, Y. Thermal expansion and lattice distortion of clathrate hydrates of cubic structures I and II. J. Supramol. Chem. 2002, 2, 453-458.
19. Matsumoto, M.; Tanaka, H. On the structure selectivity of clathrate hydrates. J. Phys. Chem. B 2011, 115, 8257-8265.
20. Jacobson, L.C.; Hujo, W.; Molinero, V. Thermodynamic stability and growth of guest-free clathrate hydrates: A low-density crystal phase of water. J. Phys. Chem. B 2009, 113, 10298-10307.
21. Wooldridge, P.J.; Richardson, H.H.; Devlin, J.P. Mobile Bjerrum defects-A criterion for ice-like crystal-growth. J. Chem. Phys. 1987, 87, 4126-4131.
22. Yagasaki, T.; Matsumoto, M.; Tanaka, H. Anomalous thermodynamic properties of ice XVI and metastable hydrates. Phys. Rev. B 2016, 93, 054118.
23. Clark, S.J.; Segall, M.D.; Pickard, C.J.; Hasnip, P.J.; Probert, M.J.; Refson, K.; Payne, M.C. First principles methods using CASTEP. Z. Kristallogr. 2005, 220, 567-570.
24. Hammer, B.; Hansen, L.B.; Norskov, J.K. Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 1999, 59, 7413-7421. [Green Version]
25. Zhang, P.; Tian, L.; Zhang, Z.P.; Shao, G.; Li, J.C. Investigation of the hydrogen bonding in ice Ih by first-principles density function methods. J. Chem. Phys. 2012, 137, 044504.
26. Zhang, P.; Han, S.H.; Yu, H.; Liu, Y. A calculating proof on hydrogen bonding in ordinary ice by the first-principles density functional theory. RSC Adv. 2013, 3, 6646-6649.
27. Yuan, Y.Z.; Zhang, P.; Yao, S.K.; Lu, Y.B.; Yang, H.Z.; Lou, H.W.; Zhao, Z.J. Computational assignments of lattice vibrations of ice Ic. RSC Adv. 2017, 7, 36801-36806. [Green Version]
28. Zhang, K.; Zhang, P.; Wang, Z.R.; Zhu, X.L.; Lu, Y.B.; Guan, C.B.; Li, Y.H. DFT simulations of the vibrational spectrum and hydrogen bonds of ice XI. Molecules 2018, 23, 1781.
29. Jiang, L.; Yao, S.K.; Zhang, K.; Wang, Z.R.; Luo, H.W.; Zhu, X.L.; Gu, Y.; Zhang, P. Exotic spectra and lattice vibrations of ice X using the DFT method. Molecules 2018, 23, 2780.
School of Space Science and Physics, Shandong University, Weihai 264209, China
*Author to whom correspondence should be addressed.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2019. This work is licensed under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
In practice, we reduced the number of possible states by restricting the search to those with a uniformly distributed dipole moment. [...]the hydrogen atoms were disordered, and the water molecules were linked by hydrogen bonds (H-bonds) to form a clathrate structure as shown in Figure 1. The geometry optimization was performed under zero pressure. Because RPBE slightly underestimates H-bond interactions, the density of the calculated crystal was 0.77 g/cm2, 5% lower than the experimental value. [...]we used these three simulated spectra as reference data for the INS, Raman, and IR experiments. Because of the wide range of spectral intensities of Raman and IR, Figure 2 presents four sections of each of the spectra separately for comparison, on appropriately adjusted intensity scales. According to our investigations, in such a local structure, when the central molecule vibrates along its angle bisector, the four connected H-bonds oscillate simultaneously to yield a four-bond mode.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer