| 查看: 2887 | 回复: 8 | ||
中科院苦行僧木虫 (小有名气)
|
[求助]
初学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
- 附件 2 : TDmPW1PW91OPtd.log
2013-10-29 16:51:01, 3.6 M
2013-10-29 16:51:13, 4.03 M
» 猜你喜欢
304求调剂
已经有5人回复
-
PDF在线压缩软件
已经有5人回复
-
物理化学论文润色/翻译怎么收费?
已经有120人回复
-
材料/化学相关专业2026级学术/专业型硕士研究生
已经有19人回复
-
求催化作用基础第三版 （甄开吉) 课件
已经有3人回复
-
26年材料/化学相关专业博士研究生招生 温州大学化学与材料工程学院(26年秋季入学)
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
304求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯TD-DFT计算荧光发射谱的问题
已经有7人回复
- 初学matlab 求解方程 大家帮帮忙啊。。。
已经有11人回复
- 关于TD-DFT优化构型并计算频率的问题
已经有12人回复
- 计算出错，求解。。。
已经有4人回复
- TD-DFT激发态优化不收敛
已经有4人回复
- 静态计算出错，求解
已经有6人回复
- MATLAB初学者，这样的式子怎么就是画不了图，总是出错!求指导!!!!!!!!!!!!!
已经有8人回复
- TD-DFT计算CD光谱问题
已经有12人回复
- 【求助】稀土元素光谱计算方法
已经有5人回复
- 【求助】含时密度泛函（TD-DFT）在gaussian的输入文件中如何表示？
已经有7人回复
- 【求助】初学有限元，弱弱问一句，如果已知边界都是位移边界条件，怎么求解？
已经有11人回复
- 【求助】关于TD-DFT模拟Uv-Vis光谱的问题
已经有45人回复
- 【求助】TD-DFT 优化求救
已经有7人回复
- 【求助】TD-DFT优化激发态构型后的频率问题
已经有11人回复
lishijunzong
铁杆木虫 (正式写手)
- 应助: 214 (大学生)
- 金币: 5034
- 红花: 10
- 帖子: 574
- 在线: 195.1小时
- 虫号: 1121020
- 注册: 2010-10-13
- 性别: GG
- 专业: 理论和计算化学
2楼2013-10-29 17:22:02
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
3楼2013-10-30 06:46:51
十五的月牙
金虫 (正式写手)
- 应助: 70 (初中生)
- 金币: 928.8
- 沙发: 4
- 帖子: 439
- 在线: 59.7小时
- 虫号: 2759595
- 注册: 2013-10-28
- 性别: GG
- 专业: 分离过程
4楼2013-10-30 09:48:23
PowerPaul
木虫 (初入文坛)
博士
- 应助: 2 (幼儿园)
- 金币: 2205
- 帖子: 42
- 在线: 45.1小时
- 虫号: 2129978
- 注册: 2012-11-16
- 性别: GG
- 专业: 理论和计算化学
5楼2013-10-30 15:24:49
6楼2013-10-31 22:06:08
7楼2013-10-31 22:14:59
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
8楼2013-11-01 06:35:54
中科院苦行僧
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1582.6
- 红花: 20
- 帖子: 247
- 在线: 64.2小时
- 虫号: 1501904
- 注册: 2011-11-20
- 性别: GG
- 专业: 理论和计算化学
|
我是分开算的 成功了 给你一个例子： 先算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文件 |
9楼2013-11-20 09:27:37
20