 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1290 | 回復(fù): 2 | |||||
三葉蟲的腳印新蟲 (初入文壇)
|
[求助]
如果想要用4*4*4個(gè)q點(diǎn)計(jì)算一個(gè)四方晶系的聲子譜，q點(diǎn)應(yīng)該怎么設(shè)？ 已有1人參與
|
|
如果想要用4*4*4個(gè)q點(diǎn)計(jì)算一個(gè)四方晶系的聲子譜，q點(diǎn)應(yīng)該怎么設(shè)？還有提示中的問(wèn)題和虛頻問(wèn)題 輸入qua-relax.in #Number of data sets ndtset 2 #Dataset1:Initial self-consistent run getwfk1 0 # Cancel default kptopt1 1 #for the ground state iscf1 7 #Self-consistent calculation, using algorithm 7 tolvrs1 1.0d-18 #SCF stopping criterion #Dataset2:Response-function calculation, with q=0 getwfk2 -1 #Use GS wave functions from dataset1 kptopt2 2 #for response functions at q=0 iscf2 7 #Self-consistent calculation, using algorithm 5 rfphon2 1 #Will consider phonon-type perturbation rfatpol2 1 6 #All the atoms will be displaced rfdir2 1 1 1 #Along all reduced coordinate axis nqpt2 1 #One wavevector is to be considered ngqpt 4 4 4 tolvrs2 1.0d-22 #SCF stopping criterion #Common input variables acell 7.2673047890E+00 7.2673047890E+00 4.6380985420E+00 angdeg 90 90 90 ntypat 2 znucl 14 8 natom 6 typat 2*1 4*2 xred 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 2.9967303949E-01 2.9967303949E-01 1.3729447097E-18 2.0032696051E-01 7.9967303949E-01 5.0000000000E-01 7.0032696051E-01 7.0032696051E-01 -1.3729447097E-18 7.9967303949E-01 2.0032696051E-01 5.0000000000E-01 ecut 80 ngkpt 10 10 10 nshiftk 1 shiftk 0.5 0.5 0.5 nstep 100 #bands and occupation tsmear 0.02 # kptopt 1 diemac 12.0 ixc 1 optforces 1 chksymbreak 0 log文件中的提示： --- !WARNING src_file: distrb2.F90 src_line: 113 message: | nkpt*nsppol (75) is not a multiple of nproc_kpt (24) The k-point parallelisation is not efficient. ... --- !WARNING src_file: mpi_setup.F90 src_line: 622 message: | Your number of spins*k-points (=75) will not distribute correctly with the current number of processors (=24). You will leave some empty. YOU ARE STRONGLY ADVICED TO ACTIVATE AUTOMATIC PARALLELIZATION! PUT "AUTOPARAL=1" IN THE INPUT FILE. ... --- !WARNING src_file: m_hdr.F90 src_line: 3980 message: | input nkpt=500 not equal disk file nkpt=75 ... --- !WARNING src_file: m_hdr.F90 src_line: 4214 message: | input kptopt= 2 /= disk file kptopt= 1 ...  --- !WARNING src_file: m_hdr.F90 src_line: 4509 message: | Restart of self-consistent calculation need translated wavefunctions. ... out文件中的虛頻。。 Phonon frequencies in cm-1 : - -1.704008E+02 1.012839E+01 1.398780E+01 1.398780E+01 4.443393E+02 - 5.406203E+02 5.406203E+02 7.389403E+02 7.902434E+02 7.902434E+02 - 8.843817E+02 8.843817E+02 9.283579E+02 1.058092E+03 1.066444E+03 - 1.180454E+03 1.180454E+03 1.298081E+03 小女子不才，初學(xué)ABINIT,望各路大牛海涵~ |
回復(fù)此樓» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有5人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有6人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有21人回復(fù)
-
26申博自薦
已經(jīng)有5人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有3人回復(fù)
-
調(diào)劑310
已經(jīng)有12人回復(fù)
-
化學(xué)0703 調(diào)劑 306分 一志愿211
已經(jīng)有6人回復(fù)
souledge
專家顧問(wèn) (著名寫手)
-
專家經(jīng)驗(yàn): +174 - 1ST強(qiáng)帖: 19
- 應(yīng)助: 303 (大學(xué)生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時(shí)
- 蟲號(hào): 941486
- 注冊(cè): 2010-01-12
- 專業(yè): 結(jié)構(gòu)陶瓷
- 管轄: 第一性原理
|
呃~現(xiàn)在可以直接用nqpt了么？記得以前老版本時(shí)的套路是先在kptopt = 1，nshiftk = 1，shiftk = (0,0,0)的情況下跑一個(gè)簡(jiǎn)單的scf，通過(guò)ngkpt得到q-mesh。 然后，再在計(jì)算聲子的輸入文件中，按照上一個(gè)計(jì)算中k點(diǎn)的量來(lái)設(shè)置ndtset，然后每個(gè)k點(diǎn)在本輸入文件中為qpt，每個(gè)qpt是一個(gè)單獨(dú)的聲子計(jì)算…… 另外，不排除一些計(jì)算在運(yùn)用asr之前有小的虛頻，所以還是需要anaddb處理之后看圖中是否存在虛頻，不過(guò)-170cm^(-1)的虛頻倒是不小…… |
2樓2016-12-22 23:47:27
Maple_autumn
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 314.8
- 紅花: 1
- 帖子: 53
- 在線: 20.3小時(shí)
- 蟲號(hào): 5692387
- 注冊(cè): 2017-02-22
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
| 首先，有個(gè)基本的問(wèn)題，你在用格點(diǎn)法取k點(diǎn)和q點(diǎn)的時(shí)候，怎么能取ngkpt 10 10 10，ngqpt 4 4 4呢？這樣的話，在你取得聲子的格點(diǎn)上你計(jì)算電子了嗎？ngqpt的參數(shù)需要能整除ngkpt。第二，acell三個(gè)值不相等，格點(diǎn)怎么能取10 10 10,4 4 4了？第三，rf計(jì)算能用iscf=7嗎？第四，聲子計(jì)算的時(shí)候似乎應(yīng)該考慮參數(shù)nline？這個(gè)我也不是很清楚。至于報(bào)錯(cuò)似乎是個(gè)并行運(yùn)算問(wèn)題，不過(guò)你先把上面的問(wèn)題改好試試？樓上說(shuō)的方法在現(xiàn)在的版本也能用來(lái)作為參考的。 |
3樓2017-03-05 17:25:26
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 三葉蟲的腳印 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +7 | 李李不服輸 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +4 | 零八# 2026-03-28 | 4/200 |
|
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|