| 查看: 1748 | 回復: 13 | ||
夏蟲語冰8新蟲 (小有名氣)
|
[求助]
知道贗勢基組Hay-Wadt VDZ (n+1) ECP和原子坐標,怎樣寫高斯的幾何優(yōu)化輸入文件 已有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 我已經(jīng)自己試過寫一些輸入文件,怎么都不對,請高手幫忙。 不勝感謝。。。 |

新蟲 (小有名氣)

木蟲 (著名寫手)
|
你試寫過的輸入文件是怎么樣的?不對是不對在哪里?高斯如何報錯的? 高斯手冊gen,genecp這章是否有看過? 不知道你現(xiàn)在會到什么程度,只能先寫個大概給你 # opt b3lyp/genecp 分子描述后空一行開始贗勢的定義 這塊內(nèi)容略復雜,可以參考大神的隨筆:http://hi.baidu.com/sobereva/item/6d3dc1b4639b52f163388e77 可以把你的輸入文件貼出來,大家?guī)湍憧纯础? |
新蟲 (小有名氣)
|
# 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 **** |

新蟲 (小有名氣)

新蟲 (小有名氣)

木蟲 (著名寫手)
|
# 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跑你本來的輸入文件也能跑。。。 |
新蟲 (小有名氣)

木蟲 (著名寫手)
新蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +4 |  ̄^ ̄゜汗 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南大,0703化學,分數(shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 南京大學化學376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 一志愿西安交通大學 學碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時間沒變,給編輯又寫郵件了,沒回復
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|