24小时热门版块排行榜    

查看: 3178  |  回复: 10

Local_Spin

铜虫 (小有名气)

[求助] 电声耦合计算中lambd.x的一些问题

各位虫子大家好,我刚刚接触pwscf,用来计算电声耦合。我有几个关于lambda.x的问题想请问一下大家:

1、在lambda.x的input文件的第一行,需要输入3个参量,分别是emax,degauss,以及smearing method。请问这三个参数有什么含义,取法有什么讲究,对于计算的结果有什么影响?

2、在lambda.x的output文件里面,分别列出了degauss从0.005到0.05的lambda及Tc的计算结果,不同的degauss对应的结果不同,那么,到底哪一个结果才是我们想要的最终结果呢?


谢谢大家!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gemucai

木虫 (正式写手)

★ ★ ★
sunyang1988: 金币+3, 谢谢指点 2012-11-22 10:44:06
电声耦合啊,苦逼的孩纸。。。
emax是比最高声子能量还高一点的一个值,取多少自己定
degauss和smearing method是在布里渊区内取点的时候用的展宽方式和大小,但具体该怎么取我还不懂
因为你布里渊区取点,k点数目是一定的,那么degauss越大越容易收敛,所以你就看degauss比较大的几个lambda和Tc,是否随着degauss的增加有收敛的趋势。如果有这种趋势,那么恭喜你。

下面是例子及其讲解:
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)

» 本帖已获得的红花(最新10朵)

2楼2012-11-21 16:33:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Local_Spin

铜虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by gemucai at 2012-11-21 16:33:19
电声耦合啊,苦逼的孩纸。。。
emax是比最高声子能量还高一点的一个值,取多少自己定
degauss和smearing method是在布里渊区内取点的时候用的展宽方式和大小,但具体该怎么取我还不懂
因为你布里渊区取点,k点数 ...

有一种茅塞顿开的感觉。然后坐等楼下高手继续介绍初始degauss的取法
3楼2012-11-21 16:57:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
2楼: Originally posted by gemucai at 2012-11-21 16:33:19
电声耦合啊,苦逼的孩纸。。。
emax是比最高声子能量还高一点的一个值,取多少自己定
degauss和smearing method是在布里渊区内取点的时候用的展宽方式和大小,但具体该怎么取我还不懂
因为你布里渊区取点,k点数 ...

您好,最近老板又想算电声子了,他自己也不是很懂,我只能自己学了,对于example里的这个lambda。in文件,我也是很头疼啊,不知道第一行的emax,degauss,以及smearing method这几个值,您现在知道怎么确定了没有,我连第一个emax都不知道去哪个文件里看声子的最高能量,还有中间的8个q点,这不就是下面那8个么,但是为什么又不一样呢,还是无所谓的,直接写下面那8个就可以了,谢谢您指导了!
4楼2013-07-03 11:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-07-08 22:37:53
引用回帖:
4楼: Originally posted by 刘仕晨 at 2013-07-03 11:31:15
您好,最近老板又想算电声子了,他自己也不是很懂,我只能自己学了,对于example里的这个lambda。in文件,我也是很头疼啊,不知道第一行的emax,degauss,以及smearing method这几个值,您现在知道怎么确定了没有, ...

是这样的,中间是k点,下面那8个是输入的文件名称。声子能量,直接看最高声子频率就是,要取以THz为单位的。
5楼2013-07-03 17:44:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
5楼: Originally posted by gemucai at 2013-07-03 17:44:01
是这样的,中间是k点,下面那8个是输入的文件名称。声子能量,直接看最高声子频率就是,要取以THz为单位的。...

这个我知道中间的是K点,下面的是前面计算的文件名称,但是我的疑惑是,就是就是取得q点,为什么中间的K点和下面的八个文件名上面的K点是不一样的,不知道您有没有注意,比如中间第5个K点是 -0.2500000 -0.2500000  0.7500000 而下面第五个文件名后面的K点是elph. 0.750000.-0.250000. 0.750000,还有后面的degauss, smearing method 以及最后一行的0.01,请问您知道这个该怎么取么,谢谢!
6楼2013-07-03 18:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-07-08 22:37:32
引用回帖:
6楼: Originally posted by 刘仕晨 at 2013-07-03 18:39:46
这个我知道中间的是K点,下面的是前面计算的文件名称,但是我的疑惑是,就是就是取得q点,为什么中间的K点和下面的八个文件名上面的K点是不一样的,不知道您有没有注意,比如中间第5个K点是 -0.2500000 -0.2500000 ...

degauss 和 smearing方法的取值我确实不知,最后一行取0.1到0.15之间,是代表电子之间的库仑排斥作用。还有就是在周期性条件下,0.75跟-0.25是一样的。
7楼2013-07-04 18:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
7楼: Originally posted by gemucai at 2013-07-04 18:11:30
degauss 和 smearing方法的取值我确实不知,最后一行取0.1到0.15之间,是代表电子之间的库仑排斥作用。还有就是在周期性条件下,0.75跟-0.25是一样的。...

谢谢,非常感谢您每次的帮助,已经很感谢了,那我去发个新帖子,看看有没有其他人能知道这方面的,有了答案,第一时间通知您,也方便以后学习的人学习交流,还是很谢谢您!
8楼2013-07-04 18:43:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

引用回帖:
7楼: Originally posted by gemucai at 2013-07-04 18:11:30
degauss 和 smearing方法的取值我确实不知,最后一行取0.1到0.15之间,是代表电子之间的库仑排斥作用。还有就是在周期性条件下,0.75跟-0.25是一样的。...

您好,我想问下
1. 您现在指导smearing方法怎么取了吗?
2. 还有中间的K点是不是就按照生成的elph. 这些文件后面跟的那些作为K点依次排列就可以了?还有那个K点的权重是根据什么得出的?
3. “0.10       ! \mu the Coloumb coefficient in the modified”  这一行是代表电子之间的库伦排斥作用,取0.1到0.15之间,这个取值有什么原则吗?还是随便取啊?
4. “           ! Allen-Dynes formula for T_c (via \omega_log)”这一行指的是什么?为什么是空的呢?
9楼2016-05-06 21:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

楼主现在知道第二个问题怎么解决了吗
10楼2016-05-09 16:00:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Local_Spin 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 317求调剂 +8 申子申申 2026-03-19 13/650 2026-03-21 00:09 by 刘国森
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[考研] 一志愿西安交通大学 学硕 354求调剂 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
信息提示
请填处理意见