计算化学公社

标题: NICS计算问题 [打印本页]

作者
Author: shanbei123 时间: 2018-3-18 19:00
标题: NICS计算问题
使用赝势基组，用Multiwfn生成nics计算输入文件，用G09运行出现401错误，是什么原因呢？输入文件如下：

#p b3lyp/genecp guess=huckel NMR geom=connectivity

b3lyp opted

0 1

C 4.52232400 -2.94758500 0.28358300

C 3.24113000 -3.44141000 0.49361100

C 2.12540000 -2.59222800 0.37990900

C 2.26893600 -1.24139900 0.05967100

C 3.58869000 -0.79228900 -0.15016400

C 4.71414700 -1.60031500 -0.04610600

.......

Bq -10.566555 -7.861814 -5.467011

Bq -10.566555 -7.861814 -5.168435

Bq -10.566555 -7.861814 -4.869859

Bq -10.566555 -7.861814 -4.571283

Bq -10.566555 -7.861814 -4.272707

Bq -10.566555 -7.861814 -3.974132

..........

.....

8000

C H O N 0

6-31G(d)

****

Pt 0

LanL2DZ

****

Pt 0

LanL2DZ

作者
Author: sobereva 时间: 2018-3-19 07:22
这叫ICSS不叫NICS

没输出文件没法说，能否运行成功和版本也有关系
把输入文件压缩后上传

作者
Author: shanbei123 时间: 2018-3-20 09:48
这是我的输入文件

作者
Author: sobereva 时间: 2018-3-20 12:11
G16 A.03下去掉guess=huckel后可以跑
G09的话，改成guess=core再试

作者
Author: shanbei123 时间: 2018-3-21 21:56
我根据NICS的计算，运行任务后一直出现l612错误，这是什么原因呢？
输入如下
#p b3lyp/genecp NMR pop=NBO6read

opted

0 1
C                4.52232400 -2.94758500 0.28358300
C                3.24113000 -3.44141000 0.49361100
C                2.12540000 -2.59222800 0.37990900
C                2.26893600 -1.24139900 0.05967100
C                3.58869000 -0.79228900 -0.15016400
C                4.71414700 -1.60031500 -0.04610600
C                -0.56898700 -1.28159000 0.30139800
Pt                0.84625100 0.14946400 0.08345000
O                0.93013400 -3.21760100 0.60658000
N                3.65747400 0.59259200 -0.48665200
C                4.70278400 1.34285300 -0.99559400
C                4.18323400 2.55517600 -1.31752600
H                5.70480200 0.95939400 -1.09999300
H                4.64716400 3.43145800 -1.74341500
N                -2.50531000 0.32209900 -0.04414400
C                -2.00276400 1.47160200 0.55616900
C                -2.87399200 2.39077700 1.18466500
N                -0.66132800 1.66446100 0.58195900......
bq 3.41010484 -2.10253784 0.17008379
bq 3.4827782130 -1.8756068701 1.1412794015
bq 3.3374314670 -2.3294688099 -0.8011118215
bq -1.61697969 -2.30019869 0.17914314
bq -1.7910789979 -2.2421851626 1.1621608721
bq -1.4428803821 -2.3582122174 -0.8038745921
bq -3.06041974 -0.70211423 -0.17099082
bq -3.3306221740 -0.6651198474 0.7911017225
bq -2.7902173060 -0.7391086126 -1.1330833625
bq -5.12275448 -0.28920043 -0.80651110
bq -5.4390276382 -0.2091224140 0.1387712921
.....
C H O N 0
6-31G(d)
****
Pt 0
LanL2DZ
****

Pt 0
LanL2DZ

52 53 54 55 56 57 58 59 60 61 61 62 63 64 65 66 67 68 69 0
S 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0
0.005 1

$NBO NCS <XYZ MO> $END

作者
Author: sobereva 时间: 2018-3-22 06:58

shanbei123 发表于 2018-3-21 21:56
我根据NICS的计算，运行任务后一直出现l612错误，这是什么原因呢？
输入如下
#p b3lyp/genecp NMR pop=NB ...

52 53 54 55 56 57 58 59 60 61 61 62 63 64 65 66 67 68 69 0
S 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0
0.005 1

意义不明，格式不对

作者
Author: shanbei123 时间: 2018-3-22 09:47
本帖最后由 shanbei123 于 2018-3-22 09:53 编辑

对于NICS计算中虚原子的基函数怎么设定呢？格式是什么呢？
(, 下载次数 Times of downloads: 81) 而且，去掉了这部分的设置，同样出现l612错误，是为什么呢？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)