| 查看: 2801 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】请问大家在做声子计算时，这个 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 |
回复此楼» 猜你喜欢
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
极限拉扯，如何抉择
已经有7人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 采用vasp ，phononpy或frophon计算声子谱的具体方法是什么？
已经有5人回复
- 今天用CASTEP计算声子谱，出现下列提示，有能看懂的大侠帮忙解释下
已经有11人回复
- 为何用晶格常数的实验值计算声子？
已经有18人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- Abinit计算声子问题
已经有10人回复
- ZnO的声子计算
已经有18人回复
- MS计算晶体声子谱的时候，横坐标有负能量值，为什么？
已经有11人回复
- castep计算Fe2O3声子谱失败，原因求教~~
已经有10人回复
- 自旋体系的声子计算应该用什么软件呢？
已经有5人回复
- 计算Sb2Te3的声子色散曲线，请教用什么软件合适,VASP+phonopy怎么样
已经有7人回复
- 【求助】castep计算声子谱，为什么不能进行振动分析
已经有9人回复
- 【求助】请教CASTEP计算声子谱后出现虚频
已经有6人回复
- 【求助】有人用有限差分方法计算过G点声子谱吗？
已经有10人回复
- 【求助】问大家一个基础概念，什么是晶体的简正频率？在声子谱中能读出简正频率吗？
已经有4人回复
- 【求助】用MS中的castep计算声子振动模频率，具体怎么算
已经有10人回复
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
qasd(金币+3,VIP+0):xiexie~ 12-28 20:15
小木虫(金币+0.5):给个红包，谢谢回帖交流
qasd(金币+3,VIP+0):xiexie~ 12-28 20:15
|
其实我到不是很懂你的意思。 因为关键的问题是你想得到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
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
TIGERYZZ(金币+3,VIP+0):谢谢指点 12-28 08:29
wyez(金币+1,VIP+0):谢谢 12-28 09:57
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
3楼2009-12-28 08:44:37
hhj0514
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2540.2
- 帖子: 364
- 在线: 110.3小时
- 虫号: 345910
- 注册: 2007-04-14
- 专业: 凝聚态物性I:结构、力学和
4楼2009-12-28 10:58:37
