24小时热门版块排行榜

返回列表

jerryeast

金虫 (初入文坛)

应助: 4 (幼儿园)
金币: 1427.7
散金: 4
帖子: 36
在线: 34.7小时
虫号: 910362
注册: 2009-11-23
性别: GG
专业: 结构化学

[求助] 表面计算，有没有人解释一下这个是什么情况

之前一直用Castep计算，最近在学QE，试计算了一个表面，结构优化后出现下面的结果，什么原因？Castep里面优化完就是表面原子略微浮动，QE这个计算一看也知道肯定不对，什么原因？

下面分别是优化前和优化后的图

Screenshot from 2012-12-17 18^%13^%00.png

Screenshot from 2012-12-17 18^%13^%14.png

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

有没有人遇到在烧制氧化锌低压压敏瓷时后样品表面出现一些小气泡的情况？已经有8人回复
纳米零价铁表面的米勒指数是什么？已经有4人回复
ZnS表面计算能带和DOS显示金属性质是怎么回事？已经有19人回复
（文献学习）一篇计算HfO2表面结构与原子性质的文章已经有26人回复
计算Cl-在Fe表面吸附问题时体系加电荷方法已经有27人回复
怎样通过实验计算固体的表面自由能？已经有6人回复
CH4在Ir（111）表面的解离吸附过渡态计算：Henkelman的PRL文章已经有10人回复
计算Cu表面能,为什么不随层数增加而减少?反而逐渐增加? 已经有15人回复
【求助】表面能的计算：几层薄膜最稳定？已经有7人回复

1楼2012-12-18 10:14:19

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

huazhorg

铁杆木虫 (正式写手)

应助: 21 (小学生)
金币: 4993.8
散金: 5
红花: 2
帖子: 381
在线: 363.2小时
虫号: 629899
注册: 2008-10-18
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
jerryeast: 金币+2, ★有帮助 2012-12-18 10:40:51

金红石110面？只计算应该不会算成这个样子，我怀疑是你晶格参数设置有误。能否贴出你的输入文件？

赞一下

回复此楼

2楼2012-12-18 10:33:09

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

jerryeast

金虫 (初入文坛)

应助: 4 (幼儿园)
金币: 1427.7
散金: 4
帖子: 36
在线: 34.7小时
虫号: 910362
注册: 2009-11-23
性别: GG
专业: 结构化学

引用回帖:

2楼: Originally posted by huazhorg at 2012-12-18 10:33:09
金红石110面？只计算应该不会算成这个样子，我怀疑是你晶格参数设置有误。能否贴出你的输入文件？

&CONTROL
                     title = P42MNM110 ,
               calculation = 'relax' ,
            restart_mode = 'from_scratch' ,
                  wfcdir = '/tmp' ,
                  prefix = '110' ,
               lkpoint_dir = .true. ,
                  disk_io = 'default' ,
               verbosity = 'default' ,
            etot_conv_thr = 1.D-5 ,
            forc_conv_thr = 1.D-3 ,
                     nstep = 100 ,
                  tstress = .false. ,
                  tprnfor = .true. ,
                  tefield = .false. ,
                  dipfield = .false. ,
                  lelfield = .false. ,
                     lorbm = .false. ,
                  lberry = .false. ,
/
&SYSTEM
                     ibrav = 14,
                        A = 5.9119254486128001 ,
                        B = 6.4354402059501000 ,
                        C = 32.1613233140812010 ,
                     cosAB = 0 ,
                     cosAC = 0 ,
                     cosBC = 0 ,
                     nat = 48,
                     ntyp = 2,
                  ecutwfc = 40 ,
                  ecutrho = 400 ,
                     nbnd = 400,
               tot_charge = 0.000000,
               occupations = 'smearing' ,
      one_atom_occupations = .false. ,
      starting_spin_angle = .false. ,
                  degauss = 1 ,
                  smearing = 'gaussian' ,
                  noncolin = .false. ,
         exxdiv_treatment = 'gygi-baldereschi' ,
/
&ELECTRONS
         electron_maxstep = 300,
         scf_must_converge = .true. ,
                  conv_thr = 1.D-6 ,
            adaptive_thr = .false. ,
               startingpot = 'atomic' ,
               startingwfc = 'atomic+random' ,
               mixing_mode = 'plain' ,
               mixing_beta = 0.5 ,
               mixing_ndim = 10,
         diagonalization = 'david' ,
            diago_thr_init = 1.D-5 ,
            diago_full_acc = .false. ,
                     tqr = .false. ,
/
&IONS
            ion_dynamics = 'bfgs' ,
            ion_positions = 'from_input' ,
               phase_space = 'full' ,
         pot_extrapolation = 'second_order' ,
         wfc_extrapolation = 'second_order' ,
         remove_rigid_rot = .false. ,
/
ATOMIC_SPECIES
O 15.99900  O.pbe-van_bm.UPF
Mn 54.93800  Mn.pbe-sp-van_mit.UPF
ATOMIC_POSITIONS crystal
O    0.250000000 0.500000000 0.000000000 0  0  0
O    0.250000000 0.000000000 0.100049369 0  0  0
O    0.250000000 0.500000000 0.200098738
O    0.250000000 0.000000000 0.300148107
O    0.250000000 0.500000000 0.077986933 0  0  0
O    0.250000000 0.000000000 0.178036302 0  0  0
O    0.250000000 0.500000000 0.278085671
O    0.250000000 0.000000000 0.378135040
O    0.000000000 0.194871127 0.038993467 0  0  0
O    0.000000000 0.694871127 0.139042835 0  0  0
O    0.000000000 0.194871127 0.239092204
O    0.000000000 0.694871127 0.339141573
O    0.000000000 0.805128873 0.038993467 0  0  0
O    0.000000000 0.305128873 0.139042835 0  0  0
O    0.000000000 0.805128873 0.239092204
O    0.000000000 0.305128873 0.339141573
O    0.750000000 0.500000000 0.000000000 0  0  0
O    0.750000000 0.000000000 0.100049369 0  0  0
O    0.750000000 0.500000000 0.200098738
O    0.750000000 0.000000000 0.300148107
O    0.750000000 0.500000000 0.077986933 0  0  0
O    0.750000000 0.000000000 0.178036302 0  0  0
O    0.750000000 0.500000000 0.278085671
O    0.750000000 0.000000000 0.378135040
O    0.500000000 0.194871127 0.038993467 0  0  0
O    0.500000000 0.694871127 0.139042835 0  0  0
O    0.500000000 0.194871127 0.239092204
O    0.500000000 0.694871127 0.339141573
O    0.500000000 0.805128873 0.038993467 0  0  0
O    0.500000000 0.305128873 0.139042835 0  0  0
O    0.500000000 0.805128873 0.239092204
O    0.500000000 0.305128873 0.339141573
Mn    0.250000000 0.000000000 0.038993467 0  0  0
Mn    0.250000000 0.500000000 0.139042835 0  0  0
Mn    0.250000000 0.000000000 0.239092204
Mn    0.250000000 0.500000000 0.339141573
Mn    0.000000000 0.500000000 0.038993467 0  0  0
Mn    0.000000000 0.000000000 0.139042835 0  0  0
Mn    0.000000000 0.500000000 0.239092204
Mn    0.000000000 0.000000000 0.339141573
Mn    0.750000000 0.000000000 0.038993467 0  0  0
Mn    0.750000000 0.500000000 0.139042835 0  0  0
Mn    0.750000000 0.000000000 0.239092204
Mn    0.750000000 0.500000000 0.339141573
Mn    0.500000000 0.500000000 0.038993467 0  0  0
Mn    0.500000000 0.000000000 0.139042835 0  0  0
Mn    0.500000000 0.500000000 0.239092204
Mn    0.500000000 0.000000000 0.339141573
K_POINTS gamma

赞一下

回复此楼

3楼2012-12-18 10:40:29

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