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

热线：400-152-6858

测试狗科研服务

预存免费试测登录

Document

当前位置：文库百科 › 文章详情

LAMMPS实操系列（三）: 大量FCC-CoCrCuFeNi高熵合金建模与结构筛选

来源：时间：2022-11-02 16:24:10 浏览：2614次

前面我们介绍了对于格点原子无序的FCC-CoCrCuFeNi高熵合金理论模拟中的建模思想，以及LAMMPS的程序实现，并详细解读我们编写的LAMMPS脚本，下面我们将继续介绍快速后处理结果的过程以及程序实现。

当我们用PBS提交计算任务完成后我们可以看到，计算任务文件夹中多了1000个FCC-CoCrCuFeNi高熵合金的结构data文件以及每个结构优化过程中对应的log文件。由于我们建模时保证了每个模型都是不一样的，因此静力学结构优化后，一定可以将这些体系自由能量从底到高排列，能量最低的那个就是我们需要找的目标结构。任意打开一个log文件，可以看到每一次优化完成之后都会有自由能的输出：

而对于每一个log文件里面会打印三次这样的信息，我们只需要抓取最后次的第三列就可以了。将每次获取的文件编号及能量写成列表，再按能量大小从小到大排序后，第一个就是我们筛选到的最优结构编号。处理程序如下：

#!/bin/bash

echo `echo -e "Number""\t""Energy_(eV)"` >>result.txt

for i in {1..1000}

do

energy=$(echo "`grep -A 1 'Energy initial' log.$i | tail -1|awk '{printf "%12f \n", $3}'`")

echo `echo -e $i"\t"$energy` >> result.txt

done

echo "print(''.join(sorted(open('result.txt'), key=lambda s: s.split()[1],reverse=1)))" > sort.py

python sort.py > sorted.txt

rm sort.py

该程序命名为ad.sh，我们通过再次计算任务文件夹中./ad.sh运行该程序，可得到排序后的文件sorted.txt。打开sorted.txt第一行就是我们要找的能量最小，也就是最优的结构。

从图中可以看到，能量最小最优的结构是620号data文件，我们找到FeNiCrCoCu_620.data文件，通过OVITO可视化软件打开，就可以得到最优结构的可视化图如下：

到此为止，我们利用LAMMPS程序以及bash脚本建立大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构的筛选就介绍完了，接下来我们将根据我们筛选出来的最优结构，做一些课题的研究，敬请各位读者关注我们的更新。

测试狗文库百科

上一篇：LAMMPS实操系列（四）: 高熵合金FCC-CoCrCuFeNi的退火模拟

下一篇：LAMMPS实操系列（二）: 大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构筛选

相关文章

【科研干货】电化学表征：循环伏安法详解（上）

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