 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2232 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[交流]
在VNL中導(dǎo)入LAMMPS軌跡文件
|
||||
|
[img][/img]簡(jiǎn)介 使用VNL您可以導(dǎo)入、可視化和分析由其他多種模擬代碼（如VASP，Quantum ESPRESSO等等）產(chǎn)生的輸出文件。也可以導(dǎo)入由LAMMPS產(chǎn)生的經(jīng)典分子動(dòng)力學(xué)軌跡文件，并且通過(guò)Movie Tool、3D Viewer和MD Analyzer等工具進(jìn)一步分析您的模擬結(jié)果。 LAMMPS是一個(gè)經(jīng)典分子動(dòng)力學(xué)代碼，輸出的分子動(dòng)力學(xué)軌跡格式可以非常靈活[Pli95]，大部分的輸出文件為文本文件。將LAMMPS軌跡導(dǎo)入VNL有兩種選擇： 導(dǎo)入為T(mén)rajectory對(duì)象，僅可以用3D Viewer可視化和拖拽入Builder； 導(dǎo)入為MDTrajectory，可以進(jìn)一步使用Movie Tool或者M(jìn)DAnalyzer進(jìn)行分析。 為了完整的使用VNL的MD分析功能，在進(jìn)行LAMMPS模擬的時(shí)候，您需要遵守一些小的規(guī)則（在本教程中會(huì)說(shuō)明）。 您想學(xué)到如何將LAMMPS模擬軌跡文件導(dǎo)入Virtual Nanolab和如何使用可視化和分析工具，以及如何建立基于輸入構(gòu)象的進(jìn)一步模擬的輸入文件。 準(zhǔn)備LAMMPS腳本 LAMMPS模擬是由輸入腳本控制的，所有的模擬參數(shù)列在輸入腳本中。輸出軌跡文件是通過(guò) dump 命令來(lái)詳列的。VNL目前能夠閱讀寫(xiě)成 atom或者 custom模式的軌跡文件。對(duì)于VNL來(lái)講，軌跡文件的擴(kuò)展名不需要是某一明確的類(lèi)型，只要和其他沒(méi)有沖突就好。 （1） 導(dǎo)入為基本的Trajectory格式 只需選定atom格式，額外的指定就不需要了。命令也會(huì)變得非常簡(jiǎn)單，例如： dump d_trajectory all atom 100 trajectory.dat 它會(huì)每100步將所有原子的坐標(biāo)寫(xiě)入trajectory.dat 文件。atom模式不寫(xiě)入速度，這樣得到的軌跡文件不能導(dǎo)入為MDTrajectory，僅可以導(dǎo)入為T(mén)rajecto ry對(duì)象。 （2）要導(dǎo)入為MDTrajectory對(duì)象格式，除了原子坐標(biāo)之外，如果要作為MDTrajectory，軌跡文件還需要寫(xiě)入原子速度，這可以通過(guò)這條命令實(shí)現(xiàn)： dump d_trajectory all custom 100 trajectory.dat id type x y z vx vy vz MDTrajectory對(duì)象還需要提取出勢(shì)能數(shù)據(jù)。這些可以從模擬過(guò)程中創(chuàng)建的log文件中的熱力學(xué)輸出中讀取。為此，LAMMPS中的熱力學(xué)輸出一定要和軌跡文件輸出頻率同步。您可以通過(guò)關(guān)鍵詞 thermo 來(lái)指定這個(gè)輸出的頻率。而且您必須確保勢(shì)能在熱力學(xué)輸出文件中列出，這可以通過(guò)使用 thermo_style 關(guān)鍵詞來(lái)選擇。典型例子如下（與上面的軌跡dump命令一致）: thermo 100 thermo_style custom step temp pe etotal press 這里pe指出勢(shì)能將被寫(xiě)入log文件。 如果滿足了以上所有要求，VNL可以將軌跡文件作為MDTrajectory來(lái)讀取。 說(shuō)明： 在LAMMPS custom格式允許你去明確地設(shè)定記錄哪些量。與VNL一致的是，至少type, index和坐標(biāo)（x,y,z,或者 xs,ys,zs）需要在軌跡文件中列出。 運(yùn)行LAMMPS LAMMPS提供格式明確的輸出文件和初始結(jié)構(gòu)，VNL首先直接讀取包含了所有計(jì)算設(shè)置的log文件�；趌og文件中的信息，才會(huì)進(jìn)一步讀取輸出文件和其他的輸入文件。所以您需要確保LAMMPS正確輸出log文件，這可以在運(yùn)行LAMMPS的命令行中加入-log選項(xiàng)來(lái)指定log文件輸出： lmp_serial -in in.settings -log log_simulation_1.lammps 注意： 您一定要使用以.lammps結(jié)尾的文件名，因?yàn)橹挥羞@種擴(kuò)展名才能被VNL識(shí)別為L(zhǎng)AMMPS的log文件。如果您是從一個(gè)獨(dú)立的文件（典型的由data.*開(kāi)始的）讀入的輸入結(jié)構(gòu)，也需要保證這個(gè)文件在您的當(dāng)前目錄中，因?yàn)閂NL要從這個(gè)文件中讀取原子量。 導(dǎo)入VNL 如果您按照先前的段落操作，在LAMMPS模擬之后，軌跡對(duì)象應(yīng)該直接出現(xiàn)在VNL的LabFloor中。軌跡對(duì)象將會(huì)按照l(shuí)og文件命名（不含擴(kuò)展名）。log_simulation_1.lammps文件導(dǎo)入后的對(duì)象將會(huì)被命名為log_simulation_1。如果您點(diǎn)擊它的圖標(biāo)，將會(huì)在右手邊的側(cè)欄中看到哪些分析插件工具可以用于您的軌跡對(duì)象。對(duì)于一個(gè)MDTrajectory您可以使用Viewer、Movie Tool和MDAnalyzer。對(duì)于常規(guī)的Trajectory，您只能使用Viewer。 在接下來(lái)的例子里，一個(gè)硅晶體在0壓力下，在200ps時(shí)間內(nèi)溫度由300K線性上升到3000K，從而融化。硅原子之間的相互作用使用Stillinger-Weber勢(shì)函數(shù)描述。最終得到的LAMMPS軌跡文件將會(huì)展現(xiàn)在LabFloor中，并使用Movie Tool進(jìn)行可視化。  注釋?zhuān)?br /> LAMMPS模擬有時(shí)被用于無(wú)量綱單元中的體系，這可以通過(guò)關(guān)鍵詞units lj和一般的原子類(lèi)型來(lái)實(shí)現(xiàn)。這些軌跡文件也可以導(dǎo)入VNL。但是因?yàn)闆](méi)有與速度和能量相關(guān)的物理單位，軌跡文件將會(huì)被按照基本的Trajectory形式導(dǎo)入。下面的例子展示了在一個(gè)由一般的Lennard-Jones原子構(gòu)成的二維系統(tǒng)中裂縫的初始和增大。 點(diǎn)擊查看動(dòng)圖  處理LAMMPS軌跡文件 如果您將LAMMPS的軌跡文件導(dǎo)入VNL LabFloor，您當(dāng)然能夠單獨(dú)提取某個(gè)結(jié)構(gòu)，并做進(jìn)一步分析或者使用VNL和ATK計(jì)算。下面的例子向您展示了如何計(jì)算已經(jīng)用LAMMPS模擬過(guò)的多晶體銅的結(jié)構(gòu)。如果您想學(xué)習(xí)在VNL-ATK中產(chǎn)生和模擬多晶體結(jié)構(gòu)，可以參考手冊(cè) 模擬多晶體銅的蠕變實(shí)驗(yàn)。 用Movie Tool打開(kāi)LAMMPS軌跡文件。切換至想要的結(jié)構(gòu)，比如最后的一個(gè)結(jié)構(gòu)，并且使用Movie Tool窗口中右下角的 按鈕將這個(gè)結(jié)構(gòu)發(fā)送給 Script Generator 。 在Script Generator中，添加 Analysis 中LocalStructure 塊。通過(guò)雙擊打開(kāi)塊，選著合適的文件名，例如Polycrystal_local_structure.nc，不要勾選 Print的框（這個(gè)結(jié)構(gòu)非常大，如果選擇將會(huì)生成很大的文本輸出），并且點(diǎn)擊OK來(lái)關(guān)閉 LocalStructure設(shè)置。 使用Job Manager來(lái)運(yùn)行計(jì)算。  在計(jì)算完成之后，當(dāng)前的結(jié)構(gòu)文件將會(huì)出現(xiàn)在LabFloor中。 局域結(jié)構(gòu)分析（LocalStructure） 在Viewer中打開(kāi)它并在Viewer的面板欄中找到Local Structure插件。點(diǎn)開(kāi)它您將獲得當(dāng)前晶體結(jié)構(gòu)的表單。在這個(gè)例子中，我們可以高亮顯示多晶體的晶界結(jié)構(gòu): 選擇除了FCC之外所有的結(jié)構(gòu)類(lèi)型，F(xiàn)CC是在硅晶粒中原始的晶體結(jié)構(gòu)。 點(diǎn)擊 Properties 插件用于打開(kāi)屬性界面。在 Atoms 標(biāo)簽里變化選擇原子的顏色，例如白色，用于區(qū)別這些原子與晶粒原子。 另外您可以點(diǎn)擊 Ball，用于關(guān)閉晶界內(nèi)部的鍵顯示。這樣可以增加多晶體的透明效果。 關(guān)閉Properties窗口將激活您作的改變，并且晶界也變得清晰可見(jiàn)。  進(jìn)一步進(jìn)行其他計(jì)算 除了當(dāng)前結(jié)構(gòu)分析，您還能夠進(jìn)行其他模擬，例如使用ATK-DFT引擎對(duì)提取的某個(gè)結(jié)構(gòu)進(jìn)行密度泛函計(jì)算：這只需要在Script Generator中加入 New Calculator 塊 ，這為您的經(jīng)典的LAMMPS模擬提供量子力學(xué)計(jì)算接口。 參考文獻(xiàn) [Pli95] Steve Plimpton. Fast parallel algorithms for short-range molecular dynamics. Journal of Computational Physics, 117(1):1 – 19, 1995. doi:10.1006/jcph.1995.1039. 鏈接：LAMMPS網(wǎng)站LAMMPS網(wǎng)站[swf][/swf] 鏈接：LAMMPS官方手冊(cè)LAMMPS官方手冊(cè) 閱讀原文：www.fermitech.com.cn 費(fèi)米wiki |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 分子動(dòng)力學(xué)模擬 |
» 猜你喜歡
數(shù)一英一271專(zhuān)碩（085401）求調(diào)劑，可跨
已經(jīng)有8人回復(fù)
-
295求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有10人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑
已經(jīng)有6人回復(fù)
-
289求調(diào)劑
已經(jīng)有12人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有10人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專(zhuān)碩
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 搶金幣啦！回帖就可以得到:
-
東北石油大學(xué)PCAFM實(shí)驗(yàn)室研究生調(diào)劑招生（化學(xué)化工、材料等相關(guān)專(zhuān)業(yè)）
+1/92
-
西南交通大學(xué)醫(yī)學(xué)院招收生物醫(yī)學(xué)工程和材料與化工專(zhuān)業(yè)碩士研究生
+1/85
-
Postdoc position in Hangzhou Normal University （陳健翔組博后招聘）
+1/83
-
中飛院航空鋰電池團(tuán)隊(duì)接受碩士生調(diào)劑
+1/74
-
國(guó)自然面上祈福
+1/71
-
青島大學(xué)固態(tài)電池研究團(tuán)隊(duì)招收材料、化學(xué)、物理等方向博士研究生
+2/46
-
環(huán)境化學(xué)、地學(xué)、環(huán)境毒理方向，擅長(zhǎng)污染物遷移轉(zhuǎn)化降解、分析檢測(cè)、計(jì)算模擬者優(yōu)先
+1/33
-
汕頭大學(xué)海洋與生命學(xué)科李升康教授團(tuán)隊(duì)招收生物類(lèi)調(diào)劑生
+1/21
-
哈爾濱工業(yè)大學(xué)-大灣區(qū)大學(xué)2026年聯(lián)合培養(yǎng)博士研究生春季招生-報(bào)名3月31號(hào)截至
+1/14
-
韓國(guó)慶熙大學(xué) Yu Taekyung 教授招2026年博士研究生
+1/10
-
三峽大學(xué)材化學(xué)院能源環(huán)保晶態(tài)材料課題組接收調(diào)劑，招收化學(xué)類(lèi)碩士
+1/10
-
南京林業(yè)大學(xué)-國(guó)家級(jí)青年人才團(tuán)隊(duì) 招2026級(jí)博士、碩士（合成化學(xué)、植物化學(xué)方向）
+1/9
-
中國(guó)地質(zhì)大學(xué)（北京)工程技術(shù)學(xué)院安老師組招收2026調(diào)劑生
+1/9
-
歡迎來(lái)英國(guó)拉夫堡大學(xué)交流-結(jié)晶，蛋白質(zhì)，AI，水處理等
+1/7
-
中國(guó)科學(xué)院大學(xué)納米科學(xué)與技術(shù)學(xué)院-羅聃課題組誠(chéng)招碩士研究生
+1/5
-
杭州電子科技大學(xué) 磁電功能材料創(chuàng)新團(tuán)隊(duì) 招收碩士調(diào)劑生
+1/4
-
蘇州高校招收08方向調(diào)劑
+1/3
-
【博士招生】天津理工大學(xué)國(guó)家杰青王鐵課題組招收2026年博士研究生
+1/2
-
香港大學(xué) 化學(xué)系 博后招聘
+1/1
-
廣東唯一石化院�！Y源與環(huán)境專(zhuān)碩招生
+1/1
shuxue2012
金蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1356
- 帖子: 169
- 在線: 100.7小時(shí)
- 蟲(chóng)號(hào): 1928755
9樓2016-07-24 20:09:02
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 085600，專(zhuān)業(yè)課化工原理，321分求調(diào)劑 +3 | 大饞小子 2026-03-28 | 3/150 |
|
|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書(shū)。 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|