| 5 | 1/1 | 返回列表 |
| 查看: 8861 | 回復: 23 | |||||||||||||||||
| 本帖產生 2 個 模擬EPI ,點擊這里進行查看 | |||||||||||||||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||||||||||||||
[交流]
關于勢函數(shù)的求助
|
|||||||||||||||||
|
相信做分子模擬親們都會用到勢函數(shù)? 但是請問大家的勢函數(shù)參數(shù)如何取得的呢? 我也相信有部分童鞋是查得文獻,同時肯定也有部分童鞋是自己擬合的。 請問是通過什么方法或軟件擬合呢?希望有詳細的步驟和例子。 感謝大家!感謝蟲子們! |
資源收集 | 分子動力學 | 分子模擬 | 模擬知識匯總 |
MS 吸附模擬 | 模型 - 模擬 - 統(tǒng)計 - 熱力學 | 我學習計算的一些帖子 | 動力學模擬 |
模擬軟件安裝精華淘 | 動力學強帖 | 我的分享 | MS學習 |
力場 | 分子動力學 |
» 搶金幣啦!回帖就可以得到:
+1/192
+1/192
+1/184
+1/81
+1/71
+1/37
+1/34
+1/21
+1/15
+1/9
+1/9
+1/8
+1/6
+1/6
+1/5
+1/5
+1/4
+1/3
+1/3
+1/2
|
非常感謝你的回答,說實話,我沒有完全聽懂你的意思,我是搞材料的,對計算化學、計算物理懂一點點,皮毛,嘿嘿。 我先介紹一下目前我的情況吧。我做的材料屬于離子氧化物,用的勢函數(shù)形式是buckingham和shel(暫時不討論),但目前我最關心的是buckingham這個二體勢,我想擬合這個勢的參數(shù)。我從文獻中已經找到了一些勢參數(shù),但是我在驗證的這些參數(shù)的時候,發(fā)現(xiàn)得到的晶格常數(shù)與實驗值相差4%,對這些通過文獻得到的參數(shù)有點不滿,所以我就想自己去擬合這些參數(shù),請問該如何辦呢?材料體系的離子是Bi、Zn、Nb、O。 |
|
MS中的Gulp模塊適合你。 看看手冊吧。我也沒用過。 晶體結構(a,b,c)4%的誤差是稍微有點大。擬合力場晶體結構方面的偏差一般能達到2%以內,密度方面小于5%, 晶格能方面小于10%。因為必須考慮實驗值也是有誤差的。晶格能,蒸發(fā)焓這些實驗值的誤差也在10%左右,刻意弄成跟實驗值完全一致沒有必要,只要跟類似的物質比,結合能定性大小順序符合即可; 另外,你還要問自己,即你要研究的問題是什么?這個問題需要多精確的晶格常數(shù)值?因為有時候你調了參數(shù),看上去與實驗值更接近了,可是一跑動力學,一升高溫度,得到的結果反而更差。比如你現(xiàn)在要預測體積膨脹系數(shù),那我建議你先用現(xiàn)在的參數(shù)跑一些模擬,得到一些結果,看看跟實驗值的誤差有多大。如果在5%以內是能夠接受的。 初學者往往會產生錯覺,以為分子模擬能真實再現(xiàn)所有原子過程并給出精確的實驗結果。但其實經典力場本身是個十分粗糙的模型,它忽略了電子的運動,把原子看成是帶電荷的球,分子間的作用通過一個簡化的數(shù)學模型(Buckingham+coulomb)來表達,沒有多體作用(可達~5%),沒有極化(可達1%),請問你覺得這在多大程度上符合真實的原子結構和相互作用?一般來說, 我們使用分子力學模型并非是為了完全精確預測實驗結果,而是來定性地研究某一性質的變化規(guī)律,或者現(xiàn)象背后的原子過程。即,關注相對值,而非絕對值。 我并不是不建議你調參數(shù),相反我認為做模擬的都要嘗試驗證參數(shù),可能的話通過擬合更多的實驗數(shù)據(jù)來改進參數(shù)。但是在調參數(shù)的時候,要認識到力場模型的局限性,以及它的誤差范圍和精度范圍。想像你拿的是一把最小刻度為cm的尺子,而兩根棍子的長度差別是幾個mm。你能得到精確的結果嗎? [ Last edited by ChemiAndy on 2011-10-26 at 11:32 ] |
|
你好!你的回答讓我茅塞頓開,非常感謝! 我研究的是材料的缺陷,是化學缺陷,不是物理缺陷,舉個例子吧:1、研究晶體中某個原子的取代,某個原子的空缺,某個原子離開自己的位置跑到間隙去了等等一系列的問題,然后計算這個過程產生的能量,再舉個例子吧:2、當我向這個材料體系中去摻雜的時候,我要摻雜的離子(非常微量、基本不影響原材料體系的晶體結構)可能會占據(jù)原體系材料的那個位置,通過模擬計算得出能量值,然后就可以判斷雜質離子可能會進入哪個位置,然后來指導我的實驗。 在這種情況下,我在選取勢參數(shù)(原材料體系的)的時候,驗證參數(shù)的合理性,應該考慮它的哪些因素?比如晶格常數(shù)、彈性模量、密度、結合能等等,你覺得哪個最重要,其實我也知道研究不同的內容,要求的不一樣。也就是說選取的勢參數(shù)不可能同時保證晶格常數(shù)、彈性模量、密度、結合能都和實驗值想符合,在這種情況下,你覺得我在判斷勢參數(shù)合理性時,最應該關注哪一個參考量呢?目前我只考慮了晶格常數(shù)這一個參考量。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學材料與化工(085600)296求調劑 +12 | 稻妻小編 2026-03-26 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 一志愿211 初試270分 求調劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
[考研] 086000生物與醫(yī)藥調劑 +3 | Feisty。 2026-03-28 | 7/350 |
|
|
[考研] 283求調劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 070300求調劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 266分,求材料冶金能源化工等調劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 339求調劑,想調回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 295求調劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 287求調劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 調劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥292求調劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 0854人工智能方向招收調劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 調劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 一志愿國科過程所081700,274求調劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境,總分308求調劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|