www.毛片在线观看_怡红院亚洲红怡院天堂麻豆_欧美jizz18欧美_国产精品免费久久久久影院

預存

登錄注冊

Document

當前位置：文庫百科 ? 文章詳情

LAMMPS實操系列（二）: 大量FCC-CoCrCuFeNi高熵合金建模與最穩定結構篩選

來源：時間：2022-11-02 16:40:23 瀏覽：5205次

上一期我們介紹了對于格點原子無序的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高熵合金。

上一篇：LAMMPS實操系列（三）: 大量FCC-CoCrCuFeNi高熵合金建模與結構篩選

下一篇：LAMMPS實操系列(一): 大量FCC-CoCrCuFeNi高熵合金建模與最穩定結構篩選

2020-05-03

手把手教你用ChemDraw 畫化學結構式：基礎篇

2021-06-19

晶體結構可視化軟件 VESTA使用教程（下篇）

2021-01-22

【科研干貨】電化學表征：循環伏安法詳解（上）

2019-10-25

【科研干貨】電化學表征：循環伏安法詳解（下）

2019-10-25

XRD的基本原理與應用

2020-11-03

熱門文章/popular

基礎理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結構）

手把手教你用ChemDraw 畫化學結構式：基礎篇

晶體結構可視化軟件 VESTA使用教程（下篇）

【科研干貨】電化學表征：循環伏安法詳解（上）

電化學實驗基礎之電化學工作站篇（二）三電極和兩電極體系的搭建和測試

【科研干貨】電化學表征：循環伏安法詳解（下）

微信掃碼分享文章

www.毛片在线观看_怡红院亚洲红怡院天堂麻豆_欧美jizz18欧美_国产精品免费久久久久影院

舉報有獎

TEL: 191-3608-6524

如：在網絡上惡意使用“測試狗”等相關關鍵詞誤導用戶點擊、惡意盜用測試狗商標、冒稱官方工作人員等情形，請您向我們舉報，經查實后，我們將給予您獎勵。

舉報內容：

200

上傳附件：

文件格式不正確，請重新上傳文件格式不正確，請重新上傳文件格式不正確，請重新上傳

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

聯系方式

姓名

電話

提交意見

意見反饋

Suggestions

您可以在此留下您寶貴的意見，您的意見或問題反饋將會成為我們不斷改進的動力。

意見類型

測試服務

網站功能

財務報賬

其他類型

意見內容

200

聯系方式

姓名

電話

提交意見

+

你好，很高興為您服務！

發送

www.毛片在线观看_怡红院亚洲红怡院天堂麻豆_欧美jizz18欧美_国产精品免费久久久久影院

項目推薦

舉報有獎

意見反饋