24小时热门版块排行榜    

查看: 2196  |  回复: 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的回帖

wyez

木虫 (小有名气)

谢谢指教,明白啦!!!谢谢
6楼2009-12-28 19:05:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

y1ding

铁杆木虫 (著名写手)

★ ★ ★ ★
TIGERYZZ(金币+3,VIP+0):谢谢指点 12-28 08:29
wyez(金币+1,VIP+0):谢谢 12-28 09:57
ANADDB utility (see later), that is, they should form a reduced set of symmetry-inequivalent wavevectors, corresponding to a regularly spaced grid. In principle, they ought not include the Gamma point, but it is recommended to have it in the set, in order for the Fourier interpolation not to introduce errors at that important point. In order to minimize the number of preliminary non-self-consistent calculations, it is advised to take a q point mesh that is adjusted to the k point mesh used for the electronic structure : all q wavevectors should connect two k point wavevectors from this grid.

Such a set of q wavevectors can be generated straightforwardly by running a GS calculation with kptopt=1, nshiftk=1, shiftk=0 0 0 (to include gamma) and taking the output kpt set file as this qpt set. One might set nstep=1 and nline=1, so only one iteration runs.
The input file ~abinit/tests/tutorespfn/Input/trf2_2.in is precisely an input file that can be used to generate such a set of k points. Copy it in the present Work_rf2 directly, as well as the accompanying ~abinit/tests/tutorespfn/Input/trf2_2.files. Examine these files, then run this calculation (it is very rapid - it won't hurt the trf2_1 job). The following k point set is obtained :

       kpt    0.00000000E+00  0.00000000E+00  0.00000000E+00
              2.50000000E-01  0.00000000E+00  0.00000000E+00
              5.00000000E-01  0.00000000E+00  0.00000000E+00
              2.50000000E-01  2.50000000E-01  0.00000000E+00
              5.00000000E-01  2.50000000E-01  0.00000000E+00
             -2.50000000E-01  2.50000000E-01  0.00000000E+00
              5.00000000E-01  5.00000000E-01  0.00000000E+00
             -2.50000000E-01  5.00000000E-01  2.50000000E-01
It is, as promised, the same as the q point set in the trf2_1.in file.
2楼2009-12-28 00:47:37
已阅   回复此楼   关注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的回帖

hhj0514

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
TIGERYZZ(金币+2,VIP+0):谢谢交流 12-28 14:05
做scf的时候,nstep=1,nline=1 自洽只算一次。
kptopt=1,nshiftk=1,shiftk=0 0 0  
ngkpt=6 6 6 (这里是你想要的多少个k点)
输出的时候有些人是看out文件,但是如果k点多于50个的时候,在kpt是写不全的。这是可以看EIG文件,利用脚本把你想要的k点提出来。

我想请问下,做scf的时候,是在那一步的时候开始啊?没有看懂你的意思啊啊?
4楼2009-12-28 10:58:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见