计算化学公社

标题: 矿物进行红外和拉曼光谱计算和一般有机物的计算有差异吗 [打印本页]

作者
Author: shuihua 时间: 2022-3-3 09:30
标题: 矿物进行红外和拉曼光谱计算和一般有机物的计算有差异吗
本帖最后由 shuihua 于 2022-3-3 09:33 编辑

我之前计算过油酸等一些有机物的光谱，gjf如下所示：
%chk=3.chk
#p opt freq=(raman,savenormalmodes,intmode) b3lyp/6-31g(d)

Title Card Required

0 1
C       4.79210002092103 -2.72647800385196 1.34080247157848
C       4.00029414385384 -1.77749960664400 0.44112735332560
C       4.03020893803100 -2.22803114657871 -0.98640308069323
C       4.39464703668700 -1.53410787705252 -2.05410038240501
C       4.90353101724918 -0.12247411805305 -2.10464510088786

………………
我在AMCSD找了一些矿物的晶格结构cif文件，请问是不是和有机物类似，选择合适的泛函和基组就可以了呢？矿物主要是一些硅酸盐，Mg2SiO4等
我没有别的计算需求，只需要知道他们的红外和拉曼光谱即可
Thanks in advance

作者
Author: shuihua 时间: 2022-3-3 09:45
http://bbs.keinsci.com/forum.php ... hlight=%BB%B7%BE%B3
是不是可以参照SOB发的这个文章呢？

作者
Author: zjxitcc 时间: 2022-3-3 09:50
本帖最后由 zjxitcc 于 2022-3-3 09:53 编辑

非也。如果你研究的是周期性体系，且它们不适合截断为孤立体系（金属氧化物、离子晶体等）时，那么你应该采用周期性体系的计算程序，例如CP2K，Quantum Espresso，Crystal等等。

实际上一般只有有机分子的分子晶体能截断成孤立体系算，大多数晶体是不适合的。强行截断的话，（即使你对边界做了合理的处理）你会观察到SCF巨难收敛，而且算出来结果也未必可信。

如果仅含Si, O, H，这都是共价单键，有可能通过截断和边界加氢饱和，变成一个孤立体系算，结果可能还合理。而再加上Mg的话，还当成孤立体系就不太合适了。

作者
Author: shuihua 时间: 2022-3-3 10:00

zjxitcc 发表于 2022-3-3 09:50
非也。如果你研究的是周期性体系，且它们不适合截断为孤立体系（金属氧化物、离子晶体等）时，那么你应该采 ...

好的，我明白了。我这就去看一看CP2K这些程序的使用，感谢您！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)