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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】Gaussian09作业中断，可否接着断点继续计算？
121/2返回列表12下一页
查看: 3226  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xmdxn

铁虫 (初入文坛)

[交流] 【求助】Gaussian09作业中断，可否接着断点继续计算？ 已有6人参与

如题。
正在用td-b3lyp优化perylene的激发态，分子是C20H12，基组是6-311+g**，跑频率分析跑了四五天还没有跑出来，恰好又遇到服务器死机。。。
哪位知道Gaussian09的作业中间断掉可不可以接着断点计算啊？小弟裸膝跪拜，不胜感激~
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-15 12:50:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人
引用回帖:
Originally posted by xmdxn at 2010-10-15 12:50:51:
如题。
正在用td-b3lyp优化perylene的激发态，分子是C20H12，基组是6-311+g**，跑频率分析跑了四五天还没有跑出来，恰好又遇到服务器死机。。。
哪位知道Gaussian09的作业中间断掉可不可以接着断点计算啊？小弟 ...

perhaps not
回复此楼
惺惺相惜情不自禁
2楼2010-10-15 13:22:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmdxn

铁虫 (初入文坛)
引用回帖:
Originally posted by manson1998 at 2010-10-15 13:22:06:


perhaps not

what a pity...
不知道大家有没有用td优化激发态几何构型的经验，做频率分析的速度是不是对基组大小很敏感？我正在考虑换小基组做下去了。。。
赞一下 回复此楼
3楼2010-10-15 16:08:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougaozhan

至尊木虫 (职业作家)

busy
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zhou2009(金币+1): 2010-10-16 15:59:58
基组6-311+g**太大了，难怪跑不过去
赞一下(2人) 回复此楼
^_^不要生气||要争气，不要看破||要突破，不要嫉妒||要欣赏，不要托延||要积极，不要心动||要行动^_^
4楼2010-10-15 17:37:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmdxn

铁虫 (初入文坛)
引用回帖:
Originally posted by gougaozhan at 2010-10-15 17:37:57:
基组6-311+g**太大了，难怪跑不过去

我现在又换了6-31+g*和6-31g*继续优化和频率分析，激发态的优化很快，一般两个小时就可以做出来，但是频率分析还是要跑很久，目前跑了两千五百分钟还是没有跑出来，不停地出现
                                                                     
**** Warning!!: The largest alpha MO coefficient is  0.28021897D+03
                                                                     
                                                                     
**** Warning!!: The smallest alpha delta epsilon is  0.95046157D-01
                                                                     
这样的循环，the largest alpha MO coefficient is 后面的值基本上在0.28D+03上下浮动，我在想这是进入了无法跑出结果来的循环么？
之前做激发态大家都是用casscf，老师让我用td，也没人有经验，用过的虫友帮帮忙解答一下呗~~
赞一下 回复此楼
5楼2010-10-16 14:27:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
频率没法接着算，死了就自认倒霉吧！
赞一下(1人) 回复此楼
Anewday,anewhour,anewminute,anewpeople.
6楼2010-10-16 15:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zhou2009(金币+3): 2010-10-16 16:07:10
激发态优化频率用的是数值梯度，计算过程总共需要计算6N＋1次SCF和TD。
N为原子个数，你看看你的输出文件里有多少次SCF Done了？用6N＋1减去这个数字就是还需要跑的次数，大致应该能估算出时间的。

提供一个脚本，可以用来检查激发态频率已经跑了多少步，还需要跑多少步。
CODE:
#!/bin/sh

if [ $# -ne 1 ]
then
        echo "Usage: $0 Gaussian output file"
        exit 1
fi

if [ ! -e $1 ]
then
        echo "$1 not existed!"
        exit 1
fi

scf_iters=`grep "SCF Done" $1 | wc -l`
natoms=`grep "NAtoms" $1 | awk '{print $2}'| head -1`
left_iters=$(($natoms*6+1-$scf_iters))
echo ""
echo "The SCF has runned for $scf_iters cycles"
echo -e "Still need \033[45m\033[5m $left_iters \033[0m cycles"
echo ""

顺便提一下，G09的频率计算绝对可以接着算，前提是保存了rwf文件，具体如何做，请参考说明书上restart关键词的说明。
引用回帖:
Originally posted by xmdxn at 2010-10-16 14:27:36:


我现在又换了6-31+g*和6-31g*继续优化和频率分析，激发态的优化很快，一般两个小时就可以做出来，但是频率分析还是要跑很久，目前跑了两千五百分钟还是没有跑出来，不停地出现
                             ...

赞一下(2人) 回复此楼
7楼2010-10-16 15:11:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmdxn

铁虫 (初入文坛)
引用回帖:
Originally posted by pwzhou at 2010-10-16 15:11:09:
激发态优化频率用的是数值梯度，计算过程总共需要计算6N＋1次SCF和TD。
N为原子个数，你看看你的输出文件里有多少次SCF Done了？用6N＋1减去这个数字就是还需要跑的次数，大致应该能估算出时间的。

提供一个脚 ...

原来是数值梯度哦。Gaussian09说明书上写的应该是说一阶导为积分的。
谢谢这位虫友，应该可以接着算，利用rwf文件可以的，似乎chk文件也可以。
用ultraedit搜索字符串也可以看出来scf done了几次O(∩_∩)O~编程序的东西我就不大灵光啦
赞一下 回复此楼
8楼2010-10-16 18:18:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zhou2009(金币+1): 2010-10-16 18:36:21
要断点续算，rwf和chk都必须要，缺一不可。
引用回帖:
Originally posted by xmdxn at 2010-10-16 18:18:47:


原来是数值梯度哦。Gaussian09说明书上写的应该是说一阶导为积分的。
谢谢这位虫友，应该可以接着算，利用rwf文件可以的，似乎chk文件也可以。
用ultraedit搜索字符串也可以看出来scf done了几次O(∩_∩)O~ ...

赞一下(2人) 回复此楼
9楼2010-10-16 18:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmdxn

铁虫 (初入文坛)
ls这位虫友，我用6-31g*重新做频率分析，终于跑出来了。我的体系是C20H12，32个原子，但是scf迭代了44次。照6N+1来算，不止44次啊。。。又困惑了，G09的TD频率分析的梯度到底是什么状况哦。。。
赞一下 回复此楼
10楼2010-10-18 15:00:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xmdxn 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[考研] 广东省 085601 329分求调剂 +13 Eddieddd 2026-04-10 13/650 2026-04-10 23:27 by 314126402
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 346，工科0854求调剂，专硕 +5 moser233 2026-04-10 6/300 2026-04-10 18:04 by 周邵阳
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 电子信息279求调剂，有书读就行 +4 wwwooden 2026-04-08 7/350 2026-04-10 11:34 by wwwooden
[考研] 生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 求调剂，262机械专硕 +6 嗯yyl 2026-04-08 6/300 2026-04-09 12:01 by zhouyuwinner
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 化学357分，考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭