| 10 | 1/1 | 返回列表 |
| 查看: 2329 | 回復(fù): 9 | |||
wyez木蟲 (小有名氣)
|
[交流]
【求助】請問大家在做聲子計算時,這個 qpt是按照什么設(shè)置的???? 已有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 |
鐵桿木蟲 (著名寫手)
|
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. |
木蟲 (正式寫手)
|
其實我到不是很懂你的意思。 因為關(guān)鍵的問題是你想得到q點,然后用來算這些q點的頻率。而這些q點你是不知道的,除非你自己寫個腳本或者是像pwscf那樣,有個專門的腳本提供這些腳本。 abinit選擇了一種比較麻煩的方式: 因為你用的這些q點,一般是用mp方法插的q點,這個其實是跟k點的插法是一樣的。 比如你想用8*8*8的q點,那么你完全可以做一個8*8*8的scf,因為在做scf的時候會先產(chǎn)生一個8*8*8的k點,我們只想要這些k點,而不是想做scf計算,通常會選非常小的ecut,讓scf跑一步。這樣大約幾秒鐘,scf就能跑完,這樣你可以再OUT或者EIG文件件中找到這些k點坐標,這些k點都是不可約的。然后把這些k點當成q點輸入進去就可以了。 不知道你是否明白了! |
木蟲 (小有名氣)
鐵桿木蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
鐵桿木蟲 (正式寫手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +13 | Micky11223 2026-03-25 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 一志愿雙一流機械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 材料學碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 312,生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 一志愿吉林大學材料與化工303分求調(diào)劑 +4 | 為學666 2026-03-24 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|