| 5 | 1/1 | 返回列表 |
| 查看: 2333 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
wyez木蟲 (小有名氣)
|
[交流]
【求助】請問大家在做聲子計(jì)算時,這個 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. |
木蟲 (正式寫手)
|
其實(shí)我到不是很懂你的意思。 因?yàn)殛P(guān)鍵的問題是你想得到q點(diǎn),然后用來算這些q點(diǎn)的頻率。而這些q點(diǎn)你是不知道的,除非你自己寫個腳本或者是像pwscf那樣,有個專門的腳本提供這些腳本。 abinit選擇了一種比較麻煩的方式: 因?yàn)槟阌玫倪@些q點(diǎn),一般是用mp方法插的q點(diǎn),這個其實(shí)是跟k點(diǎn)的插法是一樣的。 比如你想用8*8*8的q點(diǎn),那么你完全可以做一個8*8*8的scf,因?yàn)樵谧鰏cf的時候會先產(chǎn)生一個8*8*8的k點(diǎn),我們只想要這些k點(diǎn),而不是想做scf計(jì)算,通常會選非常小的ecut,讓scf跑一步。這樣大約幾秒鐘,scf就能跑完,這樣你可以再OUT或者EIG文件件中找到這些k點(diǎn)坐標(biāo),這些k點(diǎn)都是不可約的。然后把這些k點(diǎn)當(dāng)成q點(diǎn)輸入進(jìn)去就可以了。 不知道你是否明白了! |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +4 | Nadiums 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|