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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 请问， 如何得到a2f(w) 图？
101/1返回列表
查看: 2104  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +8 顾九笙要谦虚 2026-03-24 12/600 2026-03-25 20:15 by 20021109
[考研] 26考研-291分-厦门大学（085601）-柔性电子学院材料工程专业求调剂 +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[考研] 考研一志愿苏州大学初始315（英一）求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 考研调剂 +5 呼呼？~+123456 2026-03-24 5/250 2026-03-25 18:15 by xcjcqu
[考研] 0854电子信息求调剂 +7 α____ 2026-03-22 9/450 2026-03-25 13:37 by α____
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 化学调剂 +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 0854 考研调剂 招生了！AI 方向 +5 pk3725069 2026-03-19 17/850 2026-03-24 17:30 by zhouxuan..
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 263求调剂 +6 yqdszhdap－ 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap－
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 22 350 本科985求调剂，求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
信息提示
关闭
请填处理意见
关闭