| 24 | 1/1 | 返回列表 |
| 查看: 8795 | 回復: 23 | |||||||||||||||||
| 本帖產(chǎn)生 2 個 模擬EPI ,點擊這里進行查看 | |||||||||||||||||
[交流]
關于勢函數(shù)的求助
|
|||||||||||||||||
|
相信做分子模擬親們都會用到勢函數(shù)? 但是請問大家的勢函數(shù)參數(shù)如何取得的呢? 我也相信有部分童鞋是查得文獻,同時肯定也有部分童鞋是自己擬合的。 請問是通過什么方法或軟件擬合呢?希望有詳細的步驟和例子。 感謝大家!感謝蟲子們! |
資源收集 | 分子動力學 | 分子模擬 | 模擬知識匯總 |
MS 吸附模擬 | 模型 - 模擬 - 統(tǒng)計 - 熱力學 | 我學習計算的一些帖子 | 動力學模擬 |
模擬軟件安裝精華淘 | 動力學強帖 | 我的分享 | MS學習 |
力場 | 分子動力學 |
» 搶金幣啦!回帖就可以得到:
+5/1690
+1/184
+1/46
+2/40
+1/40
+1/31
+1/30
+1/28
+1/17
+1/15
+1/8
+1/8
+1/7
+1/7
+1/6
+1/6
+1/5
+1/4
+1/3
+1/2
| 內(nèi)容已刪除 |
|
非常感謝你的回答,說實話,我沒有完全聽懂你的意思,我是搞材料的,對計算化學、計算物理懂一點點,皮毛,嘿嘿。 我先介紹一下目前我的情況吧。我做的材料屬于離子氧化物,用的勢函數(shù)形式是buckingham和shel(暫時不討論),但目前我最關心的是buckingham這個二體勢,我想擬合這個勢的參數(shù)。我從文獻中已經(jīng)找到了一些勢參數(shù),但是我在驗證的這些參數(shù)的時候,發(fā)現(xiàn)得到的晶格常數(shù)與實驗值相差4%,對這些通過文獻得到的參數(shù)有點不滿,所以我就想自己去擬合這些參數(shù),請問該如何辦呢?材料體系的離子是Bi、Zn、Nb、O。 |
|
MS中的Gulp模塊適合你。 看看手冊吧。我也沒用過。 晶體結(jié)構(gòu)(a,b,c)4%的誤差是稍微有點大。擬合力場晶體結(jié)構(gòu)方面的偏差一般能達到2%以內(nèi),密度方面小于5%, 晶格能方面小于10%。因為必須考慮實驗值也是有誤差的。晶格能,蒸發(fā)焓這些實驗值的誤差也在10%左右,刻意弄成跟實驗值完全一致沒有必要,只要跟類似的物質(zhì)比,結(jié)合能定性大小順序符合即可; 另外,你還要問自己,即你要研究的問題是什么?這個問題需要多精確的晶格常數(shù)值?因為有時候你調(diào)了參數(shù),看上去與實驗值更接近了,可是一跑動力學,一升高溫度,得到的結(jié)果反而更差。比如你現(xiàn)在要預測體積膨脹系數(shù),那我建議你先用現(xiàn)在的參數(shù)跑一些模擬,得到一些結(jié)果,看看跟實驗值的誤差有多大。如果在5%以內(nèi)是能夠接受的。 初學者往往會產(chǎn)生錯覺,以為分子模擬能真實再現(xiàn)所有原子過程并給出精確的實驗結(jié)果。但其實經(jīng)典力場本身是個十分粗糙的模型,它忽略了電子的運動,把原子看成是帶電荷的球,分子間的作用通過一個簡化的數(shù)學模型(Buckingham+coulomb)來表達,沒有多體作用(可達~5%),沒有極化(可達1%),請問你覺得這在多大程度上符合真實的原子結(jié)構(gòu)和相互作用?一般來說, 我們使用分子力學模型并非是為了完全精確預測實驗結(jié)果,而是來定性地研究某一性質(zhì)的變化規(guī)律,或者現(xiàn)象背后的原子過程。即,關注相對值,而非絕對值。 我并不是不建議你調(diào)參數(shù),相反我認為做模擬的都要嘗試驗證參數(shù),可能的話通過擬合更多的實驗數(shù)據(jù)來改進參數(shù)。但是在調(diào)參數(shù)的時候,要認識到力場模型的局限性,以及它的誤差范圍和精度范圍。想像你拿的是一把最小刻度為cm的尺子,而兩根棍子的長度差別是幾個mm。你能得到精確的結(jié)果嗎? [ Last edited by ChemiAndy on 2011-10-26 at 11:32 ] |
|
你好!你的回答讓我茅塞頓開,非常感謝! 我研究的是材料的缺陷,是化學缺陷,不是物理缺陷,舉個例子吧:1、研究晶體中某個原子的取代,某個原子的空缺,某個原子離開自己的位置跑到間隙去了等等一系列的問題,然后計算這個過程產(chǎn)生的能量,再舉個例子吧:2、當我向這個材料體系中去摻雜的時候,我要摻雜的離子(非常微量、基本不影響原材料體系的晶體結(jié)構(gòu))可能會占據(jù)原體系材料的那個位置,通過模擬計算得出能量值,然后就可以判斷雜質(zhì)離子可能會進入哪個位置,然后來指導我的實驗。 在這種情況下,我在選取勢參數(shù)(原材料體系的)的時候,驗證參數(shù)的合理性,應該考慮它的哪些因素?比如晶格常數(shù)、彈性模量、密度、結(jié)合能等等,你覺得哪個最重要,其實我也知道研究不同的內(nèi)容,要求的不一樣。也就是說選取的勢參數(shù)不可能同時保證晶格常數(shù)、彈性模量、密度、結(jié)合能都和實驗值想符合,在這種情況下,你覺得我在判斷勢參數(shù)合理性時,最應該關注哪一個參考量呢?目前我只考慮了晶格常數(shù)這一個參考量。 |
|
你的問題主要是缺陷或取代造成的能量變化,當然能量信息對你來說是最重要的,不僅包括實驗測定的晶格能,而且高級量子化學計算的粒子間相互作用能也是很重要的參考。結(jié)構(gòu)雖然也很重要,但結(jié)構(gòu)再好,能量算不對,一切白搭。 你的體系是“目標晶體+離子參雜”,我不知道你的目標晶體,貌似以靜電相互作用為主。對于這樣的體系,使用于分子力學來進行估算,應該是可行的,相信你已經(jīng)看過類似的文獻了。我覺得文獻中之所以忍受4%的結(jié)構(gòu)誤差,很有可能是為了得到更好的能量數(shù)據(jù)。純屬猜測,你再仔細看看文獻中參數(shù)是怎么來的。力場是模型,參數(shù)是整個研究是否合理的關鍵。 |
| 24 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +12 | zjl的號 2026-03-16 | 17/850 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 生物學調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 化學求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿蘇州大學材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 中南大學化學學碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] A區(qū)線材料學調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|