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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 【求助】请问大家在做声子计算时，这个 qpt是按照什么设置的？？？？
101/1返回列表
查看: 2192  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2009-12-27 23:33:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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.
赞一下(6人) 回复此楼
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点提出来。
赞一下(1人) 回复此楼
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的时候，是在那一步的时候开始啊？没有看懂你的意思啊啊？
赞一下(8人) 回复此楼
4楼2009-12-28 10:58: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点输入进去就可以了。
不知道你是否明白了！
赞一下(11人) 回复此楼
5楼2009-12-28 14:50:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyez

木虫 (小有名气)
谢谢指教，明白啦！！！谢谢
赞一下 回复此楼
6楼2009-12-28 19:05:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哥是传说

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
用手册上的2-2来自动生成Q点
赞一下(1人) 回复此楼
7楼2010-08-14 09:14:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lufly89

禁虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
本帖内容被屏蔽
8楼2011-12-27 20:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lufly89

禁虫 (小有名气)
本帖内容被屏蔽
9楼2011-12-27 20:43:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
17674006楼: Originally posted by y1ding at 2009-12-28 00:47:37:
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 po ...

请问计算声子带结构的时候得到了trf2_5_band2eps.freq文件，其中包含了声子频率可以用band2eps做图，但是我想用origin做图，不知道该如何转换格式能得到origin 识别的文件做图
赞一下 回复此楼
10楼2012-05-13 17:11:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wyez 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭