| 查看: 1737 | 回复: 13 | ||
[求助]
知道赝势基组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 我已经自己试过写一些输入文件，怎么都不对，请高手帮忙。 不胜感谢。。。 |
回复此楼» 猜你喜欢
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
南大化学调剂
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有118人回复
-
材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流
已经有5人回复
-
26申博
已经有2人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
-
0703调剂
已经有2人回复
-
1
已经有19人回复
-
0703 调剂
已经有4人回复
-
西南某211招硕士、博士
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
2楼2014-05-13 11:57:26
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
gkf高: 金币+3 2014-05-13 16:12:14
感谢参与，应助指数 +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
|
# 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 **** |
4楼2014-05-13 15:58:39
5楼2014-05-13 15:59:27
6楼2014-05-13 16:06:59
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
夏虫语冰8: 金币+200, ★★★★★最佳答案, 非常感谢。 2014-05-15 14:47:54
夏虫语冰8: 金币+200, ★★★★★最佳答案, 非常感谢。 2014-05-15 14:47:54
|
# 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.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 **** Pd 0 pd-ecp 3 28 Pd ul 5 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 5 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 5 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 5 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 我把你的文件改成这样是可以运行了，但是我不知道结果是否可靠。 另外你的高斯是03还是09，我用09跑你本来的输入文件也能跑。。。 |
7楼2014-05-14 08:52:28
8楼2014-05-14 09:11:50
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
9楼2014-05-14 09:48:19
10楼2014-05-14 14:20:58
