双腿间已经湿成一片,亚洲啪啪综合AV一区,厨房的春潮A片,日本一道久久高清国产

快捷下單入口關(guān)于合作招聘新人手冊(cè) 會(huì)員中心

熱線：400-152-6858

測(cè)試狗科研服務(wù)

預(yù)存免費(fèi)試測(cè) 登錄

Document

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

實(shí)用教程丨如何用VASP+phono3py快速計(jì)算晶格熱導(dǎo)率？

來源：科學(xué)10分鐘時(shí)間：2021-08-17 21:31:06 瀏覽：3812次

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

由于該計(jì)算過程特別消耗計(jì)算資源，所以筆者用虛擬機(jī)做示范，用含兩個(gè)原子的金屬W單包為例，且在有限位移計(jì)算力常數(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)行計(jì)算。

計(jì)算的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

第五步：計(jì)算熱導(dǎo)率。

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

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

這是筆者選用不同贗勢(shì)，然后利用該方法計(jì)算的結(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：推薦一個(gè)非常好用的科研小程序“測(cè)試GO”，專注【材料測(cè)試+模擬計(jì)算】科研服務(wù)，測(cè)試狗團(tuán)隊(duì)開發(fā)，提供同步輻射、球差電鏡和常用材料表征，承接第一性原理計(jì)算、分子動(dòng)力學(xué)、有限元計(jì)算等

上一篇：研究生生活，就這？你以為 VS 實(shí)際上，太扎心了

下一篇：重磅發(fā)布：2021軟科世界大學(xué)學(xué)術(shù)排名，快來看看你的學(xué)校排第幾？

相關(guān)文章

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

2019-10-25

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

2021-01-22

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

2021-06-19

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

2019-10-25

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

2020-05-03

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

2020-08-24

舉報(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

上傳附件：