Lammps 初学者——in 文件中文解析
1.初始化
Units:
Units 命令是用来设定模拟的原子类型。
Units style(lammps 现在提供的有 style=LJ、real、metal、si、cgs、electron)
LJ 是硬球模型
Real 是真实的原子模型
Metal 是金属原子模型Si 是硅(半导体)原子模型Cgs
Electron 是电子模型
Dimension:
Dimension 命令是用来定义模拟的维度,默认情况为三维。
Dimension N (N=2,3)
Boundary:
Boundary 命令是用来设定模拟的边界条件。
Boundary x y z (x,y,z=p,s,f,m 四种类型中的一种或者两种)
P 是周期性边界条件
S 是自由边界条件(但是具有收缩性)
F 是固定边界条件
M 是具有最小值的自由边界条件(但是具有收缩性)
Atom-style:
Atom-style 定义了模拟体系中的原子属性。
Delete-atoms: Pair-style:
Pair-style 定义了相互作用力场类型,即势函数。
- 原子定义
1、
Read-data/read-restart:
Read-data 或 read-restart 为从 data 或 restart 文件中读取内容来定义原子。
2、
Lattice:
lattice 是用来原子类型构建模型晶格结构。
Lattice style scale keyword values...
Style 表示点阵类型,点阵类型有 none、sc、bcc、fcc、hcp、diamond、sq、
sq2、hex、custom。
Scale 表示构建的点阵类型的单位长度
Keyword 表示后面可以追加 0 个或者多个关键字,关键字有 origin、orient、
spacing、a1、a2、a3、basis。各关键字意思为:
Region:
Create-box:
Create-box 命令用于在 region box 命令指定的区域内创建一个模拟的盒子。
Create_box N region-ID
N 为使用该原子类型进行模拟
Region-ID 为使用区域的 ID 使用模拟域
Create-atoms:
Create-atoms 用于在所创建的模拟的盒子中填充某种原子。
Create_atoms type style args keyword values ...
- 设定: 力场系数
Pair-coeff: 模拟参数
Neighbor
Neighbor 是定义
Neighbor skin style
Skin 为
Style 可取 bin、nsq、multi
Neighbor-modify: Group: Timestep:
Timestep 是定义模拟的时间步长。
Timestep dt
Dt 为时间步长(时间单位),默认为 1ps
Reset-timestep:
Fix
Fix 为定义对部分原子或者对系综的操作。
http://lammps.sandia.gov/doc/fix.html Fix ID group-ID style args
ID=指定的操作名称编号(如 1、2、3…….) Group-ID=制定操作范围的原子(如 all……)
Style=操作内容,如何操作
Args=相关操作内容的相关操作参数
Fix-nvt:恒定 N(原子数量)、V(体积)、T(温度)对时间积分Fix-npt:恒定 N、P、T 对时间积分
Unfix:
Compute/compute-modify: 输出选项
Thermo:
Thermo 命令用于定义每隔一定时间步长输出一次采样内容。
Thermo N
每隔 N 个时间步长输出一次采样内容至屏幕或者 log 文件
Thermo_style:
Thermo_style 命令主要是用来说明采样的内容,即输出到屏幕或者 LOG
文件中的内容。(http://lammps.sandia.gov/doc/thermo_style.html)
Thermo_style style args
Style=one、multi、custom
Args=特定的参数列表One args=none Multi args=none
Custom args=list of attributes
Possible attributes=step、elaosed、elaplong、dt、cpu、tpcpu、spcpu、
Atoms、temp、press、pe、ke、otatal、enthalpy、
Evdwl、ecoul、epair、ebond、eangle、edihed、eimp、
Emol、elong、etail、
Vol、lx、ly、lz、xlo、xhi、ylo、yhi、zlo、zhi、
Xy、xz、yz、xlat、ylat、zlat、
Pxx、pyy、pzz、pxy、pyz、pxz、
Fmax、fnorm、
Cella、cellb、cellc、cellalpha、cellbeta、cellgamma、
c_ID、c_ID[I]、c_ID[I][J]、
f_ID、f_ID[I]、f_ID[I][J]、V_name
Dump:
Dump 命令主要是用来输出计算后单个原子的相关信息的。
Dump ID group-ID style N file args
ID:自己给这个 dump 命令定义的一个代号(自定义)
group-ID:那些信息需要被输出的原子群(预先定义好 group) Style:类型
N:每经过多少时间步输出信息
File: 输 出 名 称Args:相关参数(每个类型不同)
Restart:
- 运行
Run: Minimize:
例一
units metal # 单位为 lammps 中的 metel 类型
boundary p p p # 周期性边界条件
atom_style atomic # 原子模式
lattice fcc 3.61 # Cu 的晶格常数 3.61
region box block 0 4 0 4 0 4 # x,y,z 各方向上的晶胞重复单元数,也即区域大小
create_box 1 box # 将上述区域指定为模拟的盒子create_atoms 1 box # 将原子按晶格填满盒子pair_style eam # 选取 Cu 的 EAM 势作为模型pair_coeff * * Cu_u3.eam # EAM 势文件名称
run 0 # 运行 0 步,仅为启动 lammps 的热力学数据计算
variable E equal pe # 定义变量 E 为系统总势能
variable N equal atoms # 定义变量 N 为系统总原子数
print "the number of atoms & system energy now are $N $E" # 打印信息
create_atoms 1 single 2.45 2.05 2.05 # 在该位置插入一个原子
min_style sd # 能量最小化模式,sd
minimize 1.0e-12 1.0e-12 1000 1000 # 能量最小化参数,指数越大最小化程度越深
print "interstitial introduced, minimized: $N atoms, energy is $E"
fix 1 all nvt 100 100 100 drag 0.2 # nvt 系综,原子数、体积和温度保持不变;T=100K
timestep 0.005 # 步长 0.005fs run 1000 # 运行 1000 步
print "nvt performed, temperature up: $N atoms, total energy is $E"
fix 1 all nvt 100 0.0001 100 drag 0.2 # nvt 系综,温度由 100K 到 0.0001K
run 1000 # 运行 1000 步
print "nvt performed, temperature down: $N atoms, total energy is $E" compute 3 all pe/atom # 计算每个原子的势能
compute 4 all ke/atom # 计算每个原子的动能
compute 5 all coord/atom 3.0 # 计算每个原子的近邻原子数
dump 1 all custom 1 dump.atom id xs ys zs c_3 c_4 c_5 # 将信息写入 dump.atom min_style sd
minimize 1.0e-12 1.0e-12 10000 10000 # 再次能量最小化
print "the final state: $N atoms, total energy is $E" # 打印信息
~
例二
boundary p s s # 边界条件,拉伸方向
是周期性,其余是自由边界;如果是薄膜拉伸则是两个周期性,块体则是三个周期性
units metal #单位制定义为
metal
atom_style atomic #原子类型
neighbor 2.0 bin #截断半径相关的东西
neigh_modifydelay 1 check yes # 邻近原子列表更新速度
#create geometry
lattice fcc 3.61 # 定义晶胞为 fcc,晶格常数 3.61A
region box block 0 30 0 3 0 3 #定义一个长方体区域叫 box,长
30,宽和高是 3
create_box 1 box #创建了这样一个 box
create_atoms 1 box #在 box 里创建了一种原子
mass 1 63.546 #定义这种原子的质量
是 63.546
# potentials
pair_style eam # 定义势函数是 EAM
pair_coef * * Cu_u3.eam #势所需要的参数在此文件里
#define groups
region 1 block INF 1 INF INF INF INF # 定义了一个叫
1 的区域
group left region 1
#定义此区域里的原子叫 left
region 2 block 29 INF INF INF INF INF #定义了一个叫 2 的区域
group right region 2 #
定义此区域里的原子叫 right
group boundary union left right #定义 left+ right = boundary
group mobile subtract all left # 定义
mobile= all - left
# initialvelocities
velocity left set 0.0 0.0 0.0 # 设置
原子初速度为 0
computep all pressure thermo_temp # 计算应力,计算结果记为 p
variable pressx equal c_p[1] #定义变量 pressx=c_p[1],c_p[1]
的意思是 p 里第一个值
variable pressy equal c_p[2] #定义变量 pressy=c_p[2],c_p[2]
的意思是 p 里第二个值
variable pressz equal c_p[3] #定义变量 pressz=c_p[3],c_p[3]
的意思是 p 里第三个值
thermo_style custom step temp etotal press v_pressx v_pressy v_pressz vol
dump 1 all atom 1000 dump.tensile # 输出结果到
dump.tensile
minimize 1.0e-6 1.0e-6 1000 1000 # 驰豫
fix 1 left setforce 0.0 NULL NULL #固定左边的原子
fix 2 all npt 1.0 1.0 1.0 aniso 0.0 0.0 NULL NULL NULL NULL 10.0 drag 1.0
thermo 1000 #
每 1000 步输出一次结果
timestep 0.002 #
时间步长 0.002ps
run 10000
#运行 10000 步
unfix2
fix 3 all nve
fix 4 all temp/rescale 100 1.0 1.0 5 1.0
fix 5 all deform 1000 x erate 0.002 units box #设置沿 x 方向拉伸,应变率为 0.002(1/ps)
compute s all stress/atom #计算每个原子上的应力,计算结果记为 s
compute strea all reduce sum c_s[1] #将 s 里第一个值求和,计算结果记为 strea
variable stressa equal c_strea/vol # 定义变量 stressa= strea/vol, vol 是体积
compute streb all reduce sum c_s[2] variable stressb equalc_streb/vol compute strec all reduce sum c_s[3] variable stressc equal c_strec/vol
thermo_style custom step tempetotal press v_stressa v_stressb v_stressc vol run 300000
最后,有相关需求欢迎通过微信公众号联系我们。
Lammps 初学者——in 文件中文解析相关推荐
- java 解析csv 乱码_Java采用opencsv解析csv文件以及解析中文乱码问题
参考资料: opencsv文档 写下本文时最新版的opencsv版本为4.0, maven坐标如下: com.opencsv opencsv 4.0 测试用test.csv id,date,name ...
- 解决ini-parser解析ini文件中文乱码问题
解决ini-parser解析ini文件中文乱码问题 参考文章: (1)解决ini-parser解析ini文件中文乱码问题 (2)https://www.cnblogs.com/nodegis/p/95 ...
- cisco 基础配置命令中文解析 1
命令 描述 <?xml:namespace prefix = o ns = "urn:schemas-microsoft- ...
- XML解析 (JAVA解析xml文件)java+Dom4j+Xpath xml文件解析根据子节点得到父节点 查找校验xml文件中相同的节点属性值 java遍历文件夹解析XML
XML解析 (JAVA解析xml文件)java+Dom4j+Xpath xml文件解析根据子节点得到父节点 以及查找xml文件中相同的节点属性值 项目背景:这是本人实习中所碰到的项目,当时感觉很棘手, ...
- 最全的基于c++的serialport.cpp与serialport.h类文件(解析版)
最全的基于c++的serialport.cpp与serialport.h类文件(解析版) 备注 参考原文件为:serialport.cpp与serialport.h类文件源代码 作者:c344935 ...
- studiolibrary安装_初学者daz studio中文基础安装布局教程
本视频教程是关于初学者daz studio中文基础安装布局教程 ,时长:10小时,大小:800 MB,MP4高清视格式,教程使用软件:DAZ Studio,语言:英语. 使用方法1:安装软件,直接将C ...
- Web项目,网页上传excel文件并解析实战示例
最近写了一个基于poi解析excel文件的工具类,所以想在web项目中测试一下,就做了这个简单的项目.本项目主要使用了 SpringMVC+RESTful+Maven的风格.适合有一定基础的人员. 源 ...
- python读取xml文件信息_python读取xml文件方法解析
关于python读取xml文章很多,但大多文章都是贴一个xml文件,然后再贴个处理文件的代码.这样并不利于初学者的学习,希望这篇文章可以更通俗易懂的教如何使用python来读取xml文件. 什么是xm ...
- activiti7(二):画业务流程图的工具以及bpmn文件内容解析
上一篇我们学习了BPMN2.0规范 链接 ,知道了BPMN2.0规范是为了帮助我们规范的画流程图的,那么到底应该怎么画呢?这里给大家介绍一些遵循BPMN2.0规范的画图工具. 说明:业务流 ...
- Office文件的奥秘——.NET平台下不借助Office实现Word、Powerpoint等文件的解析(完)...
原文 http://www.cnblogs.com/mayswind/archive/2013/04/01/2991271.html [题外话] 这是这个系列的最后一篇文章了,为了不让自己觉得少点什么 ...
最新文章
- 自己编写linux系统,自己动手 编写Linux系统的设备驱动程序
- debugInit.c tomcat启动错误
- python爬虫文件代码大全-Python网络爬虫实战项目代码大全(长期更新,欢迎补充)...
- proxy负载均衡、读写分离
- MySQL 报错 1055
- qdu_ACM集训队3月5号组队训练
- linux 打包排除多个目录,tar打包整个目录(可排除子目录)几种方法
- cobol和java区别,COBOL语法和文法(1)
- rgmanager 介绍
- Android自动调整TextView的大小
- gettype拿不到值_这五种古董,别说是买豪车豪宅,放在现实中最多就值一顿饭钱...
- 高级操作系统——XV6进程管理
- 【Solr原理】Leader Shard选举
- 云表平台VS永中office,谁更胜一筹?
- 用python来开发webgame服务端系列
- 后台管理系统——02首页布局
- 无线射频专题《协议类,IEEE 802.11/802.11b/802.11a/802.11g/802.11n/802.11ac标准简介》
- python爬虫基础(12:app数据爬取)
- 【易语言助手插件】IDE 鼠标手势1.2 原汁原味移植浏览器手势
- js jsp 文件是什么?