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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 【讨论】菜鸟学GAMESS的心得
172/2返回列表上一页12
查看: 1393  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

loovfnd

至尊木虫 (著名写手)
结果怎么样了啊？
回复此楼
莫强求！
11楼2009-03-10 20:38:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peptide2882

木虫 (正式写手)
★ ★ ★
yuhuobuku(金币+3):欢迎参加讨论
我还没开始做过渡态呢，连顺式都没开始算，一直在优化反式，现在主要的问题是不收敛，看了上面老虎大王的建议，有那么点开窍了，我前面只知道瞎改优化参数，而不是有的放矢，老虎大王说：“收敛的标准好像是在$STATPT 中的OPTTOL变量”，我对着说明看了看，知道默认的收敛标准是0.0001（说明原文：OPTTOL = gradient convergence tolerance, in Hartree/Bohr.
      Convergence of a geometry search requires the largest component of the gradient to be less than OPTTOL, and the root mean square gradient
less than 1/3 of OPTTOL. (default=0.0001)）
而我的计算中，每次跑完50步，都基本从零点几接近0.0002了，因此我把RMS最低的一组的坐标重新拷回输入文件，又跑了50步，终于在其中发现能量最小的那组的RMS低于0.0001了，是0.000083，可是结果还显示没收敛：  
MAXIMUM GRADIENT =  0.0002357    RMS GRADIENT = 0.0000831
          HESSIAN UPDATED USING THE BFGS FORMULA
          ACTUAL ENERGY CHANGE WAS  -0.0000086726
          PREDICTED ENERGY CHANGE WAS  -0.0000068772 RATIO=  1.261
          MIN SEARCH, CORRECT HESSIAN, TRYING PURE NR STEP
          NR STEP HAS LENGTH         =   0.125972
          RADIUS OF STEP TAKEN=   0.12597  CURRENT TRUST RADIUS=   0.15972
..
..
  ***** FAILURE TO LOCATE STATIONARY POINT, TOO MANY STEPS TAKEN *****
     UPDATED HESSIAN, GEOMETRY, AND VECTORS WILL BE PUNCHED FOR RESTART
**** THE GEOMETRY SEARCH IS NOT CONVERGED! ****

  是不是还得比0.0001低到一定程度才认为收敛呢，貌似不能只看RMS，还有gradient也得符合标准，我又把最新的坐标进一步优化试，如果在不行，可能得考虑这个收敛是不是和别的因素有关。

[ Last edited by peptide2882 on 2009-3-10 at 22:08 ]
赞一下(1人) 回复此楼
终于变成木虫了
12楼2009-03-10 21:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peptide2882

木虫 (正式写手)
泪奔啊，出结果了：
MAXIMUM GRADIENT =  0.0000414    RMS GRADIENT = 0.0000109

      ***** EQUILIBRIUM GEOMETRY LOCATED *****
看来真得把每次最优的坐标当成下一步的输入文件坐标阿（或者可以直接加长搜索步数，我每次都限定50步），终于看到这行单词了，下面的工作是顺式和更难的过渡态了。明天出去找找工作，回来再继续算。
赞一下 回复此楼
终于变成木虫了
13楼2009-03-10 23:46:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
yuhuobuku(金币+3):欢迎参加讨论
呵呵。你要注意说明书里的话：
OPTTOL = gradient convergence tolerance, in Hartree/Bohr.
      Convergence of a geometry search requires the largest component of the gradient to be less than OPTTOL, and the root mean square gradient
less than 1/3 of OPTTOL. (default=0.0001)）

这个话的意思是说，最大能量梯度要小于OPTTOL，而均方根梯度(root mean square gradient, 即RMS gradient)要小于OPTTOL的1/3。注意这里有两个条件。

你在12楼中贴出来的
MAXIMUM GRADIENT =  0.0002357    RMS GRADIENT = 0.0000831
两个条件都不满足。0.0000831是RMS， 不是最大梯度。

你在13楼中贴出来的
MAXIMUM GRADIENT =  0.0000414    RMS GRADIENT = 0.0000109
则同时满足了两个条件。
赞一下(1人) 回复此楼
14楼2009-03-11 12:03:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peptide2882

木虫 (正式写手)
恩，谢谢啊，GAMESS自带的说明虽然很全，但我觉得写得不是很生动，对于一个初学者或没有量子化学基础的人来说，有时候即使看懂每个单词的意思，也领会不出它的真正意思，但在朋友们的点拨下，往往能达到一语惊醒梦中人的作用，非常感谢热心的回帖啊
赞一下 回复此楼
终于变成木虫了
15楼2009-03-11 21:49:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peptide2882

木虫 (正式写手)
最近忙着联系工作,计算也在断断续续的试, 结果一直不好, 可咋整啊, 理论上反式能量更低, 顺式能量应该更高, 可是我的结果算了好几遍, 跟理论正好相反,顺式能量反而比反式的低, 是不是我的输入文件有问题,
我是这样做的初始坐标:  反式:Chem3D默认的结构;
                                 顺式: 在反式的基础上把键手动旋转成顺式.
顺式和反式Chem3D简单能量最小化后保存为GAMESS输入文件,然后对命令作些完善,但顺反二者的命令是保持一致的， 是不是我手动调成顺式的做法不正确呢,应该怎样好呢,望高手指点,谢谢
赞一下 回复此楼
终于变成木虫了
16楼2009-03-22 21:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主
谢谢 大讨论啊
回复此楼
17楼2009-06-13 21:51:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 peptide2882 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 9/450 2026-11-11 04:54 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:59 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:32 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 01:21 by Z7qfCsdIiD1T
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:12 by Z7qfCsdIiD1T
[考博] 申请环境，微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭