版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

gloomy2004

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 1461.1
散金: 64
红花: 32
沙发: 1
帖子: 1177
在线: 167.5小时
虫号: 49848
注册: 2004-07-10

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

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

回复此楼

» 猜你喜欢

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

塔菲尔曲线如何得到？从CHI上的塔菲尔图测试直接得到的曲线准确吗？已经有20人回复
请问从这个红外谱图上可以得到什么信息，要进一步确认分子结构还需要做什么测试！已经有9人回复
响应面如何得到这种图啊已经有10人回复
求助castep计算中：晶胞总能量E 与晶胞体积V 的关系图如何得到，看内容已经有24人回复
动力学跑NPT时如何得到比体积与时间的关系图？已经有9人回复
请问如何得到矢量图？已经有11人回复
【求助】这样的散点图如何得到呢？已经有7人回复
【求助】如何得到bode图？已经有14人回复
【求助】请问红外光谱测得到的谱图怎么保存成word形式的啊！已经有10人回复
如何由PPT上的图片得到高分辨率图片已经有3人回复
如何从X射线粉末衍射图谱得到晶体的外形已经有11人回复

1楼 2012-02-05 17:14:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

★ ★
感谢参与，应助指数 +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

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

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

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

★
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

木虫 (正式写手)

应助: 16 (小学生)
金币: 1774.6
散金: 157
红花: 4
帖子: 693
在线: 201.7小时
虫号: 816478
注册: 2009-07-27
专业: 凝聚态物性 II ：电子结构

引用回帖:

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的回帖