 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2151 | 回復(fù): 6 | ||
zhaoyue_ucla鐵蟲 (初入文壇)
|
[求助]
鎢（tungsten）cohesive energy計(jì)算
|
|
本人QE剛?cè)腴T小白，W的scf計(jì)算input 如下，請各位大俠幫忙看看問題（我算的結(jié)果5.8ev跟實(shí)際結(jié)果8.5ev差距較大。如果有可能的話，有QE高手愿意語音或者視頻幫助解決問題，可以給30~50RMB/per hour的酬勞，不甚感激。 &CONTROL calculation = 'scf' , pseudo_dir = './' , prefix = 'W_exc' , lkpoint_dir = .false. , / &SYSTEM ibrav = 1, celldm(1) = 5.964, nat = 2, ntyp = 1, ecutwfc = 300 , occupations = 'smearing' , degauss = 0.06 , smearing = 'marzari-vanderbilt' , / &ELECTRONS mixing_beta = 0.7 , / ATOMIC_SPECIES W 183.84000 W.pw-mt_fhi.UPF ATOMIC_POSITIONS alat W 0.000000000 0.000000000 0.000000000 W 0.500000000 0.500000000 0.500000000 K_POINTS automatic 6 6 6 1 1 1 對ecutwfc的測試： cut(Ry) 30 50 150 300 E（Ry） -129.05979680 -148.60944687 -209.35569942 -210.14896228 350 370 -210.15691865 -210.16924485 為啥我都測到370Ry, 總能還不能完全收斂��？（變化>0.01Ry) 以及K點(diǎn)設(shè)置的時(shí)候如何確定要不要平移啊？ |
回復(fù)此樓» 猜你喜歡
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有5人回復(fù)
-
279求調(diào)劑
已經(jīng)有3人回復(fù)
-
336求調(diào)劑
已經(jīng)有3人回復(fù)
-
286求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有3人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有14人回復(fù)
-
332求調(diào)劑
已經(jīng)有3人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué)
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 急求解決-計(jì)算Fe的cohesive energy
已經(jīng)有3人回復(fù)
- 【求助】cohesive energy如何用CASTEP計(jì)算
已經(jīng)有17人回復(fù)
liqizuiyang
木蟲 (著名寫手)
- 應(yīng)助: 59 (初中生)
- 金幣: 6859.5
- 散金: 261
- 紅花: 94
- 帖子: 2093
- 在線: 708.9小時(shí)
- 蟲號: 2366686
- 注冊: 2013-03-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★
ljw4010: 金幣+2, 謝謝指導(dǎo)解答！ 2014-11-20 20:38:14
ljw4010: 金幣+2, 謝謝指導(dǎo)解答！ 2014-11-20 20:38:14
|
1. W屬副族重金屬元素，相對論效應(yīng)顯著，因此必須選用full relativistic贗勢，并且包含半芯態(tài)。W.pw-mt_fhi.UPF只是一般的模守恒贗勢，交換關(guān)聯(lián)泛函是LDA，很有可能不符合要求。 2.W屬體心立方晶系，應(yīng)該是 ibrav = 2, ATOMIC_POSITIONS crystal W 0.000000000 0.000000000 0.000000000 不需要第二個(gè)原子 另外 celldm(1)一定要是優(yōu)化后的晶格常數(shù)。 |
2樓2014-11-18 14:24:15
liqizuiyang
木蟲 (著名寫手)
- 應(yīng)助: 59 (初中生)
- 金幣: 6859.5
- 散金: 261
- 紅花: 94
- 帖子: 2093
- 在線: 708.9小時(shí)
- 蟲號: 2366686
- 注冊: 2013-03-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-11-18 14:25:36
zhaoyue_ucla
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 104.9
- 帖子: 14
- 在線: 18.7小時(shí)
- 蟲號: 3465223
- 注冊: 2014-10-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
非常感謝你的回復(fù)！ 還有問題想請教你一下 1.celldm(1)我是去wikipedia上找的晶格常數(shù)5.964,請問你說的優(yōu)化怎么進(jìn)行呀，我剛接觸這個(gè)方向不是很懂。 2.問一個(gè)小白問題，ibrav=2,nat= 1,ntyp=1,W 0 0 0和 ibrav = 1,nat = 2,ntyp=1，W 0 0 0 , W 0.5 0.5 0.5等價(jià)嗎？我用xcryden畫出來兩個(gè)結(jié)果是一樣的，但是總能量E卻相差一倍，有點(diǎn)困惑。 3.關(guān)于贗勢的選擇。我導(dǎo)師的文章上用的是siesta算的，他讓我用QE來算，導(dǎo)師文章上是這么說的：The ab initio calculations are performed using the SIESTA [22] code, which carries out self-consistent density functional theory calculations with norm-conserving pseudopotentials and a flexible linear combination of numerical atomic orbital basis. The core electrons are replaced by Troullier–Martins pseudopotentials [24] and the valence electrons are described by a double-zeta polar- ized orbital basis set [22]. For the exchange–correlation part of the energy functional, we used the local density approximation (LDA) parameterized by [4]. The conver- gence of the energy is achieved for a plane-wave cut-off of 300 Ry. Throughout the calculations, both local and volumetric relaxations are carried out using the Broyden conjugate gradient algorithm and a force tolerance of 0.01 eV A ̊ 􏰛。其實(shí)我完全不知道怎么去選這個(gè)贗勢，我就是在QE官網(wǎng)上找了個(gè)跟文章上最接近的贗勢:W.pw-mt_fhi.UPF Pseudopotential type: NORMCONS Method: Martins-Troullier Functional type: Perdew-Wang LDA non relativistic 請問我應(yīng)該選擇官網(wǎng)上給的哪一個(gè)贗勢才能算的比較準(zhǔn)呀？或者我是不是得自己做個(gè)贗勢出來？不甚感激~ |
4樓2014-11-18 15:01:43
zhaoyue_ucla
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 104.9
- 帖子: 14
- 在線: 18.7小時(shí)
- 蟲號: 3465223
- 注冊: 2014-10-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2014-11-18 15:06:30
liqizuiyang
木蟲 (著名寫手)
- 應(yīng)助: 59 (初中生)
- 金幣: 6859.5
- 散金: 261
- 紅花: 94
- 帖子: 2093
- 在線: 708.9小時(shí)
- 蟲號: 2366686
- 注冊: 2013-03-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★
ljw4010: 金幣+2, 謝謝指導(dǎo)解答！ 2014-11-20 20:38:43
ljw4010: 金幣+2, 謝謝指導(dǎo)解答！ 2014-11-20 20:38:43
|
1. 這個(gè)可以在站內(nèi)找下優(yōu)化晶格常數(shù)的教程. 2.QE有兩種指定晶體結(jié)構(gòu)的方法，第一種基于晶胞，用內(nèi)置的十四種布拉伐格子，W是體心立方，所以是2。另一種是用原胞，自己手輸基矢坐標(biāo)。體心立方的基矢不是(a,0,0)(0,a,0)(0,0,a)，你原先是按照(a,0,0)(0,a,0)(0,0,a)來理解的。 3.http://www.quantum-espresso.org/ ... 0&origin_id去這里找標(biāo)著pz/pw和scalar relativistic的贗勢。 如果計(jì)算資源允許，推薦用ELK或exciting。這兩個(gè)程序不使用贗勢，并且使用綴加平面波，算W這種重元素可能更準(zhǔn)確。 |
6樓2014-11-18 15:15:44
7樓2017-03-08 11:46:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhaoyue_ucla 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +11 | 呼呼呼。。。。 2026-03-17 | 11/550 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|