lammps教程:如何从dump文件导出带键角的data文件
大家好,我是小马老师。
本文介绍从dump文件导出data文件的方法。
在lammps模拟中,可能会遇到从前一个模拟dump文件中提取data文件的情况。
平常用的最多的方法是用ovito打开dump文件,点击“file”-“export file”的方式导出data文件。
如果模拟的体系没有键角,用这种方式是可以的,但如果模拟的是聚合物体系,用这种方法导出的data文件只有原子坐标,没有键角等信息。
显然,这样导出的data文件无法作为下一步模拟的data文件。
下面介绍一种方法,可以很方便的从dump文件中导出包含键角信息的data文件。
使用的软件还是python和ovito模块。
python ovito模块不仅可以一次性的读入dump文件,也可以分别读入包含键角信息的拓扑文件和包含模拟结果的轨迹文件。
在本例中,包含键角的拓扑文件可以用初始结构文件pp.data,轨迹文件就是模拟产生的dump.xyz文件。
读入代码为:
pipeline = import_file("pp.data")
traj_mod = LoadTrajectoryModifier()
traj_mod.source.load('dump.xyz')
pipeline.modifiers.append(traj_mod)
读入以上两个文件后,使用export_file()命令输出data文件。
export_file(pipeline, "output.data", "lammps/data", atom_style="full")
打开output.data文件可以看到,data文件包含键角信息:
在ovito中也可以看到键连接情况:
所有代码为:
from ovito.io import import_file, export_file
from ovito.modifiers import LoadTrajectoryModifier
#读入结构和轨迹文件
pipeline = import_file("pp.data")
traj_mod = LoadTrajectoryModifier()
traj_mod.source.load('dump.xyz')
pipeline.modifiers.append(traj_mod)
#输出data
export_file(pipeline, "output.data", "lammps/data", atom_style="full")
公众号:lammps加油站
lammps教程:如何从dump文件导出带键角的data文件相关推荐
- SQL Server 2008如何导出带数据的脚本文件
SQL之前的版本导出脚本的时候只能把数据表的结构导出,却不能把其中的数据一齐导出来,在平时的应用中难免有时希望获取数据导出的脚本,比如有时SQL版本不一样的时候,无法把高版本的数据库备份直接还原到低版 ...
- vue+element-ui文件导出模板及导入xlsx文件
vue文件的导出与导入 前言 一.导出模板和导入xlsx文件 二.导出JS方法,单独写示例 1.接口方法 2.页面调用导出方法 三.导入表格JS代码 前言 当我们做后台管理的时候,不免会遇到导出模版, ...
- 重复文件导出html,分享|删除重复文件的神器:dupeGuru
最近,我需要清理我父亲的文件和文件夹.有一个难题是,里面存在很多不正确的名字的重复文件.有移动硬盘的备份,同时还为同一个文件编辑了多个版本,甚至改变的目录结构,同一个文件被复制了好几次,名字改变,位置 ...
- sqlserver导出带数据的脚本文件
如何使用sqlserver导出数据库脚本 在程序员的工作中,经常会遇到数据库更换服务器的情况,那么如何有效快速安全的将数据库搬迁呢?下面我介绍一种通过导出sql脚本,然后到新服务器中执行的方法来实现. ...
- HDFS文件导出本地合并为一个文件
HDFS受限于Block大小,大文件会分割成多个块分布在多个节点,导出本地的时候直接用: hadoop fs -get 命令会在本地创建一个目录存放多个块. 要想合并为一个大文件可以这样: hadoo ...
- windows下用bat文件启动带参数的bat文件的方法
背景: 在windows下想启动zookeeper和kafka,无奈每次都要打开两个cmd,先切换到zk和kafka的目录,然后敲一大长串的命令执行,挺费事的,所以想写个脚本去完成一键启动. 使用以下 ...
- golang文件操作:打开关闭文件,带缓冲区的读文件
- lammps教程:delete_atoms删除原子后原子ID不连续的两种解决方案
大家好,我是小马老师. 使用lammps自带命令建模时,有时候会用delete_atoms命令删除部分原子,但是在删除原子后,lammps并不一定会重新分配原子的ID,造成data文件中的原子ID序号 ...
- lammps教程:lj/cut力场使用方法
lammps模拟设置中,力场的设置非常重要,力场设置合理才能准确的描述原子之间的受力. 如果两种原子之间没有键连接,可以使用pair_style命令设置这两种原子的力场. pair_style能够设置 ...
最新文章
- android获取图片缩略图,Android系获取图片和视频的缩略图
- pcb结构链表_第2章 2-1进程与PCB
- Python语言学习之字母E开头函数使用集锦:endswith使用用法之详细攻略
- 如何从复杂的货币种类中选出人民币格式的数字
- Python集合list,tuple,dict,set
- elasticsearch_script_01
- 内部类及静态内部类的实例化
- Android 4.1初识--7月13号
- 回文自动机算法+马拉车算法概述及习题【回文串系列问题】
- Maven配置ojdbc14-10.2.0.4.0.jar
- PHP能源管理系统ems,能源管理系统(EMS)方案.doc
- 发现并充分发挥你的长处—盖洛普优势测试
- Jvavscript- 跟着李南江学编程
- python小游戏课程设计报告_贪吃蛇游戏课程设计报告
- 前端html小技巧(form篇):复选框、单选按钮的使用
- 武侠小说阅读记录 2009-02-18 21:14
- 看看这些鲜为人知的宝藏Python数据科学包吧!
- 信阳市公安刑侦信息管理系统
- 运用AI3D绕转工具制作立体LOGO
- C语言菜鸟第一个代码