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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】TD-DFT 优化求救
81/1返回列表
查看: 2089  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

snoopyzhao

至尊木虫 (职业作家)

[交流] 【求助】TD-DFT 优化求救已有6人参与

我的输入文件如下：
CODE:
%chk="benz_td_631gd.chk"
%mem=12GB
%nprocs=8
#P b3lyp/6-31g(d) opt

dft opt

0,1
C
C,1,CC
C,1,CC,2,120.
C,2,CC,1,120.,3,0.,0
C,3,CC,1,120.,2,0.,0
C,4,CC,2,120.,1,0.,0
H,1,CH,2,120.,4,180.,0
H,2,CH,1,120.,7,0.,0
H,3,CH,1,120.,7,0.,0
H,4,CH,2,120.,8,0.,0
H,5,CH,3,120.,9,0.,0
H,6,CH,4,120.,10,0.,0
      Variables:
CC=1.38270049
CH=1.07316153

--Link1--
%chk="benz_td_631gd.chk"
%mem=12GB
%nprocs=8
#P b3lyp/6-31g(d) geom=check td

td-dft Excited States

0,1

--Link1--
%chk="benz_td_631gd.chk"
%mem=12GB
%nprocs=8
#P td(root=1,read) b3lyp/6-31g(d) guess=read geom=check opt freq

tddft Excited States Opt Freq

0,1

运行到最后一个任何时，出现了如下错误：
引用回帖:
******************************************
Gaussian 03:  EM64T-G03RevD.02  1-Mar-2006
                27-May-2010
******************************************
%chk="benz_td_631gd.chk"
%mem=12GB
%nprocs=8
Will use up to    8 processors via shared memory.
Default route:  MaxDisk=300GB
----------------------------------------------------------------
#P td(root=1,read) b3lyp/6-31g(d) guess=read geom=check opt freq
----------------------------------------------------------------

Warning:  this job cannot use analytic gradients
           and so will do many energy evaluations.

1/14=-1,26=3,29=20002,38=1/1,14;
2/29=3,40=1/2;
3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,74=-5/1,2,8,3;
4/5=1/1;
5/5=2,38=6/2;
8/6=1,10=2,27=40265318400/1;
9/8=1,27=40265318400,42=1,49=3/14;
6/7=2,8=2,9=2,10=2/1;
1/14=-1/14(3);
2/29=3/2;
6/7=2,8=2,9=2,10=2/1;
99//99;
2/29=3/2;
3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,74=-5/1,2,8,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
8/6=1,10=2,27=40265318400/1;
9/8=1,27=40265318400,42=1,49=4/14;
1/14=-1/14(-6);
2/29=3/2;
6/7=2,8=2,9=2,10=2/1;
99//99;
Leave Link    1 at Thu May 27 23:33:59 2010, MaxMem= 1610612736 cpu:       0.4
(Enter /home/software/Gaussian/g03/l101.exe)
-----------------------------
tddft Excited States Opt Freq
-----------------------------
Redundant internal coordinates taken from checkpoint file:
"benz_td_631gd.chk"
Charge =  0 Multiplicity = 1
C,0,-0.0120573677,0.,-0.0069613245
C,0,-0.0120573677,0.,1.3896618145
C,0,1.1974537502,0.,-0.7052728939
C,0,1.1974537502,0.,2.0879733839
C,0,2.406964868,0.,-0.0069613245
C,0,2.406964868,0.,1.3896618145
H,0,-0.9534223209,0.,-0.5504586336
H,0,-0.9534223209,0.,1.9331591236
H,0,1.1974537502,0.,-1.7922675122
H,0,1.1974537502,0.,3.1749680022
H,0,3.3483298212,0.,-0.5504586336
H,0,3.3483298212,0.,1.9331591236
Recover connectivity data from disk.
                    Isotopes and Nuclear Properties:
(Nuclear quadrupole moments (NQMom) in fm**2, nuclear magnetic moments (NMagM)
  in nuclear magnetons)

  Atom         1           2           3           4           5           6           7           8           9          10
IAtWgt=          12          12          12          12          12          12           1           1           1           1
AtmWgt=  12.0000000  12.0000000  12.0000000  12.0000000  12.0000000  12.0000000   1.0078250   1.0078250   1.0078250   1.0078250
NucSpn=           0           0           0           0           0           0           1           1           1           1
AtZEff=  -3.6000000  -3.6000000  -3.6000000  -3.6000000  -3.6000000  -3.6000000  -1.0000000  -1.0000000  -1.0000000  -1.0000000
NQMom=    0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   0.0000000
NMagM=    0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   0.0000000   2.7928460   2.7928460   2.7928460   2.7928460

  Atom        11          12
IAtWgt=           1           1
AtmWgt=   1.0078250   1.0078250
NucSpn=           1           1
AtZEff=  -1.0000000  -1.0000000
NQMom=    0.0000000   0.0000000
NMagM=    2.7928460   2.7928460
Leave Link  101 at Thu May 27 23:33:59 2010, MaxMem= 1610612736 cpu:       0.2
(Enter /home/software/Gaussian/g03/l114.exe)

NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-
NUMERICAL EIGENVECTOR FOLLOWING MINIMUM SEARCH
INITIALIZATION PASS


************************************************
** ERROR IN INITNF. NUMBER OF VARIABLES (  0) **
**   INCORRECT (SHOULD BE BETWEEN 1 AND 50)   **
************************************************


Error termination via Lnk1e in /home/software/Gaussian/g03/l114.exe at Thu May 27 23:33:59 2010.

请问是什么原因造成的？如何修改输入文件？

非常感谢！！
回复此楼

» 猜你喜欢

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

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上
★
snoopyzhao(金币+3):谢谢！回头我试一下，另外，第二步中没有变量部分啊，只有第一步中有这些变量，把它们后面加一个 D0？ 2010-05-28 13:03:25
aylayl08(金币+1):谢谢交流 2010-05-28 13:37:55
在第二步，把变量部分改成这样试一试。
CC=1.38270049 D
CH=1.07316153 D
赞一下(1人) 回复此楼
前途光明，出路难觅！
2楼2010-05-28 08:10:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Roseshu

金虫 (正式写手)
snoopyzhao(金币+2):变量太多是什么意思？ 2010-05-28 19:38:57
不是因为变量太多了吗？
赞一下 回复此楼
3楼2010-05-28 14:03:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Roseshu

金虫 (正式写手)
snoopyzhao(金币+5):多谢！那么如何从多步任务中获取上一步的结构，并且是 z-matrix 的格式呢？ 2010-05-29 11:58:13
INITNF. NUMBER OF VARIABLES SHOULD BE BETWEEN   1 AND 50
xyz坐标它可能判断不对自由度，写成z坐标，保证变量数在1-50之间
赞一下 回复此楼
4楼2010-05-29 10:16:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
非常感谢！

在多步任务如何使得结构一直使用 z-matrix 呢？
赞一下 回复此楼
5楼2010-05-29 11:59:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangraul88

铁虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你这个有个很明显的错误啊，输入命令里td后面只有root=1，没有nstats啊！指定下nstats数目试试，应该OK！
赞一下 回复此楼
6楼2012-07-16 17:07:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tudoudigua

金虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+2 2012-07-18 20:24:54
你的输入文件出现很大的错误
b3lyp/6-31g(3df,3pd) td(nstate=40,singlet)

如果需要拷贝几何，在加入这个geom=check guess=read


nstate=多少看具体情况而定
赞一下 回复此楼
人生总是在习惯那些不习惯的习惯！
7楼2012-07-17 15:19:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwenpeng

新虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
140952楼: Originally posted by snoopyzhao at 2010-05-27 23:55:10
我的输入文件如下：


%chk="benz_td_631gd.chk"
%mem=12GB
%nprocs=8
#P b3lyp/6-31g(d) opt

dft opt

0,1
C
C,1,CC
C,1,CC,2,120.
C,2,CC,1,120.,3,0.,0
C,3,CC,1,120.,2,0.,0
...

请问这个td-dft的输入文件是怎么整出来的？
赞一下 回复此楼
8楼2012-09-28 17:04:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 snoopyzhao 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭