| 查看: 4810 | 回复: 1 | |||
[求助]
QE计算声子谱电声耦合系数和临界温度
|
|
第一步，scf自洽pw.x 1 &control 2 calculation='scf' 3 restart_mode='from_scratch', 4 !pseudo_dir='directory where pseudopotentials are stored/', 5 !outdir='directory where large files are written/' 6 pseudo_dir='./' 7 outdir='./tmp' 8 prefix='PHMgB2' 9 tstress=.true. 10 tprnfor=.true. 11 / 12 &system 13 ibrav=0, celldm(1)=1.89, nat=3, ntyp=2, ecutwfc=60,ecutrho=480, 14 occupations='smearing', smearing='gaussian', degauss=0.01, 15 la2F=.true. 16 / 17 &electrons 18 diagonalization='david' 19 conv_thr=1.0d-8 20 mixing_beta=0.7 21 / 22 CELL_PARAMETERS 23 3.0699288366337654 -0.0000000004270052 0.0000000000036074 24 1.5349644186324731 2.6586363584606119 0.0000000000024496 25 -0.0000000000056698 -0.0000000000017806 3.5174059808366489 26 ATOMIC_SPECIES 27 Mg 24.3050 Mg.pz-n-vbc.UPF 28 B 10.811 B.pz-bhs.UPF 29 ATOMIC_POSITIONS {angstrom} 30 Mg -0.0000000050616800 0.0000000000000088 0.0000000000000476 31 B 0.3333333201942207 0.6666666660819658 0.5000000000000511 32 B 0.6666666548674568 0.3333333339180254 0.4999999999999015 33 K_POINTS {automatic} 34 16 16 16 0 0 0 第二步，疏的k点自洽pw.x 1 &control 2 calculation='scf' 3 restart_mode='from_scratch', 4 !pseudo_dir='directory where pseudopotentials are stored/', 5 !outdir='directory where large files are written/' 6 pseudo_dir='./' 7 outdir='./tmp' 8 prefix='PHMgB2' 9 tstress=.true. 10 tprnfor=.true. 11 / 12 &system 13 ibrav=0, celldm(1)=1.89, nat=3, ntyp=2, ecutwfc=60,ecutrho=480, 14 occupations='smearing', smearing='gaussian', degauss=0.01, 15 !la2F=.true. 16 / 17 &electrons 18 diagonalization='david' 19 conv_thr=1.0d-8 20 mixing_beta=0.7 21 / 22 CELL_PARAMETERS 23 3.0699288366337654 -0.0000000004270052 0.0000000000036074 24 1.5349644186324731 2.6586363584606119 0.0000000000024496 25 -0.0000000000056698 -0.0000000000017806 3.5174059808366489 26 ATOMIC_SPECIES 27 Mg 24.3050 Mg.pz-n-vbc.UPF 28 B 10.811 B.pz-bhs.UPF 29 ATOMIC_POSITIONS {angstrom} 30 Mg -0.0000000050616800 0.0000000000000088 0.0000000000000476 31 B 0.3333333201942207 0.6666666660819658 0.5000000000000511 32 B 0.6666666548674568 0.3333333339180254 0.4999999999999015 33 K_POINTS {automatic} 34 8 8 8 0 0 0 第三步，ph.x 1 Phonons for MgB2 2 &inputph 3 tr2_ph=1.0d-14, 4 amass(1)=24.3050, amass(2)=10.811, 5 prefix='PHMgB2', 6 ldisp=.true., trans=.true. 7 !outdir='directory where large files are written/' 8 outdir='./tmp', 9 !elph=.true. 10 electron_phonon='interpolated' 11 fildyn='MgB2.dyn', 12 nq1=4, nq2=4, nq3=4 13 / 第四步，matydn.x 1 &input 2 asr='simple', amass(1)=24.3050, amass(2)=10.811, 3 flfrc='PHMgB2.fc', flfrq='PHMgB2.freq', la2F=.true., dos=.false. 4 / 5 100 6 0.0000 0.0000 0.0000 0 7 0.0000 0.0263 0.0000 0 8 0.0000 0.0526 0.0000 0 9 0.0000 0.0789 0.0000 0 10 0.0000 0.1053 0.0000 0 .................... （1）这一步生成的.freq文件用plotband.x处理的时候出现错误，如下： /vol-th/software/espresso-5.0.2/bin/plotband.x PHMgB2.freq Reading 9 bands at 100 k-points Error reading k-point # 1 （2）另外，将生成的.freq.gp文件导入origin作声子谱画出来的图像也不对。 请大神解答是哪里出错了呢 第五步，matydn.x计算lamada 1 &input 2 asr='simple', amass(1)=24.3050, amass(2)=10.811, 3 flfrc='PHMgB2.fc', flfrq='PHMgB2lamada.freq', 4 la2F=.true., dos=.true., 5 fldos='PHMgB2.dos', 6 nk1=10, nk2=10,nk3=10, 7 ndos=50 8 / 第六步，据说是使用lamada.in文件可以直接计算出Tc，论坛里有个例子是： 10 0.12 1 ! emax (something more than highest 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 Coloumb coefficient in the modified ! Allen-Dynes formula for T_c (via \omega_log) 我想问的是： （3）emax 是从第四步的freq文件中读出来的吗？ 具体怎么读呢，感觉看不太懂呢，degauss, smearing method 的取值依据是什么呢？ （4）第三到十八行是只写高对称带点不用像第四步那样都写上吗？第十一到十八行和我们通过前面步骤计算得到的哪些elph.文件是一个意思吗？ 谢谢大家，大家帮忙看一下。 |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有25人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有3人回复
-
为什么网上很多人说本周 12号出结果
已经有6人回复
-
听说今天filecode变了
已经有49人回复
-
帮忙看看fileCode
已经有12人回复
-
FileCode能看出啥？
已经有17人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
filecode
已经有9人回复
XIHA女王1990
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 407.4
- 散金: 50
- 红花: 1
- 帖子: 165
- 在线: 113.9小时
- 虫号: 1372814
- 注册: 2011-08-19
- 专业: 色谱分析
2楼2016-12-25 18:48:54
