24小时热门版块排行榜    

查看: 2092  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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.

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

非常感谢!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 8 个回答

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
前途光明,出路难觅!
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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见