24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1803  |  回复: 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的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +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的回帖
查看全部 14 个回答

夏虫语冰8

新虫 (小有名气)

第一次上手高斯就要处理这种比较麻烦的情况,请大家不要见笑,多多帮忙啊。
学校,学习ing...
2楼2014-05-13 11:57:26
已阅   回复此楼   关注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的回帖

夏虫语冰8

新虫 (小有名气)

我自己写出来时这样的,但是高斯一直报错,
学校,学习ing...
5楼2014-05-13 15:59:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 接受任何调剂 +4 也就是栗子 2026-04-17 4/200 2026-04-17 17:57 by Equinoxhua
[考研] 320求调剂 +4 深郊akm 2026-04-17 4/200 2026-04-17 17:32 by 一切OK
[考研] 急需调剂 +7 绝不放弃22 2026-04-15 7/350 2026-04-16 22:09 by SUSE_CL
[考研] 0831生医工第一轮调剂失败求助 +14 小熊睿睿_s 2026-04-11 18/900 2026-04-16 20:41 by lpl364211
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 7/350 2026-04-16 14:31 by zhouxiaoyu
[考博] 申博自荐 +3 Linxia林夏 2026-04-13 3/150 2026-04-16 12:55 by 墨荷之露
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 求助调剂,跨调 +19 X十甫寸Y 2026-04-11 20/1000 2026-04-15 21:18 by cuisz
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 调剂结束 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
信息提示
请填处理意见