24小时热门版块排行榜

>论坛更新日志 (3807)
>导师招生 (699)
>论文投稿 (390)
>休闲灌水 (217)
>硕博家园 (211)
>文献求助 (137)
>考博 (135)
>新药研发 (98)
>公派出国 (78)
>学术会议 (73)
>考研 (67)
>虫友互识 (65)
>基金申请 (56)
>论文道贺祈福 (48)
>教师之家 (44)
>招聘信息布告栏 (39)

返回列表

当前主题已经存档。

wolfzhong

木虫 (正式写手)

应助: 23 (小学生)
金币: 3429.5
散金: 274
红花: 9
帖子: 986
在线: 751.4小时
虫号: 619989
注册: 2008-10-08
专业: 理论和计算化学

[交流] 【求助】ECP的问题

昨天尝试ECP过程中遇到一个很麻烦的问题。
我第一个例子：
#p b3lyp/gen opt(tight,gdiis) frequency=noraman scf(dsymm) Pseudo=read

Se

-4 1
Fe          0       1.495886 0.000000 0.000000
Fe          0    -1.495886 0.000000 0.000000
Se          0       0.000000 0.000000 1.656633
Se          0       0.000000 0.000000 -1.656633
Cl          0    -3.433878 2.141363 0.000000
Cl          0    -3.433878 -2.141363 0.000000
Cl          0       3.433878 -2.141363 0.000000
Cl          0       3.433878 2.141363 0.000000

Fe 0
6-311+G*
****
Cl 0
6-311+G*
****
Se 0
lanl2dz
****

Se 0
lanl2dz

可以正常运行。

但金属换成同族的Te原子时，即：
Te

-4 1
Fe          0       1.495886 0.000000 0.000000
Fe          0    -1.495886 0.000000 0.000000
Te          0       0.000000 0.000000 1.656633
Te          0       0.000000 0.000000 -1.656633
Cl          0    -3.433878 2.141363 0.000000
Cl          0    -3.433878 -2.141363 0.000000
Cl          0       3.433878 -2.141363 0.000000
Cl          0       3.433878 2.141363 0.000000

Fe 0
6-311+G*
****
Cl 0
6-311+G*
****
Te 0
sdd
****

Te 0
sdd

结果在读ECP时候给出这样提示：
Centers: 1  2
6-311+G*
****
Centers: 5  6  7  8
6-311+G*
****
Centers: 3  4
sdd
****
===================================================================================================================================
Pseudopotential Parameters
===================================================================================================================================
  Center    Atomic    Valence    Angular    Power                                                    Coordinates
  Number    Number    Electrons    Momentum    of R    Exponent       Coefficient             X          Y          Z
===================================================================================================================================
EOF while reading ECP pointer card.
Error termination via Lnk1e in /share/apps/gauss/g03/l301.exe at Sat Mar  6 11:13:35 2010.
Job cpu time:  0 days  0 hours  0 minutes 20.8 seconds.

别人文献已有报道，而我这里读ecp怎么出错。谢谢

回复此楼