计算化学公社

标题: 除了 B3LYP 外對振動光譜強度預測表現不錯的泛函 [打印本页]

作者
Author: sdd3390 时间: 2018-11-13 20:53
标题: 除了 B3LYP 外對振動光譜強度預測表現不錯的泛函
各位好

小弟最近正在做 surface enhenced raman 光譜的計算
基本想法就是把分子放在某一個基材的表面上(沒有任何共價鍵形成) 和沒放在基材表面上的振動光譜差異
分子是使用 ˊ6-31g(d) 基材是使用 3-21g* basis set

在計算振動光譜時  因為只想看基材對分子振動光譜的影響  所以把基材的所有原子都固定住 (做opt時也是如此)
目前的結果看來是 B3LYP 對某一個normal mode 的強度(intensity)有非常大的增強

但是 PBE1PBE1 的強度不變 M062X 強度反而減弱
所以想請問說還有哪一個泛函對振動光譜強度預測的表現是良好的?  因為想多進行測試以驗證計算結果

作者
Author: sobereva 时间: 2018-11-14 04:26
不要所有基底原子都固定，起码把与小分子接触的部分放开
M06-2X对于过渡金属描述不好，别用。PBE1PBE原理上可以
B3PW91虽然也可以用，但应该结果和B3LYP没显著区别。其它HF成分更低的诸如TPSSh也可以考虑

作者
Author: sdd3390 时间: 2018-11-16 20:57
感謝 sobereva 的回答
因為整個體系已經接近 250 個原子  加上分子和基材的接觸面積極大所以才在優化時將基材定住
不過用 DFT-D3 做優化前  有先用 PM7-D3 優化過結構
我本來的想法是就是看基材對 RAMAN 光譜的影響  不希望基材算出的光譜對分子的光譜有干擾
最一開始是想用 ONIOM 下去做  但是 ONIOM 好像很多人都不大推
整個體系都是有機分子  這樣使用低成分 HF 的 TPSSh OK嗎?

作者
Author: sobereva 时间: 2018-11-19 10:25

sdd3390 发表于 2018-11-16 20:57
感謝 sobereva 的回答
因為整個體系已經接近 250 個原子加上分子和基材的接觸面積極大所以才在優化時 ...

固定接触层并不能节约计算时间，反倒破坏物理意义
没有PM7-D3。PM7本身已经有了色散校正，形式用的不是DFT-D
ONIOM对这种问题严重不适合，此方法根本没法很好表现不同层之间的耦合，还会带来人为因素
TPSSh可以

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