IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. The performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases is evaluated. Based on a selected DFT computational scheme, the IR and Raman spectra of ice Ih and XI are derived and compared. When the spectra, both IR and Raman, of ice Ih and ice XI are compared, the librational vibrations are found to be the most affected by the proton ordering. The spectroscopic fingerprint of ice XI can be used to distinguish ferroelectric ice XI from ice Ih in the universe. Furthermore, the existence of only one kind of H-bond in ice Ih is demonstrated from the overlapping subspectra for different types of H-bonded pair configurations in 16 isomers of ice Ih, which provides an illustration to the historic debate on whether one or two kinds of H-bonds existed in ice.
Funding agencies: Swedish Research Council (VR); Swedish Supercomputer Center (SNIC/NSC); State Scholarship Fund of China Scholarship Council 
In this thesis H-bonded systems (natural gas hydrates, water clusters, and crystal ice) are studied by density functional theory (DFT) computations.
Natural gas hydrates (NGHs) play an important role in energy and environmental fields: NGHs are considered as a promising backup energy resource in the near-future due to their tremendous carbon content; improper exploration of NGHs could induce geological disasters and aggravate the greenhouse effect. In addition, many technologies based on gas hydrates are being applied and developed. The thermodynamic stabilities of various water cavities in different clathrate crystalline phases occupied by hydrocarbon gas molecules are studied by dispersion-corrected hybrid functionals. The Raman spectra of C-C and C-H stretching vibrations of hydrocarbon molecules in various water cavities in the solid state are derived. The trends of C-H stretching vibrational frequencies are found to follow the prediction by the “loose cage ─| tight cage” model. In addition, the trends and origins of 13C NMR chemical shifts of hydrocarbon molecules in various NGHs are presented. These theoretical results will enlarge the database of C-C and C-H stretching vibrational frequencies and 13C NMR parameters of hydrocarbon molecules in NGHs, and provide valuable information to help identify the types of clathrate phases and varieties of guest molecules included in NGHs samples taken from natural sites.
The behavior of water clusters may help to understand the properties of its liquid and solid states. The thermodynamic stabilities and IR spectra of a small-, medium-, and large-sized water cluster are studied in this work. After full optimization of (H2O)20,54,100 using the hybrid functional B3LYP, the electronic energies, zero-point energies, internal energies, enthalpies, entropies, and Gibbs free energies of the water clusters are computed. The OH stretching vibrational IR spectra of (H2O)20,54,100 are also presented and split into sub-spectra for different H-bond types based on the specific contributions from each group. It is found that the OH stretching vibrational frequencies of water are sensitive to the conformations of the H-bonds and the vibrations of the H-bonds belonging to different types are located in separated regions in the IR spectra. Thus, the spectroscopic fingerprints will reflect the H-bond topology of the water molecules in a water cluster.
Ice XI has been suggested to be involved in the process of planetary formation as a considerable electric field might be formed from the ferroelectric ice XI in space. IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. Due to the difficulty for DFT to describe non-covalent systems, the performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases are assessed. Based on the computational accuracy and cost, the IR and Raman spectra of ice Ih and XI are derived and compared. The librational vibrations are found to be the identifier which can be used to distinguish ice Ih and ice XI in the universe. In addition, the existence only one kind of H-bond in ice Ih is demonstrated from the overlapping sub-spectra for different types of H-bonded pair configurations in 16 isomers of ice Ih.
The region of water under negative pressure is an exotic land in lack of exploitation. Guest free clathrate hydrate (clathrate ice) of sII type has been recently confirmed experimentally at negative pressure. Does any other clathrate ice phase exist at negative pressure region? Since clathrate hydrate are isostructural with silica clathrate minerals and semiconductor clathrates, and crystal structure prediction by analogy with known structures and first-principles computations is an effective way to find new crystalline phases of solid materials, we are motived to look for new clathrate ice phases from silica or semiconductor clathrate materials based on first-principles computations. Borrowing the idea new clathrate frameworks of ZnO and SiC can be constructed by connecting their bubble clusters in different ways, new clathrate ice phases (sL, sL_I, sL_II, and sL_III) are generated by connecting the water bubble clusters according to different rules. Using the non-local dispersion-corrected vdW-DF2 functional, clathrate ice sL with ultralow density (0.6 g/cm3) is predicted by first-principles phase diagram computations to be stable under larger negative pressures than the sII phase. The phase diagram of water is thus extended into the lower negative pressure region.