计算化学公社

标题: 请教sob大神一个问题 [打印本页]

作者
Author: lastzealot 时间: 2015-3-29 23:56
标题: 请教sob大神一个问题
我想优化五水合铀酰 [UO2(H2O)5]2+的结构，其中铀酰是类似于CO2的直线型结构，五个水在中间以铀为中心的赤道平面对称排布（文献中的计算结果）。
首先，我将这种结构按照文献里所述做出原始结构，如图1所示，对称性非常好，每个O-U-O夹角72度，都在迟到平面上。
(, 下载次数 Times of downloads: 41)
把关键词设置为#p b3lyp/gen pseudo=read opt=(cartesian, calcfc, maxcycle=200, gediis) scf=maxcycle=1000 symm=tight
但是优化过程每次都出错，多次L508出错，后来就是L9999出错，实在不知道怎么才能算完成。
在打开未完成的.out文件后，发现赤道平面的水分子已经不在赤道平面上，与文献中差距很大。我不知道文献中是如何把这些水分子限制在赤道平面上的，请问是用内坐标的方法固定某些角度为90度吗？还是其他方法？
还请sob大神为小弟解答下，最好说的详细一些。作为一个量化计算菜鸟，想您请教。

作者
Author: sobereva 时间: 2015-3-30 05:14
你当前的关键词不应该出现508错误，那只会在使用scf=qc或xqc时才会出现。普通SCF计算不会涉及到508模块

基组和文献用的一致否？从输出文件中看，程序一开始是否已经合理判断对了点群？

作者
Author: lastzealot 时间: 2015-3-30 08:58
恩，后来我加上了scf=qc，这一点忘说了。基组和文献不一样。文献中用小核赝势RSC，我电脑运行能力有限，用了大核赝势。但是我

作者
Author: lastzealot 时间: 2015-3-30 08:59
后来就是多次出现l9999问题，不知道是怎么回事

作者
Author: sobereva 时间: 2015-3-30 11:38
一般就是几何优化不收敛，具体看输出文件的信息。解决方法见《量子化学计算中帮助几何优化收敛的常用方法》
http://sobereva.com/164

通常不推荐用SCF=QC，如果是SCF不收敛，解决方法见《解决SCF不收敛问题的方法》
http://sobereva.com/61

至于优化后不在一个平面的问题，先看点群是否判断对了。

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