24小时热门版块排行榜    

查看: 1024  |  回复: 8

xiaoluoj

金虫 (正式写手)

[求助] 带电缺陷的优化为什么看不到形变

在文献中看到如下图,研究的是Si3N4,168 个原子,一个N空位。通过优化我发现“中性态”和“—”的优化后的结构是一样的,看不到如图的结果
  PS:我是用vasp的,对于“—”的优化,我是在“中性态”优化的结果上通过设置NELECT来设置带电状态的
  请高手不吝赐教,是不是哪里错了,还是方法不对
回复此楼

» 猜你喜欢

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

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

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢指教 2012-05-21 15:24:17
xiaoluoj: 金币+5 2012-05-22 12:21:19
有没有变化不是你能看出来的,你查看POSCAR里面的坐标就知道了
嘟啦啦嘟啦啦嘟
2楼2012-05-21 08:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 贺仪 at 2012-05-21 08:18:55:
有没有变化不是你能看出来的,你查看POSCAR里面的坐标就知道了

我是通过对比坐标,没有变化
3楼2012-05-21 14:28:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by xiaoluoj at 2012-05-21 14:28:42:
我是通过对比坐标,没有变化

有这种可能,毕竟这个是周期性的体系,最稳定的结构肯能会一样,具体的没试过。
但我觉得不管你怎么优化,在你优化完后,再用同样的优化方式去优化 坐标都会有变化的,不可能完全一样,我觉得你这个优化过程有问题。
嘟啦啦嘟啦啦嘟
4楼2012-05-21 16:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 贺仪 at 2012-05-21 16:52:00:
有这种可能,毕竟这个是周期性的体系,最稳定的结构肯能会一样,具体的没试过。
但我觉得不管你怎么优化,在你优化完后,再用同样的优化方式去优化 坐标都会有变化的,不可能完全一样,我觉得你这个优化过程有问题

http://g.zhubajie.com/urllink.php?id=12363221lzve3y9b50ejjeyi
这是我的INCAR文件,麻烦指出错误,谢谢
5楼2012-05-21 17:10:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
5楼: Originally posted by xiaoluoj at 2012-05-21 17:10:14:
http://g.zhubajie.com/urllink.php?id=12363221lzve3y9b50ejjeyi
这是我的INCAR文件,麻烦指出错误,谢谢

SYSTERM=beta silicon nitride  168
ISTART=0
#ICHARG=       !Default:if ISTART=0, 2;else 0.  
               ! 0_ Calculate charge density from initial wave function            
               !  1_Read the charge density from file CHGCAR
              ! 2_Take superposition of atomic charge densities         
#INIWAV=
LWAVE=.F
PREC=high      !in order to avoid wrap around errors when ELF,we specify high  --mannual 46
#electronic relaxation
ENCUT=500       !default taken from POTCAR-file;Cut-off energy for plane wave basis set in eV;
                !it has default in potcar,but  we  should set and keep it constantly for some reason
NELM=80         !default:60;NELM gives the maximum number of electronic SC (selfconsistency) steps which may be performed
EDIFF=0.1E-05   !default:1e-4  allowed error in total energy
#EDIFFG          !default:EDIFF*10
ISPIN=2
#ionic relaxtion
IBRION=2        !default:-1,ions are not moved.determines how the ions are updated and moved;if statistic ,then default
ISIF=3          !contrals whether the stress tensor is calculated. Defalt:if IIBRION=0,0; else 2;
#EDIFFG=         !default:EDIFF*10;break condition for the ionic relaxation loop;
#
ISMEAR=0         !default:1-Femi-smearing. If we calculate supercell,set it as now;or kpoints error!
SIGMA=0.05       !default:0.2
LREAL=Auto
#LELF=.TRUE      !determines whether to creat ELFCAR
#DOS && band
#LORBIT=11        
#NBANDS =       !(NELECT+NIONS)/2
#something specail
ISYM=0          !default:1; switch symmetry on (1.2 or 3) or off(0) ;if on ,i think it reduce many calculations
#NELECT



上面是你的INCAR,设置的特别乱,你的里面没有nsw,也即没有优化的步数,体系没有优化……
嘟啦啦嘟啦啦嘟
6楼2012-05-21 22:56:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

引用回帖:
6楼: Originally posted by 贺仪 at 2012-05-21 22:56:08:
SYSTERM=beta silicon nitride  168
ISTART=0
#ICHARG=       !Default:if ISTART=0, 2;else 0.  
               ! 0_ Calculate charge density from initial wave function            
               !

,我优化时一直没有设置NSW,但为什么以前我用这个INCAR优化缺陷结构时,发现缺陷周围的原子位置还是会动的?按理说没有离子优化,其位置是不动的;另外,NSW一般设置多大?
7楼2012-05-22 12:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
7楼: Originally posted by xiaoluoj at 2012-05-22 12:26:15:
,我优化时一直没有设置NSW,但为什么以前我用这个INCAR优化缺陷结构时,发现缺陷周围的原子位置还是会动的?按理说没有离子优化,其位置是不动的;另外,NSW一般设置多大?

看你的资源,一般设置最大,让尽可能的收敛,对于初始结构不好的设置10左右然后再继续优化,到最后一步我设置的是500,但是优化收敛好了 一般200左右就可以了
嘟啦啦嘟啦啦嘟
8楼2012-05-22 14:55:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kwanqsh

金虫 (小有名气)

楼主好,我看到你的INCAR上没有设置NELECT这个参数啊 那你的带点缺陷是怎么实现的? 是在POSCAR中去掉一个原子然后设置NELECT吗
9楼2014-10-16 15:45:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoluoj 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后网上办公系统维护 +9 取名字烦人 2024-06-25 10/500 2024-06-26 10:05 by 地理学1995
[基金申请] 青年基金E02口青基去年几个函评专家? +6 他山攻玉之石 2024-06-25 8/400 2024-06-26 10:00 by 降温器
[基金申请] 今天能不能出来名单 +8 地理学1995 2024-06-25 10/500 2024-06-26 09:46 by msjy
[考博] 对象没有,还非常想读博,难以抉择 +25 pvrw0224 2024-06-23 47/2350 2024-06-26 09:32 by 我是小小小小鸟
[基金申请] 博后面上到底什么时候出结果??? +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[有机交流] 对苯乙烯磺酰氯的合成机理 25+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[基金申请] 能看出是否上会了吗 +10 articlefan 2024-06-23 15/750 2024-06-25 16:05 by 请慎重修改昵称
[基金申请] 基你太美 +13 lyfbangong 2024-06-24 14/700 2024-06-25 14:16 by 华氏911.
[教师之家] 复旦夏同学退学理由说明,哪儿可以下载? +7 苏东坡二世 2024-06-21 8/400 2024-06-25 12:42 by 药嘿1233
[有机交流] 高温酯化反应喷料 20+4 hl24678 2024-06-21 7/350 2024-06-25 09:45 by ghj-hebust
[基金申请] 今年什么时候会评啊 +7 lancet0903 2024-06-24 7/350 2024-06-24 18:48 by Pickfoot
[基金申请] 说博后基金7月出的真打电话了吗? +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[基金申请] 博后面上和特助今天出吗? +41 逗您玩 2024-06-21 77/3850 2024-06-24 09:47 by weichuny
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
[基金申请] 请问三个字的中文名字,在文中到底怎么表示? 5+3 zhongyuan764 2024-06-22 6/300 2024-06-23 17:42 by 鱼翔浅底1
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
[基金申请] 国产期刊影响因子大于12的有20多个个了 +17 babu2015 2024-06-20 22/1100 2024-06-23 13:12 by 德尚中行
[基金申请] 教育部基金 +5 m1393 2024-06-21 5/250 2024-06-21 21:13 by odes
[有机交流] 怎么萃取出锡盐内包裹的化合物 +4 硕六过 2024-06-19 5/250 2024-06-21 09:50 by 光超嘟嘟
信息提示
请填处理意见