24小时热门版块排行榜    

查看: 1845  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 267一志愿南京工业大学0817化工求调剂 +8 SUICHILD 2026-03-12 8/400 2026-03-18 20:55 by winsuccess
[考研] 0703化学调剂 +7 妮妮ninicgb 2026-03-15 11/550 2026-03-18 19:48 by macy2011
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 0703化学调剂 ,六级已过,有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
信息提示
请填处理意见