| 查看: 815 | 回复: 0 | |||
cenwanglai荣誉版主 (知名作家)
老和山猥琐派九段
|
[交流]
【求助完毕】超胞DOS计算与K-mesh及ISMEAR选取
|
|
结构优化 INCAR: ENCUT=480 EDIFF = 1E-4 EDIFFG = -0.01 ISMEAR =0 IBRION =2 ISPIN = 2 MAGMON = 2*3 64*0 30*0 ALGO =F LREAL = A KPOINTS A 0 M-K 1 1 1 0 0 0 静态计算 INCAR: #IBRION = 2 ENCUT = 480 EDIFF = 1E-4 #EDIFFG = -0.01 ISMEAR = -5 SIGMA =0.05 ISPIN = 2 MAGMOM = 2*3 64*0 30*0 KPOINTS: A 0 M-K 3 3 3 0 0 0 ****************************************** 问题： （1）结构优化后收敛，静态计算电子自洽慢很多，最后得到的力也远大于0.01 （2）结构优化后自旋mag=6.00；静态计算mag=1.x 为什么会这样呢？结构优化后，我修改了ISMEAR, SIGMA和k-mesh，其它基本没变呢。 后来发现是这样的： ISMEAR的选取与K点取样是有关系的。 对于半导体超胞结构，K点取样数可以适当减少，与之对应的ISMEAR比较适合取0。 超胞Gamma点结构优化后，加密k-mesh（如3×3×3），ISMEAR=0做静态自洽计算，得到DOSCAR。要得到高精度的DOS，有两个办法。一个是取较大的k-mesh，如 8×8×8,然后ISMEAR=-5;第二个是选取特别的k点，然后全部列出这些k点，ISMEAR = -5.对于超胞，第二种方法计算量可以接受，第一种则可能吃不消。 [ Last edited by cenwanglai on 2010-8-26 at 12:56 ] |
回复此楼» 猜你喜欢
评审有感
已经有21人回复
-
反应很差，大量原料没有反应
已经有6人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
-
上海大学实验技术岗位非升即走
已经有13人回复
-
今年审到国自然15份，谈谈感受
已经有18人回复
-
献血感触
已经有11人回复
-
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有5人回复
-
2026山东省优青
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
关于影响因子计算方法的疑问
已经有55人回复
- 【求助】在计算单点能时scf=tight的作用，非加不可吗？
已经有6人回复
- 细胞裂解
已经有13人回复
- CASTEP中GGA+U计算纳米线
已经有4人回复
- 关于VASP计算得到CHGCAR文件的数据处理以及数据单位的
已经有16人回复
- 三方晶体的晶胞体积怎么计算？应该有公式的吧？
已经有13人回复
- split_dos分割以后，如何找到对应的原子？
已经有11人回复
- 北京军区总医院生物治疗中心招聘干细胞、免疫细胞研究助研数人
已经有13人回复
- [求助]如何处理这样的审稿意见：The study is descriptive in nature
已经有12人回复
- 关于TiN超晶胞掺杂的问题
已经有11人回复
- vasp计算的半导体结构，bands带隙比dos里面的大，何故？
已经有33人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- 【求助】理论计算新手，求帮忙解决带隙及DOS的问题
已经有21人回复
- 【求助】ISMEAR的探讨
已经有18人回复
- 【求助】DOS计算中，设置ISMEAR=-5为什么还要用SIGMA？
已经有4人回复
- 【求助完毕】Xcrysden提供的k-path是reciprocal的还是cartesian格式的？
已经有7人回复
- 【重点讨论】计算能带时如何确定高对称点和k-points path?
已经有40人回复
- 【求助成功】supercell中DOS计算如何得到合理的电荷密度和波函数？
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30
