请教sob大神一个问题

lastzealot · 发表于 Post on 2015-3-29 23:56:45

我想优化五水合铀酰 [UO2(H2O)5]2+的结构，其中铀酰是类似于CO2的直线型结构，五个水在中间以铀为中心的赤道平面对称排布（文献中的计算结果）。
首先，我将这种结构按照文献里所述做出原始结构，如图1所示，对称性非常好，每个O-U-O夹角72度，都在迟到平面上。

把关键词设置为#p b3lyp/gen pseudo=read opt=(cartesian, calcfc, maxcycle=200, gediis) scf=maxcycle=1000 symm=tight
但是优化过程每次都出错，多次L508出错，后来就是L9999出错，实在不知道怎么才能算完成。
在打开未完成的.out文件后，发现赤道平面的水分子已经不在赤道平面上，与文献中差距很大。我不知道文献中是如何把这些水分子限制在赤道平面上的，请问是用内坐标的方法固定某些角度为90度吗？还是其他方法？
还请sob大神为小弟解答下，最好说的详细一些。作为一个量化计算菜鸟，想您请教。

sobereva · 发表于 Post on 2015-3-30 05:14:48

你当前的关键词不应该出现508错误，那只会在使用scf=qc或xqc时才会出现。普通SCF计算不会涉及到508模块

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

lastzealot · 发表于 Post on 2015-3-30 08:58:54

恩，后来我加上了scf=qc，这一点忘说了。基组和文献不一样。文献中用小核赝势RSC，我电脑运行能力有限，用了大核赝势。但是我

lastzealot · 发表于 Post on 2015-3-30 08:59:50

后来就是多次出现l9999问题，不知道是怎么回事

sobereva · 发表于 Post on 2015-3-30 11:38:35

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

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 请教sob大神一个问题

浏览过的版块