| 5 | 1/1 | 返回列表 |
| 查看: 2330 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
wyez木蟲 (小有名氣)
|
[交流]
【求助】請問大家在做聲子計(jì)算時(shí),這個(gè) 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 |
木蟲 (正式寫手)
|
做scf的時(shí)候,nstep=1,nline=1 自洽只算一次。 kptopt=1,nshiftk=1,shiftk=0 0 0 ngkpt=6 6 6 (這里是你想要的多少個(gè)k點(diǎn)) 輸出的時(shí)候有些人是看out文件,但是如果k點(diǎn)多于50個(gè)的時(shí)候,在kpt是寫不全的。這是可以看EIG文件,利用腳本把你想要的k點(diǎn)提出來。 我想請問下,做scf的時(shí)候,是在那一步的時(shí)候開始。繘]有看懂你的意思啊。 |
鐵桿木蟲 (著名寫手)
|
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)你是不知道的,除非你自己寫個(gè)腳本或者是像pwscf那樣,有個(gè)專門的腳本提供這些腳本。 abinit選擇了一種比較麻煩的方式: 因?yàn)槟阌玫倪@些q點(diǎn),一般是用mp方法插的q點(diǎn),這個(gè)其實(shí)是跟k點(diǎn)的插法是一樣的。 比如你想用8*8*8的q點(diǎn),那么你完全可以做一個(gè)8*8*8的scf,因?yàn)樵谧鰏cf的時(shí)候會先產(chǎn)生一個(gè)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ā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +6 | minht 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 348求調(diào)劑 +5 | 小懶蟲不懶了 2026-03-28 | 5/250 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|