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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 初学TD-DFT计算一直出错，求解
91/1返回列表
查看: 2679  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

中科院苦行僧

木虫 (小有名气)

[求助] 初学TD-DFT计算一直出错，求解

最近刚开始学者做TDDFT的计算，算了很多还一直出错，很急，各位大神，给小子支支招：

TD-B3LYP

我是在基态B3LYP构型的基础上做TD计算的:输入文件：

%nprocshared=8
%mem=4000MW
%chk=TDb3lypoptb.chk
#p b3lyp/6-311++g(d,p) td(nstate=6,singlets,root=1) pop=full opt freq

结果：
           (ROBERT BURNS 'TO A LOUSE')
Leave Link 9999 at Mon Oct 28 18:40:57 2013, MaxMem= 4194304000 cpu:       1.0
Job cpu time:  0 days  2 hours 48 minutes  4.3 seconds.
File lengths (MBytes):  RWF=    107 Int=      0 D2E=      0 Chk=     17 Scr=      1
Normal termination of Gaussian 09 at Mon Oct 28 18:40:57 2013.
(Enter /public/grid_app//g09/l1.exe)
Link1:  Proceeding to internal job step number  2.
----------------------------------------------------------------------
#P Geom=AllCheck Guess=TCheck SCRF=Check GenChk RTD-B3LYP(FC)/6-311++G
(d,p) Freq
----------------------------------------------------------------------
QPErr --- A syntax error was detected in the input line.
Check SCRF=Check GenChk RTD-B3LYP(FC)/6-
                         '
Last state="GCL"
TCursr= 1046 LCursr=   48
Error termination via Lnk1e in /public/grid_app//g09/l1.exe at Mon Oct 28 18:40:57 2013.
Job cpu time:  0 days  0 hours  0 minutes  0.4 seconds.
File lengths (MBytes):  RWF=    107 Int=      0 D2E=      0 Chk=     17 Scr=      1


log文件件附件1

TD-mPW1PW91

同样在基态mPWpw91的构型上优化S1太的构型

输入文件：

%nprocshared=8
%mem=4000MW
%chk=TDmPW1PW91OPtd.chk
#p mpw1pw91/6-311++g(d,p) td(nstate=6,singlets,root=1) pop=full opt freq

错误：

BE LIKE A POSTAGE STAMP.  STICK TO ONE THING UNTIL
YOU GET THERE -- JOSH BILLINGS
Leave Link 9999 at Mon Oct 28 18:22:01 2013, MaxMem= 4194304000 cpu:       0.8
Job cpu time:  0 days 15 hours 20 minutes  3.7 seconds.
File lengths (MBytes):  RWF=    139 Int=      0 D2E=      0 Chk=     17 Scr=      1
Normal termination of Gaussian 09 at Mon Oct 28 18:22:01 2013.
(Enter /public/grid_app//g09/l1.exe)
Link1:  Proceeding to internal job step number  2.
----------------------------------------------------------------------
#P Geom=AllCheck Guess=TCheck SCRF=Check GenChk RTD-mPW1PW91(FC)/6-311
++G(d,p) Freq
----------------------------------------------------------------------
QPErr --- A syntax error was detected in the input line.
Check SCRF=Check GenChk RTD-mPW1PW91(FC)
                         '
Last state="GCL"
TCursr= 1046 LCursr=   48
Error termination via Lnk1e in /public/grid_app//g09/l1.exe at Mon Oct 28 18:22:01 2013.
Job cpu time:  0 days  0 hours  0 minutes  0.3 seconds.
File lengths (MBytes):  RWF=    139 Int=      0 D2E=      0 Chk=     17 Scr=      1
log文件件附件2
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : TDb3lypoptd.log
    • 2013-10-29 16:51:01, 3.6 M
  • 附件 2 : TDmPW1PW91OPtd.log
    • 2013-10-29 16:51:13, 4.03 M

» 猜你喜欢

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

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
出错这一行写的有问题，程序不认了。这样还不如分开计算就行，再拿你算好的继续算吧
赞一下 回复此楼
2楼2013-10-29 17:22:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
gkf高: 金币+2 2013-10-30 09:56:41
这个问题我觉得是高斯的bug，即在计算激发态时不能选择opt+freq这种形式，必须分开计算。
一种方式是先做opt，再另写文件读入opt的chk文件算freq。
另一种方式是在同一个输入文件里先做opt，再用--Link1--分隔后写freq的输入文件。
赞一下(2人) 回复此楼
3楼2013-10-30 06:46:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

十五的月牙

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
反正我一般都是先优化 没怎么出过错
赞一下 回复此楼
4楼2013-10-30 09:48:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

PowerPaul

木虫 (初入文坛)

博士

【答案】应助回帖

感谢参与，应助指数 +1
分开算吧。我之前也遇到过类似的问题，最后还是要分开算……
赞一下 回复此楼
Youcandoitwell.
5楼2013-10-30 15:24:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pq33714

金虫 (正式写手)
引用回帖:
5楼: Originally posted by PowerPaul at 2013-10-30 15:24:49
分开算吧。我之前也遇到过类似的问题，最后还是要分开算……

怎么分开算，我是新手，不知道优化后怎么构建输入文件进行下一步了，求指教，谢谢！
赞一下 回复此楼
6楼2013-10-31 22:06:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pq33714

金虫 (正式写手)
引用回帖:
3楼: Originally posted by 枪下游魂 at 2013-10-30 06:46:51
这个问题我觉得是高斯的bug，即在计算激发态时不能选择opt+freq这种形式，必须分开计算。
一种方式是先做opt，再另写文件读入opt的chk文件算freq。
另一种方式是在同一个输入文件里先做opt，再用--Link1--分隔后写 ...

用什么软件读入chk？
赞一下 回复此楼
7楼2013-10-31 22:14:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by pq33714 at 2013-10-31 22:14:59
用什么软件读入chk？...

在输入文件里用%chk，如果这个chk文件存在，会先读取里面的数据再覆盖。
赞一下 回复此楼
8楼2013-11-01 06:35:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

中科院苦行僧

木虫 (小有名气)
引用回帖:
6楼: Originally posted by pq33714 at 2013-10-31 22:06:08
怎么分开算，我是新手，不知道优化后怎么构建输入文件进行下一步了，求指教，谢谢！...

我是分开算的 成功了 给你一个例子：

先算opt，构型坐标是基态的B3LYP构型，直接拿来用：
%nprocshared=8
%mem=2GB
%chk=ofscisS1-opt-tdb3lyp-normal.chk
#p opt td=(singlets,root=1) rb3lyp/6-311++g(d,p) nosymm

Title Card Required

0 1
C                 -0.75352300    1.25521900    0.00019900
C                  0.65389300    1.27483400    0.00067300
C                  1.24110700    2.54686800   -0.00010000
C                  0.51186600    3.72316300   -0.00114700
C                 -0.87979000    3.66066100   -0.00149800
C                 -1.51327000    2.41916900   -0.00088600
H                 -1.25063400    0.29102400    0.00071200
H                  1.04159900    4.66781500   -0.00165900
H                 -1.45968600    4.57598800   -0.00234500
H                 -2.59506200    2.35770800   -0.00121400
F                  2.59451700    2.65688600    0.00028200
C                  1.37927900   -0.00180500    0.00183500
H                  0.71980600   -0.86684600    0.00247900
C                  2.69847200   -0.22896200    0.00202900
H                  3.43477800    0.56163100    0.00128300
H                  3.06645500   -1.24845400    0.00279600

算完了用GV打开，在此构型上写写freq的输入文件：


%nprocshared=8
%mem=2GB
%chk=ofscisS1-freq-tdb3lyp-normal.chk
#p freq=noraman rb3lyp/6-311++g(d,p) nosymm

Title Card Required

0 1
C                 -0.77498000    1.22967800    0.00022800
C                  0.65928200    1.22086800    0.00058800
C                  1.23473500    2.53609700   -0.00019200
C                  0.50743800    3.75060400   -0.00122200
C                 -0.89655900    3.70115100   -0.00152500
C                 -1.51835000    2.43681700   -0.00079200
H                 -1.29810200    0.28217900    0.00077000
H                  1.06185500    4.67949100   -0.00174600
H                 -1.48326200    4.60938000   -0.00230700
H                 -2.60075100    2.37858400   -0.00101000
F                  2.56779400    2.65735300    0.00005100
C                  1.38001300   -0.00486700    0.00162900
H                  0.73765800   -0.88244900    0.00207700
C                  2.73711700   -0.23882800    0.00214600
H                  3.47693300    0.54774300    0.00179800
H                  3.09898600   -1.25890100    0.00294700

两步计算完就OK了，实验室的师兄说计算频率时不是用TD，而是用基态的方法，这个也不知道为什么~~~

我准备验证下用能不能用TD算S1的频率。

关于CHk文件的读入，我给个例子：
%nprocshared=8
%mem=1000MW
%chk=ofscisS1-opt-tdb3ltp.chk
#p opt td(singlets,root=1) B3lyp/6-311++g(d,p) nosymm geom=checkpoint

Title Card Required

0 1
  这里的构型坐标是直接由geom=checkpoint给定的，还可以用其他的关键字，你可以查查相关的关键字。关键是提交计算的时候，你要把gif文件盒chk文件提交到一个任务上，且chk文献的名称必须是与属于文件的%CHk文件的名称相同，计算时服务器会先读这个CHk文件，计算时会产生新的同名CHk文件覆盖原有CHk文件
赞一下 回复此楼
beayardstickogquality
9楼2013-11-20 09:27:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员PzX7SI 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭