| 查看: 643 | 回复: 6 | ||
llg1987金虫 (知名作家)
鱼馆的鲢鱼，其他家族的小罗
|
[求助]
【求助】计算Mg3Ce的DOS的问题
|
|
本人计算Mg3Ce，出现了一些问题，想请教一下专家。 结构优化 incar如下 SYSTEM = mg3ce ENCUT = 350 ISTART=0;ICHARG=2 ISMEAR=1;SIGMA=0.2 NSW=60;IBRION=2 ISIF=3;ISPIN=2 EDIFF=1E-5;EDIFFG=-1E-3 PREC=Accurate kpoints 选的是9*9*9，g和m-p都试过 静态自洽 incar如下 SYSTEM = mg3ce ENCUT = 350 ISTART=0;ICHARG=2 ISMEAR=-5;ISPIN=2 PREC=Accurate 其余同上（将c*拷贝过来，将contcar改为poscar） 非自洽 incar如下 SYSTEM = Mg3Ce ENCUT = 350 ISTART = 1;ICHARG = 11 ISMEAR = -5;SIGMA = 0.2 EDIFF = 1E-5 EDIFFG = -1E-3 EMIN = -3 EMAX = 16 NEDOS = 4000 LORBIT = 11 PREC = Accurate GGA = 91 VOSKOWN=1 其余同上（将c*拷贝过来，将contcar改为poscar） 计算的结果导出来后是如下图所示 自旋非自旋，G和M-P都试过了，但是结果还是一样啊 感觉那个图为什么到0.2左右就没有了呢？ 请高手支招啊，谢谢了先。  |
回复此楼» 本帖@通知
» 猜你喜欢
A期刊撤稿
已经有5人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
-
26申博自荐
已经有7人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有4人回复
-
带资进组求博导收留
已经有9人回复
-
求助大佬们，伤口沾上了乙腈
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 看看别人仅用VASP做DOS计算发好文章（两篇PRL）
已经有30人回复
- vasp计算的半导体结构，bands带隙比dos里面的大，何故？
已经有33人回复
- 虫友帮忙！band和dos计算得到的gap值差别很大，无法对应！
已经有11人回复
- vasp 计算表面DOS与文献对应不上 请求大家帮助.......
已经有20人回复
- 【求助】理论计算新手，求帮忙解决带隙及DOS的问题
已经有21人回复
- 【求助成功】VASP算出来的能带和DOS图对应不上，哪里出了问题？
已经有20人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助成功】supercell中DOS计算如何得到合理的电荷密度和波函数？
已经有3人回复
- 【求助】计算dos时很慢，而且WAVECAR很大
已经有21人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2011-10-13 16:05:43
llg1987
金虫 (知名作家)
鱼馆的鲢鱼，其他家族的小罗
- 应助: 7 (幼儿园)
- 贵宾: 2.352
- 金币: 18895.6
- 散金: 11592
- 红花: 61
- 沙发: 29
- 帖子: 6245
- 在线: 5228.4小时
- 虫号: 1089455
- 注册: 2010-09-04
- 性别: GG
- 专业: 金属材料的微观结构
- 管辖: 电脑使用交流
3楼2011-10-13 16:09:03
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★
youzhizhe(金币+2): 谢谢交流。已计入专家考核存档。 2011-10-13 21:10:50
llg1987(金币+50): 感谢侯老师提供的帮助 2011-10-14 07:39:43
youzhizhe(金币+2): 谢谢交流。已计入专家考核存档。 2011-10-13 21:10:50
llg1987(金币+50): 感谢侯老师提供的帮助 2011-10-14 07:39:43
|
你设置40000是对NEDOS，这个相当于定义态密度数据中多少个能量点。 NBANDS用来设置计算多个能带或对每个k点计算多少个本征值。默认值就是:总的电子数/2 + m (这个数的值，程序默认给的很小) 当把NBANDS的设置的很大，比如: 总电子数/2 + 100~200，Fermi level以上100～200个能带将会被计算，相应的态密度也会被计算。 此时可以不用设置EMIN和EMAX。 |
4楼2011-10-13 16:19:52
llg1987
金虫 (知名作家)
鱼馆的鲢鱼，其他家族的小罗
- 应助: 7 (幼儿园)
- 贵宾: 2.352
- 金币: 18895.6
- 散金: 11592
- 红花: 61
- 沙发: 29
- 帖子: 6245
- 在线: 5228.4小时
- 虫号: 1089455
- 注册: 2010-09-04
- 性别: GG
- 专业: 金属材料的微观结构
- 管辖: 电脑使用交流
5楼2011-10-13 17:11:31
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
6楼2011-10-13 20:33:37
liuhuimeiwoo
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1469.1
- 散金: 10
- 红花: 3
- 帖子: 235
- 在线: 40.2小时
- 虫号: 911880
- 注册: 2009-11-24
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
7楼2011-10-13 23:00:31