lammps问题-关于NVT与NVE系综产生的不同
相同的材料构型,在不同的系综下产生如下不同的结果:
1.NVE系综下:
2.NVT系综如下:
产生的不同结果如上,原脚本如下:
1.
shell mkdir initial-comp
shell mkdir config-comp
shell mkdir restart-comp
#--------------------------DEFINITION-------------------
units metal
dimension 3
boundary s s s
atom_style atomic
read_data Ti.data
region 1 block INF INF INF INF INF 27 units box
region 2 block INF INF INF INF 319 INF units box
group bottom region 1
group upper region 2
group mobile subtract all bottom upper
group boundary union bottom upper
#---------- Define Interatomic Potential ---------------------
pair_style eam/fs
pair_coeff * * Ti.eam.fs Ti
#-----------------------------setting------------------------
compute eng all pe/atom
compute 1 all stress/atom NULL
compute v all voronoi/atom
variable e1 equal “step”
variable e2 equal “pe”
variable e3 equal “temp”
#------------------------------------------------CUT------------------------------------------------------
delete_atoms overlap 2 all all compress yes
#----------------------------EQUILIBRATION---------------------
timestep 0.002
velocity mobile create 300.0 12345
min_style cg
minimize 1e-18 1e-18 100000 100000
fix 2 boundary setforce 0.0 0.0 0.0
fix 1 all nve
fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
dump 1 all custom 10000 ./initial-comp/relax.*.dump id type x y z
fix def2 all print 100 “${e1} ${e2} ${e3}” file initial-pe-comp screen no
thermo 2000
thermo_style custom step lx ly lz press pe temp
run 50000
undump 1
unfix def2
variable tmp equal “lz”
variable L0 equal ${tmp}
print “Initial Length, L0: ${L0}”
#-------------------------COMPRESSION---------------------------------
reset_timestep 0
timestep 0.002
dump 2 all custom 500 ./config-comp/comp..dump id type x y z c_1[] c_v[]
restart 500000 ./restart-comp/compression..equil
variable strain equal “-(lz - v_L0)/v_L0”
variable e1 equal “v_strain”
variable s1 equal “pzz/10000”
variable s2 equal “pxx/10000”
variable s3 equal “pyy/10000”
#variable s2 equal “4v_s1lx/(ly*PI)”
fix def1 all print 1000 “${e1} ${s1} ${s2} ${s3}” file Ti1000-comp screen no
thermo 100000
thermo_style custom step temp ke pe press
velocity upper set 0 0 -0.15
velocity mobile ramp vz 0.0 -0.15 z 29 317 sum yes
run 200000
print “ALL DONE!”
shell mkdir initial-comp
shell mkdir config-comp
shell mkdir restart-comp
#--------------------------DEFINITION-------------------
units metal
dimension 3
boundary s s s
atom_style atomic
read_data Ti.data
region 1 block INF INF INF INF INF 27 units box
region 2 block INF INF INF INF 319 INF units box
group bottom region 1
group upper region 2
group mobile subtract all bottom upper
group boundary union bottom upper
#---------- Define Interatomic Potential ---------------------
pair_style eam/fs
pair_coeff * * Ti.eam.fs Ti
#-----------------------------setting------------------------
compute eng all pe/atom
compute 1 all stress/atom NULL
compute v all voronoi/atom
variable e1 equal “step”
variable e2 equal “pe”
variable e3 equal “temp”
#------------------------------------------------CUT------------------------------------------------------
delete_atoms overlap 2 all all compress yes
#----------------------------EQUILIBRATION---------------------
timestep 0.002
velocity mobile create 300.0 12345
min_style cg
minimize 1e-18 1e-18 100000 100000
fix 2 boundary setforce 0.0 0.0 0.0
fix 1 all nvt temp 300 300 100
dump 1 all custom 10000 ./initial-comp/relax.*.dump id type x y z
fix def2 all print 100 “${e1} ${e2} ${e3}” file initial-pe-comp screen no
thermo 2000
thermo_style custom step lx ly lz press pe temp
run 50000
undump 1
unfix def2
variable tmp equal “lz”
variable L0 equal ${tmp}
print “Initial Length, L0: ${L0}”
#-------------------------COMPRESSION---------------------------------
reset_timestep 0
timestep 0.002
dump 2 all custom 500 ./config-comp/comp..dump id type x y z c_1[] c_v[]
restart 500000 ./restart-comp/compression..equil
variable strain equal “-(lz - v_L0)/v_L0”
variable e1 equal “v_strain”
variable s1 equal “pzz/10000”
variable s2 equal “pxx/10000”
variable s3 equal “pyy/10000”
#variable s2 equal “4v_s1lx/(ly*PI)”
fix def1 all print 1000 “${e1} ${s1} ${s2} ${s3}” file Ti1000-comp screen no
thermo 100000
thermo_style custom step temp ke pe press
velocity upper set 0 0 -0.15
velocity mobile ramp vz 0.0 -0.15 z 29 317 sum yes
run 200000
print “ALL DONE!”
lammps问题-关于NVT与NVE系综产生的不同相关推荐
- 了解lammps中NVE/NVT/NPT三种系综的特性和区别
了解lammps中NVE/NVT/NPT三种系综的特性和区别: 系综是lammps模拟中的一个重要环节,系综种类及对应参数选择的直接决定了模拟结果程度的好坏.然而,由于各种体系.特定条件.系综种类和参 ...
- lammps一对一课程学习大纲
对于初学者来说,自学lammps入门会有些困难. 大多数情况是,看视频学习觉得都能听懂,自己写in文件的时候却又无从下手. 一个小小的错误可能也会困扰你好几天,如果有人指导,可能只需要几分钟就能解决. ...
- lammps案例:团簇融化过程模拟案例学习
分享一个团簇融化的lammps模拟案例,由Carsten Svaneborg博士编写. 模拟单位为LJ,2D模拟体系,周期性边界条件,力场为最常见的LJ力场. 模拟思路: (1)设置一个4040的模拟 ...
- lammps案例:Al冲击过程模拟案例代码
大家好,我是小马老师. 本文介绍lammps模拟冲击的一个案例. 冲击是指一个物体以较快的速度冲向另一个物体,这个过程会有温度的升高,因此,在lammps冲击模拟中,冲击过程一般采用nve系综. 冲击 ...
- Lammps命令与in文件
一.各类文件 1 in :程序写入 2 log :写入状态信息(if the switch is used?) 3 screen : 决定结果的是否进行屏幕输出 4 var name : 定义一个变 ...
- 分子动力学系综小结 (转)
系综(ensemble) 是指在一定的宏观条件下(约束条件),大量性质和结构完全相同的.处于各种运动状态的.各自独立的系统的集合.全称为统计系综.系综是用统计方法描述热力学系统的统计规律性时引入的一个 ...
- 避免lammps弛豫出错的一个小技巧
大家好,我是小马老师. 本文介绍避免lammps弛豫出错的一个方法. 在lammps模拟中,原子的移动速度不能太大,如果在一个timestep内移动距离过大,容易造成模拟出错. "bond ...
- lammps笔记-入门
1 输入文件 输入脚本文件(命令文件?)要包含五个部分: # 1) Initialization # 2) System definition # 3) Simulation settings # 4 ...
- lammps案例:流体绕过障碍物模拟代码讲解
大家好,我是小马老师. 本文介绍lammps官方案例中的一个流体案例,在模型中设置了两个球形障碍物,模拟流体原子绕过障碍物的过程. 本文对案例代码进行注释,以供初学者参考. 本例lammps in文件 ...
最新文章
- C/C++中static关键字用法汇总
- Android上传队列使用
- PowerDesigner建立UML序列图
- 根据时间格式字符串取出时分秒各自的数值
- Java 的Tuple(类似.net等的元组)
- 12月29日二周五次【Python基础语法】
- centos 修改密码_centos7忘记root密码怎么改
- Java基础(1):Java简介和开发环境配置
- 远程桌面服务器连接失败,Windows服务器远程桌面连接失败是什么原因
- Semantic UI 之 下拉菜单 dropdown
- linux BT面板的安装
- 第三次作业_03功能测试
- 求职过程中展示您的社交技巧
- Python可视化--条形图
- Easy_RL 01:强化学习基础
- 彩屏显示入门:我要五彩斑斓的黑 | 用Arduino玩转ESP32与掌控板系列
- DedeCMS找后台目录漏洞
- vue 根据字符串生成表单_vue自定义表单生成器,可根据json参数动态生成表单
- 第一课 request传参
- 基于QT实现的图元拖曳、定点滚轮旋转和缩放
热门文章
- 关于投篮的数学建模模型_球员投篮稳定性的数学模型研究
- 车联网CAN远程云网关给EPEC控制器远程升级
- 迈克耳孙干涉、光的衍射、艾里斑、瑞利判据和光栅(大学物理笔记)
- ION-DTN AMS协议介绍
- 揭秘IT领域:挑战,机遇与未来
- 用高效联接推进数字化转型,华为云会议“会聚”服务力量
- 举例说明应用计算机技术求解问题的步骤,CCCF专栏 | 裘宗燕:计算机问题求解的三类方法...
- 《软件测试技术大全》一书的目录
- 2020年,突然迷上了无线电,想做一名火腿 HAM
- 论文精度笔记(六):《YOLOv7: Trainable bag-of-freebies sets new state-of-the-art for real-time object detecto》