版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

identation

金蟲 (著名寫手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】關(guān)于QHA模塊的使用已有4人參與

以AlAs為例，想請(qǐng)教一下具體計(jì)算不同溫度的熱力學(xué)性質(zhì)怎么實(shí)現(xiàn)。

以下是espresso-4.2中給出的用QHA計(jì)算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 < &input
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 < 5 500 5
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 < PHDOS.out
0.0001
3 15 3
500 10
EOF

具體實(shí)施的時(shí)候，我的理解是這樣的，不知道對(duì)不對(duì)，想請(qǐng)高人指教一下，感激不盡！
（1）先構(gòu)建alas.scf.in文件，進(jìn)行pw.x的自洽計(jì)算；
（2）構(gòu)建alas.ph.in，根據(jù)（1）的相關(guān)輸出進(jìn)行ph.x計(jì)算，得到力常數(shù)和動(dòng)力學(xué)矩陣；
（3）用QHA計(jì)算，但這里QHA計(jì)算的in文件的寫法是怎樣的呢？我的理解是如下的，但不知對(duì)否（里面的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 < &input
amass(1)=26.98,
amass(2)=74.922,
asr='crystal',
flfrc=$FC_file,
flfrq='frequency'
&end
EOF

cat > Temperature < 5 500 5
EOF

cat >T_Debye.in < PHDOS.out
0.0001
3 15 3
500 10
EOF

回復(fù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

COMSOL Multiphysics中文使用手冊(cè)和化學(xué)反應(yīng)工程模塊教程已經(jīng)有233人回復(fù)
關(guān)于用超胞方法計(jì)算聲子的一些經(jīng)驗(yàn) 已經(jīng)有21人回復(fù)
如何通過聲子計(jì)算求解熱導(dǎo)率已經(jīng)有5人回復(fù)
Castep模塊CPU和內(nèi)存使用率已經(jīng)有7人回復(fù)
【求助】siesta中transiesta模塊的收斂標(biāo)準(zhǔn)是啥~~ 已經(jīng)有20人回復(fù)
【求助】siesta中transiesta模塊的收斂問題已經(jīng)有13人回復(fù)

1樓 2010-05-23 15:49:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時(shí)
蟲號(hào): 932458
注冊(cè): 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

EOF (end of file) 作為標(biāo)識(shí)

cat >T_Debye.in < PHDOS.out
0.0001
3 15 3
500 10
EOF

是將EOF前面的內(nèi)容寫入 T_Debye.in

qha.x < in.qha
qha.x 會(huì)讀取in.qha 設(shè)定的參數(shù)和
matdyn.init Temperature T_Debye.in三個(gè)文件

[ Last edited by lono75 on 2010-5-23 at 20:59 ]

贊一下

回復(fù)此樓

由天命,成人事！

2樓2010-05-23 20:40:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時(shí)
蟲號(hào): 932458
注冊(cè): 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
wuli8(金幣+1):謝謝 2010-05-23 20:45:08
identation(金幣+3): 2010-05-24 10:16:57

EOF (end of file) 作為標(biāo)識(shí)

cat > Temperature < 5 500 5
EOF

是將EOF前面的內(nèi)容寫入Temperature

贊一下(1人)

回復(fù)此樓

由天命,成人事！

3樓2010-05-23 20:43:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by lono75 at 2010-05-23 20:43:24:
EOF (end of file) 作為標(biāo)識(shí)

cat > Temperature < 5 500 5
EOF

是將EOF前面的內(nèi)容寫入Temperature

那QHA.in文件的格式是怎樣的呢？

贊一下

回復(fù)此樓

4樓2010-05-24 10:17:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時(shí)
蟲號(hào): 932458
注冊(cè): 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

identation(金幣+5): 2010-05-24 22:56:14
identation(金幣+18): 2010-05-25 15:35:17

# Parameters required for QHA calculations
# Total Phonon DOS file
# output file for C_V, S, Internal energy
# T_start, T_end, T_step for QHA calculations
#
cat >fqha.in < PHDOS.out
$Sysname.QHA.out
5 500 5
EOF

贊一下

回復(fù)此樓

由天命,成人事！

5樓2010-05-24 13:20:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

goldenfisher

金蟲 (著名寫手)

1ST強(qiáng)帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時(shí)
蟲號(hào): 117764
注冊(cè): 2005-11-26
性別: GG
專業(yè): 原子和分子物理

identation(金幣+10): 2010-05-24 22:56:28

引用回帖:

Originally posted by identation at 2010-05-24 10:17:21:

那QHA.in文件的格式是怎樣的呢？

格式就是前面說說的那樣的，注意的是flfrc的文件名應(yīng)該是你前面ph.x算出來的文件名。
&input
amass(1)=26.98,
amass(2)=74.922,
asr='crystal',
flfrc=$FC_file,
flfrq='frequency'
&end

贊一下

回復(fù)此樓

6樓2010-05-24 22:31:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by goldenfisher at 2010-05-24 22:31:42:

格式就是前面說說的那樣的，注意的是flfrc的文件名應(yīng)該是你前面ph.x算出來的文件名。
&input
amass(1)=26.98,
amass(2)=74.922,
asr='crystal',
flfrc=$FC_file,
flfrq='frequ ...

用這個(gè)輸入文件在哪里指定溫度呢？

贊一下

回復(fù)此樓

7樓2010-05-24 22:56:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時(shí)
蟲號(hào): 932458
注冊(cè): 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
cenwanglai(金幣+1):謝謝參與! 2010-05-25 09:23:28
zzy870720z(金幣+1):謝謝 2010-05-25 14:29:00
identation(金幣+4): 2010-05-25 15:35:33

引用回帖:

Originally posted by identation at 2010-05-24 22:56:04:

用這個(gè)輸入文件在哪里指定溫度呢？

上面已經(jīng)清楚咯
# T_start, T_end, T_step for QHA calculations
5 500 5

贊一下(2人)

回復(fù)此樓

由天命,成人事！

8樓2010-05-25 07:01:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

goldenfisher

金蟲 (著名寫手)

1ST強(qiáng)帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時(shí)
蟲號(hào): 117764
注冊(cè): 2005-11-26
性別: GG
專業(yè): 原子和分子物理

★
zzy870720z(金幣+1):謝謝交流 2010-05-25 14:29:15
identation(金幣+10): 2010-05-25 15:35:39
identation(金幣+5):又有問題請(qǐng)教，謝謝 2010-06-01 08:42:44

是的，樓上的是對(duì)的，我前面寫錯(cuò)了，那個(gè)應(yīng)該是聲子的計(jì)算文件。
PHDOS.out
0.0001
3 15 3
500 10
這樣寫就可以了。你也可以直接運(yùn)行qha.x，然后按照提示一步一步來。

贊一下(1人)

回復(fù)此樓

9樓2010-05-25 12:56:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

又有問題了

圖2示出了XXXX直流電導(dǎo)率與溫度倒數(shù)（180─300K）的關(guān)系點(diǎn)。XXXX在180K的直流電導(dǎo)率在10-4S/m左右，在該溫度XXXX呈現(xiàn)電阻特性。XXXX的電導(dǎo)率隨溫度上升而提高，并在室溫時(shí)達(dá)到10-2S/m。計(jì)算得到的激活能為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+離子的兩個(gè)自由電子，其長程跳躍可能對(duì)STMO的電導(dǎo)率產(chǎn)生顯著的貢獻(xiàn)。

請(qǐng)教如下問題
1、這里的激活能Ea=0.237eV（±0.001 eV）是帶隙么？
2、PWSCF的GGA+U計(jì)算只指定ldau=.true.和Hubbard_U(i)是不是就足夠了？至于starting_ns_eigenvalue(m,ispin,I)需要指定么？

[ Last edited by identation on 2010-6-1 at 08:45 ]

贊一下

回復(fù)此樓

10樓2010-06-01 08:42:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 identation 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求材料，環(huán)境專業(yè)調(diào)劑 +3	18567500178 2026-03-18	3/150	2026-03-23 23:50 by 熱情沙漠
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +3	17501029541 2026-03-23	4/200	2026-03-23 23:47 by Txy@872106
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 材料292調(diào)劑 +4	橘頌思美人 2026-03-23	4/200	2026-03-23 23:16 by peike
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 一志愿中國石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +4	石能偉 2026-03-17	4/200	2026-03-23 17:51 by 17862566385
[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	3/150	2026-03-23 10:32 by Iveryant
[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx

24小時(shí)熱門版塊排行榜

[交流] 【求助】關(guān)于QHA模塊的使用 已有4人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

又有問題了

[交流] 【求助】關(guān)于QHA模塊的使用已有4人參與