| 9 | 1/1 | 返回列表 |
| 查看: 2864 | 回復(fù): 8 | |||
| 本帖產(chǎn)生 1 個(gè) QC強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | |||
baiyang2000銀蟲 (正式寫手)
|
[求助]
求助Gaussian 09如何進(jìn)行單線態(tài)的開殼層計(jì)算 已有3人參與
|
||
| 請(qǐng)問(wèn)如何用Gaussian 09用DFT方法進(jìn)行開殼層的單線態(tài)計(jì)算(open-shell singlet state),我使用了UB3LYP,但是對(duì)于單線態(tài),加U不加U的結(jié)果是一樣,加了沒(méi)有作用,請(qǐng)問(wèn)這種情況該如何設(shè)置關(guān)鍵詞? |

銀蟲 (正式寫手)
專家顧問(wèn) (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |

木蟲 (小有名氣)
|
對(duì)于開殼和閉殼計(jì)算,是有可能結(jié)果一樣的,不過(guò)兩種計(jì)算的基本原理不一樣。 閉殼是假設(shè)每個(gè)軌道上的兩個(gè)電子(alpha和beta)除了自旋不一樣,其波函數(shù)是一樣的,計(jì)算過(guò)程中只需要算一次(alpha),beta的波函數(shù)等于alpha(這個(gè)可以在計(jì)算完結(jié)果中看軌道分布得到驗(yàn)證)。 但如果是開殼計(jì)算,alpha和beta是分開算的,同一軌道的兩個(gè)電子的波函數(shù)可能不相等。不過(guò)開殼計(jì)算一般是針對(duì)磁性材料 對(duì)于你的情況,存在這種可能,如果你計(jì)算的分子是沒(méi)有磁性,有可能幾遍alpha和beta分開算,你得到波函數(shù)也可能是一樣的。 如果你想驗(yàn)證這一點(diǎn),你可以在計(jì)算中前面加p,后面加gfprint pop=full,你可以在out文件中看到分子軌道的參數(shù)。如果閉殼應(yīng)該只有Molecular Orbital Coefficients這個(gè)數(shù)據(jù)塊, 如果是開殼則有Alpha Molecular Orbital Coefficients 和Beta Molecular Orbital Coefficients,你的情況應(yīng)該是這兩個(gè)數(shù)據(jù)塊是一樣的(或者Density Matrix數(shù)據(jù)塊也是如此)。 |
新蟲 (正式寫手)
專家顧問(wèn) (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |

新蟲 (初入文壇)
版主 (著名寫手)
|
需要用UDFT, 因?yàn)槭莖pen-shell 然后因?yàn)間aussian 當(dāng)設(shè)置singlet 用udft算還是會(huì)默認(rèn)強(qiáng)制paired up。如果要破壞成對(duì)電子 可以用 GUESS=(mix,always) 把HOMO和LUMO 混合起來(lái) 破壞spatial symmetry 就有可能得到open shell singlet 如果不行 就試試 先算一個(gè) triplet出來(lái) 然后用生成的軌道信息和geometry作為初猜 再用上面的命令算一遍 成功的幾率會(huì)大一點(diǎn) |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +4 | 皮皮皮123456 2026-03-22 | 4/200 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|