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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 请问， 如何得到a2f(w) 图？
101/1返回列表
查看: 2334  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 22:18 by xxniao123
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭