Dmol3算铀掺杂石墨烯结构优化出现问题

卡开发发

我做了测试，取了6x6的表面（不知道什么原因，小于6x6超胞我测试了能带结构，无论k取多大给出的结果都有点问题），做了a=b=14.74~14.83A的测试，拟合抛物线内插得到晶格常数的能量极小点为a=b=14.793，真空层我取了30A。

在默认Fine精度的基准上修改了k、smearing，k为11x11x1，这个在先前的帖子中做了测试，smearing降低到0.001Ha(约0.03eV，其实也可以不用)。

截断半径取全局为5.9A，我做了测试，如果只有U取5.9，C取4.5，可能导致力不收敛，以往的经验来说也不推荐取不同的截断半径，有时候可能造成scf不收敛，有时候可能造成结构弛豫不收敛等问题。因为截断半径和k都很大，可能计算速度稍慢些，12核心服务器预计0.8~1.5m做一次scf。

另外，你的初始结构O原子指向比较糟糕。

我把以及优化好的U-graphite和O@U-graphite的结构做成cif（避免版本差异）放在这，你可以参考这个结构做一些后续处理。

lmwnbpes

卡开发发 发表于 2015-1-22 08:44
我做了测试，取了6x6的表面（不知道什么原因，小于6x6超胞我测试了能带结构，无论k取多大给出的结果都有点 ...

非常感谢您的详细解答，学到很多。
我想问一下  我尝试把Max.step size从0.3  调到0.1 ，就收敛了，这个怎么理解

鹃子

是说每次跳.3 正好错过了能量最低点吗? 猜的

卡开发发

lmwnbpes 发表于 2015-1-26 09:30
非常感谢您的详细解答，学到很多。
我想问一下  我尝试把Max.step size从0.3  调到0.1 ，就收敛了，这 ...

结构振荡比较严重的情况可以限制分配的下一步原子移动，这样每次计算能量和力之后分配出各个原子的移动即便大于0.1A也按照0.1A进行移动，有时候比较管用。不过一般不建议这个值设定的太小。