24小时热门版块排行榜    

查看: 2199  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wyez

木虫 (小有名气)

[交流] 【求助】请问大家在做声子计算时,这个 qpt是按照什么设置的???? 已有3人参与

nqpt   1            # One qpt for each dataset (only 0 or 1 allowed)
                         # This is the default for all datasets and must
                         #  be explicitly turned off for dataset 1.

     qpt2   0.00000000E+00  0.00000000E+00  0.00000000E+00
     qpt3   0.00000000E+00  0.00000000E+00  0.00000000E+00
     qpt4   2.50000000E-01  0.00000000E+00  0.00000000E+00
     qpt5   5.00000000E-01  0.00000000E+00  0.00000000E+00
     qpt6   2.50000000E-01  2.50000000E-01  0.00000000E+00
     qpt7   5.00000000E-01  2.50000000E-01  0.00000000E+00
     qpt8  -2.50000000E-01  2.50000000E-01  0.00000000E+00
     qpt9   5.00000000E-01  5.00000000E-01  0.00000000E+00
     qpt10 -2.50000000E-01  5.00000000E-01  2.50000000E-01
回复此楼

» 猜你喜欢

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

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

xbelinda

银虫 (小有名气)

★ ★
wyez(金币+2,VIP+0):谢谢帮助!我在研究下 12-28 09:56
做scf的时候,nstep=1,nline=1 自洽只算一次。
kptopt=1,nshiftk=1,shiftk=0 0 0  
ngkpt=6 6 6 (这里是你想要的多少个k点)
输出的时候有些人是看out文件,但是如果k点多于50个的时候,在kpt是写不全的。这是可以看EIG文件,利用脚本把你想要的k点提出来。
3楼2009-12-28 08:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
qasd(金币+3,VIP+0):xiexie~ 12-28 20:15
引用回帖:
Originally posted by hhj0514 at 2009-12-28 10:58:
做scf的时候,nstep=1,nline=1 自洽只算一次。
kptopt=1,nshiftk=1,shiftk=0 0 0  
ngkpt=6 6 6 (这里是你想要的多少个k点)
输出的时候有些人是看out文件,但是如果k点多于50个的时候,在kpt是写不全的。这是 ...

其实我到不是很懂你的意思。
因为关键的问题是你想得到q点,然后用来算这些q点的频率。而这些q点你是不知道的,除非你自己写个脚本或者是像pwscf那样,有个专门的脚本提供这些脚本。
abinit选择了一种比较麻烦的方式:
因为你用的这些q点,一般是用mp方法插的q点,这个其实是跟k点的插法是一样的。
比如你想用8*8*8的q点,那么你完全可以做一个8*8*8的scf,因为在做scf的时候会先产生一个8*8*8的k点,我们只想要这些k点,而不是想做scf计算,通常会选非常小的ecut,让scf跑一步。这样大约几秒钟,scf就能跑完,这样你可以再OUT或者EIG文件件中找到这些k点坐标,这些k点都是不可约的。然后把这些k点当成q点输入进去就可以了。
不知道你是否明白了!
5楼2009-12-28 14:50:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wyez 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见