Next Article in Journal
Paramylon, a Potent Immunomodulator from WZSL Mutant of Euglena gracilis
Previous Article in Journal
Bio-Polyethylene-Based Composites Reinforced with Alkali and Palmitoyl Chloride-Treated Coffee Silverskin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Analysis of Exotic Molecular and Atomic Vibrations in Ice XV

School of Space Science and Physics, Shandong University, Weihai 264209, China
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(17), 3115; https://doi.org/10.3390/molecules24173115
Submission received: 19 August 2019 / Revised: 26 August 2019 / Accepted: 26 August 2019 / Published: 27 August 2019
(This article belongs to the Section Molecular Structure)

Abstract

:
It is always difficult to assign the peaks of a vibrational spectrum in the far-infrared region. The two distinct peaks seen in many ice phases are still a mystery to date. The normal modes of ice XV were calculated using the CASTEP code based on first-principles density functional theory. On the basis of vibrational modes analysis, we divided the translational modes into three categories: four-bond vibrations, which have the highest energy levels; two-bond vibrations, which have medium levels of energy; and relative vibrations between two sublattices, which have the lowest energy. Whale et al. found that some intramolecular stretching modes include the isolated vibration of only one O–H bond, whereas the others do not vibrate in ice XV. We verified this phenomenon in this study and attributed it to local tetrahedral deformation. Analysis of normal modes, especially in the translation and stretching band of ice XV, clarified the physical insights of the vibrational spectrum and can be used with other ice phases.

1. Introduction

More than 17 ice phases have been experimentally identified at various temperatures and pressures [1]. These phases can be divided into two main categories on the basis of hydrogen long-range order or disorder. Among ice phases, seven hydrogen-ordered phases are currently known, ices II, VIII, IX, XI, XII, XIV, and XV. Most of them can be converted into hydrogen-disordered ice phase, in which the water molecules exhibit orientation disorder after heating. There are six known counterparts: Ih/XI, III/IX, V/XIII, VI/XV, VII/VIII, and XII/XIV [2]. The only exception is ice II, for which a hydrogen-disordered counterpart has not been observed. Upon heating, it transforms to either ice Ill or ice V depending on the pressure. The high-pressure phase of ice VI was first discovered by Bridgman [3]. The crystal structure of ice VI was solved by Kamb via X-ray diffraction at ambient pressure at the temperature of liquid nitrogen [4], and confirmed by Kuhs et al. via neutron diffraction [5,6]. Considering only oxygen lattices, the primitive cell of ice VI contains 10 water molecules with space group P42/nmc, which is a hydrogen-disordered phase. Ice VI consists of two interpenetrating yet independent hydrogen-bonded networks, and has thus been described as a “self-clathrate” [4].
In 2009, Salzmann et al. found a new hydrogen-ordered phase by doping HCl in ice VI at a lower temperature, which is called ice XV. This hydrogen-ordered structure is antiferroelectric ( P 1 ¯ ) and can be stable from 0.8 GPa to 1.5 GPa at temperatures below 130 K [7]. Whale et al. measured the Raman spectrum and performed simulations with the first-principles density functional theory (DFT) [8]. Salzmann et al. later presented another hydrogen-ordered structure of VI, the ferromagnetic structure with space group Cc [9,10]. In 2018, Gasser et al. demonstrated the existence of one more XV phase, which they called the ice beta-XV phase [11].
For theoretical investigations of vibrational modes, the most interesting results were those of Whale et al., who found exotic “rigid network” molecular vibrations in the translation band and “isolated stretching” modes of atomic vibrations in the stretching band [8]. We also found “rigid network” modes in ice VIII and VII [2] and “isolated stretching” modes in ice XIV [12]. According to our series of investigations of ice phases, we identified two intrinsic categories of molecule vibrations in the translational band that constitute two sharp peaks [2,12,13,14,15]. In this paper, we confirm the presence of cluster vibrations in a much lower region of the translation band. We review the normal modes in ice XV and aim to clarify the physical insights of these exotic vibration modes.

2. Computational Methods

We first constructed the oxygen structure of ice VI. Then, we used a self-compiled program to find two possible hydrogen-ordered structures. We selected the antiferromagnetic one that is the same as reported by Whale for this work. The density is 1.32 g/cm3, in which the three orthorhombic constants of a primitive cell are a = 6.23 Å, b = 6.22 Å, and c = 5.83 Å. Using the CASTEP [16] code based on first-principles density functional theory (DFT), we calculated the phonon density of states (PDOS) of antiferroelectric ice XV. The energy threshold and self-consistent field tolerance were set as 1 × 10−9 eV/atom. The energy cutoff was 750 eV, and the k-point was 2 × 2 × 2. As the electron density has a large distribution in ice lattice, we chose generalized gradient approximation (GGA) for this work. To test the exchange-correlation (XC) functional, we tried RPBE, PBE, PW91, WC, and PBESOL separately for geometric optimization. The results showed that the densities are 1.21, 1.40, 1.41, 1.55, and 1.56 g/cm3, respectively. As the length of intermolecular hydrogen bond (H-bond) is relative to volume change and the vibrational frequency is very sensitive to H-bond length, the best selection is RPBE [17], which produces a slight redshift in the translation band, but greater accuracy in the intramolecular vibrations. The hydrostatic pressure was set to 0.9 and 1.4 GPa for comparison with the k mesh was set as 4 × 4 × 4. We found the H-bond frequencies of 1.4 GPa have a big redshift than 0.9 GPa (please see Supplementary File S1). Thus, we selected the simulation data of 0.9 GPa for discussions hereinafter. The optimized dimensions are a = 6.39 Å, b = 6.44 Å, and c = 6.02 Å.

3. Results and Discussion

Figure 1 shows the computing spectra of Raman scattering, infrared (IR) absorption, and PDOS, which are divided into four parts according to four separate vibrational regions. Because of the wide scale of the intensity in different regions, we adjusted their proportions for comparison. A primitive cell contains 10 molecules; thus, there are 10 × 3 × 3 − 3 = 87 optical normal modes.
Table 1 shows the details of PDOS, normal modes, and Raman scattering peaks. Column 1 lists the main peaks of PDOS. We identified 34 distinguishable peaks, of which 11 are major. The PDOS curve integrates phonon signals throughout the Brillouin zone (BZ), so it could correspond to an inelastic neutron scattering (INS) spectrum. Because reports of INS data are lacking, we present PDOS herein for reference. The 87 normal modes are shown in column 2, and the IR and Raman active intensities are shown in the next two columns. Sixteen Raman peaks are measured in the experiment: 114, 158, 186, 220, 230, 300, 312, 387, 451, 499, 548, 3224, 3335, 3345, 3396, and 3425 cm−1 [8]. However, the peak at 387 cm−1 did not match any normal mode. We doubt that this is a false peak because a gap should be present between the translation and libration bands. There is a slight underestimation of the hydrogen bond (H-bond) by RPBE, but it is more accurate than Whale’s work by PBE [8].
In the translational band, the Raman peaks are from optic modes near the BZ that are detectable subject to selection rules. On the basis of harmonic approximation, 27 optic normal modes are present between 73 and 285 cm−1. Note that the peaks are matched with normal modes by the calculated intensities.
Our previous study of ice Ic revealed two kinds of intrinsic translational modes in the ice lattice [13]. In the strong mode, the molecule vibrates together with four connecting H-bonds, which is called the four-bond mode. Two perpendicular degenerate weak modes involve only two oscillating H-bonds, which is called the two-bond mode. Treating the two H-bonded molecules as spring, a simple harmonic oscillator model yielded a strength ratio of 2 .
Although the tetrahedral structure of ice XV presents deformation under pressure to some extent, these two kinds of modes can also be distinguished in the translation region. For the mode at 285 cm−1 shown in Figure 2, typical molecules vibrate along its HOH angular bisector, which yields a strong vibrational mode (please see Supplementary File S2). Note that some small vibrations of other molecules can keep the mass center of a primitive cell static. The mode at 206 cm−1 is a two-bond mode, in which a molecule vibrates toward its two neighbors (please see Supplementary File S3). Whale et al. found that the three lowest modes at 66, 78, and 85 cm−1 present vibrations of the two rigid ice XV networks. In this study, the three modes are at 73, 106, and 107 cm−1, respectively. They actually represent three orthorhombic vibrations in 3D space that originate from the interactions of two sublattices without H-bonds, and thus present the lowest vibrational energy. Figure 2 illustrates an example of the mode at 107 cm−1 (please see the dynamic process from Supplementary File S4). We also found this phenomenon in ice VII and VIII [1]; however, the most common case in other ice phases is that many cluster vibration or skeleton deformation modes are present in the much lower region. Considering the translational band of ice, we can divide normal modes into three categories: four-bond vibrations, which possess the highest energy; two-bond vibrations, with medium energy levels; and cluster vibrations, which have the lowest energy. Figure 3 shows a diagram of fitted curves of three kinds of modes in ice XV compared with the inserted PDOS. Note that the data for normal modes are in the BZ center, whereas the real PDOS covers the entire first BZ. Thus, some disagreements occur as a result of the dispersion of phonons.
As shown in Figure 2, the libration band contains 30 normal modes from 477 to 936 cm−1. There are three types of vibration modes: rocking, or rotation of the whole molecule around an axis perpendicular to the molecular plane; twisting, or rotation around an axis coincident with the HOH angle bisector; and wagging, or rotation of the molecule around an axis in the molecular plane, perpendicular to the bisector of the HOH angle. For an optic mode, the mass center of the primitive cell remains static while the molecules within vibrate. Figure 2 shows that the main vibrating mode at 477 cm−1 is rocking, whereas the mode at 913 cm−1 is twisting. For the highest frequency at 936 cm−1, the vibrating mode of molecules is a mixture of wagging, rocking, and twisting.
The bending band has 10 normal modes from 1609 to 1709 cm−1, as shown in Figure 4. The vibration pattern of each molecule is bending, whereas the collective vibrations can be distinguished as in-phase and out-of-phase. For the highest-frequency mode at 1709 cm−1, the vibrating mode of crystal lattice is in-phase. The lowest mode at 1609 cm−1 is out-of-phase. The mode at 1684 cm−1 shown in Figure 4 has the greatest intensity of IR absorption, but has no Raman activity. In contrast, the mode at 1655 cm−1 has the greatest intensity of Raman activity, but no IR activity in this band. Because the intensities are very small compared with the stretching band, it is always difficult to detect signals for IR and Raman experiments.
There are 20 normal modes in the stretching band from 3278 to 3485 cm−1, including intramolecular symmetrical and asymmetrical stretching vibrations. For the lowest frequency at 3281 cm−1, all molecules show symmetric stretching. For the frequency at 3461 cm−1, the main vibrations show asymmetric stretching. The mode at 3280 cm−1 has the greatest intensity of Raman activity, but no IR activity, whereas the opposite is found in the mode at 3278 cm−1.
Whale et al. reported the normal modes in this band; some vibrational modes include isolated vibration of only one O–H bond, whereas the other does not vibrate [8]. A typical example is shown in Figure 4 of the mode at 3367 cm−1 (please see Supplementary File S5). We also found this vibration pattern in other ice phases, such as ice XIV and VII. We regard this “isolated stretching” mode as a special case of asymmetric stretching and attribute it to local tetrahedral deformation by pressure and a hydrogen-disordered lattice.

4. Conclusions

It is traditionally difficult to assign the peaks of a vibrational spectrum in the fingerprint region. We analyze the normal modes of ice XV via DFT simulation. On the basis of our series of investigations in the translational band of ice, we divide normal modes into three categories: four-bond vibrations that possess the highest energies; two-bond vibrations that possess medium energies; and cluster vibrations that have the lowest energy. In ice XV, VII, and VIII, the cluster modes are relative vibrations between two sublattices. Thus, we can analyze the translational modes in the far IR region of ice. This study presents physical insights into the vibrational spectrum and can be used with other ice phases.
In addition, we reviewed the “isolated stretching” mode first reported by Whale et al. As we also found this phenomenon in other ice phases, we do not regard this mode as a third type in addition to symmetric and asymmetric stretching, and instead attribute it to a special asymmetric stretching due to local tetrahedral deformation by pressure and a hydrogen-disordered lattice.

Supplementary Materials

The following are available online. File S1: The results of calculations at 1.4 GPa and 0.9 GPa, the k mesh was set as 4 × 4 × 4. Video S2: Dynamic process of four-bond vibration at 285 cm−1. Video S3: Dynamic process of two-bond vibration at 206 cm−1. Video S4: Dynamic process of relative vibration between two sublattices at 107 cm−1. Video S5: Dynamic process of isolated stretching at 3366 cm−1.

Author Contributions

X.-L.Q. drafted the manuscript. X.-L.Z. and J.-W.C. assisted structural modeling and data processing. L.J., Y.G., and X.-C.W. participated in the discussion of results. All authors gave final approval for publication. P.Z. conducted simulations and data analysis and edited the manuscript.

Funding

This research received no external funding.

Acknowledgments

The numerical calculations were conducted on the supercomputing system in the Supercomputing Center, Shandong University, Weihai.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Gu, Y.; Zhu, X.L.; Jiang, L. Comparative analysis of hydrogen bond vibrations in ice VIII and VII. J. Phys. Chem. C 2019, 123, 14880–14883. [Google Scholar] [CrossRef]
  2. Moberg, D.R.; Sharp, P.J.; Paesani, F. Molecular-level interpretation of vibrational spectra of ordered ice phases. J. Phys. Chem. B 2018, 122, 10572–10581. [Google Scholar] [CrossRef] [PubMed]
  3. Bridgman, P.W. Water, in the Liquid and Five Solid Forms, under Pressure. Proc. Am. Acad. Arts Sci. 1912, 47, 441. [Google Scholar] [CrossRef]
  4. Kamb, B. Structure of ice VI. Science 1965, 150, 205. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  6. Kuhs, W.F.; Ahsbahs, H.; Londono, D.; Finney, J.L. In-situ crystal growth and neutron four-circle diffractometry under high pressure. Phys. B 1989, 684, 156–157. [Google Scholar] [CrossRef]
  7. 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. [Google Scholar] [CrossRef] [PubMed]
  8. Whale, T.F.; Clark, S.J.; Finney, J.L. DFT-assisted interpretation of the Raman spectra of hydrogen-ordered ice XV. J. Raman Spectrosc. 2013, 44, 290–298. [Google Scholar] [CrossRef]
  9. Komatsu, K.; Noritake, F.; Machida, S. Partially ordered state of ice XV. Sci. Rep. 2016, 6, 28920. [Google Scholar] [CrossRef] [PubMed]
  10. Nanda, K.D.; Beran, G.J.O. What governs the proton ordering in ice XV? J. Phys. Chem. Lett. 2013, 4, 3165–3169. [Google Scholar] [CrossRef]
  11. Gasser, T.M.; Thoeny, A.; Plaga, L. Experimental evidence for a second hydrogen ordered phase of ice VI. Chem. Sci. 2018, 9, 4224–4234. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, K.; Zhang, P.; Wang, Z.R.; Zhu, X.L.; Lu, Y.B.; Guan, C.B.; Li, Y. DFT simulations of the vibrational spectrum and hydrogen bonds of ice XIV. Molecules 2018, 23, 1781. [Google Scholar] [CrossRef] [PubMed]
  13. Yuan, Z.Y.; Zhang, P.; Yao, S.K.; Lu, Y.B.; Yang, H.Z.; Luo, H.W.; Zhao, Z.J. Computational assignments of lattice vibrations of ice Ic. RSC Adv. 2017, 7, 36801. [Google Scholar] [CrossRef]
  14. Wang, Z.R.; Zhu, X.L.; Jiang, L.; Zhang, K.; Luo, H.W.; Gu, Y.; Zhang, P. Investigations of the hydrogen bonds and vibrational spectra of clathrate ice XVI. Materials 2019, 12, 246. [Google Scholar] [CrossRef]
  15. Zhu, X.L.; Yuan, Z.Y.; Jiang, L.; Zhang, K.; Wang, Z.R.; Luo, H.W.; Gu, Y.; Cao, J.W.; Qin, X.L.; Zhang, P. Computational analysis of vibrational spectrum and hydrogen bonds of ice XVII. New J. Phys. 2019, 21, 043054. [Google Scholar] [CrossRef]
  16. Clark, S.J.; Segall, M.D.; Pickard, C.J.; Hasnip, P.J.; Probert, M.I.J.; Refson, K.; Payne, M.C. First principles methods using CASTEP. Z. Kristallogr. 2005, 220, 567–570. [Google Scholar] [CrossRef] [Green Version]
  17. 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. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Computational spectra of Raman, IR, and phonon density of states (PDOS) of ice XV. Four vibrational bands: (a) intermolecular translation band; (b) intermolecular libration band; (c) intramolecular bending band; and (d) intramolecular stretching band.
Figure 1. Computational spectra of Raman, IR, and phonon density of states (PDOS) of ice XV. Four vibrational bands: (a) intermolecular translation band; (b) intermolecular libration band; (c) intramolecular bending band; and (d) intramolecular stretching band.
Molecules 24 03115 g001
Figure 2. Top view of six normal modes in translational band (107, 206, and 285 cm−1) and libration band (477, 913, and 936 cm−1). Typical molecular vibrations are shown in gold.
Figure 2. Top view of six normal modes in translational band (107, 206, and 285 cm−1) and libration band (477, 913, and 936 cm−1). Typical molecular vibrations are shown in gold.
Molecules 24 03115 g002
Figure 3. Fitted curves of four-bond modes (blue), two-bond modes (red), and non-H-bond modes (green) of ice XV in translational band. Inset PDOS curve is shown for comparison.
Figure 3. Fitted curves of four-bond modes (blue), two-bond modes (red), and non-H-bond modes (green) of ice XV in translational band. Inset PDOS curve is shown for comparison.
Molecules 24 03115 g003
Figure 4. Top view of six normal modes in bending band (1609, 1656, and 1709 cm−1) and stretching band (3278, 3366, and 3461 cm−1), respectively. Mode at 3366 cm−1 is a typical isolated stretching vibration.
Figure 4. Top view of six normal modes in bending band (1609, 1656, and 1709 cm−1) and stretching band (3278, 3366, and 3461 cm−1), respectively. Mode at 3366 cm−1 is a typical isolated stretching vibration.
Molecules 24 03115 g004
Table 1. Comparison of calculated results with Raman data. The normal modes (0.9 GPa) with IR/Raman intensities are compared with experimental Raman peaks. N means Raman inactive. PDOS, phonon density of states.
Table 1. Comparison of calculated results with Raman data. The normal modes (0.9 GPa) with IR/Raman intensities are compared with experimental Raman peaks. N means Raman inactive. PDOS, phonon density of states.
PDOSNormal ModesIR IntensityRaman IntensityRaman Exp. [8/11]
73N0.12
106N0.24
107N0.15
113136N0.22114/113
1381.65N
130148N0.72
1501.14N
1581.60N
162N0.27
165N1.20
1671.39N
1750.92N158/159
1591810.24N
182N0.57
184N0.40
1781851.22N
192N1.55
2010.22N186/185
206N0.33
221N0.32220/215
2172250.37N
2310.33N
232N0.04230/
2340.13N
243236N0.30300/
2682750.64N
285N1.01312/303
477N1.81451/450
4894841.09N
5065080.80N
529531N7.56499/495
5509.21N
566N3.06548/544
5792.93N
606611N1.43
6184.52N
621N2.42
63163110.29N
635N3.31
6491.36N
65216.64N
657662N4.75
67115.25N
679677N1.89
6858.19N
702N3.94
7227185.70N
731N1.79
7477485.20N
768775N3.73
815817N2.33
86411.87N
867N1.51
894892N1.08
9132.73N
9269283.25N
936N1.97
162116095.37N
16301629N1.94
16491648N3.68
1655N5.27
16565.46N
16603.70N
16845.57N
16921692N2.71
17052.37N
17061709N2.91
327863.85N
32963280N3777.673224/3222
332557.33N
33183326N552.403335/
333833397.64N
3341N346.173345/3349
335457.53N
33603357N288.67
3366N278.88
336750.96N
3391340227.87N
34113403N164.023396/
342927.91N
34323432N394.533425/3420
3442N154.45
344839.83N
345334529.24N
3461N100.59
34693472N383.12
348512.40N

Share and Cite

MDPI and ACS Style

Qin, X.-L.; Zhu, X.-L.; Cao, J.-W.; Jiang, L.; Gu, Y.; Wang, X.-C.; Zhang, P. Computational Analysis of Exotic Molecular and Atomic Vibrations in Ice XV. Molecules 2019, 24, 3115. https://doi.org/10.3390/molecules24173115

AMA Style

Qin X-L, Zhu X-L, Cao J-W, Jiang L, Gu Y, Wang X-C, Zhang P. Computational Analysis of Exotic Molecular and Atomic Vibrations in Ice XV. Molecules. 2019; 24(17):3115. https://doi.org/10.3390/molecules24173115

Chicago/Turabian Style

Qin, Xiao-Ling, Xu-Liang Zhu, Jing-Wen Cao, Lu Jiang, Yue Gu, Xue-Chun Wang, and Peng Zhang. 2019. "Computational Analysis of Exotic Molecular and Atomic Vibrations in Ice XV" Molecules 24, no. 17: 3115. https://doi.org/10.3390/molecules24173115

APA Style

Qin, X. -L., Zhu, X. -L., Cao, J. -W., Jiang, L., Gu, Y., Wang, X. -C., & Zhang, P. (2019). Computational Analysis of Exotic Molecular and Atomic Vibrations in Ice XV. Molecules, 24(17), 3115. https://doi.org/10.3390/molecules24173115

Article Metrics

Back to TopTop