av影院在线观看-久久久久久国产精品视频-亚洲综合第一页-日日夜夜天天-快射视频网-国产精品国产-欧美三级午夜理伦三级中视频-国产日韩未满十八禁止观看-亚洲一区国产-亚洲成av人片-久草综合视频-天天射影院-91精品一区二区-国产又粗又硬-黄色91免费-亚洲精品无码永久在线观看-中文字幕人妻高清乱码

預(yù)存

登錄注冊(cè)

Document

當(dāng)前位置：文庫(kù)百科 ? 文章詳情

LAMMPS實(shí)操系列（三）: 大量FCC-CoCrCuFeNi高熵合金建模與結(jié)構(gòu)篩選

來(lái)源：時(shí)間：2022-11-02 16:24:10 瀏覽：4218次

前面我們介紹了對(duì)于格點(diǎn)原子無(wú)序的FCC-CoCrCuFeNi高熵合金理論模擬中的建模思想，以及LAMMPS的程序?qū)崿F(xiàn)，并詳細(xì)解讀我們編寫的LAMMPS腳本，下面我們將繼續(xù)介紹快速后處理結(jié)果的過(guò)程以及程序?qū)崿F(xiàn)。

當(dāng)我們用PBS提交計(jì)算任務(wù)完成后我們可以看到，計(jì)算任務(wù)文件夾中多了1000個(gè)FCC-CoCrCuFeNi高熵合金的結(jié)構(gòu)data文件以及每個(gè)結(jié)構(gòu)優(yōu)化過(guò)程中對(duì)應(yīng)的log文件。由于我們建模時(shí)保證了每個(gè)模型都是不一樣的，因此靜力學(xué)結(jié)構(gòu)優(yōu)化后，一定可以將這些體系自由能量從底到高排列，能量最低的那個(gè)就是我們需要找的目標(biāo)結(jié)構(gòu)。任意打開一個(gè)log文件，可以看到每一次優(yōu)化完成之后都會(huì)有自由能的輸出：

而對(duì)于每一個(gè)log文件里面會(huì)打印三次這樣的信息，我們只需要抓取最后次的第三列就可以了。將每次獲取的文件編號(hào)及能量寫成列表，再按能量大小從小到大排序后，第一個(gè)就是我們篩選到的最優(yōu)結(jié)構(gòu)編號(hào)。處理程序如下：

#!/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，我們通過(guò)再次計(jì)算任務(wù)文件夾中./ad.sh運(yùn)行該程序，可得到排序后的文件sorted.txt。打開sorted.txt第一行就是我們要找的能量最小，也就是最優(yōu)的結(jié)構(gòu)。

從圖中可以看到，能量最小最優(yōu)的結(jié)構(gòu)是620號(hào)data文件，我們找到FeNiCrCoCu_620.data文件，通過(guò)OVITO可視化軟件打開，就可以得到最優(yōu)結(jié)構(gòu)的可視化圖如下：

到此為止，我們利用LAMMPS程序以及bash腳本建立大量FCC-CoCrCuFeNi高熵合金建模與最穩(wěn)定結(jié)構(gòu)的篩選就介紹完了，接下來(lái)我們將根據(jù)我們篩選出來(lái)的最優(yōu)結(jié)構(gòu)，做一些課題的研究，敬請(qǐng)各位讀者關(guān)注我們的更新。

測(cè)試狗文庫(kù)百科

上一篇：LAMMPS實(shí)操系列（四）: 高熵合金FCC-CoCrCuFeNi的退火模擬

下一篇：LAMMPS實(shí)操系列（二）: 大量FCC-CoCrCuFeNi高熵合金建模與最穩(wěn)定結(jié)構(gòu)篩選

相關(guān)文章

基礎(chǔ)理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結(jié)構(gòu)）

2020-05-03

手把手教你用ChemDraw 畫化學(xué)結(jié)構(gòu)式：基礎(chǔ)篇

2021-06-19

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程（下篇）

2021-01-22

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（上）

2019-10-25

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（下）

2019-10-25

XRD的基本原理與應(yīng)用

2020-11-03

舉報(bào)有獎(jiǎng)

TEL: 191-3608-6524

如：在網(wǎng)絡(luò)上惡意使用“測(cè)試狗”等相關(guān)關(guān)鍵詞誤導(dǎo)用戶點(diǎn)擊、惡意盜用測(cè)試狗商標(biāo)、冒稱官方工作人員等情形，請(qǐng)您向我們舉報(bào)，經(jīng)查實(shí)后，我們將給予您獎(jiǎng)勵(lì)。

舉報(bào)內(nèi)容：

200

上傳附件：

文件格式不正確，請(qǐng)重新上傳文件格式不正確，請(qǐng)重新上傳文件格式不正確，請(qǐng)重新上傳

文件格式：jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar

聯(lián)系方式

姓名

電話

提交意見

意見反饋

Suggestions

您可以在此留下您寶貴的意見，您的意見或問(wèn)題反饋將會(huì)成為我們不斷改進(jìn)的動(dòng)力。

意見類型

測(cè)試服務(wù)

網(wǎng)站功能

財(cái)務(wù)報(bào)賬

其他類型

意見內(nèi)容

200

聯(lián)系方式

姓名

電話

提交意見

+

你好，很高興為您服務(wù)！

發(fā)送

項(xiàng)目推薦

舉報(bào)有獎(jiǎng)

意見反饋