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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】Gaussian09作业中断，可否接着断点继续计算？
121/2返回列表12下一页
查看: 3187  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿华南师大 070300（化学）304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 354求调剂 +5 Tyoumou 2026-03-18 8/400 2026-03-21 00:35 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
信息提示
关闭
请填处理意见
关闭