国产主播欧美精品,在线视频cao,中文字幕免费一区二区,日韩三级影视

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

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

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

    第一步：通過以下命令產生隨機位移，筆者這里為了節省算力超胞設的是111。


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

    第二步：將所有產生的images用批量腳本提交任務進行計算。

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

    任務批量提交腳本：


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

    第三步：提取三階力常數。


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

    第四步：提取二階力常數。


    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里作圖就可以了。

    這是筆者選用不同贗勢，然后利用該方法計算的結果：

    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”，專注【材料測試+模擬計算】科研服務，測試狗團隊開發，提供同步輻射、球差電鏡和常用材料表征，承接第一性原理計算、分子動力學、有限元計算等
    評論 / 文明上網理性發言
    12條評論
    全部評論 / 我的評論
    最熱 /  最新
    全部 3小時前 四川
    文字是人類用符號記錄表達信息以傳之久遠的方式和工具?，F代文字大多是記錄語言的工具。人類往往先有口頭的語言后產生書面文字，很多小語種，有語言但沒有文字。文字的不同體現了國家和民族的書面表達的方式和思維不同。文字使人類進入有歷史記錄的文明社會。
    點贊12
    回復
    全部
    查看更多評論
    相關文章

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

    2020-05-03

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

    2021-01-22

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

    2021-06-19

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

    2019-10-25

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

    2019-10-25

    Zeta電位的基本理論、測試方法和應用

    2020-08-24

    熱門文章/popular

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

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

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

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

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

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

    微信掃碼分享文章
    主站蜘蛛池模板： 牟定县| 巩留县| 荆门市| 新密市| 永和县| 淅川县| 红原县| 高平市| 苗栗县| 双城市| 平遥县| 永年县| 同仁县| 休宁县| 荃湾区| 义马市| 静海县| 吉水县| 达日县| 新巴尔虎左旗| 嘉禾县| 崇文区| 武汉市| 宝兴县| 永春县| 比如县| 无为县| 衢州市| 临洮县| 稷山县| 海安县| 闻喜县| 朝阳区| 且末县| 景谷| 阳原县| 定西市| 商河县| 四子王旗| 杭州市| 军事|