| 5 | 1/1 | 返回列表 |
| 查看: 8869 | 回復(fù): 23 | |||||||||||||||||
| 本帖產(chǎn)生 2 個(gè) 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | |||||||||||||||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||||||||||||||
[交流]
關(guān)于勢函數(shù)的求助
|
|||||||||||||||||
|
相信做分子模擬親們都會用到勢函數(shù)? 但是請問大家的勢函數(shù)參數(shù)如何取得的呢? 我也相信有部分童鞋是查得文獻(xiàn),同時(shí)肯定也有部分童鞋是自己擬合的。 請問是通過什么方法或軟件擬合呢?希望有詳細(xì)的步驟和例子。 感謝大家!感謝蟲子們! |
資源收集 | 分子動力學(xué) | 分子模擬 | 模擬知識匯總 |
MS 吸附模擬 | 模型 - 模擬 - 統(tǒng)計(jì) - 熱力學(xué) | 我學(xué)習(xí)計(jì)算的一些帖子 | 動力學(xué)模擬 |
模擬軟件安裝精華淘 | 動力學(xué)強(qiáng)帖 | 我的分享 | MS學(xué)習(xí) |
力場 | 分子動力學(xué) |
» 搶金幣啦!回帖就可以得到:
+1/287
+1/184
+1/178
+1/50
+1/39
+1/37
+1/23
+1/16
+1/10
+1/10
+1/10
+1/10
+1/10
+1/9
+1/9
+1/6
+1/6
+1/3
+1/3
+1/2
|
非常感謝你的回答,說實(shí)話,我沒有完全聽懂你的意思,我是搞材料的,對計(jì)算化學(xué)、計(jì)算物理懂一點(diǎn)點(diǎn),皮毛,嘿嘿。 我先介紹一下目前我的情況吧。我做的材料屬于離子氧化物,用的勢函數(shù)形式是buckingham和shel(暫時(shí)不討論),但目前我最關(guān)心的是buckingham這個(gè)二體勢,我想擬合這個(gè)勢的參數(shù)。我從文獻(xiàn)中已經(jīng)找到了一些勢參數(shù),但是我在驗(yàn)證的這些參數(shù)的時(shí)候,發(fā)現(xiàn)得到的晶格常數(shù)與實(shí)驗(yàn)值相差4%,對這些通過文獻(xiàn)得到的參數(shù)有點(diǎn)不滿,所以我就想自己去擬合這些參數(shù),請問該如何辦呢?材料體系的離子是Bi、Zn、Nb、O。 |
|
MS中的Gulp模塊適合你。 看看手冊吧。我也沒用過。 晶體結(jié)構(gòu)(a,b,c)4%的誤差是稍微有點(diǎn)大。擬合力場晶體結(jié)構(gòu)方面的偏差一般能達(dá)到2%以內(nèi),密度方面小于5%, 晶格能方面小于10%。因?yàn)楸仨毧紤]實(shí)驗(yàn)值也是有誤差的。晶格能,蒸發(fā)焓這些實(shí)驗(yàn)值的誤差也在10%左右,刻意弄成跟實(shí)驗(yàn)值完全一致沒有必要,只要跟類似的物質(zhì)比,結(jié)合能定性大小順序符合即可; 另外,你還要問自己,即你要研究的問題是什么?這個(gè)問題需要多精確的晶格常數(shù)值?因?yàn)橛袝r(shí)候你調(diào)了參數(shù),看上去與實(shí)驗(yàn)值更接近了,可是一跑動力學(xué),一升高溫度,得到的結(jié)果反而更差。比如你現(xiàn)在要預(yù)測體積膨脹系數(shù),那我建議你先用現(xiàn)在的參數(shù)跑一些模擬,得到一些結(jié)果,看看跟實(shí)驗(yàn)值的誤差有多大。如果在5%以內(nèi)是能夠接受的。 初學(xué)者往往會產(chǎn)生錯(cuò)覺,以為分子模擬能真實(shí)再現(xiàn)所有原子過程并給出精確的實(shí)驗(yàn)結(jié)果。但其實(shí)經(jīng)典力場本身是個(gè)十分粗糙的模型,它忽略了電子的運(yùn)動,把原子看成是帶電荷的球,分子間的作用通過一個(gè)簡化的數(shù)學(xué)模型(Buckingham+coulomb)來表達(dá),沒有多體作用(可達(dá)~5%),沒有極化(可達(dá)1%),請問你覺得這在多大程度上符合真實(shí)的原子結(jié)構(gòu)和相互作用?一般來說, 我們使用分子力學(xué)模型并非是為了完全精確預(yù)測實(shí)驗(yàn)結(jié)果,而是來定性地研究某一性質(zhì)的變化規(guī)律,或者現(xiàn)象背后的原子過程。即,關(guān)注相對值,而非絕對值。 我并不是不建議你調(diào)參數(shù),相反我認(rèn)為做模擬的都要嘗試驗(yàn)證參數(shù),可能的話通過擬合更多的實(shí)驗(yàn)數(shù)據(jù)來改進(jìn)參數(shù)。但是在調(diào)參數(shù)的時(shí)候,要認(rèn)識到力場模型的局限性,以及它的誤差范圍和精度范圍。想像你拿的是一把最小刻度為cm的尺子,而兩根棍子的長度差別是幾個(gè)mm。你能得到精確的結(jié)果嗎? [ Last edited by ChemiAndy on 2011-10-26 at 11:32 ] |
|
你好!你的回答讓我茅塞頓開,非常感謝! 我研究的是材料的缺陷,是化學(xué)缺陷,不是物理缺陷,舉個(gè)例子吧:1、研究晶體中某個(gè)原子的取代,某個(gè)原子的空缺,某個(gè)原子離開自己的位置跑到間隙去了等等一系列的問題,然后計(jì)算這個(gè)過程產(chǎn)生的能量,再舉個(gè)例子吧:2、當(dāng)我向這個(gè)材料體系中去摻雜的時(shí)候,我要摻雜的離子(非常微量、基本不影響原材料體系的晶體結(jié)構(gòu))可能會占據(jù)原體系材料的那個(gè)位置,通過模擬計(jì)算得出能量值,然后就可以判斷雜質(zhì)離子可能會進(jìn)入哪個(gè)位置,然后來指導(dǎo)我的實(shí)驗(yàn)。 在這種情況下,我在選取勢參數(shù)(原材料體系的)的時(shí)候,驗(yàn)證參數(shù)的合理性,應(yīng)該考慮它的哪些因素?比如晶格常數(shù)、彈性模量、密度、結(jié)合能等等,你覺得哪個(gè)最重要,其實(shí)我也知道研究不同的內(nèi)容,要求的不一樣。也就是說選取的勢參數(shù)不可能同時(shí)保證晶格常數(shù)、彈性模量、密度、結(jié)合能都和實(shí)驗(yàn)值想符合,在這種情況下,你覺得我在判斷勢參數(shù)合理性時(shí),最應(yīng)該關(guān)注哪一個(gè)參考量呢?目前我只考慮了晶格常數(shù)這一個(gè)參考量。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 11ggg 2026-03-30 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +6 | 恬淡ye 2026-03-31 | 7/350 |
|
|
[考研] 289求調(diào)劑 +6 | BrightLL 2026-03-29 | 6/300 |
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 研s. 2026-03-25 | 4/200 |
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +3 | 鐵憨憨123425 2026-03-31 | 3/150 |
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
[考研] 求收留 +6 | 1943443204 2026-03-28 | 6/300 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我!。 +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考博] 26申博 +3 | 加油沖。 2026-03-26 | 3/150 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|