24小时热门版块排行榜

>论坛更新日志 (3016)
>虫友互识 (732)
>文献求助 (167)
>休闲灌水 (134)
>导师招生 (133)
>考博 (124)
>教师之家 (80)
>找工作 (73)
>硕博家园 (72)
>博后之家 (71)
>论文投稿 (70)
>论文道贺祈福 (65)
>基金申请 (61)
>考研 (55)
>公派出国 (49)
>招聘信息布告栏 (31)

返回列表

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

[求助] 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.文件是一个意思吗？

谢谢大家，大家帮忙看一下。

回复此楼

» 猜你喜欢

1楼 2016-03-29 15:27:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

XIHA女王1990

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 407.4
散金: 50
红花: 1
帖子: 165
在线: 113.9小时
虫号: 1372814
注册: 2011-08-19
专业: 色谱分析

请问楼主临街温度算出来了吗

发自小木虫IOS客户端

赞一下

回复此楼

2楼2016-12-25 18:48:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主菲顶顶的主题更新

返回列表