24小时热门版块排行榜    

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

夏虫语冰8

新虫 (小有名气)

[求助] 知道赝势基组Hay-Wadt VDZ (n+1) ECP和原子坐标,怎样写高斯的几何优化输入文件 已有1人参与

赝势基组为Hay-Wadt VDZ (n+1) ECP
#  Hay-Wadt VDZ (n+1) ECP  EMSL  Basis Set Exchange Library  5/12/14 8:35 PM
# Elements                             References
# --------                             ----------
# K - Au: P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 299 (1985).
#         P. J. Hay (Private Communication).
#   


BASIS "ao basis" PRINT
#BASIS SET: (5s,6p,4d) -> [3s,3p,2d]
Pd    S
      2.7870000             -1.6102393              1.3540775        
      1.9650000              1.8489842             -1.6780848        
      0.6243000              0.6037492             -0.8559381        
      0.1496000              0.0000000              1.0200299        
Pd    S
      0.0436000              1.0000000        
Pd    P
      5.9990000             -0.1034910        
      1.4430000              0.7456952        
      0.5264000              0.3656494        
Pd    P
      0.7368000              0.0763285        
      0.0899000              0.9740065        
Pd    P
      0.0262000              1.0000000        
Pd    D
      6.0910000              0.0376146        
      1.7190000              0.5200479        
      0.6056000              0.5706071        
Pd    D
      0.1883000              1.0000000        
END
# Elements                             References
# --------                             ----------
# K - Au: P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 299 (1985).
#         P. J. Hay (Private Communication).
#   


ECP
Pd nelec 28
Pd ul
0    598.3336444             -0.0563177        
1    162.4298290            -20.1288036        
2     51.5714771           -105.8197923        
2     16.4888260            -42.5733345        
2      5.8287656             -3.6165086        
Pd S
0     73.3806304              3.0003651        
1     14.7550438             32.4350093        
2     17.8350204            459.0830383        
2     12.7111477           -868.0629029        
2      9.3292063            514.4726098        
Pd P
0     55.6689877              4.9593099        
1     64.2337771             21.1711029        
2     17.6254952            605.0560092        
2     11.9058155           -726.9641846        
2      8.5100832            396.3274883        
Pd D
0     49.9994728              3.0508745        
1     39.7477547             22.2506580        
2     11.4321366            674.8357698        
2      9.1790080          -1040.8554048        
2      7.5624429            505.9375147        
END
b3lyp The density functional hybrid model Becke3LYP  was used together with the valence triple-fbasis set 6-311+G(d,p) Pd was treated by  Hay-Wadt VDZ (n+1)(ECP)

原子坐标为
2 3
Pd                -0.47468421   -0.55468477   -0.00857500
N                 -1.68368998    0.15122248   -0.00857500
N                 -1.49131094   -1.22207451   -0.00857500
N                  0.79386306    0.73216134    0.00756931
N                  1.33215486   -0.57102727    0.01400269
H                  0.98591474    1.71354265    0.01019095
H                  2.16077145   -1.13075032    0.02423222
H                 -2.28857332    0.94753657   -0.00857500
H                 -1.92092591   -2.12508668   -0.00857500

我已经自己试过写一些输入文件,怎么都不对,请高手帮忙。
不胜感谢。。。
回复此楼
学校,学习ing...
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

引用回帖:
8楼: Originally posted by 夏虫语冰8 at 2014-05-14 09:11:50
我知道可以跑,但是算出来的out文件提示错误,我看了下最终的优化结构也是不对的。...

是的,因为你pd没有指定基组,你开始没讲清楚我以为你是直接就报错了。。。
你先用我这个跑跑看吧,如果不对的话应该是pd的基组定义那块还要修改,ECP那块应该就这样了。
9楼2014-05-14 09:48:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

夏虫语冰8

新虫 (小有名气)

第一次上手高斯就要处理这种比较麻烦的情况,请大家不要见笑,多多帮忙啊。
学校,学习ing...
2楼2014-05-13 11:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
gkf高: 金币+3 2014-05-13 16:12:14
你试写过的输入文件是怎么样的?不对是不对在哪里?高斯如何报错的?
高斯手册gen,genecp这章是否有看过?
不知道你现在会到什么程度,只能先写个大概给你

# opt b3lyp/genecp

分子描述后空一行开始赝势的定义

这块内容略复杂,可以参考大神的随笔:http://hi.baidu.com/sobereva/item/6d3dc1b4639b52f163388e77

可以把你的输入文件贴出来,大家帮你看看。
3楼2014-05-13 13:13:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夏虫语冰8

新虫 (小有名气)

# opt B3LYP GEN scf=qc

title

2 3
Pd                -0.47468421   -0.55468477   -0.00857500
N                 -1.68368998    0.15122248   -0.00857500
N                 -1.49131094   -1.22207451   -0.00857500
N                  0.79386306    0.73216134    0.00756931
N                  1.33215486   -0.57102727    0.01400269
H                  0.98591474    1.71354265    0.01019095
H                  2.16077145   -1.13075032    0.02423222
H                 -2.28857332    0.94753657   -0.00857500
H                 -1.92092591   -2.12508668   -0.00857500

H,N 0
6-31g**
****

Pd     0
S   3   1.00
      2.7870000             -1.6102393        
      1.9650000              1.8489842        
      0.6243000              0.6037492        
S   4   1.00
      2.7870000              1.3540775        
      1.9650000             -1.6780848        
      0.6243000             -0.8559381        
      0.1496000              1.0200299        
S   1   1.00
      0.0436000              1.0000000        
P   3   1.00
      5.9990000             -0.1034910        
      1.4430000              0.7456952        
      0.5264000              0.3656494        
P   2   1.00
      0.7368000              0.0763285        
      0.0899000              0.9740065        
P   1   1.00
      0.0262000              1.0000000        
D   3   1.00
      6.0910000              0.0376146        
      1.7190000              0.5200479        
      0.6056000              0.5706071        
D   1   1.00
      0.1883000              1.0000000        
****
学校,学习ing...
4楼2014-05-13 15:58:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 11/550 2026-03-19 08:11 by 脱颖而出
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 材料专硕274一志愿陕西师范大学求调剂 +6 薛云鹏 2026-03-13 6/300 2026-03-18 14:14 by 脱颖而出
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见