 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1001 | 回復(fù): 0 | |||
[交流]
分子動力學(xué)模擬心得
|
|
適合新手1.看完《分子模擬從算法到應(yīng)用》那本書的第四章，不用全看完，但是至少要對分子動力學(xué)模擬過程有一個了解。2.試著按照書的過程做個Ar的NVE，其實Ar和離子晶體以及其它的任何材料的差別僅僅是勢函數(shù)的問題，雖然由勢函數(shù)帶來了一些問題，但是這些都不是本質(zhì)問題。3.從初始化的原子數(shù)，原子位置，初始速度，時間步長，初始溫度等等這些初始化結(jié)束了以后，選擇一個簡單的積分算法，如6階的Gear預(yù)測校正，不要控溫控壓，就是一個簡單的NVE，不要考慮任何的提高效率的鄰位算法，因為這個時候我們可以選擇5×5×5的超原胞，總共的原子數(shù)也就500個，不需要考慮鄰位算法。4.開始循環(huán)計算：預(yù)測----計算原子的力和能量—校正5.輸出能量這樣最簡單的NVE就編成了，總共也就1000多行，是個很小的程序。自己先試著感受一下。其實當(dāng)這樣的小程序完成以后，你會覺得分子動力學(xué)編程也很簡單，那么接下來的復(fù)雜的分子動力學(xué)也不會是什么問題。做完了這一些，你需要知道的是那些是和材料無關(guān)的東西，那么就盡量的分離，開始使用一個個的函數(shù)。例如，原子的位置是和具體的材料相關(guān)，但是初始速度卻和材料一點關(guān)系都沒有，同樣的數(shù)值積分中的預(yù)測和校正也是和材料無關(guān)的，以后的控溫和控壓算法也是和材料無關(guān)的。當(dāng)規(guī)模大了以后，鄰位算法也是和材料無關(guān)的，像這些和材料無關(guān)的部分最好自己做成小函數(shù)，選擇調(diào)用。以后換材料的時候程序也不會有太大的改動。編程的一些小技巧：1.選擇用intel編譯器，個人喜歡用10.1或者9.1的版本，打開優(yōu)化選擇，類似的/QaxS /QxS /Qipo /Qprec-div-等等，以后可以使用Openmp的并行計算（具體的可以參考intel編譯器的幫助手冊）2.盡量的簡化計算，例如2×a就要寫成a+a，在計算機中，加減是一個數(shù)量的計算，乘是一個，除是一個。計算量逐漸增加，所以選擇a+a來替代2×a 會減少一些計算量3.選擇數(shù)據(jù)來替代結(jié)構(gòu)體，結(jié)構(gòu)體看的比較方便，但是計算效率要低，尤其是在編譯過程中的矢量化的時候，而數(shù)組則可以很好的矢量化，也更適合并行。結(jié)構(gòu)體的計算是先找到結(jié)構(gòu)體指針，然后再找里面的參數(shù)，當(dāng)并行的時候，多線程同時找結(jié)構(gòu)體指針，會很大的減低速度。4.除法盡量用乘法來替代，有了程序以后，自己可以仔細的分析里面的計算消耗，這樣可以更好的優(yōu)化計算，當(dāng)然最主要的消耗是在力和能量的計算 版權(quán)聲明：本文為CSDN博主「hdpai2018」的原創(chuàng)文章，遵循 CC 4.0 BY-SA 版權(quán)協(xié)議，轉(zhuǎn)載請附上原文出處鏈接及本聲明。原文鏈接：https://blog.csdn.net/hdpai2018/article/details/84250685 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 猜你喜歡
德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復(fù)
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試
已經(jīng)有11人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有206人回復(fù)
-
2026-CJ開始申報了
已經(jīng)有1人回復(fù)
-
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 瘋不覺KG 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+5
|
15779376950 2026-03-01 | 7/350 |
|
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +7 | 土土小蟲 2026-03-03 | 7/350 |
|
|
|
[考研] 中國科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑，有科研和獲獎經(jīng)歷 +9 | wsxw 2026-03-02 | 10/500 |
|
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +10 | 簡單點0 2026-03-02 | 14/700 |
|
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
|
[考研] 26考研報考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
|
[考研] 281電子信息求調(diào)劑 +3 | jhtfeybgj 2026-03-02 | 6/300 |
|
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +6 | 愛擦汗的可樂冰 2026-02-28 | 7/350 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|