| 5 | 1/1 | 返回列表 |
| 查看: 2342 | 回復(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 |
|
其實(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í)候會(huì)先產(chǎn)生一個(gè)8*8*8的k點(diǎn),我們只想要這些k點(diǎn),而不是想做scf計(jì)算,通常會(huì)選非常小的ecut,讓scf跑一步。這樣大約幾秒鐘,scf就能跑完,這樣你可以再OUT或者EIG文件件中找到這些k點(diǎn)坐標(biāo),這些k點(diǎn)都是不可約的。然后把這些k點(diǎn)當(dāng)成q點(diǎn)輸入進(jìn)去就可以了。 不知道你是否明白了! |
鐵桿木蟲 (著名寫手)
|
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. |
木蟲 (正式寫手)
|
做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í)候開始。繘]有看懂你的意思啊啊? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +10 | 美味可樂雞翅 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
[基金申請] 面上5B能上會(huì)嗎? +8 | redcom 2026-03-29 | 8/400 |
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +7 | 里子木yy 2026-03-29 | 7/350 |
|
|
[考研] 286求調(diào)劑 +5 | 丟掉懶惰 2026-03-27 | 8/400 |
|
|
[考研] 291求調(diào)劑 +12 | Y-cap 2026-03-29 | 15/750 |
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 289求調(diào)劑 +16 | 新時(shí)代材料 2026-03-27 | 16/800 |
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|