快捷下单入口 关于 合作 招聘 新人手册 会员中心

    热线:400-152-6858

    测试狗科研服务

    预存 免费试测 登录
    Document
    当前位置:文库百科文章详情
    LAMMPS实操系列(二): 大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构筛选
    来源: 时间:2022-11-02 16:40:23 浏览:2081次

    上一期我们介绍了对于格点原子无序的FCC-CoCrCuFeNi高熵合金理论模拟中的建模思想,以及LAMMPS的程序实现。本期文章,我们会详细解读上一期贴出的脚本,以便让读者更深的了解LAMMPS程序的运行过程。下面就是程序实现过程的详细解读与分析:

    units      metal      #定义单位系统为metal

    boundary      p p p        #定义三个方向都为周期性边界条件

    atom_style    atomic       #定义原子类型为atomic

    timestep      0.001        #设定时间步长为1 fs

    #构建模型

    lattice       fcc 3.56       #定义fcc晶格常数

    region        box block 0 10 0 10 0 10   #定义模拟盒子大小

    create_box    5 box        #定义模拟盒子中最多可以包含5种原子

    label       loop       #定义脚本循环起始标签为loop

    variable      a loop 1000      #定义跟踪循环次数变量a为1至1000     

    log        log.$a       #定义记录标准输出文件的文件名以次数结尾

    create_atoms   1 box   #现在模拟盒子中将格点原子用第1种原子Fe代替

    set    type 1 type/ratio 2 0.2 7603  #将1/5的Fe原子换成第2种原子Ni

    set    type 1 type/ratio 3 0.5 7603  #将1/2的Fe原子换成第3种原子Cr

    set    type 1 type/ratio 4 0.5 7603  #将1/2的Fe原子换成第4种原子Co

    set    type 3 type/ratio 5 0.5 7603  #将1/2的Cr原子换成第5种原子Cu

    #注意set进行的是依次操作,这样下来Ni、Cr、Co、Cu各占1/5

    #设置原子质量

    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     #选用eam合金势

    pair_coeff * * FeNiCrCoCu-with-ZBL.eam.alloy Fe Ni Cr Co Cu

    #选择文献提供的Fe Ni Cr Co Cu势函数[1]

    min_style     cg     #选用共轭梯度法对模拟盒子进行结构优化

    minimize      1.0e-6 1.0e-7 10000 10000

    #能量和力的容差为分别为1.0e-6 和1.0e-7,最大迭代次数为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      #写出第a次循环中结构优化的模型

    group ALL type 1 2 3 4 5         #将模拟盒子中的原型设为group ALL

    delete_atoms group ALL        #删除模拟盒子中的所有原子

    shell ./mseed.sh         #跟换下一轮建模计算的随机随机数

    next      a        #跳到下轮计算的跟中变量符a

    jump FeNiCrCoCu.in loop         

    #跳到FeNiCrCoCu.in loop文件中的label loop开始下一轮计算

    在以上LAMMPS脚本中,需要特别注意的是,命令“set    type 1 type/ratio 2 0.2 7603”执行的结果与随机数7603有关,相同的随机数通过该命令生成的结果是一样的因此我们在生成下一个模型时需要提前对随机数进行更改。我们的更改方式是通过外部shell脚本mseed.sh进行更改的。调用命令为shell ./mseed.sh。mseed.sh中的内容如下:

    #!/bin/bash

    a=$RANDOM

    sed -i "13c set    type 1 type/ratio 2 0.2 $a" FeNiCrCoCu.in

    sed -i "14c set    type 1 type/ratio 3 0.5 $a" FeNiCrCoCu.in

    sed -i "15c set    type 1 type/ratio 4 0.5 $a" FeNiCrCoCu.in

    sed -i "16c set    type 3 type/ratio 5 0.5 $a" FeNiCrCoCu.in

    该脚本中,变量RANDOM是Linux shell中的一个特殊变量,每一次对其进行取值都会任意取到0-32767之间的任意整数,且在同一个shell进程取值不会重复,取完为止。后面四行通过sed命令将拿到的随机用来修改FeNiCrCoCu.in文件,这样就可以保证在循环运行FeNiCrCoCu.in脚本时,生成的结构都是不同的。

    最后提交任务时只需要在任务文件夹中放入FeNiCrCoCu.in、mseed.sh、FeNiCrCoCu-with-ZBL.eam.alloy用PBS脚本进行提交就可以了。其中FeNiCrCoCu-with-ZBL.eam.alloy为eam势函数文件,该势函数来自于O.R. Deluigi等人的文章[1],势文件可以直接从网上下载[2]。


    下一期文章我们将继续介绍脚本运行完成之后的后处理,筛选出最稳定­结构的FCC-CoCrCuFeNi高熵合金。

    评论 / 文明上网理性发言
    12条评论
    全部评论 / 我的评论
    最热 /  最新
    全部 3小时前 四川
    文字是人类用符号记录表达信息以传之久远的方式和工具。现代文字大多是记录语言的工具。人类往往先有口头的语言后产生书面文字,很多小语种,有语言但没有文字。文字的不同体现了国家和民族的书面表达的方式和思维不同。文字使人类进入有历史记录的文明社会。
    点赞12
    回复
    全部
    查看更多评论
    相关文章

    【科研干货】电化学表征:循环伏安法详解(上)

    2019-10-25

    晶体结构可视化软件 VESTA使用教程(下篇)

    2021-01-22

    手把手教你用ChemDraw 画化学结构式:基础篇

    2021-06-19

    【科研干货】电化学表征:循环伏安法详解(下)

    2019-10-25

    基础理论丨一文了解XPS(概念、定性定量分析、分析方法、谱线结构)

    2020-05-03

    Zeta电位的基本理论、测试方法和应用

    2020-08-24

    热门文章/popular

    【科研干货】电化学表征:循环伏安法详解(上)

    晶体结构可视化软件 VESTA使用教程(下篇)

    手把手教你用ChemDraw 画化学结构式:基础篇

    电化学实验基础之电化学工作站篇 (二)三电极和两电极体系的搭建 和测试

    【科研干货】电化学表征:循环伏安法详解(下)

    基础理论丨一文了解XPS(概念、定性定量分析、分析方法、谱线结构)

    微信扫码分享文章

    意见反馈

    有奖举报

    商务合作

    ...

    更多

    公众号

    关注我们 了解更多

    小程序

    随时预约 掌握进度

    举报有奖

    TEL: 191-3608-6524

    如:在网络上恶意使用“测试狗”等相关关键词误导用户点击、恶意盗用测试狗商标、冒称官方工作人员等情形,请您向我们举报,经查实后,我们将给予您奖励。

    举报内容:

    200

    上传附件:
    文件格式不正确,请重新上传文件格式不正确,请重新上传文件格式不正确,请重新上传
    文件格式:jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar
    联系方式
    姓名
    电话
    提交意见

    意见反馈

    Suggestions

    您可以在此留下您宝贵的意见,您的意见或问题反馈将会成为我们不断改进的动力。

    意见类型
    测试服务
    网站功能
    财务报账
    其他类型
    意见内容

    200

    联系方式
    姓名
    电话
    提交意见

    收起

    01

    专属信用额度,先测后付0元下单

    02

    下单享高额积分,万千好礼免费兑

    200
    200元无门槛优惠券
    立即激活 立即下单
    已使用
    已作废

    全流程在线可视化,便捷高效触手可及

    如下单过程中有任何疑问或需要帮助,请随时咨询专属顾问~
    9
    9折无门槛优惠券

    支付一笔订单后可领取

    立即领取 立即下单
    已使用
    已作废

    免费测+惊喜盲盒+高额福利,多重福利大放送

    新人免费测

    双双开盲盒(100%中奖)

    邀请人专享

    受邀人专享

    创建/加入团队,解锁定制化权益

    01

    1500元团队专属优惠券

    02

    万元大额信用额度,享先测后付

    03

    团队成员统一开票报销;

    04

    专业工程师课题专属服务

    领取成功,请下单
    请您支付一笔订单后才可以领取优惠券
    Document
    关于我们 新手帮助 测试干货 商务合作 基金查询 相关资质 模拟计算 现场测试 服务项目 科研绘图 同步辐射 电池行业

    联系方式/contact

    400-152-6858

    工作时间/work time

    09:00-18:00

    测试狗公众号

    关注我们 了解更多

    测试狗小程序

    随时预约 掌握进度

    蜀公网安备51010602000648号

    蜀ICP备17005822号-1

    成都世纪美扬科技有限公司

    Copyright@测试狗·科研服务