应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 请问， 如何得到a2f(w) 图？
101/1返回列表
查看: 2105  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gloomy2004

金虫 (著名写手)

[求助] 请问， 如何得到a2f(w) 图？

请问从pwscf计算电声相互作用过程中， 应该有a2F(w)-w关系得到， 在哪个文件里面呢？我查看了几个文件， 格式都不对啊， 不是二维数据， 需要数据后续处理才能得到这个关系图吗？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-02-05 17:14:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
ben_ladeng(金币+2, 专家考核): 感谢 2012-02-06 00:26:37
没关心过超导方面，不过从PWscf的手册来看，PWscf可以做的工作参见example 07。
得到的结果是lamda和e-ph coupling系数。
然后，根据Eliashberg公式：

求得不同声子频率处的a2F。
具体的算法偶没详细看～关于超导的计算和理论，可以参考书：Theory of Superconductivity-From weak to strong coupling, A.S.Alexandrov[著], Institute of Physics Publishing 2003年出版。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
2楼2012-02-05 22:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

gloomy2004(金币+10): ★★★很有帮助 我也认为中间过程可以产生用来做这个图的文件，就是没找到是哪个文件？谢谢您。 2012-02-05 22:10:47
gloomy2004(金币+10): ★★★★★最佳答案 看来源代码， 我得学习看了。 2012-02-05 22:22:23
ben_ladeng(专家考核): 2012-02-06 00:26:47
不过看PWscf的代码，似乎lamda.f90中写着：
CODE:
  ! INPUT from standard input:
  !    emax  degaussq  ngaussq
  !    nks
  !    q(1,1)    q(2,1)    q(3,1)    wk(1)
  !      ...       ...       ...      ...
  !    q(1,nks)  q(2,nks)  q(3,nks)  wk(nks)
  !    filelph(1)
  !     ...
  !    filelph(nks)
  !
  ! emax (THz)    : alpha2F is plotted from 0 to "emax" in "nex" steps
  ! degaussq (THz): gaussian smearing for sum over q
  !                 NB: not the same used in phonon !
  ! ngaussq       : 0 for simple gaussian, 1 for Methfessel-Paxton etc.
  ! nks           : number of q-points used in the sum
  ! q, wk         : q-points and weights
  ! filelph       : output files from phonon, one for each q-point
  !                 May contain "nsig" calculations done with different
  !                 broadenings for the sum over k - all of them are used
  !
  ! OUTPUT in xmgr-readable format: files 'lambda.dat' and 'alpha2F.dat'
  !

没准可以直接做出Eliashberg方程的曲线～
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
3楼2012-02-05 22:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
ben_ladeng(专家考核): 2012-02-06 00:27:00
liliangfang(金币+1): 多谢提示 2012-02-06 09:05:33
PH中，这个不算看源代码，只是代码注释中提到了输入参数而已～这是再PWscf手册中提到的，没有认真看PWscf那只有几页的手册呢吧～
赞一下(2人) 回复此楼
思想重于技巧，内涵重于表象
4楼2012-02-05 23:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
3楼: Originally posted by souledge at 2012-02-05 22:04:47:
不过看PWscf的代码，似乎lamda.f90中写着：
[code]
  ! INPUT from standard input:
  !    emax  degaussq  ngaussq
  !    nks
  !    q(1,1)    q(2,1)    q(3,1)    wk(1)
  !      ...      ...

刚才找了老长时间，就是没找到lambda.f90文件。然后我的输出结果里也真心没有Eliashberg计算结果'alpha2F.dat'文件。这是怎么回事呢？
赞一下 回复此楼
5楼2012-02-20 10:17:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
另外，我的PWSCF是v4.1版本，会不会是版本问题？
赞一下 回复此楼
6楼2012-02-20 10:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
franch(金币+2): 鼓励交流。。呵呵 2012-02-20 16:02:59
问题解决了，原来matdyn.x有集成lambda.x的功能。相对于每个高斯展宽的a2f(w)都算出来了。
赞一下 回复此楼
7楼2012-02-20 10:30:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyongkai

木虫 (职业作家)
引用回帖:
7楼: Originally posted by gemucai at 2012-02-20 10:30:23
问题解决了，原来matdyn.x有集成lambda.x的功能。相对于每个高斯展宽的a2f(w)都算出来了。

您好，我利用kpoints.x产生的lambda.in的输入文件与ph.x计算得到的输出文件不相匹配是怎么回事呢？
赞一下 回复此楼
凌云
8楼2013-04-02 21:53:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
8楼: Originally posted by weiyongkai at 2013-04-02 21:53:46
您好，我利用kpoints.x产生的lambda.in的输入文件与ph.x计算得到的输出文件不相匹配是怎么回事呢？...

ph.x的输出文件不能直接用到lambda.in的，所以肯定不匹配
赞一下 回复此楼
9楼2013-04-03 10:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyongkai

木虫 (职业作家)
引用回帖:
5楼: Originally posted by gemucai at 2012-02-20 10:17:06
刚才找了老长时间，就是没找到lambda.f90文件。然后我的输出结果里也真心没有Eliashberg计算结果'alpha2F.dat'文件。这是怎么回事呢？...

QE5.0中有的，大侠
赞一下 回复此楼
凌云
10楼2013-06-21 14:04:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gloomy2004 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械学硕总分317求调剂！！！！ +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 各位老师您好：本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 求调剂 一志愿 本科 北科大 化学 343 +4 13831862839 2026-03-24 5/250 2026-03-25 09:47 by 无际的草原
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境） +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 材料调剂 +6 匹克i 2026-03-23 6/300 2026-03-24 21:09 by greychen00
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 一志愿中南化学（0703）总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
信息提示
关闭
请填处理意见
关闭