应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】新手求助高斯几何优化的问题
141/2返回列表12下一页
查看: 741  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

huabangou

[交流] 【求助】新手求助高斯几何优化的问题

本人是新手，对磷酸镧进行几何优化，建了很多结构都不成功，主要是不收敛，能量没有出现震荡，希望高手给予指点．输入文件如下：
%chk=24conventional cell.chk
%mem=300MW
%nproc=1

# B3LYP/LanL2DZ opt test

Title Card Required
０１
La                 0.00000000    1.49648495    2.12686250
La                 0.00000000    6.46932505    6.38058750
La                 3.93761000    5.47938995    2.12686250
La                 3.93761000    2.48642005    6.38058750
P                  5.28300573    0.00000000    4.25372500
P                  5.28300573    7.96581000    4.25372500
P                  6.52982427    3.98290500    8.50745000
P                  6.52982427    3.98290500    0.00000000
O                  0.47592026    2.75509384    4.03856895
O                  0.47592026    5.21071616    4.46888105
O                  2.21588310    4.18761771    3.01245189
O                  2.21588310    3.77819229    5.49499811
O                  5.65933690    4.18761771    1.24127311
O                  5.65933690    3.77819229    7.26617689
O                  3.46168974    1.22781116    8.29229395
O                  3.46168974    6.73799884    0.21515605
O                  1.72172690    0.20471271    1.24127311
O                  1.72172690    7.76109729    7.26617689
O                  4.41353026    1.22781116    4.46888105
O                  4.41353026    6.73799884    4.03856895
O                  7.39929974    5.21071616    8.29229395
O                  7.39929974    2.75509384    0.21515605
O                  6.15349310    7.76109729    5.49499811
O                  6.15349310    0.20471271    3.01245189
回复此楼

» 猜你喜欢
1楼2009-06-17 15:11:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

S07111072

铁杆木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+2,VIP+0):谢谢交流 6-17 18:03
楼主优化时是怎么不收敛啊，如果构型没有问题，可以改变波函数的初始猜测，如Guess=indo
赞一下(7人) 回复此楼
2楼2009-06-17 16:20:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+4,VIP+0):谢谢交流 6-17 18:03
首先你给出的输入文件格式有问题，不知道你在计算中是否就是用这样的文件
其次你搭的结构也有问题，哪有那么多氧原子堆在一块儿的

另外你给的信息太少了，
计算出错了吗？
如果出错了，贴出输出文件最后的出错信息
如果没有出错，需要贴出不收敛的信息，
在你的输出文件中查找“Threshhold” 字段，找到最后一个
将其下面的四行贴出来，
顺便看看最后几个“Threshhold” 字段下面的情况

看你的chk文件的名字，好像还用了cell
赞一下(7人) 回复此楼
据权威报纸调查，84%的青少年犯罪者有玩电子游戏的经历，故应该禁止经营电子游戏。据我们所知，100%的青少年犯罪者有穿鞋的经历，故制鞋厂应该停业。
3楼2009-06-17 16:59:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
Guess=indo

有什么用啊?
赞一下(1人) 回复此楼
4楼2009-06-17 21:02:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

orangedp6597

银虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+2,VIP+0):谢谢交流 6-18 07:56
你的分子中有金属，应该用使用赝势的输入文件，具体参照高斯说明。
赞一下(7人) 回复此楼
5楼2009-06-18 01:49:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+2,VIP+0):欢迎交流 6-18 21:25
引用回帖:
Originally posted by orangedp6597 at 2009-6-18 01:49:
你的分子中有金属，应该用使用赝势的输入文件，具体参照高斯说明。

他用的就是赝势啊，只不过是全赝势基组。他的输入格式是对的
赞一下(11人) 回复此楼
6楼2009-06-18 13:07:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★
yuhuobuku(金币+3,VIP+0):欢迎交流 6-18 21:26
引用回帖:
Originally posted by erylingjet at 2009-6-17 21:02:
Guess=indo
有什么用啊?

Guess=indo的意思是：使用Gaussian98的默认初始猜测：对第一行元素是INDO，第二行元素是CNDO，第三行及以后的元素是Huckel。
如果不写Guess=indo，则使用Gaussian03的默认初始猜测：全部使用Huckel初始猜测。
Guess=indo的用处：使比较差的初始构型，能比较顺利地通过初始猜测。
想不到楼主的输入文件写的一塌糊涂，却居然知道Guess=indo！

[ Last edited by yjcmwgk on 2009-6-18 at 13:11 ]
赞一下(10人) 回复此楼
7楼2009-06-18 13:10:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+2,VIP+0):欢迎交流 6-18 21:27
引用回帖:
Originally posted by huabangou at 2009-6-17 15:11:
本人是新手，对磷酸镧进行几何优化，建了很多结构都不成功，主要是不收敛，能量没有出现震荡，希望高手给予指点．输入文件如下：
%chk=24conventional cell.chk
%mem=300MW
%nproc=1

# B3LYP/LanL2DZ opt t ...

你的计算根本上就行不通。LanL2DZ方法，处理La，误差相当大。对于稀土元素来说，LanL2DZ仅仅对Y还算有一定准确度（参看ChemPhysChem 2008, 9, 781-792和ChemEurJ2007，13，9503-9514）。但是对于La的准确度就太低了（参看JPCA，2007，111，392-400）建议采取SDD（如果你的资源比较多）或CEP-121G（如果你的资源比较少）
赞一下(11人) 回复此楼
8楼2009-06-18 13:30:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaochao85

铜虫 (正式写手)
★
yuhuobuku(金币+1,VIP+0):欢迎交流 6-18 21:27
SDD 确实比较准确
赞一下(10人) 回复此楼
9楼2009-06-18 13:59:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huabangou

非常感谢大家的回帖,我也知道问题,所以才拿出来让大家帮着看看,太感谢你们了
下面是最后三个Threshold的值,还请大家帮忙分析:
     Item                                Value     Threshold  Converged?
Maximum Force                0.093072     0.000450     NO
RMS     Force                    0.022360     0.000300     NO
Maximum Displacement    0.074187     0.001800     NO
RMS     Displacement        0.016424     0.001200     NO
Predicted change in Energy=-5.518068D-02

         Item                             Value     Threshold  Converged?
Maximum Force                 0.082378     0.000450     NO
RMS     Force                     0.020155     0.000300     NO
Maximum Displacement     0.010794     0.001800     NO
RMS     Displacement         0.002636     0.001200     NO
Predicted change in Energy=-2.249731D-05

         Item                              Value     Threshold  Converged?
Maximum Force                 0.084392     0.000450     NO
RMS     Force                     0.020264     0.000300     NO
Maximum Displacement     0.199430     0.001800     NO
RMS     Displacement         0.062802     0.001200     NO
Predicted change in Energy=-3.685631D-03

[ Last edited by huabangou on 2009-6-18 at 15:42 ]
赞一下 回复此楼
10楼2009-06-18 15:19:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huabangou 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭