| 查看: 1375 | 回复: 5 | |||
[交流]
在VNL中导入LAMMPS轨迹文件
|
|||
|
简介 使用VNL您可以导入、可视化和分析由其他多种模拟代码（如VASP，Quantum ESPRESSO等等）产生的输出文件。也可以导入由LAMMPS产生的经典分子动力学轨迹文件，并且通过Movie Tool、3D Viewer和MD Analyzer等工具进一步分析您的模拟结果。 LAMMPS是一个经典分子动力学代码，输出的分子动力学轨迹格式可以非常灵活[Pli95]，大部分的输出文件为文本文件。将LAMMPS轨迹导入VNL有两种选择： 导入为Trajectory对象，仅可以用3D Viewer可视化和拖拽入Builder； 导入为MDTrajectory，可以进一步使用Movie Tool或者MDAnalyzer进行分析。 为了完整的使用VNL的MD分析功能，在进行LAMMPS模拟的时候，您需要遵守一些小的规则（在本教程中会说明）。 您想学到如何将LAMMPS模拟轨迹文件导入Virtual Nanolab和如何使用可视化和分析工具，以及如何建立基于输入构象的进一步模拟的输入文件。 准备LAMMPS脚本 LAMMPS模拟是由输入脚本控制的，所有的模拟参数列在输入脚本中。输出轨迹文件是通过 dump 命令来详列的。VNL目前能够阅读写成 atom或者 custom模式的轨迹文件。对于VNL来讲，轨迹文件的扩展名不需要是某一明确的类型，只要和其他没有冲突就好。 （1） 导入为基本的Trajectory格式 只需选定atom格式，额外的指定就不需要了。命令也会变得非常简单，例如： dump d_trajectory all atom 100 trajectory.dat 它会每100步将所有原子的坐标写入trajectory.dat 文件。atom模式不写入速度，这样得到的轨迹文件不能导入为MDTrajectory，仅可以导入为Trajecto ry对象。 （2）要导入为MDTrajectory对象格式，除了原子坐标之外，如果要作为MDTrajectory，轨迹文件还需要写入原子速度，这可以通过这条命令实现： dump d_trajectory all custom 100 trajectory.dat id type x y z vx vy vz MDTrajectory对象还需要提取出势能数据。这些可以从模拟过程中创建的log文件中的热力学输出中读取。为此，LAMMPS中的热力学输出一定要和轨迹文件输出频率同步。您可以通过关键词 thermo 来指定这个输出的频率。而且您必须确保势能在热力学输出文件中列出，这可以通过使用 thermo_style 关键词来选择。典型例子如下（与上面的轨迹dump命令一致）: thermo 100 thermo_style custom step temp pe etotal press 这里pe指出势能将被写入log文件。 如果满足了以上所有要求，VNL可以将轨迹文件作为MDTrajectory来读取。 说明： 在LAMMPS custom格式允许你去明确地设定记录哪些量。与VNL一致的是，至少type, index和坐标（x,y,z,或者 xs,ys,zs）需要在轨迹文件中列出。 运行LAMMPS LAMMPS提供格式明确的输出文件和初始结构，VNL首先直接读取包含了所有计算设置的log文件。基于log文件中的信息，才会进一步读取输出文件和其他的输入文件。所以您需要确保LAMMPS正确输出log文件，这可以在运行LAMMPS的命令行中加入-log选项来指定log文件输出： lmp_serial -in in.settings -log log_simulation_1.lammps 注意： 您一定要使用以.lammps结尾的文件名，因为只有这种扩展名才能被VNL识别为LAMMPS的log文件。如果您是从一个独立的文件（典型的由data.*开始的）读入的输入结构，也需要保证这个文件在您的当前目录中，因为VNL要从这个文件中读取原子量。 导入VNL 如果您按照先前的段落操作，在LAMMPS模拟之后，轨迹对象应该直接出现在VNL的LabFloor中。轨迹对象将会按照log文件命名（不含扩展名）。log_simulation_1.lammps文件导入后的对象将会被命名为log_simulation_1。如果您点击它的图标，将会在右手边的侧栏中看到哪些分析插件工具可以用于您的轨迹对象。对于一个MDTrajectory您可以使用Viewer、Movie Tool和MDAnalyzer。对于常规的Trajectory，您只能使用Viewer。 在接下来的例子里，一个硅晶体在0压力下，在200ps时间内温度由300K线性上升到3000K，从而融化。硅原子之间的相互作用使用Stillinger-Weber势函数描述。最终得到的LAMMPS轨迹文件将会展现在LabFloor中，并使用Movie Tool进行可视化。 [img][/img]  注释： LAMMPS模拟有时被用于无量纲单元中的体系，这可以通过关键词units lj和一般的原子类型来实现。这些轨迹文件也可以导入VNL。但是因为没有与速度和能量相关的物理单位，轨迹文件将会被按照基本的Trajectory形式导入。下面的例子展示了在一个由一般的Lennard-Jones原子构成的二维系统中裂缝的初始和增大。 点击查看动图 处理LAMMPS轨迹文件 如果您将LAMMPS的轨迹文件导入VNL LabFloor，您当然能够单独提取某个结构，并做进一步分析或者使用VNL和ATK计算。下面的例子向您展示了如何计算已经用LAMMPS模拟过的多晶体铜的结构。如果您想学习在VNL-ATK中产生和模拟多晶体结构，可以参考手册 模拟多晶体铜的蠕变实验。 用Movie Tool打开LAMMPS轨迹文件。切换至想要的结构，比如最后的一个结构，并且使用Movie Tool窗口中右下角的 按钮将这个结构发送给 Script Generator 。 在Script Generator中，添加 Analysis 中LocalStructure 块。通过双击打开块，选着合适的文件名，例如Polycrystal_local_structure.nc，不要勾选 Print的框（这个结构非常大，如果选择将会生成很大的文本输出），并且点击OK来关闭 LocalStructure设置。 使用Job Manager来运行计算。 在计算完成之后，当前的结构文件将会出现在LabFloor中。 局域结构分析（LocalStructure） 在Viewer中打开它并在Viewer的面板栏中找到Local Structure插件。点开它您将获得当前晶体结构的表单。在这个例子中，我们可以高亮显示多晶体的晶界结构: 选择除了FCC之外所有的结构类型，FCC是在硅晶粒中原始的晶体结构。 点击 Properties 插件用于打开属性界面。在 Atoms 标签里变化选择原子的颜色，例如白色，用于区别这些原子与晶粒原子。 另外您可以点击 Ball，用于关闭晶界内部的键显示。这样可以增加多晶体的透明效果。 关闭Properties窗口将激活您作的改变，并且晶界也变得清晰可见。 进一步进行其他计算 除了当前结构分析，您还能够进行其他模拟，例如使用ATK-DFT引擎对提取的某个结构进行密度泛函计算：这只需要在Script Generator中加入 New Calculator 块 ，这为您的经典的LAMMPS模拟提供量子力学计算接口。 参考文献 [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网站 链接：LAMMPS官方手册 |
回复此楼» 猜你喜欢
0805 316求调剂
已经有6人回复
-
085601求调剂总分293英一数二
已经有3人回复
-
考研调剂
已经有3人回复
-
08工学调剂
已经有17人回复
-
340求调剂
已经有4人回复
-
311求调剂
已经有3人回复
-
食品专硕 一志愿双一流 328
已经有4人回复
-
材料调剂
已经有6人回复
-
300分，材料，求调剂，英一数二
已经有5人回复
-
求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
南京医科大学生殖医学与子代健康国重实验室-董飞宏课题组-2026年博士招生-9月入学
+1/85
-
济南大学林秀娟教授招收材料类博士生
+1/82
-
西南交通大学医学院招收生物医学工程和材料与化工专业硕士研究生
+1/80
-
祈福2026国基面上项目
+5/70
-
招收2026年材料与化工专硕调剂，山东省属重点大学，过国家线即可调剂，名额充足
+2/32
-
从专科一路到211硕士，到TOP5博士，也想进高校体制
+1/29
-
年薪约32-60万！成会明院士团队诚聘材料理论计算博士后（长期有效）
+1/21
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/17
-
丽水学院 材料与能源工程学院 陈江飞老师课题组硕士招生----材料与化工专业
+1/15
-
三峡大学 计算力学课题组招收调剂（力学/水利/土木/材料/矿业）
+1/9
-
中国民航飞行学院先进连接团队招收航空宇航、航空机械专业硕士研究生（学硕专硕均可）
+1/7
-
中国人民大学物理学系蔡鹏课题组招收优秀硕士博士研究生
+1/7
-
中国中医科学院中医药健康产业研究所（江西中医药健康产业研究院）2026年高层次人才等
+1/6
-
中国人民大学物理学系蔡鹏课题组招收优秀硕士博士研究生
+1/6
-
调剂 22408 340分
+1/4
-
天津工业大学双一流高校招博士、硕士、博士后
+1/3
-
上海师范大学有机化学专业研究生招生
+1/3
-
【科研干货】纳米抗体凭什么成为实验室“新宠”？附优质制备服务商参考
+1/1
-
上海交通大学AIMS-Lab招收AI for Science【智能制造】方向2026级博士生
+1/1
-
环境工程专业 2026年硕士研究生调剂名额（有学硕1名，专硕2名）
+1/1
简单回复
psylhh2楼
2016-06-10 22:30
回复
sxj9420(金币+1): 谢谢参与
2016-06-10 22:40
回复
sxj9420(金币+1): 谢谢参与
tzynew4楼
2016-06-10 22:44
回复
sxj9420(金币+1): 谢谢参与
tableman5楼
2016-06-10 22:53
回复
sxj9420(金币+1): 谢谢参与
nono20096楼
2016-06-12 08:39
回复
sxj9420(金币+1): 谢谢参与
·
