24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

172292350

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 26
帖子: 14
在线: 75.1小时
虫号: 2660004
注册: 2013-09-16
专业: 凝聚态物性 II ：电子结构

[交流] 求助QE计算电声耦合及Tc的问题已有3人参与

请问在第四步matydn.x 与第六步lamada.x中输入文件中的K-POINTS 与q-points怎么获得，对一般体系的话。有辅助工具么? 非常感谢~

回复此楼

» 收录本帖的淘帖专辑推荐

第一性

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有8人回复
求推荐博导已经有4人回复
为什么nbs上溴没有产物点出现呢已经有5人回复
存款400万可以在学校里躺平吗已经有34人回复
求助一下有机合成大神已经有4人回复
最失望的一年已经有16人回复
求推荐英文EI期刊已经有5人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
疑惑？已经有5人回复

1楼 2016-06-19 10:21:53

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

wlandy

铁杆木虫 (正式写手)

知道自己不知道

应助: 0 (幼儿园)
金币: 7953.4
红花: 3
帖子: 886
在线: 83.9小时
虫号: 2923551
注册: 2014-01-09
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

首先说明：q点和k点的选取方式是一样的。
方法如下：
step1，可以在qe中设置对应的q点。例如 Al 3 3 3 ,计算模式设置为 relax
step2，pw.x 运行
step3，在out文件中找到如下信息：
   number of k points=    4  Marzari-Vanderbilt smearing, width (Ry)=  0.0500
                     cart. coord. in units 2pi/alat
      k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0740741
      k( 2) = (  -0.3333333 0.3333333  -0.3333333), wk = 0.5925926
      k( 3) = ( 0.0000000 0.6666667 0.0000000), wk = 0.4444444
      k( 4) = ( 0.6666667  -0.0000000 0.6666667), wk = 0.8888889
step4，gamma点是1，k2 就是 8， k3 就是 6， k4就是 12
step5，对照example里的Al lambda.in
10  0.12  1 ! emax (something more than highest phonon mode in THz), degauss, smearing method
4       ! Number of q-points for which EPC is calculated,
0.000000000 0.000000000 0.000000000 1.00  ! the first q-point, use kpoints.x program to calculate
-0.333333333 0.333333333  -0.333333333 8.0 ! q-points and their weight
0.000000000 0.666666667 0.000000000 6.0
0.666666667  -0.000000000 0.666666667  12.0
elph_dir/elph.inp_lambda.1 ! elph output file names,
elph_dir/elph.inp_lambda.2 ! in the same order as the q-points before
elph_dir/elph.inp_lambda.3
elph_dir/elph.inp_lambda.4
0.10                   ! \mu the Coloumb coefficient in the modified
                     ! Allen-Dynes formula for T_c (via \omega_log)