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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 【讨论】菜鸟学GAMESS的心得
北京石油化工学院2026年研究生招生接收调剂公告
172/2返回列表上一页12
查看: 1134  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 专硕310求调剂 +6 捞捞我…. 2026-04-04 7/350 2026-04-08 17:02 by kookay8
[考研] 085600材料与化工301分求调剂院校 +27 刺痛jk 2026-04-06 28/1400 2026-04-08 16:16 by luoyongfeng
[考研] 269求调剂 +5 啊啊我我 2026-04-07 5/250 2026-04-08 16:06 by sunhuadong
[考研] 308求调剂 +15 墨墨漠 2026-04-06 15/750 2026-04-08 15:14 by screening
[考研] 298求调剂 +6 钉叮咚冬瓜 2026-04-07 8/400 2026-04-08 10:51 by 中飞院空管学院�
[考研] 求调剂！生物与医药专硕 +9 逆转陆先生 2026-04-01 10/500 2026-04-08 10:37 by maolC
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 一志愿武理车辆专硕总分 281 求调剂 +4 上岸研究生. 2026-04-02 4/200 2026-04-07 09:52 by 加油向未来啊
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 工科08-机械专硕-求调剂 +3 雷欧飞踢 2026-04-02 3/150 2026-04-05 18:49 by 蓝云思雨
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
[考研] 085601，一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
信息提示
关闭
请填处理意见
关闭