24小时热门版块排行榜    

查看: 808  |  回复: 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%!很不解。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
leoaqu: 金币+10, ★★★很有帮助, 谢谢 2013-01-13 18:59:21
第一,不清楚Rutile TiO2是不是在晶格原子位置上是完全占据的还是分数占据的晶体。
第二,LZ使用PSLib提供的KJPAW赝势,那就应该看到,Ti的建议截断能是50Ry,O的建议截断能是47Ry,那么就是说,计算时的截断能最小也得是50Ry。
这是我也使用同样的赝势,在LZ的输入文件的基础上,修改了截断能的设置:
CODE:
&system
   ibrav=6 ,
   celldm(1)=8.73 ,
   celldm(3)=0.643 ,
   nat=6 ,
   ntyp=2 ,
   occupations='fixed' ,
   ecutwfc = 55 ,
   ecutrho = 550 ,
/

以及得到的相应结果:
引用回帖:
CELL_PARAMETERS (alat=  8.73000000)
   1.004793478   0.000000138   0.000000000
   0.000000138   1.004793478   0.000000000
   0.000000000   0.000000000   0.641970536

与实验结果的偏差非常的小,小于1%。
思想重于技巧,内涵重于表象
2楼2013-01-13 15:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leoaqu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +21 不知道叫什么! 2026-03-15 23/1150 2026-03-22 00:12 by BruceLiu320
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿武理材料工程348求调剂 +3  ̄^ ̄゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
信息提示
请填处理意见