分享一个团簇融化的lammps模拟案例,由Carsten Svaneborg博士编写。

模拟单位为LJ,2D模拟体系,周期性边界条件,力场为最常见的LJ力场。

模拟思路:
(1)设置一个4040的模拟盒子,在盒子中央设置一个2020的区域,在中间区域内随机填充原子。
(2)能量最小化消除随机原子的重叠。
(3)设置LJ力场
(4)设置nve系综,使用langevin控温方式进行控温。这一步也可以直接用nvt系综代替。
(5)运行并输出轨迹文件。

案例比较简单,需要注意的是:采用随机方法生成原子,原子之间容易重叠,一定要在正式模拟之前消除原子重叠,否则原子距离过近导致排斥力过大,容易产生“lost atoms”错误。
消除重叠可使用“delete_atoms verlap”和“minimize”命令。
案例代码已注释,仅供参考。

#原子数量
variable npart  equal 500
#LJ单位温度
variable tempstart   equal 0.5
variable tempstop    equal 0.5
#模型基本设置
units       lj
dimension 2
atom_style  atomic
boundary        p p p
neighbor        0.5   bin
neigh_modify    every 1 delay 0 check yes
#模拟盒子设置
region box block -20 20 -20 20 -0.1 0.1
create_box 1 box
#2维计算设置
fix 2d  all enforce2d
#在box中央设置box2
region box2 block -10 10 -10 10 -0.1 0.1
#在box2中随机生成原子
create_atoms 1 random ${npart} 324523 box2
#摩尔质量
mass   *    1
#LJ力场设置
pair_style lj/cut 2.5
pair_coeff * * 1.0 1.0
#能量最小化
minimize 1e-4 1e-4 1000 1000
reset_timestep 0
#保存原子轨迹
dump        img all  atom  1000 dump.xyz
#nve系综,langevin控温
fix integrator all nve
fix dynamics   all langevin ${tempstart} ${tempstop} 1.0 2352325
#设置屏幕输出热力学信息
thermo_style custom step temp  ke pe
thermo      100
#模拟步长设置
timestep    0.01
#运行20000步
run     20000

更多lammps案例,关注微信公众号:lammps加油站

lammps案例:团簇融化过程模拟案例学习相关推荐

  1. lammps案例:Al冲击过程模拟案例代码

    大家好,我是小马老师. 本文介绍lammps模拟冲击的一个案例. 冲击是指一个物体以较快的速度冲向另一个物体,这个过程会有温度的升高,因此,在lammps冲击模拟中,冲击过程一般采用nve系综. 冲击 ...

  2. lammps计算团簇的方法

    大家好,我是小马老师. 本文介绍在lammps中如何进行团簇的分析. 在一个group内,原子之间的距离小于一定距离时,可以认为这组原子时一个团簇. 在lammps中,提供了团簇计算命令: compu ...

  3. Lammps的团簇(cluster)之Kmean识别—Matlab分析

    关注 M r . m a t e r i a l   , \color{Violet} \rm Mr.material\ , Mr.material ,

  4. 【案例教程】山洪径流过程模拟及洪水危险性评价技术

    GIS水文分析(ArcHydro.Spatial Anlysist等模块)是流域水文模拟建模的重要工具,能够自动提取及计算流域边界.河网水系.流向.汇流时间和其它流域特征参数.美国陆军工程兵团开发的开 ...

  5. airpak模拟案例,Airpak模拟教程-体育馆通风模拟案例-CFD数值模拟教程airpak.pdf

    Airpak模拟教程-体育馆通风模拟案例-CFD数值模拟教程airpak.pdf Airpak 模拟案例模拟案例 CFD 模拟案例模拟案例 体育馆通风案例教程体育馆通风案例教程 1 Airpak 软件 ...

  6. 用计算机来模拟二维生命,用并行计算机模拟二维氢团簇的光电离

    摘要: 随着超短超强激光技术的发展,团簇与超短超强激光的相互作用一直是近些年来的研究热点之一.作为介于原子分子与固体之间的特殊物质形态,原子和分子团簇的特殊性质使得激光能量能够非常有效的转化到离子与电 ...

  7. aspen软件对电脑的要求_学会使用Aspen Plus自带案例,让你的模拟事半功倍

    原标题:学会使用Aspen Plus自带案例,让你的模拟事半功倍 前言 我们科普化工的小伙伴们在前面给大家介绍了很多关于Aspen Plus的基础教程和应用案例,通过这些内容可以让大家对Aspen P ...

  8. airpak模拟案例,Airpak模拟教程-体育馆通风模拟案例-CFD数值模拟教程airpak

    Airpak模拟教程-体育馆通风模拟案例-CFD数值模拟教程airpak Airpak 模拟案例模拟案例 CFD 模拟案例模拟案例 体育馆通风案例教程体育馆通风案例教程 1 Airpak 软件中文学习 ...

  9. 高斯移动热源,熔池流场分布,comsol模拟案例

    高斯移动热源,熔池流场分布,comsol模拟案例. ,里面涉及到流体传热,结构网格剖分,熔化凝固,流场. ID:69260633198808345小仙女的小宠物店

最新文章

  1. 理解卷积神经网络的局限
  2. 禁用windows更新完成后的重启提示
  3. VC中使用Unicode的一些列问题
  4. alert点击完确定关闭浏览器窗口_关务实操技巧:单一窗口缓存4大清理方法,解决部分退单问题...
  5. C语言中extern修饰符的用法
  6. mysql数据库基础的简单操作指南
  7. 1.6编程基础之一维数组_04数组逆序重放
  8. 【精】【火】关于CSDN博客与博主的第二个博客之间的区别与联系
  9. 电脑音响怎么插_BMW宝马5系G38改原厂全套哈曼卡顿音响+无钥匙进入,厚街宝马原厂改装中心...
  10. lodash是否是数组_lodash源码解析:for家族
  11. 转 XMLHttpRequest().readyState的五种状态详解
  12. 转换整形数字为16进制字符串
  13. 阶段1 语言基础+高级_1-3-Java语言高级_06-File类与IO流_05 IO字符流_7_字符输出流的续写和换行...
  14. 从Gaussian 09的Hartree-Fock计算的fchk文件转换到Columbus的mocoef文件的脚本
  15. iphone6出信号 无服务器,苹果6手机无服务(非常见故障点)维修
  16. 社会新现象:“蜗婚”离婚不分家
  17. 【问题解决】xlwings处理excel复制粘贴时数字自动变成科学计数法
  18. 虚拟机的Nginx安装配置
  19. 一篇文章教会大家制作小程序,利用小程序创业。
  20. 2012年9月 《UDID罗生门》

热门文章

  1. 火狗工房(注意哦,不是热狗)
  2. 3大功能全面升级!Cinema 4D S24好用在哪?
  3. Shader视差贴图
  4. 让Project把周六和周日也计算工期
  5. workbook加载文件路径_EXCEL VBA 引用文件路径问题,请来帮忙,谢谢!
  6. vertx编程需注意的点
  7. 动态web技术(二) --- CGI
  8. 为什么软件测试计划重要性,谈谈“测试计划”这点事
  9. 【MySQL】mysql | 年月日时分秒 | 日期格式化
  10. 四、GStreamer基础