| 10 | 1/1 | 返回列表 |
| 查看: 2332 | 回復(fù): 9 | |||
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)去就可以了。 不知道你是否明白了! |
木蟲 (小有名氣)
鐵桿木蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
鐵桿木蟲 (正式寫手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +4 | 小小機(jī)靈蟲 2026-03-29 | 10/500 |
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運(yùn)來上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。。 +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|