| 8 | 1/1 | 返回列表 |
| 查看: 2539 | 回復(fù): 7 | |||
pleasant1789木蟲 (初入文壇)
|
[求助]
G09中計算AIM計算過渡金屬原子和氫原子間作用的問題
|
|
以前從沒進(jìn)行過AIM計算,這次投的文章里審稿人提出用AIM計算一下,應(yīng)該是想用AIM計算證實(shí)過渡金屬原子和氫原子間存在相互作用。 從小木蟲上現(xiàn)有的與AIM計算相關(guān)的文章里看,http://www.gaoyang168.com/bbs/viewthread.php?tid=1767755 這個帖子里給出了在Gaussian中計算AIM的方法,用AIM=all關(guān)鍵字實(shí)現(xiàn),然而我計算的時候出錯,出錯信息為“ Error termination in NtrErr: NtrErr Called from FileIO.” 換用較低的基組,則可以進(jìn)行到輸出“ I. ATTRACTORS”的部分,然后“FATAL ERROR: DUPLICATE CRITICAL POINTS FOUND: 50 41 ”609錯誤結(jié)束 計算的體系是IMes配位的Ni與另外一個分子之間的配位,重原子數(shù)34。對Ni原子使用了lanl2dz的ECP。 搜解決原因的時候看到各種說法 1. 計算軟件不合適?有人提到Gaussian的AIM有問題,故應(yīng)該用專業(yè)的AIM軟件計算,然而免費(fèi)的AIMPac需要編譯(這個太復(fù)雜,搞不定);AIMall免費(fèi)的只能算10個以下的原子,AIM2000也是商業(yè)軟件,故都不可行,平時并不用AIM分析所以也不可能因?yàn)檫@個去購買商業(yè)軟件。 2.“芳香環(huán)片段高斯不可以做,會出錯,高斯不可以算環(huán)狀結(jié)構(gòu),可以使用aimall做 ”出自 http://www.gaoyang168.com/bbs/viewthread.php?tid=1781656 是否為真?如果是如此那就真是沒辦法了,是否有相關(guān)的文獻(xiàn)說明這點(diǎn)呢。 3.也有提出ECP不能解決,但是http://www.gaoyang168.com/bbs/viewthread.php?tid=2511503 提到“所有AIM軟件都能支持ECP計算產(chǎn)生的wfn文件”,而且G09 的B01版本其實(shí)是可以產(chǎn)生可用的wfx文件 4.有人提出高斯中計算AIM不可信? 我的問題總結(jié)來說就是 1. 能否用Gaussian09計算這個34個重原子含過渡金屬Ni,且使用ECP的體系? 2. 如果不能,是否有可用的免費(fèi)軟件且不太麻煩的可以進(jìn)行計算? 3.具體如何操作和讀取相關(guān)信息? 問題非,嵥楹吐闊,望哪位同仁能夠不吝回答,非常感謝! |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
無論高斯自帶的AIM模塊能不能正常分析某個體系,我都不建議用,這個模塊年久失修,速度、可靠性、穩(wěn)定性都比較低。 如果不需要做盆積分,使用Multiwfn做AIM分析是最佳選擇,操作與AIM2000一樣方便,而速度比AIMALL還要快得多。 對于過渡金屬使用贗勢的情況,對于LANL2這樣小核贗勢,直接用各種AIM軟件去分析價層性質(zhì)都是可以的,例如Journal of Computational Chemistry, Vol. 18, No. 3, 416]429 (1996)文中給出了確切的結(jié)論(SC代表small core): For transition metals, the electron density obtained from a SC-ECP calculation is good enough to reproduce not only the correct topology of the metal]ligand bonds but the values of the local properties as well 在G09 B01及之后版本中支持wfx格式,會自動將事先擬合好的內(nèi)核密度寫進(jìn)去,使得分析結(jié)果和全電子計算十分一致,見J. Phys. Chem. A 2011, 115, 12879–12894。目前只有Multiwfn和AIMALL支持wfx格式下的拓?fù)浞治。而至于G09自帶的AIM模塊,我雖然沒試過,但是估計它自己都不支持內(nèi)核密度的處理。 |
木蟲 (初入文壇)
禁蟲 (文壇精英)
密度泛函·小卒
木蟲 (初入文壇)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
用Multiwfn做AIM分析的例子在手冊第四章就有,也可以參考專題帖子 http://www.gaoyang168.com/bbs/viewthread.php?tid=3890042&fpage=1 目前引用Multiwfn時引JCC上的介紹文章的doi號或者引網(wǎng)址都可以,等那文章正式期號出現(xiàn)后就建議直接引文章了。 之前我的回復(fù)用不著評分 |
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |

|
你好啊,我也最近在研究AIM問題,我看了一篇文獻(xiàn)感覺跟您問的問題有點(diǎn)相像,他們也用AIM分析了過渡金屬與氫原子形成的建的電子密度還有拉普拉量希望對你有所幫助。 http://g.zhubajie.com/ |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)工程321分求調(diào)劑 +18 | 大米飯! 2026-03-15 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 8/400 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +20 | 不知道叫什么! 2026-03-15 | 22/1100 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 298求調(diào)劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +3 | 沿岸有貝殼6 2026-03-21 | 3/150 |
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|