24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1638  |  回复: 7

luojin7653

银虫 (正式写手)

[求助] 请教一个分子优化的问题 四个判据

为了获得分子的相关参数(包括,键长,键角,分子能量,分子轨道能量),
加freq 进行计算。

发现高斯程序正常终止,计算结果中都不存在虚频。
能说明是极小值而非鞍点吗?

但去查看输出文件中的四个判据,有的并不全为yes。

请问下,这样不全为yes的情况对我的计算得到相关参数影响会很大吗?需要重新优化为全yes不?


新手请教大家,很多问题不明白 谢谢了。
回复此楼
今晚打老虎
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)

我也遇到这种情况   不过你可以试着把初始结构改变以下  看看优化出来的一样不
2楼2012-01-06 07:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onlyyou2002

铁虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
luojin7653(金币+3): 有帮助 2012-01-06 09:48:59
luojin7653(金币+2): 2012-01-06 09:49:37
如果四个盘踞不全是yes,一般情况下程序不应该正常终止。优化的结果也不应该是global minum。你是不是优化的次数设少了?试着加入下面这个命令再重新再算一次
opt=(maxcycle=200)
人生如梦
3楼2012-01-06 08:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luojin7653

银虫 (正式写手)

引用回帖:
3楼: Originally posted by onlyyou2002 at 2012-01-06 08:46:39:
如果四个盘踞不全是yes,一般情况下程序不应该正常终止。优化的结果也不应该是global minum。你是不是优化的次数设少了?试着加入下面这个命令再重新再算一次
opt=(maxcycle=200)

嗯 好的
我的直接用的就是 opt freq
我再试试
今晚打老虎
4楼2012-01-06 09:49:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanan027

木虫 (著名写手)

★ ★
gmy1990(金币+2): 2012-01-06 21:55:03
是不是正常收敛你看一下输出里有没有Normal Termination就知道了啊。
如果前面两项值小于判据的百分之一的话它也可以收敛的(不管后两项达不达到收敛标准),结果可以接受。
Timeandtidewaitfornoman.Comeon,baby!
5楼2012-01-06 14:00:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luojin7653

银虫 (正式写手)

引用回帖:
5楼: Originally posted by chuanan027 at 2012-01-06 14:00:21:
是不是正常收敛你看一下输出里有没有Normal Termination就知道了啊。
如果前面两项值小于判据的百分之一的话它也可以收敛的(不管后两项达不达到收敛标准),结果可以接受。

Item               Value     Threshold  Converged?
Maximum Force            0.000093     0.000450     YES
RMS     Force            0.000023     0.000300     YES
Maximum Displacement     0.003988     0.001800     NO
RMS     Displacement     0.000871     0.001200     YES
Predicted change in Energy=-1.317597D-07
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad


Test job not archived.

n您看我这个结果合适吗?后面有显示的是 test job not archived.


另外的结果显示optimization completed 但是同样未收敛
         Item               Value     Threshold  Converged?
Maximum Force            0.000004     0.000450     YES
RMS     Force            0.000001     0.000300     YES
Maximum Displacement     0.005318     0.001800     NO
RMS     Displacement     0.001233     0.001200     NO
Predicted change in Energy=-2.737867D-08
Optimization completed on the basis of negligible forces.
    -- Stationary point found.

另外我加了maxcycle=200后,结果计算结果是一模一样,连4个判据的 value也一样
今晚打老虎
6楼2012-01-06 14:45:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanan027

木虫 (著名写手)

引用回帖:
6楼: Originally posted by luojin7653 at 2012-01-06 14:45:18:
Item               Value     Threshold  Converged?
Maximum Force            0.000093     0.000450     YES
RMS     Force            0.000023     0.000300     YES
Maximum Displacement     0.00 ...

照下面那个收敛情况就可以了,具体在哪份资料上有的我记不太清楚了,就会说前面两项值达收敛标准的百分之一也认为优化结束,因为对某些体系后面两项可能很难达到收敛标准。你上面那个收敛情况应该是算频率时算的吧,一般我都不看那里的。看正常结束,优化中间体的话看频率没有虚频就可以了。
Timeandtidewaitfornoman.Comeon,baby!
7楼2012-01-06 22:32:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onlyyou2002

铁虫 (初入文坛)

【答案】应助回帖

luojin7653(金币+10): ★★★很有帮助 谢谢 我去试试, 另外我并不是需要计算频率. 只是计算电子亲和能 2012-01-07 16:23:36
luojin7653(金币+5): 2012-01-08 16:30:40
你可以再加上这个试试
opt=(tight, maxcycle=200)
这样的收敛要求更严格
人生如梦
8楼2012-01-06 23:55:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luojin7653 的主题更新
信息提示
请填处理意见