版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

leoaqu

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 1221.8
散金: 2
紅花: 1
帖子: 336
在線: 47.3小時(shí)
蟲(chóng)號(hào): 592997
注冊(cè): 2008-09-05
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] vc-relax 優(yōu)化rutile TiO2結(jié)構(gòu)問(wèn)題

剛開(kāi)始用QE，就想先從一個(gè)簡(jiǎn)單的體系開(kāi)始學(xué)習(xí)，就先學(xué)著用vc-relax優(yōu)化TiO2結(jié)構(gòu)，發(fā)現(xiàn)最后的晶格常數(shù)的優(yōu)化結(jié)果嚴(yán)重比正常值變小。小了大概15%�。�！我換過(guò)不同的potential，（GGA-PBE, PAW-PBE 等）結(jié)果都大同小異。很困惑，求高人指點(diǎn)。下面是input file。（雖然過(guò)程中我也嘗試著改過(guò)好多參數(shù)，但結(jié)果幾乎都沒(méi)太大改進(jìn)。所以就貼出下面的一個(gè)input file)
&control
calculation='vc-relax',
pseudo_dir = '~/pslib',
outdir='./',
prefix='TiO2',
&wf_collect=.true.,
wfcdir='./',
tprnfor = .true.
&restart==================
&restart_mode='from_scratch'
  &restart_mode='restart'

/
&system
ibrav=6, celldm(1)=8.73, celldm(3)=0.643,
nat=6, ntyp=2,
occupations='smearing', smearing='mp', degauss=0.05
ecutwfc=15.0,
/

&electrons
diagonalization='david'
mixing_beta=0.7
conv_thr=1.0d-5
/
&IONS
ion_dynamics='bfgs',
/

&cell
cell_dynamics='bfgs',
cell_dofree=xyz,
cell_factor=3.0D0
/

ATOMIC_SPECIES
Ti 47.88  Ti.pbe-spn-kjpaw_psl.0.2.3.UPF
O 16.0 O.pbe-n-kjpaw_psl.0.1.UPF

ATOMIC_POSITIONS crystal
Ti 0.0 0.0 0.0
Ti 0.5 0.5 0.5
O  0.3053 0.3053 0.0
O  0.6947 0.6947 0.0
O  0.8053 0.1947 0.5
O  0.1947 0.8053 0.5

K_POINTS automatic
6 6 6 0 0 0

優(yōu)化的結(jié)果大概是：
CELL_PARAMETERS (alat=  8.73000000)
0.858842467 0.000000000 0.000000000
0.000000000 0.858842467 0.000000000
0.000000000 0.000000000 0.569949709

ATOMIC_POSITIONS (crystal)
Ti    0.000000000 0.000000000 0.000000000
Ti    0.500000000 0.500000000 0.500000000
O       0.298652764 0.298652764 0.000000000
O       0.701347236 0.701347236 0.000000000
O       0.798569657 0.201430343 0.500000000
O       0.201430343 0.798569657 0.500000000
End final coordinates

其中alat是他應(yīng)該的理想值，但優(yōu)化的結(jié)果只有理想值的 86%！很不解。

回復(fù)此樓

» 猜你喜歡

【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有3人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠(chéng)求調(diào)劑已經(jīng)有14人回復(fù)
085410 一志愿211 22408分?jǐn)?shù)359求調(diào)劑已經(jīng)有4人回復(fù)
271求調(diào)劑已經(jīng)有19人回復(fù)
385分生物學(xué)（071000）求調(diào)劑已經(jīng)有3人回復(fù)
一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑已經(jīng)有3人回復(fù)
318求調(diào)劑，計(jì)算材料方向已經(jīng)有8人回復(fù)
291求調(diào)劑已經(jīng)有25人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

TiO2（101）晶面的優(yōu)化已經(jīng)有13人回復(fù)
VASP一次性優(yōu)化銳鈦礦TiO2晶格常數(shù)時(shí)出錯(cuò)�。� 已經(jīng)有6人回復(fù)

1樓 2013-01-12 11:51:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

souledge

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +174
1ST強(qiáng)帖: 19
應(yīng)助: 303 (大學(xué)生)
金幣: 11306.4
散金: 1123
紅花: 108
帖子: 2814
在線: 916.1小時(shí)
蟲(chóng)號(hào): 941486
注冊(cè): 2010-01-12
專(zhuān)業(yè): 結(jié)構(gòu)陶瓷
管轄: 第一性原理

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
leoaqu: 金幣+10, ★★★很有幫助, 謝謝 2013-01-13 18:59:21

第一，不清楚Rutile TiO2是不是在晶格原子位置上是完全占據(jù)的還是分?jǐn)?shù)占據(jù)的晶體。
第二，LZ使用PSLib提供的KJPAW贗勢(shì)，那就應(yīng)該看到，Ti的建議截?cái)嗄苁?0Ry，O的建議截?cái)嗄苁?7Ry，那么就是說(shuō)，計(jì)算時(shí)的截?cái)嗄茏钚∫驳檬?0Ry。
這是我也使用同樣的贗勢(shì)，在LZ的輸入文件的基礎(chǔ)上，修改了截?cái)嗄艿脑O(shè)置：

CODE:

&system

   ibrav=6 ,

   celldm(1)=8.73 ,

   celldm(3)=0.643 ,

   nat=6 ,

   ntyp=2 ,

   occupations='fixed' ,

   ecutwfc = 55 ,

   ecutrho = 550 ,

/

以及得到的相應(yīng)結(jié)果：

引用回帖:

CELL_PARAMETERS (alat= 8.73000000)
1.004793478 0.000000138 0.000000000
0.000000138 1.004793478 0.000000000
0.000000000 0.000000000 0.641970536

與實(shí)驗(yàn)結(jié)果的偏差非常的小，小于1%。

贊一下(2人)

回復(fù)此樓

思想重于技巧，內(nèi)涵重于表象

2樓2013-01-13 15:03:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 leoaqu 的主題更新