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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】Displacement优化不收敛怎如何处理?
71/1返回列表
查看: 998  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

loovfnd

至尊木虫 (著名写手)

[交流] 【求助】Displacement优化不收敛怎如何处理?

后二个指标早己收敛,而前面二个总是在变化, 很长时间都没有收敛的迹象, 请高手指点,拜托各位!


Maximum Force            0.017061     0.002500              NO
RMS         Force            0.003698     0.001667              NO
Maximum Displacement     0.004802     0.010000         YES
RMS         Displacement     0.001647     0.006667         YES
回复此楼

» 猜你喜欢
莫强求！
1楼 2009-11-30 20:14:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才
★ ★ ★
loovfnd(金币+3,VIP+0):xx 12-1 09:10
引用回帖:
Originally posted by loovfnd at 2009-11-30 20:14:
后二个指标早己收敛,而前面二个总是在变化, 很长时间都没有收敛的迹象, 请高手指点,拜托各位!


Maximum Force            0.017061     0.002500              NO
RMS         Force          ...

前两个在变化，那就是你的分子构型进入了一个不稳定区，一般这个区也没有最优化构型，  你可以iop（1/8=n)   n>3 试试  
跳到其他亚稳区域 进行计算
赞一下(2人) 回复此楼
我想做个科学家，可我又喜欢写诗；我想做一个诗人，可我却读的理工科..
2楼2009-11-30 21:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)
引用回帖:
Originally posted by jiewei at 2009-11-30 21:08:

前两个在变化，那就是你的分子构型进入了一个不稳定区，一般这个区也没有最优化构型，  你可以iop（1/8=n)   n>3 试试  
跳到其他亚稳区域 进行计算

我用过iop（1/8=1)，这里为什么要大于3？

[ Last edited by yjcmwgk on 2009-12-1 at 22:02 ]
赞一下 回复此楼
莫强求！
3楼2009-11-30 22:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才
★ ★ ★ ★ ★ ★
yjcmwgk(金币+3,VIP+0): 12-1 22:02
loovfnd(金币+3,VIP+0):xiexie 12-7 10:52
引用回帖:
Originally posted by loovfnd at 2009-11-30 22:49:
我用过iop（1/8=1)，这里为什么要大于3？  

如果n=1  那是步长小，那分子的优化就是在不稳定周围的一个小区间优化。
而显然你现在分子处于的这个区间是没有最稳定点的。
所以你要跳出去啊，就要加大步长 。

[ Last edited by yjcmwgk on 2009-12-9 at 09:44 ]
赞一下(11人) 回复此楼
我想做个科学家，可我又喜欢写诗；我想做一个诗人，可我却读的理工科..
4楼2009-12-01 10:12:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)
★ ★ ★ ★
loovfnd(金币+4,VIP+0):谢谢您的帮助 12-9 10:17
有可能是点群设置过高
opt=gdiis也有帮助
赞一下(1人) 回复此楼
5楼2009-12-09 09:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)
引用回帖:
Originally posted by Jasminer at 2009-12-9 09:52:
有可能是点群设置过高
opt=gdiis也有帮助

我使用的是D5h对称性坐标，我想要的是D2h的，使用symm=pg=D2h不哈用。

有什么方法可以直接改成D2h的，gv里我试了好像不行。
赞一下 回复此楼
莫强求！
6楼2009-12-09 10:16:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 12-9 20:18
只能手动改，很麻烦。
几个替代方法：
1. nosymm优化几步，对称性就被破坏了，再用gv收成D2h点群；
2. 算一个频率，沿相应的震动方式改变结构。慢些，但优点是如果你的displacement收敛得很好的话，一般可以得到一个破坏高点群的虚频，直接确认实际点群。有时算出来的是实频，但不能说明高点群是正确的。
赞一下(11人) 回复此楼
7楼2009-12-09 19:44:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 loovfnd 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭