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

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

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

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


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

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

              計算的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

              第五步:計算熱導率。


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

              打開out.txt,下拉到底就可以看到計算的熱導率了,直接復制到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ù),測試狗團隊開發(fā),提供同步輻射、球差電鏡和常用材料表征,承接第一性原理計算、分子動力學、有限元計算等
              評論 / 文明上網(wǎng)理性發(fā)言
              12條評論
              全部評論 / 我的評論
              最熱 /  最新
              全部 3小時前 四川
              文字是人類用符號記錄表達信息以傳之久遠的方式和工具?,F(xiàn)代文字大多是記錄語言的工具。人類往往先有口頭的語言后產(chǎn)生書面文字,很多小語種,有語言但沒有文字。文字的不同體現(xiàn)了國家和民族的書面表達的方式和思維不同。文字使人類進入有歷史記錄的文明社會。
              點贊12
              回復
              全部
              查看更多評論
              相關(guān)文章

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

              2020-05-03

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

              2021-06-19

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

              2021-01-22

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

              2019-10-25

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

              2019-10-25

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

              2020-11-03

              熱門文章/popular

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

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

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

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

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

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

              微信掃碼分享文章
              亚洲精品久久无码av片俺去也| 亚洲精品456人成在线| 精品亚洲AV无码一区二区| 亚洲成人中文字幕| 亚洲成a人一区二区三区| 亚洲欧洲AV无码专区| 亚洲影视一区二区| 久久精品国产亚洲AV无码娇色| 国产v亚洲v天堂无码网站| 国产亚洲精品a在线无码| 亚洲综合av永久无码精品一区二区| 亚洲色欲久久久久综合网| 亚洲国产高清在线一区二区三区| 四虎亚洲国产成人久久精品| 99亚洲乱人伦aⅴ精品| 女bbbbxxxx另类亚洲| 美国毛片亚洲社区在线观看| 精品国产亚洲一区二区三区在线观看 | 国产精品亚洲综合网站| 人人狠狠综合久久亚洲| 久久精品国产亚洲AV| 四虎精品亚洲一区二区三区| 亚洲男女内射在线播放| 久久久久无码专区亚洲av| 亚洲性猛交XXXX| 亚洲av永久无码精品古装片 | 国产亚洲一区二区三区在线不卡| 亚洲国产成人精品女人久久久| 亚洲成aⅴ人片久青草影院| 亚洲精品国产va在线观看蜜芽| 国产乱辈通伦影片在线播放亚洲 | 亚洲日韩在线视频| 色婷五月综激情亚洲综合| 亚洲色大成网站www永久网站| 亚洲欧美日韩中文高清www777 | 亚洲私人无码综合久久网| AV激情亚洲男人的天堂国语| 亚洲国产成人影院播放| 亚洲精品国产成人片| 精品亚洲国产成AV人片传媒| 亚洲一区二区三区四区视频|
              +

              你好,很高興為您服務(wù)!

              發(fā)送