| 3 | 1/1 | 返回列表 |
| 查看: 1353 | 回復(fù): 2 | |||||
snowmami新蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】wien2k-FAQ實(shí)在是看不懂,哪位高手幫我翻譯下?
|
|
下面的文字對(duì)我很重要,但實(shí)在看不懂,很著急,難者不會(huì),會(huì)者不難,各位高手能否幫忙翻譯下? WIEN2k-FAQ: Setup of antiferromagnetic calculations? ©2001 by P. Blaha, K. Schwarz and J. Luitz Back to: FAQ Index Setup of antiferromagnetic calculations requires some tricks which will be described below: 1. You need to construct a unit cell which allows for the desired AF ordering. E.g. for bcc Cr you have to put "P" lattice and 2 NON-equivalent atoms at (0,0,0) and (0.5,0.5,0.5). In order to make nn working properly, you should also specify different labels in the "name(3:10)"-field for the two AF-atoms (Eg. "Cr1" and "Cr2" .2. Run init_lapw. nn and/or sgroup determines equivalent atoms. Accept eventually the new struct_nn file (You may have to use the same procedure as described in faq:supercells). 3. Edit case.inst and flip the spin of the second AF atom. This means you must invert the spin-up and dn occupation of this atom. Atoms which carry no spin (e.g. an O between two AF metal atoms) should be made non-spinpolarized by setting equal numbers of up and dn-electrons. For AFM NiO case.inst looks like: Ni Ar 3 5 3, 2,2.0 N 3, 2,2.0 N 3,-3,3.0 N 3,-3,1.0 N 4,-1,1.0 N 4,-1,1.0 N Ni Ar 3 5 3, 2,2.0 N 3, 2,2.0 N 3,-3,1.0 N <=== spin flipped 3,-3,3.0 N <=== spin flipped 4,-1,1.0 N 4,-1,1.0 N O He 3 5 2,-1,1.0 N 2,-1,1.0 N 2, 1,1.0 N 2, 1,1.0 N 2,-2,1.0 N <=== equal occup for up/dn 2,-2,1.0 N <=== equal occup for up/dn **** End of Input **** End of Input This case.inst file can also be prepared by the new: instgen_lapw -ask , where you can specify for each atom spin-up, dn or non-magnetic configurations. 4. Continue with init_lapw (rerunning lstart) and select "spin-polarized" calculations. 5. You may now run an scf cycle using runsp_lapw; or 6. You may try to select "AFM-calculations", generate "case.inclmcopy" and use "runafm_lapw". For simplest use of "afminput" you should provide the nonmagnetic (or ferromagnetic) supergroup in "case.struct_supergroup". Please read the usersguide for more details. (help_lapw and search for AFM). If you succeed to generate proper input files, you can do the scf-cycle by "runafm_lapw", which will save at least half the computer time, since only one spin-direction is calculated and the other one is generated by symmetry. However, since these programs are little tested, I recommend to check the final results of "runafm_lapw" by 7. "runsp_lapw". If the charge distance (grep IS case.scf) remains small and the magnetic moments do not change, runafm_lapw was successful, otherwise continue with runsp_lapw until convergence.Peter Blaha, [ Last edited by snowmami on 2010-12-8 at 21:59 ] |
wien2k常見(jiàn)問(wèn)題集錦 | wien2k大雜燴 |
» 搶金幣啦!回帖就可以得到:
+1/183
+1/82
+1/77
+1/47
+1/26
+1/23
+1/22
+1/16
+1/13
+1/9
+1/7
+1/7
+1/5
+1/4
+1/4
+1/3
+1/3
+1/2
+1/2
+1/2
金蟲(chóng) (著名寫(xiě)手)
|
進(jìn)行反鐵磁計(jì)算的步驟: 1、創(chuàng)建一個(gè)AF需要的單胞。例如:對(duì)于體心Cr,應(yīng)該選取晶格對(duì)稱(chēng)性P(取消對(duì)稱(chēng)性),2個(gè)不等價(jià)原子位置(0,0,0)和(0.5,0.5,0.5).為了使nn選取合適,應(yīng)該給兩個(gè)AF原子不同的標(biāo)記(例如:Cr1&Cr2) 2、運(yùn)行init_lapw。nn或者sgroup需要等價(jià)原子,接受新結(jié)構(gòu)文件struct_nn 3、編輯case.inst并且翻轉(zhuǎn)第二個(gè)AF原子自旋。這意味著翻轉(zhuǎn)這個(gè)原子占據(jù)態(tài)的上下自旋。對(duì)于不帶自旋的原子(例如O原子),應(yīng)該使用非自旋極化(上下自旋電子數(shù)相等)。 對(duì)于AFM的NIO 的case,inst為:略 ***這個(gè)case.inst文件也可以用instgen_lapw -ask命令產(chǎn)生,這個(gè)命令可以標(biāo)記任何原子為上自旋,下自旋,或者非磁。(我在計(jì)算AFM時(shí)都是使用這個(gè)命令) 4、繼續(xù)init_lapw并且選擇自旋極化計(jì)算,完成初始化。 5、初始化完成后運(yùn)行runsp_lapw或者進(jìn)行第6操作 6、選擇AFM-calculations,生成case.inclmcopy文件,使用runafm_lapw。對(duì)于afminput最簡(jiǎn)單的使用,你應(yīng)該用case.struct_supergroup提供非磁或鐵磁supergroup。詳細(xì)閱讀userguide(help_lapw and search for AFM) 成功生成輸入文件后,可以運(yùn)行自恰循環(huán)通過(guò)runafm_lapw。因?yàn)橹挥?jì)算上自旋,其他自旋方向的計(jì)算都是通過(guò)對(duì)稱(chēng)性完成的,所以可以減少一般的計(jì)算時(shí)間。 7、如果charge distance 仍然很小,磁矩沒(méi)有改變,運(yùn)行runsp_lapw,再runafm_lapw。如果還不成功繼續(xù)runsp_lapw直到收斂。 |
新蟲(chóng) (正式寫(xiě)手)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +6 | @taotao 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 307求調(diào)劑 +8 | 超級(jí)伊昂大王 2026-03-24 | 9/450 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|