基底四层原子，固定地层上面放开，发现优化后上层原子完全变形

lastzealot

patience2014 发表于 2015-5-11 17:47
补充一句，楼主所用的表面不是四层吧。

是四层，从上到下依次是 O Si O Si

卡开发发

lastzealot 发表于 2015-5-11 19:28
你说的dft-d2是泛函吗？

D2是色散校正，DMol3似乎是6.0还是6.1版本以后可以做GGA+Grimme06（即D2），版本低的话就是没法做-D而已（测试下来D2比TS靠谱，站内有相关的帖子聊过这个问题）。

g09算的话B3LYP算其他的体系的我也见过，有些处理就是固定了边缘的原子。DMol3没对周期体系使用B3LYP，如果文献只用了B3LYP或PBE、PW91而不考虑色散校正也被认可的话投这类期刊审稿人不问马马虎虎也就将就那样了，从科研本身角度来说，这样的体系还是加上-D靠谱。

如果硬要是杂化泛函来处理的话应该CP2K值得尝试，但速度能多快也不好说。数值基组的方法的程序支持杂化泛函并且公开的并不多见。你这个尺寸的基底MP2基本不用考虑

sobereva

lastzealot 发表于 2015-5-11 19:48
是四层，从上到下依次是 O Si O Si

对于SiO2，三四层不是指有三四层原子，而是说SiO2有三四层，例如此图（TiO2的例子）

你的基底太薄，底层原子完全体现不出体相状态。

lastzealot

卡开发发 发表于 2015-5-11 21:22
D2是色散校正，DMol3似乎是6.0还是6.1版本以后可以做GGA+Grimme06（即D2），版本低的话就是没法做-D而已 ...

非常感谢老师耐心指导。
我用高斯试一试您说的cp2k泛函，并且固定边界原子，看看能否得到较理想的结果。

lastzealot

sobereva 发表于 2015-5-11 22:29
对于SiO2，三四层不是指有三四层原子，而是说SiO2有三四层，例如此图（TiO2的例子）

你的基底太薄， ...

原来如此，看来我的基底只有一层半二氧化硅。

卡开发发

lastzealot 发表于 2015-5-11 22:32
非常感谢老师耐心指导。
我用高斯试一试您说的cp2k泛函，并且固定边界原子，看看能否得到较理想的结果。

cp2k是个程序，一种利用平面波与原子轨道相互变换的方法。色散校正应当被考虑。

xjamesmao

sobereva 发表于 2015-5-11 22:29
对于SiO2，三四层不是指有三四层原子，而是说SiO2有三四层，例如此图（TiO2的例子）

你的基底太薄， ...

正好请教一下各位方家。我最近在做一个类似的体系，SiO2表面和一个小分子（类苯环）的吸附，需要优化结构和计算它们之间的作用力。对于SiO2部分，我建了一个六层Si的结构，差不多有160个原子（35个Si，89个O，剩下的是H）。我的计划是固定SiO2部分的Si和O的坐标，优化其余的原子，然后计算作用力。对SiO2非顶层的原子打算用小一些的基组。哪位有经验的，能不能建议一下，顶层和非顶层都分别用什么基组比较可行？或者哪些类似的文献可以参考？
谢谢大家了！

sobereva

xjamesmao 发表于 2015-5-13 11:58
正好请教一下各位方家。我最近在做一个类似的体系，SiO2表面和一个小分子（类苯环）的吸附，需要优化结构 ...

对于底层SiO2最便宜能用的是3-21G。顶层的就用6-31G*，不算大，精度也没问题。

nkallwar

sobereva 发表于 2015-5-13 12:09
对于底层SiO2最便宜能用的是3-21G。顶层的就用6-31G*，不算大，精度也没问题。

Sob老师，PM6这种半经验法和B3lyp/3-21G 哪个更便宜些？

看到有些帖子说PM6在有些体系结果比低级别的DFT还好些？

lastzealot

sobereva 发表于 2015-5-13 12:09
对于底层SiO2最便宜能用的是3-21G。顶层的就用6-31G*，不算大，精度也没问题。

请问sob老师，同一种元素、不同的原子可以使用不同的基组吗？

sobereva

lastzealot 发表于 2015-5-13 13:12
请问sob老师，同一种元素、不同的原子可以使用不同的基组吗？

可以

sobereva

nkallwar 发表于 2015-5-13 12:36
Sob老师，PM6这种半经验法和B3lyp/3-21G 哪个更便宜些？

看到有些帖子说PM6在有些体系结果比低级别 ...

PM6比B3LYP/3-21G便宜一两个数量级。
你可以用PM6做初步优化。半经验方法有时候结果比DFT结合低级别基组好，但是稳妥性没那么高。