| 查看: 3217 | 回復(fù): 17 | |||
identation金蟲 (著名寫手)
|
[交流]
【求助】關(guān)于QHA模塊的使用 已有4人參與
|
|
以AlAs為例,想請教一下具體計算不同溫度的熱力學(xué)性質(zhì)怎么實現(xiàn)。 以下是espresso-4.2中給出的用QHA計算alas例子中的Edit_Me文件 ############################################################################## # Optional parameters, any information specific for the system studied # SysInfo='AlAs' # Mandatory parameters # Specify SystemName and Force Constants matrix Sysname='AlAs' FC_file='alas444.fc' # # Specify lattice type (used to create ttrinp file). It should be the same as in scf.in file # Specify atoms in the unit cell as they specified in scf.in file # Specify atomic masses for these atoms in the same order as in scf.in # Specify the frequency step (delta_e) as well, but 0.75 is a good choice ibrav=2 atoms="Al As " mass="26.98 74.922 " delta_e=0.75 # Edit ONLY amass parameters # Please do not change flfrq='frequency' line # leave asr (acoustic sum rule) and flfrc lines cat >matdyn.init < amass(1)=26.98, amass(2)=74.922, asr='crystal', flfrc=$FC_file, flfrq='frequency' &end EOF # # In most cases there is no need to edit files listed below, but if you like ... # # Temperature range for thermodynamic properties # T_start, T_end, T_step for QHA calculations cat > Temperature < EOF # Debye Temperature calculations # Phonon DOS filename (total phonon DOS, not projected), leave it as PHDOS.out # accuracy (limited 1.d-5, more accuracy is not required ) # Low_Temp_start, Low_Temp_end, and Low_Temp_step for Low Temperature limit, up to 15-30K # Hihg temperature and T_step for HT limit cat >T_Debye.in < 0.0001 3 15 3 500 10 EOF 具體實施的時候,我的理解是這樣的,不知道對不對,想請高人指教一下,感激不盡! (1)先構(gòu)建alas.scf.in文件,進(jìn)行pw.x的自洽計算; (2)構(gòu)建alas.ph.in,根據(jù)(1)的相關(guān)輸出進(jìn)行ph.x計算,得到力常數(shù)和動力學(xué)矩陣; (3)用QHA計算,但這里QHA計算的in文件的寫法是怎樣的呢?我的理解是如下的,但不知對否(里面的EOF和cat > Temperature < alas.qha.in文件 SysInfo='AlAs' Sysname='AlAs' FC_file='alas444.fc' ibrav=2 atoms="Al As " mass="26.98 74.922 " delta_e=0.75 cat >matdyn.init < amass(1)=26.98, amass(2)=74.922, asr='crystal', flfrc=$FC_file, flfrq='frequency' &end EOF cat > Temperature < EOF cat >T_Debye.in < 0.0001 3 15 3 500 10 EOF |
鐵桿木蟲 (著名寫手)

鐵桿木蟲 (著名寫手)

金蟲 (著名寫手)
鐵桿木蟲 (著名寫手)

金蟲 (著名寫手)
金蟲 (著名寫手)
鐵桿木蟲 (著名寫手)

金蟲 (著名寫手)
金蟲 (著名寫手)
|
圖2示出了XXXX直流電導(dǎo)率與溫度倒數(shù)(180─300K)的關(guān)系點。XXXX在180K的直流電導(dǎo)率在10-4S/m左右,在該溫度XXXX呈現(xiàn)電阻特性。XXXX的電導(dǎo)率隨溫度上升而提高,并在室溫時達(dá)到10-2S/m。計算得到的激活能為Ea=0.237eV(±0.001 eV)(用Arrhenius關(guān)系σ=σ0exp(-Ea/kBT)擬合得到,σ0是指前因子、Ea是激活能、kB是Boltzmann常數(shù)、T是溫度),這是由自由電荷載流子(如電子、特定種類的缺陷、自由離子)的長程遷移所貢獻(xiàn)的。如此低的激活能(<0.3 eV)通常表明,一般與鈣鈦礦型化合物高溫?zé)Y(jié)有關(guān)的氧空位不參與長程電導(dǎo)。帶電載流子,特別是局域于Mn4+離子的兩個自由電子,其長程跳躍可能對STMO的電導(dǎo)率產(chǎn)生顯著的貢獻(xiàn)。 請教如下問題 1、這里的激活能Ea=0.237eV(±0.001 eV)是帶隙么? 2、PWSCF的GGA+U計算只指定ldau=.true.和Hubbard_U(i)是不是就足夠了?至于starting_ns_eigenvalue(m,ispin,I)需要指定么? [ Last edited by identation on 2010-6-1 at 08:45 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|