| 2 | 1/1 | 返回列表 |
| 查看: 1152 | 回復(fù): 1 | |||
shuaizhaoz新蟲(chóng) (初入文壇)
|
[求助]
Multiwfn畫DFT的電子密度圖時(shí),格點(diǎn)(grid point)處的值是如何得到的? 已有1人參與
|
|
新手使用multiwfn后,有一些疑問(wèn)希望能夠搞清楚,grid point(default: 200,200)是可調(diào)的參數(shù),那么很有可能是通過(guò)插值的辦法來(lái)得到,那么fchk文件中是否有電子密度的參考格點(diǎn)的值? 如果有,要如何在fchk中查看呢? 如果沒(méi)有,是如何計(jì)算得到的呢? 我采用的泛函是B3LYP,基組是6-31g,發(fā)現(xiàn)DFT積分格點(diǎn)是一個(gè)可調(diào)的參數(shù),例如int=ultrafine或superfine等,但是不太清楚具體空間的格點(diǎn)數(shù)量/密度,不像VASP可以通過(guò)NGX,NGY,NGZ直接指定,所以有些困惑,希望大神來(lái)解答一下。 |
版主 (著名寫手)
|
高斯的積分格點(diǎn)可以直接用數(shù)字來(lái)控制. 不過(guò)并不是通過(guò)笛卡爾的形式而是極坐標(biāo)形式。也就是定義徑向(radial)和 角度(angular)的格點(diǎn)。 用int 關(guān)鍵詞來(lái)控制。 int(Grid=99302) 就代表了 徑向有99個(gè) 而角度是302個(gè)格點(diǎn)?偢顸c(diǎn)數(shù)=99*302 個(gè) 也就是差不多3萬(wàn)個(gè) 。 而默認(rèn)的int=ultrafine 實(shí)際上是 int(Grid=99590). 高斯的有一些固定的格點(diǎn)級(jí)別可以用關(guān)鍵詞來(lái)調(diào)用 從低到高 Coarse(33,110)< SG1 (50,194)< Fine (99,302)<UltraFine (99,590)< SuperFine (250,974), 它對(duì)H和He原子是(175,974) 不想用上述的默認(rèn)設(shè)置好的格點(diǎn)關(guān)鍵詞,就用int(Grid=xxxxxx) 來(lái)自己設(shè)置。 除此之外,高斯在做freq的時(shí)候會(huì)用到CPHF,這里也需要格點(diǎn)。CPHF的格點(diǎn)在進(jìn)行普通freq的時(shí)候會(huì)默認(rèn)比Int的格點(diǎn)低兩個(gè)級(jí)別比如默認(rèn) int=ultrafine 的時(shí)候 對(duì)于freq的CPHF來(lái)說(shuō)格點(diǎn)就是 SG1 級(jí)別的 ,某些特殊任務(wù)下的cphf會(huì)高一級(jí)達(dá)到fine 比如 計(jì)算非諧振之類的。 而設(shè)置CPHF格點(diǎn)使用 CPHF(Grid=xxxxxx) |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|