LAMMPS实操系列(一): 大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构筛选
高熵合金(High-entropy alloys,HEA),是由五种或五种以上等量或大约等量金属形成的合金。由于高熵合金可能具有许多优秀的物理化学性质,因此在材料科学领域受到广泛关注。一般的合金中,金属成分一般只有一到两种,例如铁基、镍基合金,它们一般以铁或镍为基础,再加入一些量的元素来提升其性能。以往的观念中,若合金中加的金属种类越多,会使其材质脆化,但高熵合金和以往的合金不同,有多种金属却不会脆化,是一种新的金属材料。
本章我们将介绍利用LAMMPS结合bash脚本来批量建立具有格点原子无序的FCC-CoCrCuFeNi高熵合金,并进行结构优化并统计能量,得出能量最小结构的过程。
我们的建模过程是先建立10x10x10的FCC-Fe的超胞,然后按比例替换铁原子,最后得到组份相等的高熵合金CoCrCuFeNi,由于高熵合金的格点原子构型是无序的,因此我们需要建立大量的模型,筛选出能量最小的一个。大量建模中需要注意以下几个方面:
(1)我们建模过程中按比例替换Fe原子时,需要给定随机数组,相同的随机数组出来的构型是一样的,因此我们每建立一个新的模型都会给一组和前面不同的随机数。
(2)我们建成千上万的模型不可能手动提交任务因此需要写循环完成。
解决以上问题,我们运用了LAMMPS程序中的shell脚本,并结合LAMMPS的循环命令成功完成。我们的LAMMPS脚本如下:
units metal
boundary p p p
atom_style atomic
timestep 0.001
#构建模型
lattice fcc 3.56
region box block 0 10 0 10 0 10
create_box 5 box
label loop
variable a loop 1000
log log.$a
create_atoms 1 box
set type 1 type/ratio 2 0.2 12499
set type 1 type/ratio 3 0.5 12499
set type 1 type/ratio 4 0.5 12499
set type 3 type/ratio 5 0.5 12499
#设置原子质量
mass 1 55.8 #Fe
mass 2 58.7 #Ni
mass 3 52 #Cr
mass 4 59 #Co
mass 5 63 #Cu
#设置势函数
pair_style eam/alloy
pair_coeff * * FeNiCrCoCu-with-ZBL.eam.alloy Fe Ni Cr Co Cu
#cg和sd交替进行能量最小化
min_style cg
minimize 1.0e-6 1.0e-7 10000 10000
run 0
min_style sd
minimize 1.0e-6 1.0e-7 10000 10000
run 0
min_style cg
minimize 1.0e-6 1.0e-7 10000 10000
run 0
write_data FeNiCrCoCu_$a.data
group ALL type 1 2 3 4 5
delete_atoms group ALL
#为建立新模型调用shell脚本更换随机数
shell ./mseed.sh
next a
jump FeNiCrCoCu.in loop
通过运行该脚本,我们获得了1000个模型如图1所示。同时,通过共轭梯度和快速下降方法交替对每个模型进行优化,获得其结构优化的能量。如图2所示,我们通过排序得出了第620号模型的能量是最低的,因此我们选用620号模型作为我们下一步研究工作的模型。图3为我们筛选出来的620号FCC-CoCrCuFeNi高熵合金模型。
图1. 随机建模生成的1000个FCC-CoCrCuFeNi高熵合金模型及其序号
图2. 按能量大小排序后的各序号对应的FCC-CoCrCuFeNi高熵合金模型
图3. 筛选出的620号FCC-CoCrCuFeNi高熵合金模型
本章,我们介绍了对于格点原子无序的FCC-CoCrCuFeNi高熵合金理论模拟中的建模思想,以及LAMMPS的程序实现,下一章我们将详细解读我们本章贴出的脚本,以及结果快速后处理的程序实现。
测试狗文库百科
LAMMPS实操系列(一): 大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构筛选相关推荐
- LAMMPS实操系列(二): 大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构筛选
上一期我们介绍了对于格点原子无序的FCC-CoCrCuFeNi高熵合金理论模拟中的建模思想,以及LAMMPS的程序实现.本期文章,我们会详细解读上一期贴出的脚本,以便让读者更深的了解LAMMPS程序的 ...
- LAMMPS实操系列(三): 大量FCC-CoCrCuFeNi高熵合金建模与结构筛选
前面我们介绍了对于格点原子无序的FCC-CoCrCuFeNi高熵合金理论模拟中的建模思想,以及LAMMPS的程序实现,并详细解读我们编写的LAMMPS脚本,下面我们将继续介绍快速后处理结果的过程以及程 ...
- LAMMPS实操系列(四): 高熵合金FCC-CoCrCuFeNi的退火模拟
我们前几期文章内容详细介绍了通过LAMMPS程序进行大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构的筛选,并最终赛选出了我们所建1000个模型中能量最小也就是最稳定的结构.本章我们将继续通过 ...
- lammps教程:高熵合金建模备用方法
大家好,我是小马老师. 在之前推文中,已经介绍过高熵合金的建模方法: 高熵合金建模方法 有些粉丝朋友反映这个命令不能用,主要原因是set type/ratio命令是lammps的一个新命令,只能用于2 ...
- lammps构建高熵合金模型+结构优化初步筛选能量最小的结构
文章目录 前言 一.ubuntu并行安装Lammps 安装教程 二.单个高熵合金建模 参考资料 操作步骤 三.构建1000个高熵合金模型筛选能量最小的结构 参考资料 操作步骤 1.构建高熵合金模型 2 ...
- lammps 案例:高熵合金FeNiCrCoCu建模及拉伸模拟
高熵合金力学性能模拟包括拉伸.剪切.摩擦等模拟,本文以高熵合金FeNiCrCoCu为例介绍高熵合金建模及拉伸过程模拟. FeNiCrCoCu原子数量比为1:1:1:1:1 高熵合金的拉伸与单晶拉伸过程 ...
- lammps模拟技巧:高熵合金势函数设置三种方法
大家好,我是小马老师. 本文继续介绍lammps模拟中的势函数设置方法. 在前面的文章中,介绍过高熵合金的建模:高熵合金建模 高熵合金包含的原子数较多,势函数的设置相对复杂,在平时一对一的咨询中,这种 ...
- 【直播回顾及资料下载】小程序云应用入门实操系列课程第一讲...
直播详情:小程序云应用入门实操系列课程第一讲:https://yq.aliyun.com/articles/698244 直播时间:2019年4月16日 20:30 直播专家: 白宦成 - Linux ...
- 不同网段的局域网怎么互通_华为实操系列 | 交换机在局域网中是怎么应用的,看完你肯定懂了!...
编辑 | 排版 | 制图 | 测试 | ©瑞哥 此文用时0小时48分钟,原创不易,坚持更不易,希望我的每一份劳动成果都可以得到大家的一个[在看] 交换机在局域网中的应用 交换机在局域网中的应用分为两种 ...
最新文章
- 中文分词工具jieba分词器的使用
- Verilog初级教程(7)Verilog模块例化以及悬空端口的处理
- java项目校内网的描述_[导入]校内网的错误信息
- python新增练习(二)
- SPOJ GCDEX (数论)
- lsof命令_lsof命令的部分说明
- 有1到100共100个数, 从1开始, 每隔1, 2, 3... 个数拿走一个数, 最后剩下几?(约瑟夫环)...
- win2012 定时自动备份mysql_SQL SERVER 2012数据库自动备份的方法
- 南无阿弥陀佛是什么意思呢
- 关于判断单峰数组的几种方法
- Photoshop 更改图片颜色
- 洛谷P1338(末日的传说)
- windows 安装metis_图划分软件Metis的使用(win10+vs2017)
- 左右植树java_Plant 模拟植树活动,编写一个java应用程序 联合开发网 - pudn.com
- http,https,TCP,UDP,apache 的ab,压测工具
- 欧特克开发者训练营(Autodesk DevCamp 2013)的视频教程可以下载
- android5.1设置ip,android5.1 imx6q接并口摄像头(IPU1-CSI1)怎么配置
- you-get使用教程
- 再见,VLOOKUP!2021 Excel数据分析已成气候!
- Ubuntu离线安装Nvidia显卡驱动