| 8 | 1/1 | 返回列表 |
| 查看: 2538 | 回復(fù): 7 | |||
pleasant1789木蟲 (初入文壇)
|
[求助]
G09中計算AIM計算過渡金屬原子和氫原子間作用的問題
|
|
以前從沒進(jìn)行過AIM計算,這次投的文章里審稿人提出用AIM計算一下,應(yīng)該是想用AIM計算證實過渡金屬原子和氫原子間存在相互作用。 從小木蟲上現(xiàn)有的與AIM計算相關(guān)的文章里看,http://www.gaoyang168.com/bbs/viewthread.php?tid=1767755 這個帖子里給出了在Gaussian中計算AIM的方法,用AIM=all關(guān)鍵字實現(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è)軟件。 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版本其實是可以產(chǎn)生可用的wfx文件 4.有人提出高斯中計算AIM不可信? 我的問題總結(jié)來說就是 1. 能否用Gaussian09計算這個34個重原子含過渡金屬Ni,且使用ECP的體系? 2. 如果不能,是否有可用的免費(fèi)軟件且不太麻煩的可以進(jìn)行計算? 3.具體如何操作和讀取相關(guān)信息? 問題非,嵥楹吐闊奈煌誓軌虿涣呋卮穑浅8兄x! |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
無論高斯自帶的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)驗: +458 |

|
你好啊,我也最近在研究AIM問題,我看了一篇文獻(xiàn)感覺跟您問的問題有點(diǎn)相像,他們也用AIM分析了過渡金屬與氫原子形成的建的電子密度還有拉普拉量希望對你有所幫助。 http://g.zhubajie.com/ |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|