| 4 | 1/1 | 返回列表 |
| 查看: 1029 | 回復(fù): 3 | |||
qzhaosdu金蟲(chóng) (著名寫(xiě)手)
|
[交流]
【討論】關(guān)于G03里guess=only的一個(gè)疑惑。
|
|
G03里GUESS關(guān)鍵詞里說(shuō)到guess=alter的用二重態(tài)NH2舉例。 首先,用了guess=only得出初始猜測(cè),說(shuō)明書(shū)里是這樣的 Initial guess orbital symmetries. Alpha Orbitals: Occupied (A1) (A1) (B2) (B1) (A1) Virtual (A1) (B2) Beta Orbitals: Occupied (A1) (A1) (B2) (A1) Virtual (B1) (A1) (B2) 但是我算了一下卻是 Orbital symmetries: Alpha Orbitals: Occupied (A1) (A1) (B2) (A1) (B1) Virtual (A1) (B2) Beta Orbitals: Occupied (A1) (A1) (B2) (A1) Virtual (B1) (A1) (B2) 也就是說(shuō)說(shuō)明書(shū)里的占據(jù)軌道里4和5是反著的。 但是用了OPT UHF去做的時(shí)候,程序在算第一步SCF的時(shí)候卻給的是和說(shuō)明書(shū)一樣的結(jié)果。 正常來(lái)說(shuō),Alpha和Beta的對(duì)稱性應(yīng)該是一一對(duì)應(yīng)相同的吧?(我也不知道這樣說(shuō)對(duì)不對(duì)?) ************************** 另外我剛看了下G09的說(shuō)明書(shū),GUESS里給的是正常的結(jié)果,也就是說(shuō)G03里寫(xiě)錯(cuò)了? ************************** 另外還有,我算了一個(gè)Ni+...C2H2的C2V結(jié)構(gòu),guess=only給的猜測(cè)是2-B2,但是b3lyp/6-31+g*算的時(shí)候給的初始猜測(cè)卻是2-B1. ************************** 所以說(shuō),G03里這個(gè)guess=only這個(gè)關(guān)鍵詞是不是有問(wèn)題?(起碼第一個(gè)SCF是有錯(cuò)誤的) [ Last edited by qzhaosdu on 2010-12-9 at 17:50 ] |
» 搶金幣啦!回帖就可以得到:
+5/480
+1/89
+1/85
+1/78
+1/67
+1/56
+1/40
+1/37
+1/30
+2/26
+1/17
+1/13
+1/9
+1/7
+1/5
+1/5
+1/4
+1/3
+1/1
+1/1
金蟲(chóng) (正式寫(xiě)手)
|
第一個(gè)問(wèn)題:不同版本的程序所用的初始猜測(cè)方法可能有所不同(即使方法相同也可能參數(shù)或算法有所不同)。不同的方法給出的初始猜測(cè)不一定相同,但只要方法合適,優(yōu)化后的結(jié)果應(yīng)該是一致的。 第二個(gè)問(wèn)題:“正常來(lái)說(shuō)占據(jù)軌道和非占據(jù)軌道的對(duì)稱性應(yīng)該是一一對(duì)應(yīng)相同的吧?”這一理解是錯(cuò)誤的!首先因?yàn)檎紦?jù)軌道的數(shù)目與空軌道的數(shù)目一般是不同的,故不可能建立一一對(duì)應(yīng)關(guān)系。其次,即使二者的數(shù)目相同,也沒(méi)有任何物理或幾何原因能在二者之間建立這種關(guān)系。 第三個(gè)問(wèn)題:有些開(kāi)殼層體系特別是自由基,初始猜測(cè)的對(duì)稱性常常是與基組有關(guān)的。原因在于你使用了單Slater行列式近似,而這種近似對(duì)于有些開(kāi)殼層體系并不是好的近似。 |
鐵桿木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (正式寫(xiě)手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +3 | 蔡大美女 2026-03-13 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 314求調(diào)劑 +8 | 無(wú)懈可擊的巨人 2026-03-12 | 8/400 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 一志愿中國(guó)海洋大學(xué),生物學(xué),301分,求調(diào)劑 +3 | 1孫悟空 2026-03-17 | 3/150 |
|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評(píng)估A類院;瘜W(xué)工程與工藝,誠(chéng)求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|