 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6679 | 回復: 33 | ||||||||||||||||
cenwanglai榮譽版主 (知名作家)
老和山猥瑣派九段
|
[交流]
【活動】vasp計算分子動力學討論合集及專貼（籌劃中） 已有26人參與
|
|||||||||||||||
|
vasp大部分用來做結(jié)構(gòu)弛豫和電子結(jié)構(gòu)計算。它在分子動力學及模擬退火方面的計算目前關(guān)注度不夠。這可能是因為vasp的操作界面不好，以及功能有限。但它本身作為一個成熟軟件，其分子動力學的功能，以及vasp DFT計算結(jié)果在動力學上的挖掘應用應該值得關(guān)注。 我會用一周的時間整理版內(nèi)的相關(guān)帖子。同時希望各位專家、版主、潛伏在版內(nèi)的牛牛們跟帖分享你的收藏、學識和智慧，給大家普及提高一下vasp在分析動力學等非結(jié)構(gòu)弛豫和電子結(jié)構(gòu)結(jié)算方面的應用心得體會，經(jīng)驗和理論。 對于有價值的貢獻，金幣應該不會少吧？  [ Last edited by cenwanglai on 2010-11-5 at 16:57 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 仿真建模與計算 | 研究生資料集 | VASP | VASP使用技巧 |
| 計算-vasp | 第一性原理和電化學 | first principle | VASP計算資料收集byZF |
| vasp 學習貼 | vasp從結(jié)構(gòu)優(yōu)化到光學性質(zhì) | 分子模擬專題 | 科研 |
| 分子模擬 |
» 本帖已獲得的紅花（最新10朵）
海風的說» 猜你喜歡
332求調(diào)劑
已經(jīng)有5人回復
-
085701求調(diào)劑初試286分
已經(jīng)有3人回復
-
材料工程294求調(diào)劑
已經(jīng)有4人回復
-
本科新能源科學與工程，一志愿華理能動285求調(diào)劑
已經(jīng)有10人回復
-
070305高分子化學與物理 304分求調(diào)劑
已經(jīng)有5人回復
-
083000學碩274求調(diào)劑
已經(jīng)有8人回復
-
一志愿北京化工大學材料與化工（085600）296求調(diào)劑
已經(jīng)有10人回復
-
求調(diào)劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過
已經(jīng)有11人回復
-
071000生物學求調(diào)劑，初試成績343
已經(jīng)有7人回復
-
0703化學調(diào)劑，求導師收
已經(jīng)有8人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- vasp分子動力學NVE，能量出現(xiàn)震蕩峰是為什么？
已經(jīng)有9人回復
- VASP分子動力學優(yōu)化的INCAR如下，怎么CONTCAR和POSCAR沒有什么變化��？
已經(jīng)有13人回復
- 【求助】VASP做分子動力學模擬NVT系綜下為啥內(nèi)能竟然是升高的？
已經(jīng)有8人回復
- CASTEP和Vasp分子動力學計算，哪個運算的快些？
已經(jīng)有3人回復
- 【任務征集】VASP跑分子動力學數(shù)據(jù)處理方法
已經(jīng)有20人回復
- 【求助】關(guān)于vasp進行分子動力學模擬的求助
已經(jīng)有12人回復
- 【求助】請教vasp計算恒溫第一原理分子動力學
已經(jīng)有9人回復
- 【求助】第一原理的分子動力學的材料學計算，軟件選擇VASP還是MS？或者其他
已經(jīng)有17人回復
- 【求助】vasp分子動力學模擬求助
已經(jīng)有15人回復
- 【原創(chuàng)】vasp的分子動力學模擬
已經(jīng)有44人回復
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
8樓2010-11-05 16:30:56
cenwanglai
榮譽版主 (知名作家)
老和山猥瑣派九段
- 1ST強帖: 5
- 應助: 46 (小學生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 仿真模擬
|
占位 http://www.gaoyang168.com/html/201006/2086568.html 作者: goldenfisher 收錄日期: 2010-06-04 發(fā)布日期: 2010-05-25 很長一段時間以來，一直有人在問BOMD與CPMD究竟有什么不同？這里我就簡單說一些二者的區(qū)別。 實際上，我想大部分人理解和接觸的第一原理分子動力學方法以CPMD居多。CPMD，就是Car和Parrinello兩個人作出的基于密度泛函的分子動力學方法，其特點是在引入電子虛擬質(zhì)量，將電子運動耦合到了運動方程中，每一步分子動力學計算后，對電子結(jié)構(gòu)的計算就不再需要自洽場迭代這一過程，因此可以大大節(jié)省計算資源，在計算機還不是那么好的80，90年代是彌足珍貴的。CPMD的建立開創(chuàng)了第一原理分子動力學方法的新時代，使其真正開始了實用化進程。 而BOMD，顧名思義，就是Born-Oppenheimer分子動力學。Born-Oppenheimer近似，也就是絕熱近似，指的是將電子和離子的求解分離開來，只處理離子的動力學部分，而認為電子可以快速跟上電子的運動。其特點是每一步分子動力學計算之后都需要對電子結(jié)構(gòu)進行自洽場迭代，使電子達到基態(tài)。正式由于這個自洽場迭代過程需要的計算量巨大，致使其一直沒有達到廣泛運用，直到90年代中后期，計算機技術(shù)的發(fā)展，才使BOMD開始逐步被人們重視。 CPMD和BOMD各有優(yōu)劣。CPMD雖然計算速度更快，不需要進行自洽場迭代計算，計算量小，但是由于計算中并沒有使體系真正達到基態(tài)，而只是盡量靠近基態(tài)，因此其準確性對電子的虛擬質(zhì)量這個參數(shù)的選取依賴程度很大，一旦計算參數(shù)不對，得到的體系很可能遠遠偏離真實的勢能面，得不到正確的動力學軌跡。同時，為了保證其盡量靠近基態(tài)，分子動力學時間步長一般選得較小。而BOMD雖然計算量大，但是由于每一步都保證系統(tǒng)達到基態(tài)，因此其分子動力學步長可以取得較大，一般1fs到5fs都有可能。綜合來說，如果能結(jié)合二者的優(yōu)勢，第一原理分子動力學計算效率將大大提高，這也成為近年來其發(fā)展的重要方向。 關(guān)于這部分近年來的發(fā)展，可以參考以下兩篇文獻，以及其參考文獻和引用他們的文獻。 http://link.aps.org/doi/10.1103/PhysRevLett.98.066401 http://iopscience.iop.org/0295-5075/88/2/20001/?rss=2.0 這部分實際上也是一個研究熱點，歡迎大家熱辣討論，并長期關(guān)注。 作者：beyondstar 頂，嚴重支持�。。。。。。。。。。。�！ 作者：zsjan 寫得非常好�。� 作者：valenhou001 另外一個CPMD處理金屬體系有是不大好的。 作者：vasp001 CPMD 用的勢場是 on the fly ,而BOMD 用的勢場是嚴格的第一性原理的勢場，即從能量對粒子位置求一次偏倒得到。 CPMD 處理金屬體系時，引入了熱庫的概念。來表示電子和離子的耦合。即電子在一定的頻率范圍以上認為和離子進行有效的耦合，即在一定的頻率薄層范圍內(nèi)，都是認為電子和離子進行有效耦合的。 作者：enola 謝謝goldenfisher 作者：freer007 將得好，謝謝樓主分享。 作者：xirainbow 牛！；） 作者：goldenfisher 目前大部分CPMD程序只能進行Gamma點的計算，而BOMD可以進行多k點計算。這也是因為CPMD一般處理流體等無對稱性的對象，所以Gamma點可以增加計算速度，且不關(guān)心體系對稱性。 作者：quantumfang 可能不能評價   作者：oulihui666 頂，好貼！  [ Last edited by cenwanglai on 2011-1-11 at 20:31 ] |
2樓2010-11-04 10:54:34
10樓2010-11-05 16:56:29
lihb734
鐵桿木蟲 (職業(yè)作家)
站在計算化學入門的門檻上
- 應助: 132 (高中生)
- 金幣: 7495.4
- 散金: 9378
- 紅花: 74
- 帖子: 3107
- 在線: 1177.3小時
- 蟲號: 448403
- 注冊: 2007-11-01
- 性別: GG
- 專業(yè): 理論和計算化學
13樓2010-11-06 20:47:50
cenwanglai
榮譽版主 (知名作家)
老和山猥瑣派九段
- 1ST強帖: 5
- 應助: 46 (小學生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 仿真模擬
|
占位 Vasp做分子動力學的好處，由于vasp是近些年開發(fā)的比較成熟的軟件，在做電子scf速度方面有較好的優(yōu)勢。 缺點：可選系綜太少。 盡管如此，對于大多數(shù)有關(guān)分子動力學的任務還是可以勝任的。 主要使用的系綜是 NVT 和 NVE。 下面我將對主要參數(shù)進行介紹！ 一般做分子動力學的時候都需要較多原子，一般都超過100個。 當原子數(shù)多的時候，k點實際就需要較少了。有的時候用一個k點就行，不過這都需要嚴格的測試。通常超過200個原子的時候，用一個k點，即Gamma點就可以了。 INCAR： EDIFF 一般來說，用1E-4 或者1E-5都可以，這個參數(shù)只是對第一個離子步的自洽影響大一些，對于長時間的分子動力學的模擬，精度小一點也無所謂，但不能太小。 IBRION=0 分子動力學模擬 IALGO=48 一般用48，對于原子數(shù)較多，這個優(yōu)化方式較好。 NSW=1000 多少個時間步長。 POTIM=3 時間步長,單位fs, 通常1到3. ISIF=2 計算外界的壓力. NBLOCK= 1 多少個時間步長，寫一次CONTCAR，CHG和CHGCAR，PCDAT. KBLOCK=50 NBLOCK*KBLOCK 個步長寫一次 XDATCAR. ISMEAR=-1 費米迪拉克分布. SIGMA =0.05 單位:電子伏 NELMIN=8 一般用6到8, 最小的電子scf數(shù).太少的話,收斂的不好. LREAL=A APACO=10 徑向分布函數(shù)距離, 單位是埃. NPACO=200 徑向分布函數(shù)插的點數(shù). LCHARG=F 盡量不寫電荷密度,否則CHG文件太大. TEBEG=300 初始溫度. TEEND=300 終態(tài)溫度。 不設的話，等于TEBEG. SMASS -3 NVE ensemble;-1 用來做模擬退火。大于0 NVT 系綜。 作者：chuzhaonan 這個要認真學習 作者：wuchenwf 感謝分享 分子動力學模擬的參數(shù)可選擇性也不小，不要完全被某個INCAR所局限，要明白每個參數(shù)是什么意思 作者：study163 這個要看看，不錯 作者：spur 感謝分享。問一下樓主：有辦法得到體系在某一個方向上(e.g. x, y, z)的徑向分布函數(shù)嗎？另外樓主能不能給大家講講SMASS=1，2，3的區(qū)別？ [ Last edited by spur on 2009-12-21 at 13:00 ] 作者：xiaowandouer 樓主太無私了，感動ing。。。。。。。。 收藏起來，以備不時之需 作者：wuchenwf QUOTE: Originally posted by spur at 2009-12-21 12:58: 感謝分享。問一下樓主：有辦法得到體系在某一個方向上(e.g. x, y, z)的徑向分布函數(shù)嗎？另外樓主能不能給大家講講SMASS=1，2，3的區(qū)別？ [ Last edited by spur on 2009-12-21 at 13:00 ] 這個還真沒有聽說過，有一個問題是徑向分布函數(shù)其本身就是一維的物理量 作者：aylayl08  作者：xbelinda QUOTE: Originally posted by spur at 2009-12-21 12:58: 感謝分享。問一下樓主：有辦法得到體系在某一個方向上(e.g. x, y, z)的徑向分布函數(shù)嗎？另外樓主能不能給大家講講SMASS=1，2，3的區(qū)別？ [ Last edited by spur on 2009-12-21 at 13:00 ] SMASS=1,2,3 是沒有區(qū)別的。都是NVT ensemble。 SMASS只要是大于0就是NVT系綜。 另一個問題，我也不會。沒用到過。 作者：gavinliu7390 it is worth to learn. 作者：csfn 樓主能否考慮簡單介紹下vasp做了md后的主要輸出文件（主要結(jié)果）分析  [ Last edited by cenwanglai on 2011-2-24 at 10:51 ] 學習了，但是要指出你的一點小錯誤。 NBLOCK= 1 指的是多少個離子步（時間步長只是對于MD而言，對于其它的計算，NBLOCK也是起控制作用的）寫一次XDATCAR KBLOCK=50 NBLOCK*KBLOCK 個離子步寫一次PCDAT. CONTCAR是每個離子步之后都會寫出來的，但是會用新的把老的覆蓋 CHG是在每10個離子步寫一次，不會覆蓋 CHGCAR是在任務正常結(jié)束之后才寫的。 [ Last edited by cenwanglai on 2011-2-24 at 10:55 ] 原文鏈接： 感謝xbelinda http://www.gaoyang168.com/bbs/viewthread.php?tid=1744011&fpage=0&view=&highlight=&page=1 [ Last edited by cenwanglai on 2011-2-24 at 10:57 ] |
3樓2010-11-04 10:54:46
cenwanglai
榮譽版主 (知名作家)
老和山猥瑣派九段
- 1ST強帖: 5
- 應助: 46 (小學生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 仿真模擬
4樓2010-11-04 10:54:59
mingxing.s
銀蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1184.9
- 紅花: 1
- 帖子: 256
- 在線: 20.7小時
- 蟲號: 1126217
- 注冊: 2010-10-19
- 性別: GG
- 專業(yè): 材料物理化學
5樓2010-11-04 11:48:03
6樓2010-11-04 13:00:57
stractor
金蟲 (著名寫手)
- 1ST強帖: 6
- 應助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時
- 蟲號: 394087
- 注冊: 2007-06-06
- 專業(yè): 半導體材料
7樓2010-11-04 18:05:44
| 支持， |
9樓2010-11-05 16:47:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cenwanglai 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程294求調(diào)劑 +4 | mikasa的圍巾 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿華北電力大學能動專碩，293，求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +8 | 辰沐5211314 2026-03-26 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 0703化學求調(diào)劑，各位老師看看我�。�！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 化學308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿山東大學藥學學碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|