 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 4819 | 回復: 2 | ||
[求助]
ADCH電荷與 Hirshfeld
|
|
看了sob等提出的ADCH 電荷的文章以及MultiWfn，想試著用一下。 有幾個問題： 1。ADCH電荷的計算中，需要提供wfn文件。對于有贗勢的原子，wfn中好像沒有相關(guān)信息，對這樣的體系，那么ADCH電荷還可以用嗎？ 2。MultiWfn自帶了算wfn的過程，請問計算使用的方法與基組與算整個分子時用的計算方法不同，這個影響很小，是嗎？另外計算原子時我看很多就用閉殼層算的，不是原子實際上的基態(tài)，這個對電荷分析過程也是影響很小，是嗎？ 3。 G09算Hirshfeld電荷，如果給關(guān)鍵詞HirshfeldEE，就會直接給出QQ，QD以及DD相互作用矩陣以及總的相互作用能，這是不是說明G09中的Hirshfeld方法中，已經(jīng)包含了原子Dipole的相互作用？ 參考文獻 [1] Lu, T.; Chen, F. W. Atomic Dipole Moment Corrected Hirshfeld Population Method, Journal of Theoretical & Computational Chemistry 2012, 11, 163-183. [2] 盧天, 陳. 原子電荷計算方法的對比, 物理化學學報 2012, 28, 1-18. [3] Lu, T.; Chen, F. W. Multiwfn: A Multifunctional Wavefunction Analyzer, J. Comput. Chem. 2012, 33, 580-592. |
回復此樓» 猜你喜歡
深圳理工大學招聘AI4Science與分子模擬方向科研人員
已經(jīng)有0人回復
-
全獎博士 英國利物浦大學 × 臺灣清華大學 聯(lián)合培養(yǎng)
已經(jīng)有0人回復
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有50人回復
-
尋博導
已經(jīng)有13人回復
-
基金正文30頁指的是報告正文還是整個申請書
已經(jīng)有9人回復
-
哈爾濱工業(yè)大學（深圳）綠色化學團隊招收2026年秋季入學博士生
已經(jīng)有11人回復
-
澳大利亞昆士蘭科技大學PhD博士獎學金
已經(jīng)有0人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 用什么軟件做電荷分布比較準確一點，請高手指點！
已經(jīng)有7人回復
- 關(guān)于VDD中原子的WFN文件
已經(jīng)有7人回復
- Hirshfeld 電荷的定義是什么
已經(jīng)有3人回復
- 【求助成功】Mulliken和Hirshfeld 電荷分析的區(qū)別
已經(jīng)有7人回復
2樓2013-04-02 11:44:27
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應助: 226 (大學生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學
【答案】應助回帖
★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
zhou2009: 金幣+6 2013-04-02 17:05:36
感謝參與，應助指數(shù) +1
zhou2009: 金幣+6 2013-04-02 17:05:36
|
不對稱的問題是原子的波函數(shù)文件對應的密度沒有球?qū)ΨQ化所導致的。所以我特別強調(diào)要對原子波函數(shù)先進行球?qū)ΨQ化處理（此文后半部分有具體說明：《使用Multiwfn作電子密度差圖》http://hi.baidu.com/sobereva/item/96eaeb5dd1b12309e7c4a5c1） 球?qū)ΨQ化是個難搞的問題。雖然用Multiwfn調(diào)用Gaussian計算原子波函數(shù)文件會自動做球?qū)ΨQ化，但是這前提是得到的組態(tài)必須正確（諸如B3LYP/6-31G*肯定沒問題），然而當某些方法和基組相組合時，得到的組態(tài)可能不對，因此Multiwfn處理后并沒有真正將密度球?qū)ΨQ化。這時得手動去修改wfn文件，然后對它作個密度的等值面圖，看看是否很像球形。 球?qū)ΨQ化問題主要是對過渡元素而言，對于主族元素處理起來很簡單。 所以結(jié)果不對稱不是Hirshfeld或者ADCH方法的問題，而是具體運算上的問題。 以后可能會給Multiwfn增加個功能，直接用s型高斯函數(shù)自動擬合密度，這樣就總能確保球?qū)ΨQ了。其實高斯計算Hirshfeld電荷時用的就是類似這樣的事先擬合好的密度，但是卻沒有說清楚密度是什么級別下得到的，如果和當前計算條件相差較多，那么Hirshfeld電荷并不準確。 1 贗勢可以用。將自行提供的贗勢下計算（并且自行將之球?qū)ΨQ化好了）的原子波函數(shù)文件放到atomwfn文件夾下即可，Multiwfn會自動使用它。也可以直接讓Multiwfn在相應贗勢下自動計算（不過用使用贗勢時不敢保證能正確自動地球?qū)ΨQ化）。 2 方法、基組不一定非得嚴格一致，大概在一個水平就行了，影響不大。有的的確不是自動在基態(tài)算的，諸如碳，Multiwfn自動用的是sp3組態(tài)算的，因為這樣密度直接就是滿足球?qū)ΨQ的。實際上，以原子的什么狀態(tài)為參考態(tài)，這有很大任意性。 3 G09并沒有把原子偶極矩等效地體現(xiàn)到原子電荷里。確實G09會輸出原子偶極矩之類的以及一些相互作用數(shù)據(jù)，但我曾驗證過，G09輸出的原子偶極矩是錯的。至于目前最新的版本是否修正了我沒檢驗過。 |
3樓2013-04-02 15:53:53
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yanrding 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研]
 295復試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 0703化學 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|