Generalized Energy-Based Fragmentation Approach for Accurate Binding Energies and Raman Spectra of Methane Hydrate Clusters

Methane hydrates(MHs)play important roles in the fields of chemistry,energy,environmental sciences,etc.In this work,we employ the generalized energy-based fragmentation(GEBF)approach to compute the binding energies and Raman spectra of various MH clusters.For the GEBF binding energies of various MH clusters,we first evaluated the various functionals of density functional theory(DFT),and compared them with the results of explicitly correlated combined coupled-cluster singles and doubles with noniterative triples corrections[CCSD(T)(F12^(*))]method.Our results show that the two best functionals are B3PW91-D3 and B97D,with mean absolute errors of only 0.27 and 0.47 kcal/mol,respectively.Then we employed GEBF-B3PW91-D3 to obtain the structures and Raman spectra of MH clusters with mono-and double-cages.Our results show that the B3PW91-D3 functional can well reproduce the experimental C-H stretching Raman spectra of methane in MH crystals,with errors less than 3 cm^(-1).As the size of the water cages increased,the C-H stretching Raman spectra exhibited a redshift,which is also in agreement with the experimental“loose cage-tight cage”model.In addition,the Raman spectra are only slightly affected by the neighboring environment(cages)of methane.The blueshifts of C-H stretching frequencies are no larger than 3 cm^(-1) for CH_(4) from monocages to doublecages.The Raman spectra of the MH clusters could be combined with the experimental Raman spectra to investigate the structures of methane hydrates in the ocean bottom or in the interior of interstellar icy bodies.Based on the B3PW91-D3 or B97D functional and machine learning models,molecular dynamics simulations could be applied to the nucleation and growth mechanisms,and the phase transitions of methane hydrates.

The relative energies of polypeptide conformers predicted by linear scaling second-order M?ller-Plesset perturbation theory

We describe an implementation of the cluster-in-molecule （CIM） resolution of the identity （RI） approximation second-order Moller-Plesset perturbation theory （CIM-RI-MP2）, with the purpose of extending RI-MP2 calculations to very large systems. For typical conformers of several large polypeptides, we calculated their conformational energy differences with the CIM-RI-MP2 and the generalized energy-based fragmentation MP2 （GEBF-MP2） methods, and compared these results with the density functional theory （DFT） results obtained with several popular functionals. Our calculations show that the conformational energy differences obtained with CIM-RI-MP2 and GEBF-MP2 are very close to each other. In comparison with the GEBF-MP2 and CIM-RI-MP2 relative energies, we found that the DFT functionals （CAM-B3LYP-D3, LC-ωPBE-D3, M05-2X, M06-2X and coB97XD） can give quite accurate conformational energy differences for structurally similar conformers, but provide less-accurate results for structurally very different conformers.