預(yù)存
              Document
              當(dāng)前位置:文庫百科 ? 文章詳情
              實用教程丨如何用VASP+phono3py快速計算晶格熱導(dǎo)率?
              來源:科學(xué)10分鐘 時間:2021-08-17 21:31:06 瀏覽:6051次

              在上期的文章:理論到實踐:VASP+Phono3py計算PhononLifetime中,筆者介紹了怎么用VASP+phono3py計算聲子壽命和頻率的關(guān)系。本文繼續(xù)介紹利用VASP+phono3py計算晶格熱導(dǎo)率。晶格熱到率的計算公式如下[1].

              由于該計算過程特別消耗計算資源,所以筆者用虛擬機(jī)做示范,用含兩個原子的金屬W單包為例,且在有限位移計算力常數(shù)的過程中沒有擴(kuò)包,讀者要在服務(wù)器上嘗試建議擴(kuò)成3x3x3的超胞,保證聲子沒有虛頻。

              第一步:通過以下命令產(chǎn)生隨機(jī)位移,筆者這里為了節(jié)省算力超胞設(shè)的是111。


              phono3py  -d --dim="1 1 1" -c POSCAR

              第二步:將所有產(chǎn)生的images用批量腳本提交任務(wù)進(jìn)行計算。

              計算的INCAR:


              PREC = Accurate    IBRION = -1    NELMIN = 5     ENCUT = 500     EDIFF = 1.000000e-08    ISMEAR = 0     SIGMA = 5.000000e-02     IALGO = 38     LREAL = .FALSE.   ADDGRID = .TRUE.     LWAVE = .FALSE.    LCHARG = .FALSE.      NPAR = 4

              KPOINTS:


              Automatic generation0M-P5 5 5 0 0 0

              任務(wù)批量提交腳本:


              #!/bin/bashfor i in {00001..00005}domkdir $i      cp POSCAR-$i ./$i/POSCAR      cp INCAR KPOINTS POTCAR ./$i      cd ./$impirun -np 12 vaspcd ../done

              第三步:提取三階力常數(shù)。


              phono3py --cf3 {00001..00005}/vasprun.xml

              第四步:提取二階力常數(shù)。


              phono3py --dim="1 1 1" --sym-fc -c POSCAR

              第五步:計算熱導(dǎo)率。


              phono3py  --fc3  --fc2  --dim="1 1 1"  --mesh="19 19 19"   -c POSCAR   --br | tee out.txt

              打開out.txt,下拉到底就可以看到計算的熱導(dǎo)率了,直接復(fù)制到originlab里作圖就可以了。

              這是筆者選用不同贗勢,然后利用該方法計算的結(jié)果:

              Reference:

              [1] Togo, Atsushi; Chaput, Laurent; Tanaka, Isao (2015). Distributions of phonon lifetimes in Brillouin zones. Physical Review B, 91(9), 094306–. doi:10.1103/physrevb.91.094306. 

              PS:推薦一個非常好用的科研小程序“測試GO”,專注【材料測試+模擬計算】科研服務(wù),測試狗團(tuán)隊開發(fā),提供同步輻射、球差電鏡和常用材料表征,承接第一性原理計算、分子動力學(xué)、有限元計算等
              評論 / 文明上網(wǎng)理性發(fā)言
              12條評論
              全部評論 / 我的評論
              最熱 /  最新
              全部 3小時前 四川
              文字是人類用符號記錄表達(dá)信息以傳之久遠(yuǎn)的方式和工具。現(xiàn)代文字大多是記錄語言的工具。人類往往先有口頭的語言后產(chǎn)生書面文字,很多小語種,有語言但沒有文字。文字的不同體現(xiàn)了國家和民族的書面表達(dá)的方式和思維不同。文字使人類進(jìn)入有歷史記錄的文明社會。
              點贊12
              回復(fù)
              全部
              查看更多評論
              相關(guān)文章

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

              2020-05-03

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

              2021-01-22

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

              2021-06-19

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

              2019-10-25

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

              2019-10-25

              Zeta電位的基本理論、測試方法和應(yīng)用

              2020-08-24

              熱門文章/popular

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

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

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

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

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

              電化學(xué)實驗基礎(chǔ)之電化學(xué)工作站篇 (二)三電極和兩電極體系的搭建 和測試

              微信掃碼分享文章
              久久激情亚洲精品无码?V| 亚洲国产成人久久综合一区| 国产成人不卡亚洲精品91| 亚洲毛片在线免费观看| 亚洲熟女乱综合一区二区| 亚洲中文字幕无码亚洲成A人片| 亚洲激情在线视频| 亚洲性久久久影院| 亚洲人AV永久一区二区三区久久| 亚洲色偷偷色噜噜狠狠99网| 国产亚洲精品bv在线观看| 亚洲小说图区综合在线| 亚洲日韩国产一区二区三区在线| 亚洲愉拍一区二区三区| 亚洲私人无码综合久久网| 亚洲另类无码专区首页| 亚洲国产精品18久久久久久| 亚洲hairy多毛pics大全| 久久亚洲精品11p| 国产亚洲成在线播放va| 亚洲AV无码乱码在线观看性色扶| 国产亚洲日韩在线a不卡| 亚洲精品A在线观看| 亚洲七七久久精品中文国产| 国产精品亚洲不卡一区二区三区| 中文字幕专区在线亚洲| 亚洲人成色7777在线观看| 人人狠狠综合久久亚洲88| 久久久久亚洲AV成人无码网站 | 国产精品xxxx国产喷水亚洲国产精品无码久久一区 | 在线91精品亚洲网站精品成人| 亚洲av综合av一区二区三区 | 久久久久亚洲精品无码网址色欲| 在线观看免费亚洲| 国产亚洲情侣一区二区无码AV| 亚洲一区精品无码| 亚洲人成中文字幕在线观看| 亚洲av无码一区二区三区网站| 亚洲高清免费在线观看| 色噜噜亚洲男人的天堂| 激情无码亚洲一区二区三区 |