大家好,我是小马老师。

本文介绍从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文件相关推荐

  1. SQL Server 2008如何导出带数据的脚本文件

    SQL之前的版本导出脚本的时候只能把数据表的结构导出,却不能把其中的数据一齐导出来,在平时的应用中难免有时希望获取数据导出的脚本,比如有时SQL版本不一样的时候,无法把高版本的数据库备份直接还原到低版 ...

  2. vue+element-ui文件导出模板及导入xlsx文件

    vue文件的导出与导入 前言 一.导出模板和导入xlsx文件 二.导出JS方法,单独写示例 1.接口方法 2.页面调用导出方法 三.导入表格JS代码 前言 当我们做后台管理的时候,不免会遇到导出模版, ...

  3. 重复文件导出html,分享|删除重复文件的神器:dupeGuru

    最近,我需要清理我父亲的文件和文件夹.有一个难题是,里面存在很多不正确的名字的重复文件.有移动硬盘的备份,同时还为同一个文件编辑了多个版本,甚至改变的目录结构,同一个文件被复制了好几次,名字改变,位置 ...

  4. sqlserver导出带数据的脚本文件

    如何使用sqlserver导出数据库脚本 在程序员的工作中,经常会遇到数据库更换服务器的情况,那么如何有效快速安全的将数据库搬迁呢?下面我介绍一种通过导出sql脚本,然后到新服务器中执行的方法来实现. ...

  5. HDFS文件导出本地合并为一个文件

    HDFS受限于Block大小,大文件会分割成多个块分布在多个节点,导出本地的时候直接用: hadoop fs -get 命令会在本地创建一个目录存放多个块. 要想合并为一个大文件可以这样: hadoo ...

  6. windows下用bat文件启动带参数的bat文件的方法

    背景: 在windows下想启动zookeeper和kafka,无奈每次都要打开两个cmd,先切换到zk和kafka的目录,然后敲一大长串的命令执行,挺费事的,所以想写个脚本去完成一键启动. 使用以下 ...

  7. golang文件操作:打开关闭文件,带缓冲区的读文件

  8. lammps教程:delete_atoms删除原子后原子ID不连续的两种解决方案

    大家好,我是小马老师. 使用lammps自带命令建模时,有时候会用delete_atoms命令删除部分原子,但是在删除原子后,lammps并不一定会重新分配原子的ID,造成data文件中的原子ID序号 ...

  9. lammps教程:lj/cut力场使用方法

    lammps模拟设置中,力场的设置非常重要,力场设置合理才能准确的描述原子之间的受力. 如果两种原子之间没有键连接,可以使用pair_style命令设置这两种原子的力场. pair_style能够设置 ...

最新文章

  1. Open vSwitch 安装
  2. Windows 2003不同域之间迁移密码方法笔记截图
  3. python伪装浏览器https_python伪造HTTP-REFERER
  4. ThhinkPHP5隐藏入口文件index.php
  5. Python桌面图形程序美化的方法论
  6. Rosenbrock函数到底什么用
  7. usr/bin/expect方式免密码登录和发送文件脚本
  8. 关于 win7 远程桌面的多用户问题
  9. ddos工具linux,DDoS常用工具大全
  10. python怎么改变字体大小_python IDLE 背景以及字体大小的修改方法
  11. IPTABLES常用配置
  12. 未能找到文件“\bin\roslyn\csc.exe”
  13. WPF 利用键盘钩子来捕获键盘,做一些不为人知的事情...完整实例
  14. Overloud Gem Modula for Mac - 三合一经典调制插件
  15. 西门子系统C语言编程,SINUMERIK西门子系统加工中心编程各代码.doc
  16. 基于RV1126平台imx291分析 --- v4l2_pipeline_pm_use
  17. 【国家统计局】统计用区划和城乡划分代码2019
  18. Win10环境VMware开WinXP虚拟机CPU占用100%
  19. PCB正片和负片的个人理解
  20. Ubuntu上无法登录网页版本微信的解决方案

热门文章

  1. XnView Standard 1.91.1
  2. python代码---银行卡冻结时间算法
  3. 俄版百度44.7G源码泄露!遭前员工背刺,涉及搜索地图打车电商等所有主要服务...
  4. 酷派大神目标年出货千万台 不降价
  5. android滑动卡片开源,CardSlideView 一行代码实现ViewPager卡片效果
  6. Arduino与掌控板通过软串口通信将数据上传阿里云
  7. TPLINK TL-WR710N 迷你无线路由器桥接模式无法上网问题解决
  8. 守望先锋服务器不稳定,《守望先锋》国服频现高延迟掉线 手把手教你如何分析网络情况...
  9. linux学习笔记(10) - vi和vim文字处理器
  10. Zipf's law