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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » Abinit 中计算 PDOS 和PJPDOS
232/3返回列表上一页123下一页
查看: 1735  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

zxzj05

荣誉版主 (著名写手)
★ ★
fegg7502(金币+2,VIP+0):鼓励讨论
的确是files问题，非常感谢watermall  的指点。
但是，我又碰到另外一个问题，原因还该不清楚，
log文件最后一句话是：
simpson_int : number of points in function must be >=6

out文件最后是：
Calculation of the interatomic forces
-begin at tcpu      0.180  and twall      0.180 sec
       Homogeneous q point set in the B.Z.
Grid q points  :        8
  1)   0.00000000E+00  0.00000000E+00  0.00000000E+00
  2)   5.00000000E-01  0.00000000E+00  0.00000000E+00
  3)   0.00000000E+00  5.00000000E-01  0.00000000E+00
  4)   5.00000000E-01  5.00000000E-01  0.00000000E+00
  5)   0.00000000E+00  0.00000000E+00  5.00000000E-01
  6)   5.00000000E-01  0.00000000E+00  5.00000000E-01
  7)   0.00000000E+00  5.00000000E-01  5.00000000E-01
  8)   5.00000000E-01  5.00000000E-01  5.00000000E-01
symkpt : not enough symmetry to change the number of k points.
symkpt : the number of k-points, thanks to the symmetries,
is reduced to 29791 .

换到其他结构体系，仍然出现类似问题：
60)  -2.50000000E-01  5.00000000E-01 -2.50000000E-01
61)   0.00000000E+00 -2.50000000E-01 -2.50000000E-01
62)   2.50000000E-01 -2.50000000E-01 -2.50000000E-01
63)   5.00000000E-01 -2.50000000E-01 -2.50000000E-01
64)  -2.50000000E-01 -2.50000000E-01 -2.50000000E-01
symkpt : the number of k-points, thanks to the symmetries,
is reduced to    27 .
symkpt : the number of k-points, thanks to the symmetries,
is reduced to  9261 .

期间，增大了ng2qpt  也没有效果。

文件中都没有error出现，还请watermall出手相助，帮忙解决这个问题？  
非常感谢！
赞一下(10人) 回复此楼
储氢家族欢迎储氢研究者！
11楼2008-12-12 19:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)
★ ★ ★
spur(金币+3,VIP+0):^_^ 3Q！！欢迎常来！
我觉得没什么问题啊
这段是abinit对BZ的K点根据对称性进行优化的操作
你给的第一个由于对称性低，所以没有什么优化“not enough symmetry to change the number of k points”
你给的第二个对称性比较好，所以“the number of k-points, thanks to the symmetries,is reduced to    27 “

我想初始的结构要建好，让abinit能找出所有的对称性，这能对之后的计算能加快速度，当然，如果你不想要对称性那就另当别论了。

就你的问题而言，abinit没有出错，呵呵
赞一下(10人) 回复此楼
12楼2008-12-13 11:15:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
★ ★
fegg7502(金币+2,VIP+0):鼓励讨论
感谢watermall  的解答，
也不知道为什么，
前面就无缘无故终止，也的确没“error”出现，
今天重新计算，
没有改动任何参数，
竟然成功完成计算！
在次感谢楼主的帮助！

我对计算结果还是存在一些疑惑想请教一下楼主：
附上计算
#
# Phonon density of states and projected DOS generated by anaddb
# Gaussian method with smearing =   6.19920938E-01meV, nqibz =9261
#
# omega     PHDOS    IPHDOS   PJDOS[1]  IPJDOS[1] ...
#
-1.67505420E+02  1.12940848-135  1.05021589-137  1.11500250-135  1.03682004-137  1.44059770-137  1.33958494-139
……………………………………………………………………………………………………………………………………………………

omega指的是能量or频率
PHDOS PJDOS[1]分别是总台密度和分波态密度
IPHDOS  IPJDOS[1]分别是总台密度和分波态密度的积分
我的问题在分波态密度为什么只有两个， 1.11500250-135   和  1.44059770-137 ，
而我计算的体系是包括s p d轨道的。
这点我不太理解，还请楼主帮忙解答，不胜感激！
赞一下(10人) 回复此楼
储氢家族欢迎储氢研究者！
13楼2008-12-13 22:51:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)
这个帖子是讨论声子态密度的
你说的是电子态密度，不是一回事~
赞一下 回复此楼
14楼2008-12-14 09:39:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
★ ★
fegg7502(金币+2,VIP+0):鼓励讨论
我贴的确的是声子态密度，计算的PDOS负值区也的确有声子谱光学支存在，这可能是软模（或称虚频）。
BDF分别是总的PDOS和 PJDOS[1]  PJDOS[2]，
我不太明白PJDOS[1]  PJDOS[2]分别是不是指s p 轨道，
如果是这样，那么d轨道似乎就不存在了，而计算体系是包括s p d轨道的。
还请watermall 帮忙解疑！不胜感激！
赞一下(10人) 回复此楼
储氢家族欢迎储氢研究者！
15楼2008-12-14 17:03:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)
★ ★
fegg7502(金币+2,VIP+0):thanks
引用回帖:
Originally posted by zxzj05 at 2008-12-14 17:03:
我贴的确的是声子态密度，计算的PDOS负值区也的确有声子谱光学支存在，这可能是软模（或称虚频）。
BDF分别是总的PDOS和 PJDOS[1]  PJDOS[2]，
我不太明白PJDOS[1]  PJDOS[2]分别是不是指s p 轨道，
如果是这样 ...

pjdos1, pjdos2 表示原子1， 原子2
不是什么s p d之类的，声子哪来的s ，p ，d之类的呢
赞一下(10人) 回复此楼
16楼2008-12-14 19:33:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
哦……………………对对对……
糗大了！
赞一下 回复此楼
储氢家族欢迎储氢研究者！
17楼2008-12-14 21:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
感谢watermall的指导帮助，
我还有个关于声子谱计算的问题，
在文献 PHYSICAL REVIEW B 70, 060101(R) (2004)
中有关于声子能量损失谱的计算，
不知道具体是怎样的原理和操作，
请watermall继续指点一下！
非常感谢！
赞一下 回复此楼
储氢家族欢迎储氢研究者！
18楼2008-12-15 17:20:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yy08

木虫 (著名写手)

PHD
呵呵  好东西  值得学习一下
赞一下 回复此楼
科研并快乐着
19楼2008-12-16 11:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)
引用回帖:
Originally posted by zxzj05 at 2008-12-15 17:20:
感谢watermall的指导帮助，
我还有个关于声子谱计算的问题，
在文献 PHYSICAL REVIEW B 70, 060101(R) (2004)
中有关于声子能量损失谱的计算，
不知道具体是怎样的原理和操作，
请watermall继续指点一下！
...

我看了下这篇文献，但也不知道作者是如何计算的

好像文章主要是计算了中子非弹性散射谱，其中的双声子过程我不知道abinit怎么算。好像abinit目前的教程中没有说如何算双声子过程的啊。
赞一下 回复此楼
20楼2008-12-17 23:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 watermall 的主题更新
232/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭