| 查看: 928 | 回复: 1 | |||
leoaqu金虫 (正式写手)
|
[求助]
vc-relax 优化rutile TiO2结构问题
|
|
刚开始用QE， 就想先从一个简单的体系开始学习，就先学着用vc-relax优化TiO2结构，发现最后的晶格常数的优化结果严重比正常值变小。小了大概15%！！！ 我换过不同的potential， （GGA-PBE, PAW-PBE 等）结果都大同小异。很困惑，求高人指点。下面是input file。（虽然过程中我也尝试着改过好多参数， 但结果几乎都没太大改进。所以就贴出下面的一个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 优化的结果大概是： 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是他应该的理想值，但优化的结果只有理想值的 86%！很不解。 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- TiO2（101）晶面的优化
已经有13人回复
- VASP一次性优化锐钛矿TiO2晶格常数时出错！！
已经有6人回复
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
2楼2013-01-13 15:03:17
投票:
