24小时热门版块排行榜    

查看: 2685  |  回复: 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

银虫 (小有名气)

★ ★ ★ ★
小木虫(金币+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的回帖
查看全部 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金未中,这种答复是模板吗? +4 zhaosm1982 2026-08-27 5/250 2026-08-27 12:17 by zhaosm1982
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +5 Kittylucky 2026-08-27 5/250 2026-08-27 12:10 by 姜太公烤鱼
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[基金申请] 国自然面上复盘~欢迎讨论 +7 晴天加油 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见