转载自:https://www.cnblogs.com/sysu/p/10817315.html  和   https://www.cnblogs.com/panscience/p/4953940.html

lammps做分子动力学模拟时,需要一个输入文件(input script),也就是in文件,以及关于体系的原子坐标之类的信息文件(data file)和势文件(potential file)。lammps在执行计算的时候,从这个in文件中读入命令,所以对LAMMPS的使用最主要的就是对in文件的编写和使用。

下面是in文件书写方面的注意事项以及对4个模块的详细解释,原文来源于百度文库:

http://wenku.baidu.com/view/152b56906bec0975f465e26a.html?re=view

由RCWDliwei于2014年4月13日上传。在转载并重新编排之前对原文作者表示致谢。

   1lammps命令分类整理

Initialization

atom_modify, atom_style, boundary,dimension, newton, processors, units

Atom definition

create_atoms, create_box, lattice, read_data, read_restart, region, replicate

Force fields

angle_coeff, angle_style, bond_coeff, bond_style, dielectric, dihedral_coeff,

dihedral_style,improper_coeff, improper_style, kspace_modify, kspace_style,

pair_coeff, pair_modify, pair_style, pair_write,special_bonds

Settings

communicate, dipole, group, mass, min_modify, min_style, neigh_modify,

neighbor, reset_timestep, run_style, set, shape, timestep, velocity

Fixes

fix, fix_modify, unfix

Computes

compute, compute_modify, uncompute

Output

dump, dump_modify, restart, thermo, thermo_modify, thermo_style, undump, write_restart

Actions

delete_atoms, delete_bonds, displace_atoms, displace_box, minimize, run, temper

Miscellaneous

clear, echo, if, include, jump, label, log, next, print, shell, variable

 2、命令书写注意事项

  1. 每一非空行都被认为是一条命令(大小写敏感,但极少有命令或参数大写的)。
  2. 各命令的顺序可能会对计算产生影响,但大部分情况下不会有影响。
  3. 每行后的“&” 表示续行(类似fortran);“#”表示注释(类似bash)。
  4. 每行命令中的不同字段由空格或者制表符分隔开来,每个字段可以由字母、数字、下划线、或标点符号构成。
  5. 每行命令中第一个字段表示命令名,之后的字段都是相关的参数。
  6. 很多命令都是在需要修改默认值的情况下才特别设置的。

 3in文件4个部分的详细解释

       3.1   Initialization

这一部分包含了关于计算体系最基本的信息,例如:

  • units: 单位系统(units style),lammps现在提供包括lj、real、metal、si和cgs几种单位系统;
  • dimension: 定义了两维或者三维模拟(默认是三维);
  • boundary: 定义了分子动力学体系使用的边界条件,例如周期性边界条件或者自由边界条件等;
  • atom_style: 定义模拟体系中的原子属性,这一命令与力场设置的参数中的原子类型(atom type)不同;
  • pair_style: 相互作用力场类型,例如范德化势或者硬球势等;
  • bond_style: 键合相互作用势类型;
  • angle_style: 键角作用势类型;
  • dihedral_style: 二面角作用势类型;
  • improper_style: 混合作用势类型;

    还有其它一些参数设置,例如newton,
    processors,boundary, atom_modify等,可查阅lammps官网。

3.2   Atom
definition

lammps提供3种定义原子方式

  1. 通过read_data或read_restart命令从data或restart文件读入,这些文件可以包含分子拓扑结构信息,这一方法在续算上也很有用。
  2. 按照晶格的方式创建原子,这种方式不包含分子拓扑信息,可能会用到如下的一些命令:lattice, region, create_box,
    create_atoms。
  3. 对已经设置好的原子可以用replicate命令复制后生成一个更大规模的计算体系(注意它与周期性边界条件是两个不同的概念)。
  4. 相应命令:

lattice ** 晶格参数
region **  选择一个区域
create_box **  创建一个盒子
create_atoms ** 创建原子
group, **  给原子分组
delete_atoms,**  删除某些原子

delete_bonds, 
displace_atoms,* 移动某些原子
replicate* 将已经构建的模型在三维重复复制放大
read_data***这个命令就是用来读取其他软件建的模型的
read_restart* 以之前运行的结果作为这次计算的初始构型  

     3.3
  Settings

原子或分子的拓扑信息定义好后,就需要制定一系列的设置,例如力场系数、模拟参数、输出选项等

      力场系数可以通过这样的一些命令来定义:pair_coeff,
bond_coeff, angle_coeff, dihedral_coeff, improper_coeff, kspace_style,
dielectric, special_bonds等。实际上力场系数也可以在关于体系的原子坐标之类的信息的文件(data file)中制定,这样具体参考read_data命令的相关介绍。

模拟参数可以由如下这样一些命令来设置:neighbor,
neigh_modify, group, timestep, reset_timestep, run_style, min_style, min_modify等。模拟过程中通过compute,
compute_modify, variable等一些命令来制定。

      输出选项可以由thermo, dump, restart等一些命令来设置。

----------------------------------------Seting---------------------------------------------

neighbor,
**  截断半径相关的设置
neigh_modify, **
timestep,** 模拟的步长
compute, *** !定义计算,一般我们需要的输出数据就需要这个命令来定义
uncompute *  取消某个compute设置
variable**这个命令可以多次使用,定义一个变量

reset_timestep, ** 步长计数器清零,有多个run命令的需注意这个
run_style, 配套run的,不一定需要
compute_modify *配套compute的

----------------------------------------Minimize-------------------------------------------
结构优化,注意和热力学平衡区分开来
minimize,**
min_style, **
min_modify **

     

     3.4   Run a simulation

通常run命令被设置在in文件的最后,使用run命令来开始一个分子动力学模拟的过程。另外,使用minimize命令来实施能量最小化计算,使用temper命令来进行复制品交换采样模拟。

关于LAMMPS计算前后的处理问题,计算前的原子初始形态文件的生成,由read_data读入一个data文件,这个文件包括体系中各个原子的xyz坐标等等相关参数,或者由其他软件生成并修改后符合LAMMPS的输入文件格式生成。而计算后的输出,因为LAMMPS不支持图形输出,需要借助第三方可视化软件实现,例如VMD。

LAMMPS的输出文件主要可以分为三种:

一种是log.lammps,这里面记录了整个计算过程屏幕上显示的所有信息,可由thermo、thermo_modify等命令控制;

另一种是输出应力、能量、原子位置、速度等等信息,由dump命令控制输出文件;

第三种是断点续算的restart文件输出信息,由write_restart命令控制。

------------------------------------------Output------------------------------------------
dump, ** 输出数据,输出的为单个原子的信息,可视化就用它
undump,* 取消某个dump设置
thermo, ** 输出全局数据,比如温度,压强,长度,总能量等
thermo_style,* 配套thermo的
thermo_modify,*
(后面的都是输出一些数据类型,我一般没怎么用,具体需要的可以单独自己了解)
dump image, 
dump_modify, 
dump movie, 
restart, *输出restart文件,对应前面read-restart,data文件只有原子坐标,restart含有部分计算信息
write_data, *将当前模型输出为一个data文件,对应之前的read-data,
write_dump, 
write_restart*
------------------------------------------Run a
simulation---------------------------------------- 
run,  ***这个命令出现之后,模拟才开始进行计算,前面的都是在设置计算过程

change_box, *
neb prd, 
rerun, 
temper

-------------------------------------------Others--------------------------------------------
这些命令主要可以实现循环以及逻辑判断,通过他们,可以向更复杂的模拟过程前进
clear,*   log, *   echo,   if, *  include,  jump, *   label, *  
next, *  print, *   shell,

in文件注意事项及详细解释(转载)的更多相关文章

  1. in文件注意事项及详细解释

    lammps做分子动力学模拟时,需要一个输入文件(input script),也就是in文件,以及关于体系的原子坐标之类的信息文件(data file)和势文件(potential file).lam ...

  2. 对单片机的modbus RTU的详细解释(转载)

    Modbus 一个工业上常用的通讯协议.一种通讯约定.Modbus协议包括RTU.ASCII.TCP.其中MODBUS-RTU最常用,比较简单,在单片机上很容易实现.虽然RTU比较简单,但是看协议资料 ...

  3. 移动web 的viewport设置注意事项的详细解释 六一快乐=-_-_

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  4. C语言scanf函数详细解释(转载)

    原文地址:https://blog.csdn.net/21aspnet/article/details/174326 scanf 函数名: scanf 功 能: 执行格式化输入 用 法: int sc ...

  5. SpringBoot-02:SpringBoot中的POM文件详细解释

    ------------吾亦无他,唯手熟尔,谦卑若愚,好学若饥------------- 我把pom文件,以及它的详细解释发出来 <?xml version="1.0" en ...

  6. 关于 /proc/sys/net/ipv4/下 文件的详细解释

    关于 /proc/sys/net/ipv4/下 文件的详细解释: 1) /proc/sys/net/ipv4/ip_forward  该文件表示是否打开IP转发.      0,禁止      1,转 ...

  7. .htaccess语法之RewriteCond与RewriteRule指令格式详细解释

    htaccess语法之RewriteCond与RewriteRule指令格式详细解释 (2012-11-09 18:09:08) 转载▼ 标签:  htaccess it 分类: 网络 上文htacc ...

  8. tar命令的详细解释

    tar命令的详细解释 标签: linuxfileoutputbashinputshell 2010-05-04 12:11 235881人阅读 评论(12) 收藏 举报  分类: linux/unix ...

  9. Linux学习笔记15——GDB 命令详细解释【转】

    GDB 命令详细解释 Linux中包含有一个很有用的调试工具--gdb(GNU Debuger),它可以用来调试C和C++程序,功能不亚于Windows下的许多图形界面的调试工具. 和所有常用的调试工 ...

随机推荐

  1. macOS 遇到 svnadmin无法使用的情况

    首先,我在网上查了很多资料,大多数开发者都是表明了mac内置svn,然后直接使用svnadmin来创建仓库,但是我用命令行敲入svnadmin却显示找不到svnadmin.最开始,有个资料说用whic ...

  2. FPAG_Microblaze_PWM_定时器

    由于Xilinx底层库的定时器没有PWM例程,调试过程中费了不少劲.生产PWM需要两个定时器同时工作,一个控制频率,一个控制占空比,位数可通过硬件设置. #include "xtmrctr_ ...

  3. [草稿]基于 Virtuoso 环境比较便捷的项目文件及权限管理方案

    https://www.cnblogs.com/yeungchie/ 假设如下情况: 1 项目名称 Project_01 2 包含 4 名研发用户,user01 和 user02 为前端工程师,use ...

  4. 5.15 牛客挑战赛40 B 小V的序列 关于随机均摊分析 二进制

    LINK:小V的序列 考试的时候 没想到正解 于是自闭. 题意很简单 就是 给出一个序列a 每次询问一个x 问序列中是否存在y 使得x^y的二进制位位1的个数<=3. 容易想到 暴力枚举. 第一 ...

  5. IDEA、maven创建webapp项目

      maven官方入门指南:http://maven.apache.org/guides/getting-started/index.html 推荐跟着官方文档学习噢~   IDEA.maven创建w ...

  6. 搭建Redis主从复制的集群

    在主从复制模式的集群里,主节点一般是一个,从节点一般是两个或多个,写入主节点的数据会被复制到从节点上,这样一旦主节点出现故障,应用系统能切换到从节点去读写数据,这样能提升系统的可用性.而且如果再采用主 ...

  7. Centos8最小化部署安装OpenStack Ussuri

    #!/bin/bash #Centos8最小化部署安装OpenStack Ussuri #共两台主机,分别是一台控制节点,一台计算节点 #.控制节点内存4096M.双网卡,分别为eth0:10.0.0 ...

  8. JAVA学习:HashMap 和 ConcurrentHashMap

     一.最基本的HashMap 和 ConcurrentHashMap 1.HashMap的结构和底层原理:由数组和链表组成,数组里面每个地方都存了Key-Value这样的实例,在Java7叫Entry ...

  9. cxuan读者的外包面试之旅

    最开始是打算面试外包公司刷经验,等经验差不多了,再去甲方尝试,可惜不太顺利. 一.迈思(面试了30分钟) 自我介绍 左连接(a 表左连接 b 表,a 表全部数据出来,b 表没有的数据为空) a 表左连 ...

  10. GIT pull 如何解决 fatal: refusing to merge unrelated histories

    在Github新建一个仓库,写了Readme.md,然后把本地一个已有内容的仓库上传. 先pull,因为两个仓库不同,发现refusing to merge unrelated histories,无 ...