關(guān)于Δρ計(jì)算的一些注意點(diǎn)

一個(gè)蟲友學(xué)作Δρ，選擇了HCl來練手，結(jié)果作出來了，發(fā)現(xiàn)與文獻(xiàn)上看到的大不一樣，問我為什么。
因?yàn)椴恢�道他是如何具體操作的，便擬寫了幾點(diǎn)讓他對(duì)照自查，不想寫下來還不少，雖然是簡(jiǎn)單一個(gè)HCl例子，但還是涉及到方方面面，麻雀雖小，肝膽具全！
現(xiàn)貼在這里供學(xué)作Δρ蟲友參考。

據(jù)我的經(jīng)驗(yàn)基本原則是：
1、HCl、H、Cl的坐標(biāo)都取HCl 中的，并用-1和nosymm將它們固定在某一個(gè)坐標(biāo)平面。
2、對(duì)HCl、H、Cl統(tǒng)一自定義網(wǎng)格點(diǎn)數(shù)、網(wǎng)格的起始點(diǎn)和步長(zhǎng)，不能讓程序自選。
3、HCl分子取向的主軸是z還是x？Cl呢？指的是Cl的那個(gè)單電子軌道的取向要與HCl一致。Cl的那個(gè)單電子軌道的取向可以在multiwfn中看到，如果需要可以在這里可以改變?nèi)∠颉?br /> 深入的考慮是：
4、關(guān)于Cl的計(jì)算，如果用ROHF，單電子將放在孤對(duì)電子之上，是HOMO，能級(jí)比H還高。如果用HF計(jì)算，單電子將放在孤對(duì)電子之下，能級(jí)比H低得多。Cl用HF計(jì)算的總能量也比ROHF的低得多。故Cl應(yīng)用HF計(jì)算的數(shù)據(jù)。
5、對(duì)于單原子Cl的計(jì)算，還不能加nosymm，因?yàn)榧恿�，不作�?duì)稱性操作，Cl的單電子取向不在某一個(gè)主軸方向，不能與HCl對(duì)應(yīng)。這時(shí)是先算Cl不加nosymm，計(jì)算自動(dòng)會(huì)將Cl放在坐標(biāo)原點(diǎn)，單電子取向在主軸z方向。再安排HCl在主軸z方向，Cl放在坐標(biāo)原點(diǎn)，并在輸入坐標(biāo)數(shù)據(jù)中在Cl和H之后加-1固定，并加nosymm以免改變了坐標(biāo)系的主軸取向和挪動(dòng)了電荷重心。
6、HCl的Δρ，既可作全電子的，又可以只作相關(guān)軌道的，如果只看成鍵，后者更明晰一些。因?yàn)镠Cl的成鍵，就只是Cl的那個(gè)單電子和H之間的，Cl上的其他電子沒有參加與H的成鍵，MO的安排充分體現(xiàn)了這點(diǎn)。這時(shí)可取HCl的成鍵MO與Cl的那個(gè)單電子MO、H的即成。如果是作全電子的，Δρ還會(huì)包含Cl因吸引了電荷過來，使其上的孤對(duì)電子能級(jí)被抬升，也會(huì)反映到Δρ的圖形上。

[ Last edited by yjcmwgk on 2011-2-10 at 10:45 ] 返回小木蟲查看更多