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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】Gaussian09作业中断，可否接着断点继续计算？
121/2返回列表12下一页
查看: 3163  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料专业298分有科研经历 硕士研究生调剂自荐信 +6 zyf上岸 2026-03-01 6/300 2026-03-01 12:43 by liqiongjy
[考研] 化工299分求调剂 一志愿985落榜 +3 嘻嘻(*^ω^*) 2026-03-01 3/150 2026-03-01 12:35 by houyaoxu
[考研] 材料学调剂 +7 提神豆沙包 2026-02-28 9/450 2026-03-01 11:50 by 王伟要上岸啊
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 材料类求调剂 +8 wana_kiko 2026-02-28 8/400 2026-03-01 11:44 by 王伟要上岸啊
[考研] 求调剂 +5 repeatt?t 2026-02-28 5/250 2026-03-01 11:43 by 王伟要上岸啊
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[考研] 291分工科求调剂 +7 science饿饿 2026-03-01 8/400 2026-03-01 10:43 by sunny81
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 10/500 2026-03-01 10:02 by 科研狗111
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 化工专硕348，一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考博] 26申博 +4 想申博！ 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
信息提示
关闭
请填处理意见
关闭