 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 3722 | 回復(fù): 2 | ||
xiaoqiu007新蟲 (小有名氣)
|
[求助]
5.0.2版PWSCF電聲耦合常數(shù)的計(jì)算過程解讀
|
|
跟上木蟲上發(fā)的這個(gè)熱帖：http://www.gaoyang168.com/bbs/viewthread.php?tid=6086352&fpage=1，也請教幾個(gè)關(guān)于5.0.2版PWSCF電聲耦合常數(shù)的計(jì)算問題，希望能得到廣大同仁的指點(diǎn)，我也會(huì)即時(shí)匯報(bào)我的計(jì)算進(jìn)展。 5.0.2版計(jì)算電聲耦合常數(shù)分為以下幾個(gè)步驟： 1. make a self-consistent calculation for Al using a dense grid of k-points. The dense grid must contain all k and k+q grid points used in the subsequent electron-phonon calculation and must be dense enough to produce accurate el-phon coefficients (in particular the double-delta integral at Ef is very critical). 問題1：第一步就涉及q點(diǎn)和k點(diǎn)設(shè)置問題，例子中是設(shè)置成4 ：2 ：1的關(guān)系，即：16x16x16, 8x8x8, 4x4x4。請問q點(diǎn)的設(shè)置依據(jù)是什么？一般應(yīng)設(shè)置成多少合適？ 什么是 the double-delta integral at Ef ？ 2. make aself-consistent calculation for Al using a grid of k-points that is suitable for good self-consistency and phonon calculation.This example uses a (888) MP grid. 3. make the phonon and electron-phonon calculation for the grid of q-points. Specify elph=.true., and the name of a file where the derivative of the potential is stored "fildvscf".This example uses a (444) MP grid of q-points (nq1=4, nq2=4, nq3=4). The output contains the results for the el-ph coefficient at each q-point lambda(q), gamma(q), and the double-delta integral at several values of the gaussian broadening (set in filePH/elphon.f90). These are useful for convergence testing. (input=al.elph.in,output=al.elph.out) The results are also written into outputfiles "a2Fq2r.*", one per value of the gaussian broadening, for further use. 問題2：第三步中（1）5.0.2新版在ph.x步驟設(shè)置 electron_phonon='interpolated', 以前的老版本是設(shè)成：elph=.true.,請問另個(gè)參數(shù)的區(qū)別在哪里？ （2）convergence testing是算完之后測試？ gaussian broadening 對計(jì)算結(jié)果有什么影響？ 4. Bring to r-space both force constants and el-phon coefficients using "q2r.x" Output in files "a2Fmatdyn.*", one per value of the gaussian broadening 問題3：第四步q2r.x，為什么要轉(zhuǎn)換到實(shí)空間？ 與gaussian broadening 有什么關(guān)系？這一步是必須的嗎？ 5. Calculate gamma on selected lines using "matdyn.x" (dos=.false.) 問題4： 這一步是在計(jì)算聲子色散曲線，它是計(jì)算電聲耦合常數(shù)必須的步驟嗎？ 6. Calculate lambda coefficient (in file"lambda"  and the a^2F(omega)function using "matdyn.x" (dos=.true.) 問題5： 這一步感覺是在計(jì)算聲子態(tài)密度，請問該步能同時(shí)給出 lambda coefficient和a^2F(omega) function 嗎？ 7. Calculate lambda coefficient (in "lambda.out" and Tc using "lambda.x"問題6： 計(jì)算仍在進(jìn)行中，請問有過計(jì)算經(jīng)歷的同仁，新版通過這一步可以直接算出Tc嗎？ 綜上，感覺第5步，第6步似乎不是計(jì)算電聲耦合常數(shù)必須的步驟，不知理解是否正確。 完整的例子文件如下所示： # #########SCF at dense k-mesh, good enough forelectronic DOS 【多密算密】 ####### # cat >al.scf.fit.in << EOF &control calculation='scf' restart_mode='from_scratch', prefix='al', pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 2, celldm(1) =7.5, nat= 1, ntyp= 1, ecutwfc =15.0, occupations='smearing',smearing='methfessel-paxton', degauss=0.05, la2F= .true., / &electrons conv_thr = 1.0d-8 mixing_beta = 0.7 / ATOMIC_SPECIES Al 26.98 Al.pz-vbc.UPF ATOMIC_POSITIONS Al 0.00 0.00 0.00 K_POINTS{automatic} 16 16 16 0 0 0 EOF $ECHO " running the scf calculation with densek-point grid...\c" $PW_COMMAND < al.scf.fit.in > al.scf.fit.out # ########## SCFat k-mesh good enough for phonons【多好算好！】###### # cat > al.scf.in<< EOF &control calculation='scf' restart_mode='from_scratch', prefix='al', pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 2, celldm(1) =7.5, nat= 1, ntyp= 1, ecutwfc =15.0, occupations='smearing',smearing='methfessel-paxton', degauss=0.05 / &electrons conv_thr = 1.0d-8 mixing_beta = 0.7 / ATOMIC_SPECIES Al 26.98 Al.pz-vbc.UPF #### Pseudopotential type: NORMCONS ATOMIC_POSITIONS Al 0.00 0.00 0.00 K_POINTS{automatic} 8 8 8 00 0 EOF $ECHO " running the scf calculation...\c" $PW_COMMAND <al.scf.in > al.scf.out ############################# cat >al.elph.in << EOF Electron-phonon coefficients for Al &inputph tr2_ph=1.0d-10, prefix='al', fildvscf='aldv', amass(1)=26.98, outdir='$TMP_DIR/', fildyn='al.dyn', electron_phonon='interpolated', #與elph=.true.,的區(qū)別何在？ trans=.true., ldisp=.true. nq1=4, nq2=4, nq3=4 / EOF $ECHO " running the el-ph calculation...\c" $PH_COMMAND <al.elph.in > al.elph.out # # q2r and matdyn # cat > q2r.in<< EOF &input zasr='simple', fildyn='al.dyn', flfrc='Al444.fc', la2F=.true. / EOF $ECHO " running q2r...\c" $Q2R_COMMAND < q2r.in > q2r.out # ###########聲子譜：聲子色散曲線#### cat >matdyn.in.freq << EOF &input asr='simple', amass(1)=26.98, flfrc='Al444.fc', flfrq='Al444.freq', la2F=.true., dos=.false. / 19 0.000 0.0 0.0 0.0 0.125 0.0 0.0 0.0 0.250 0.0 0.0 0.0 0.375 0.0 0.0 0.0 0.500 0.0 0.0 0.0 0.750 0.0 0.0 0.0 1.000 0.0 0.0 0.0 0.825 0.125 0.125 0.0 0.750 0.250 0.250 0.0 0.625 0.375 0.375 0.0 0.500 0.500 0.500 0.0 0.325 0.325 0.325 0.0 0.250 0.250 0.250 0.0 0.125 0.125 0.125 0.0 0.000 0.000 0.000 0.0 0.125 0.125 0.000 0.0 0.250 0.250 0.000 0.0 0.325 0.325 0.000 0.0 0.500 0.500 0.000 0.0 EOF $ECHO " running matdyn for frequencycalculation...\c" $MATDYN_COMMAND< matdyn.in.freq > matdyn.out.freq check_failure $? ###########聲子態(tài)密度計(jì)算#### cat >matdyn.in.dos << EOF &input asr='simple', amass(1)=26.98, flfrc='Al444.fc', flfrq='Al444.freq', la2F=.true., dos=.true. fldos='phonon.dos', nk1=10, nk2=10, nk3=10,ndos=50 / EOF $ECHO " running matdyn for a2F(omega) calculation...\c" $MATDYN_COMMAND< matdyn.in.dos > matdyn.out.dos check_failure $? $ECHO " done" # # ############Eliashberg函數(shù)alpha^2 F和電聲耦合常數(shù)lamda計(jì)算######## cat > lambda.in<< EOF 10 0.12 1 ! emax (something more thanhighest phonon mode in THz), degauss, smearing method 8 ! Number of q-points for which EPC is calculated, 0.0000000 0.0000000 0.0000000 1.00 ! the first q-point, use kpoints.x program to calculate -0.2500000 -0.2500000 0.2500000 8.00 ! q-points and their weight -0.5000000 -0.5000000 0.5000000 4.00 ! 0.0000000 0.0000000 0.5000000 6.00 ! 4th q-point, qx,qy,qz -0.2500000 -0.2500000 0.7500000 24.00 ! -0.5000000 -0.5000000 1.0000000 12.00 ! 0.0000000 0.0000000 1.0000000 3.00 ! -0.5000000 0.0000000 1.0000000 6.00 ! the last q-point elph. 0.000000.0.000000. 0.000000 ! elph output file names, elph.-0.250000.0.250000.-0.250000 ! in the same order as the q-points before elph.0.500000.-0.500000. 0.500000 elph. 0.000000.0.500000. 0.000000 elph.0.750000.-0.250000. 0.750000 elph. 0.500000.0.000000. 0.500000 elph.0.000000.-1.000000. 0.000000 elph.-0.500000.-1.000000.0.000000 0.10 ! \mu the Coloumbcoefficient in the modified ! Allen-Dynes formulafor T_c (via \omega_log) EOF $ECHO " running lambda.x for lambdacalculation...\c" $LAMBDA_COMMAND< lambda.in > lambda.out |
回復(fù)此樓» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 學(xué)習(xí)PWSCF(一)：能帶計(jì)算步驟與討論
已經(jīng)有28人回復(fù)
- pwscf是如何計(jì)算debye溫度的
已經(jīng)有7人回復(fù)
- pwscf并行計(jì)算時(shí)總是隔一段時(shí)間重新開始
已經(jīng)有5人回復(fù)
- 怎么利用PWSCF計(jì)算差分電荷密度
已經(jīng)有25人回復(fù)
- 請教，關(guān)于通過PWSCF用wannier90計(jì)算介電極化和hopping integral的方法
已經(jīng)有8人回復(fù)
- PWscf計(jì)算homo lumo或者valence conduction band，問題已解決
已經(jīng)有14人回復(fù)
- 【求助】求詳細(xì)的PWscf并行編譯過程（越詳細(xì)越好，包括各種環(huán)境和參數(shù)的設(shè)置，及軟件?/a>
已經(jīng)有29人回復(fù)
- 【求助】關(guān)于PWSCF介電常數(shù)和聲子譜的計(jì)算
已經(jīng)有7人回復(fù)
- 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計(jì)算
已經(jīng)有23人回復(fù)
- 【求助】請教PWSCF晶體建模和自洽計(jì)算，不收斂
已經(jīng)有21人回復(fù)
xiaoqiu007
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 356.5
- 散金: 7
- 帖子: 100
- 在線: 20.8小時(shí)
- 蟲號: 707167
- 注冊: 2009-02-23
- 專業(yè): 金屬功能材料
2樓2013-09-19 09:54:24
nishihao
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1.5
- 帖子: 42
- 在線: 55.8小時(shí)
- 蟲號: 2080714
- 注冊: 2012-10-23
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
3樓2015-04-29 09:00:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiaoqiu007 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|