| 24小時(shí)熱門(mén)版塊排行榜 |
| 2 | 1/1 | 返回列表 |
| 查看: 814 | 回復(fù): 1 | |||||
leoaqu金蟲(chóng) (正式寫(xiě)手)
|
[求助]
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%!很不解。 |
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +174 |
|
第一,不清楚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è)置: 以及得到的相應(yīng)結(jié)果: 與實(shí)驗(yàn)結(jié)果的偏差非常的小,小于1%。 |

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +8 | 一樣YWY 2026-04-02 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +24 | 材料學(xué)257求調(diào)劑 2026-03-28 | 25/1250 |
|
|
[考研] 一志愿9初試366 本雙非求調(diào)劑 +4 | 運(yùn)氣來(lái)得若有似?/a> 2026-04-02 | 4/200 |
|
|
[考研] 354求調(diào)劑 +4 | lxb598 2026-03-31 | 5/250 |
|
|
[考研] 292求調(diào)劑 +17 | 木蟲(chóng)er12138 2026-04-01 | 17/850 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +21 | 荷包蛋hyj 2026-03-26 | 21/1050 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠(yuǎn). 2026-03-28 | 6/300 |
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 4/200 |
|
|
[考研] 一志愿西交大080500材料學(xué)碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
[考研] 315求調(diào)劑 +6 | akie... 2026-03-28 | 7/350 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|