24小时热门版块排行榜    

查看: 1996  |  回复: 9

732749120

金虫 (初入文坛)

[求助] NTO出错

在计算完激发态的fchk文件的基础上怎么算NTO?我的输入文件如下:
%Chk=state148.chk
%mem=2500MW
%nprocl=2
# td(read)=(nstates=800,root=600) b3lyp/lanl2dz Guess=Read Pop=(SaveNTO,NTO) density=Transition=148
可是老出错,请问是什么原因,怎么改?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
最好能贴出出错的地方。
你要读取的是fchk,但是你用%chk,是可以这样用的?我是没这样用过
2楼2013-10-26 12:19:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

732749120

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2013-10-26 12:19:38
最好能贴出出错的地方。
你要读取的是fchk,但是你用%chk,是可以这样用的?我是没这样用过

我看的是G09help,
SaveNaturalTransitionOrbitals
Save the generated orbitals in the checkpoint file, replacing the canonical ones if the density was read-in from there. SaveNTO is a synonym for this option. If you want to visualize the orbitals, you need to write them back to the checkpoint file. It is a good idea to do so with a copy of the checkpoint file for each state. After the initial excited state calculation (using %Chk ex.chk), use a technique like the following to generate visualization data for each state:

$ cp ex.chk staten.chk                     Copy original checkpoint file.
$ g09 <<END                                Run a Guess=Only job just to generate orbitals.
%Chk=staten                                Use the checkpoint file copy.
# Geom=AllCheck ChkBas Guess=(Read,Only)
  Density=Transition=n Pop=SaveNTO         Repeat process for each state of interest.
END                                        End of Gaussian 09 input file.


NATURAL ORBITAL-RELATED OPTIONS
我按那步骤做的,可是老出错,
3楼2013-10-27 17:06:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

732749120

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2013-10-26 12:19:38
最好能贴出出错的地方。
你要读取的是fchk,但是你用%chk,是可以这样用的?我是没这样用过

我用g09打开fchk文件做成这样的输出文件:
%Chk=state148.chk
%mem=2500MW
%nprocl=2
# td(read)=(nstates=800,root=600) b3lyp/lanl2dz Guess=Read Pop=(SaveNTO,NTO) density=Transition=148
结果老出错,错误如下:
End of file in ZSymb.
Error termination via Lnk1e in /scratch/g09/l101.exe at Sun Oct 27 16:36:19 2013.
4楼2013-10-27 17:24:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 732749120 at 2013-10-27 17:24:12
我用g09打开fchk文件做成这样的输出文件:
%Chk=state148.chk
%mem=2500MW
%nprocl=2
# td(read)=(nstates=800,root=600) b3lyp/lanl2dz Guess=Read Pop=(SaveNTO,NTO) density=Transition=148
结果老出错, ...

L101错误说明出错在读取分子说明或标题上,看了一下你的输入,似乎没有对分子说明做定义?

尝试加geom=allcheck试试
另外输入完左右东西最后要空一行。
5楼2013-10-28 07:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

732749120

金虫 (初入文坛)

引用回帖:
5楼: Originally posted by 枪下游魂 at 2013-10-28 07:15:15
L101错误说明出错在读取分子说明或标题上,看了一下你的输入,似乎没有对分子说明做定义?

尝试加geom=allcheck试试
另外输入完左右东西最后要空一行。...

我又改了改输入文件如下:
%chk=staten.chk
%mem=2500MW
# Geom=allcheck Chkbas Guess=(read,only) density=transition=148 Pop=(NTO,saveNTO)
结果又错了,错误如下:
Transition densities not found by NTOOrb.
Error termination via Lnk1e in /scratch/g09/l601.exe
这是怎么回事?
6楼2013-10-28 09:06:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by 732749120 at 2013-10-28 09:06:01
我又改了改输入文件如下:
%chk=staten.chk
%mem=2500MW
# Geom=allcheck Chkbas Guess=(read,only) density=transition=148 Pop=(NTO,saveNTO)
结果又错了,错误如下:
Transition densities not found by N ...

出错信息提示没有找到跃迁密度。
输入文件看不出你有算激发态,只有基态,哪来第1~第148的跃迁密度?
可尝试下density=alltransition看是什么结果。

另外既然你是关心第600个激发态,为何跃迁密度只用到第148的?望指教~
7楼2013-10-28 09:50:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

732749120

金虫 (初入文坛)

引用回帖:
7楼: Originally posted by 枪下游魂 at 2013-10-28 09:50:14
出错信息提示没有找到跃迁密度。
输入文件看不出你有算激发态,只有基态,哪来第1~第148的跃迁密度?
可尝试下density=alltransition看是什么结果。

另外既然你是关心第600个激发态,为何跃迁密度只用到第148 ...

我是想看看某个激发态下的某个轨道,如
Excited State 148:      Singlet-A      2.5266 eV  490.71 nm  f=0.0205  <S**2>=0.000
     271 ->303        -0.16944
     272 ->304        -0.20987
     273 ->303         0.29919
     274 ->303        -0.14255
     275 ->304         0.15035
     280 ->304        -0.16007
     283 ->305         0.14589
     289 ->307         0.10030
     290 ->307        -0.16766
     290 ->308        -0.12994
     293 ->307        -0.10140
     294 ->309         0.12905
     298 ->310         0.11517
     301 ->314        -0.10997
我是想看看273 ->303  到底是哪个轨道到哪个轨道
8楼2013-10-28 10:06:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

732749120

金虫 (初入文坛)

引用回帖:
7楼: Originally posted by 枪下游魂 at 2013-10-28 09:50:14
出错信息提示没有找到跃迁密度。
输入文件看不出你有算激发态,只有基态,哪来第1~第148的跃迁密度?
可尝试下density=alltransition看是什么结果。

另外既然你是关心第600个激发态,为何跃迁密度只用到第148 ...

你有什么好的方法吗?望指教
9楼2013-10-28 10:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
732749120: 金币+10, 有帮助 2013-10-29 17:08:18
引用回帖:
9楼: Originally posted by 732749120 at 2013-10-28 10:08:30
你有什么好的方法吗?望指教...

我觉得你可以先找一个小分子,用你理解的方法,跑优化,算下td,用density=transition=1,尝试能不能正常跑完。
等能正常跑完了,再来算你的目标体系。
这样过一遍,至少能保证解决你在命令行输入上所遇到的问题。
10楼2013-10-28 11:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 732749120 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 6/300 2026-08-28 01:46 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +9 Kittylucky 2026-08-27 9/450 2026-08-28 01:08 by jurkat.1640
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见