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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 请问， 如何得到a2f(w) 图？
101/1返回列表
查看: 2155  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿沪9，生物学326求调剂 +8 刘墨墨 2026-04-15 8/400 2026-04-16 10:29 by 圆心602
[考研] 0854求调剂 +17 门路摸摸 2026-04-15 19/950 2026-04-16 10:24 by dick_runner
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 273求调剂 +6 白居不易. 2026-04-09 8/400 2026-04-15 22:02 by wooluyong
[考研] 296求调剂 +12 汪！？！ 2026-04-09 13/650 2026-04-15 20:01 by 学员JpLReM
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成！
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
信息提示
关闭
请填处理意见
关闭