| 查看: 1744 | 回复: 4 | ||
那年的小艾金虫 (初入文坛)
|
[求助]
关于态密度计算的一些问题
|
|
本人先做了晶体的结构优化后，将INCAR CONTCAR POTCAR KPOINTS考入新建的DOS文件夹中并将CONTCAR改成了POSCAR，准备进行态密度的计算。请问INCAR中的设置哪些需要修改的（体系中有Eu,要考虑自旋） 新手一枚，望各位多加指点。 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算态密度的时候设置的NPAR是什么意思
已经有17人回复
- vasp计算态密度结果出问题
已经有11人回复
- VASP 计算DOS，态密度中费米能的作用
已经有3人回复
- 能带结构和态密度问题
已经有19人回复
- 怎么根据态密度图积分算区间的电子数
已经有8人回复
- 关于太阳能电池态密度和载流子复合问题
已经有8人回复
- 用vasp如何计算能带图和态密度图呢？
已经有11人回复
- VASP计算态密度出现问题
已经有12人回复
- abinit能否处理一个位置被两种元素占据的固溶体的能带和态密度
已经有10人回复
- 【求助】关于MS软件中能带和态密度的计算问题
已经有8人回复
- 【求助】计算AlN的态密度请教
已经有4人回复
- 用MS的castep模块计算能带结构和态密度，为什么很长时间计算还不结束
已经有6人回复
- 【求助】呼叫zhangguangping...关于态密度的计算
已经有11人回复
- 【求助】MS中局域态密度的计算
已经有3人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】关于vasp态密度计算
已经有10人回复
- 【求助】关于态密度的计算的基础问题
已经有10人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】MS５.０ Castep 计算不出部分态密度（PDOS）怎么办？（站内有人提问MS4.3，�
已经有19人回复
- 【求助】怎样由能带结构和态密度看计算的结构是不是可能存在超导现象
已经有12人回复
- 【求助】态密度计算究竟是自洽计算，还是非自洽计算？
已经有24人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
identation: 金币+5, 谢谢指教 2013-10-28 17:43:26
那年的小艾: 金币+30, ★★★很有帮助, 非常感谢，很详细的解答 2013-10-28 20:33:50
感谢参与，应助指数 +1
identation: 金币+5, 谢谢指教 2013-10-28 17:43:26
那年的小艾: 金币+30, ★★★很有帮助, 非常感谢，很详细的解答 2013-10-28 20:33:50
|
VASP中电子态密度计算的流程 主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算 以Al－FCC为例子 第一步 结构优化 输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步 静态自洽计算 INCAR： PREC = Medium，ISTART = 0，ICHARG = 2，ISMEAR = -5 输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步计算是在结构优化的结果上进行的，所以开始第二步的时候，将第一步中的输入文件INCAR, POTCAR, POSCAR, KPOINT 以及 C* 文件放入静态自洽计算中去，并且将CONTCAR 拷贝到 POSCAR中，然后运行VASP。计算结果中的Fermi能是准确的，需要记下。 第三步 非自洽计算 INCAR： PREC = Medium，ICHARG = 11，ISMEAR = -5，LORBIT = 10 或者 11（这时可不设RWIGS)，ISTART = 1 在第二步自洽计算的基础上进行，修改输入文件INCAR, POTCAR, POSCAR, KPOINT。 INCAR文件 System=Al ISTART=1 ISMEAR=-5 SIGMA=0.2 ICHARG=11 RWIGS=1.402 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 1.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 21 21 21 0.0 0.0 0.0 VASP中文手册：在进行能带和DOS计算时，ISMEAR 不能使用N阶MP方法。因为MP方法在空轨道上有负的占据，所以求得的能带和DOS是不正确的。但是从其它地方看到---“提示：在计算能带结构时，采ISMEAR = 0或1对结果的影响非常小，可以认为是一样的。但是不能采用ISMEAR = -5 或-4。”ISMEAR到底多少？ 计算能带：ICHARG = 11 导体的话，用ISMEAR=1； 半导体或绝缘体，用ISMEAR=0 。 计算 DOS： ICHARG = 11 ISMEAR = -5 计算的时候，金属可用0、1，非金属不要大过0，体材料可用－4、－5（面的话就用－1、0吧？） 参考http://210.45.66.120/doc/vasp5.0/node224.html 设置完成后进行计算，计算完后，得到包含了态密度值的DOSCAR文件，采用split_dos对态密度文件DOSCAR进行分割，得到总态密度DOS0，各个原子的分波态密度DOS1，DOS2……。另外在运行split_dos程序对DOSCAR文件分割时，要保证当前目录下有对应的OUTCAR和POSCAR文件。 分割后的DOS0，DOS1…等文件的能量值是以费米能级作为能量参考零点。DOS0的第一列数据是能量值，单位为eV；第二列数据是总态密度的值，单位 State/eV.unit cell；第三列数据是总态密度的积分值，也就是电子数，单位为electrons。DOS1是第一个原子的分波态密度值，其中的第一列数据是能量值，单位为eV；第二、三、四列数据分别对应于s、p、d态的分波态密度值，单位为State/eV.atom。其他的DOS文件与DOS1类似。 |
2楼2013-10-28 16:20:37
3楼2013-11-13 12:46:03
4楼2013-11-13 15:17:36
5楼2013-11-13 17:01:35
