| 5 | 1/1 | 返回列表 |
| 查看: 2268 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
xiaoxiao45至尊木蟲 (正式寫手)
|
[求助]
berkeleygw 輸入文件問(wèn)題 已有1人參與
|
||
|
有哪位大神指導(dǎo)一下berkeleygw中的WFN.in文件lattice vector怎么計(jì)算的嗎?為什么kgrid.x運(yùn)行后bravais lattice和crystal 的space group總是對(duì)不上?謝謝了 另外知道晶體的point group: R3c a=b=c= 5.667 alfa=beta=gama=59.17 原子坐標(biāo) 1 Bi,2 Fe,3 O 1 0.00000000 0.00000000 0.00000000 1 0.50000000 0.50000000 0.50000000 2 0.22120000 0.22120000 0.22120000 2 0.72120000 0.72120000 0.72120000 3 0.53500000 0.94000000 0.39400000 3 0.94000000 0.39400000 0.53500000 3 0.03500000 0.89400000 0.44000000 3 0.39400000 0.53500000 0.94000000 3 0.89400000 0.44000000 0.03500000 3 0.44000000 0.03500000 0.89400000 |
至尊木蟲 (正式寫手)
|
我的數(shù)據(jù)是用的是cta Materialia 59 (2011) 4229–4234 原子分?jǐn)?shù)坐標(biāo)是按照Wyckoff position設(shè)置的,對(duì)稱性為R3c,然后我把所有的位置數(shù)據(jù)拷貝到WFN.in中,只是不知道lattice vector (WFN.in 文件需要輸入的)怎么計(jì)算, 然后用你所說(shuō)的方法先算scf(pw.x),進(jìn)而把計(jì)算結(jié)果帶到WFN.in中計(jì)算,得到的WFN.log文件,文件中有布拉菲對(duì)稱和晶體對(duì)稱性,結(jié)果得不到R3c的結(jié)果,因?yàn)槲腋杏X如果得不到所需要的結(jié)果,那么后面得到約化的K-point是不是也不對(duì)。 那么接來(lái)的wfn_q和wfn_co就無(wú)法計(jì)算了。 所以,這幾天一直糾結(jié)著這問(wèn)題,頭發(fā)又掉不少,不知道liqizuiyang遇到過(guò)此類問(wèn)題沒(méi)? 求幫忙指條明路 |
木蟲 (著名寫手)
至尊木蟲 (正式寫手)
木蟲 (著名寫手)
|
如果用QE算DFT部分,有個(gè)簡(jiǎn)便方法。 先進(jìn)行scf計(jì)算,輸出文件中這一部分就是基矢坐標(biāo): crystal axes: (cart. coord. in units of alat) a(1) = ( 4.582581 4.582581 0.000000 ) a(2) = ( 0.000000 4.582581 4.582581 ) a(3) = ( 4.582581 0.000000 4.582581 ) 這一部分就是同一坐標(biāo)系下的原子坐標(biāo): Cartesian axes site n. atom positions (alat units) 1 O tau( 1) = ( 2.3756976 2.2068836 2.3756976 ) 2 O tau( 2) = ( 4.6669876 4.4981749 2.3756989 ) 3 O tau( 3) = ( 2.3756989 4.4981749 4.6669876 ) 4 O tau( 4) = ( 4.6669876 6.7894636 4.6669876 ) 5 O tau( 5) = ( 4.4981749 2.3756989 4.4981749 ) 6 O tau( 6) = ( 6.7894636 4.6669876 4.4981749 ) 7 O tau( 7) = ( 4.4981749 4.6669876 6.7894636 ) 8 O tau( 8) = ( 6.7894648 6.9582788 6.7894648 ) 9 Mg tau( 9) = ( 2.2912906 0.0000000 2.2912906 ) 10 Mg tau( 10) = ( 4.5825812 2.2912906 2.2912906 ) 11 Mg tau( 11) = ( 2.2912906 2.2912906 4.5825812 ) 12 Mg tau( 12) = ( 4.5825812 4.5825812 4.5825812 ) 13 Ca tau( 13) = ( 0.0000000 0.0000000 0.0000000 ) 14 Ca tau( 14) = ( 2.2912906 2.2912906 0.0000000 ) 15 Ca tau( 15) = ( 0.0000000 2.2912906 2.2912906 ) 16 Ca tau( 16) = ( 2.2912906 4.5825812 2.2912906 ) |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 349求調(diào)劑 +9 | zwjjjjjj 2026-03-31 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +13 | jjjjjjjjjjl 2026-03-27 | 13/650 |
|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +9 | 雙馬尾痞老板2 2026-04-01 | 9/450 |
|
|
[考研] 271求調(diào)劑 +15 | 勒布朗@ 2026-03-31 | 20/1000 |
|
|
[基金申請(qǐng)] esi高被引論文是不是能對(duì)中標(biāo)有所加分和幫助呢 +4 | redcom 2026-04-01 | 4/200 |
|
|
[考研] 08生物與醫(yī)藥專碩初試346找調(diào)劑 +6 | dianeeee 2026-04-01 | 7/350 |
|
|
[考研] 350求調(diào)劑 +7 | 阿佳~ 2026-03-31 | 7/350 |
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
[考研] 考研材料工程351分調(diào)劑 +5 | 整個(gè)好的 2026-03-31 | 5/250 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-31 | 4/200 |
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|