24小时热门版块排行榜    

查看: 1740  |  回复: 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

新虫 (小有名气)

引用回帖:
7楼: Originally posted by 枪下游魂 at 2014-05-14 08:52:28
# opt b3lyp/gen pseudo=read scf=qc

title

2 3
Pd                -0.47468421   -0.55468477   -0.00857500
N                 -1.68368998    0.15122248   -0.00857500
N                 -1.4913109 ...

我知道可以跑,但是算出来的out文件提示错误,我看了下最终的优化结构也是不对的。
学校,学习ing...
8楼2014-05-14 09:11:50
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +13 rare12345 2026-03-18 13/650 2026-03-19 17:30 by jean5056
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 12/600 2026-03-19 16:17 by 脱颖而出
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 317求调剂 +3 申子申申 2026-03-19 6/300 2026-03-19 14:16 by 申子申申
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 321求调剂(食品/专硕) +3 mxcz321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见