矿物进行红外和拉曼光谱计算和一般有机物的计算有差异吗

shuihua

我之前计算过油酸等一些有机物的光谱，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

shuihua

http://bbs.keinsci.com/forum.php ... hlight=%BB%B7%BE%B3
是不是可以参照SOB发的这个文章呢？

zjxitcc

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

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

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

shuihua

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

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