Vasp 石墨烯能带计算

这里暂时不做结构优化，使用数据库结构

首先做自洽计算产生正确的基态电子密度

INCAR

SYSTEM = C
PREC = accutat
ISTART = 0
ICHARG = 2
EDIFF = 0.1E-03
NSW = 0
IBRION = 2
ISIF = 2

POSCAR

C
2.45944
1.000 0.000 0.000
0.500 0.866 0.000
0.000 0.000 15.000
2
Cartesian
0.000000 0.000000 0.000000
0.500000 0.288667 0.000000

KPOINTS

Auto
0
Monkhorst Pack
5 5 5
0 0 0

POTCAR文件使用GGA泛函，从赝势文件拷过来。

执行计算脚本

第二步，固定电子密度，得到选取给定点的能量

修改INCAR

SYSTEM = C
PREC = accutat
ISTART = 1
ICHARG = 11
EDIFF = 0.1E-03
NSW = 0
IBRION = 2
ISIF = 2

修改KPOINTS

Auto
40
Line-mode
cart
0.0 0.0 0.0
1.000000 0.0 0.0

再次提交任务计算，计算的结果已经写在vasprun.xml文档中

后续借助pymatgen工具绘制能带图

导入用到的python模块

import matplotlib.pyplot as plt
from pymatgen.io.vasp.outputs import Vasprun
from pymatgen.electronic_structure.plotter import BSDOSPlotter,
BSPlotter,BSPlotterProjected,DosPlotter

matplotlib.pyplot用于画图，

pymatgen.io.vasp.outputs.Vasprun读取vasp的输出文件vasprun.xml。

pymatgen.electronic_structure.plotter是专门用来绘制Dos和Band图的模块。

读入vasp的计算结果

bs_vasprun = Vasprun("./vasprun.xml",parse_projected_eigen=True)
bs_data = bs_vasprun.get_band_structure(line_mode=True)

第一行用来解析vasprun.xml文件生成指定的对象bs_vasprun

第二行用来获取能带数据

由于后面的KPOINTS不是对称取，所以暂时没有解决画图问题

这个坑未完待续

Vasp 石墨烯能带计算相关推荐

VASP载流子有效质量计算
VASP载流子有效质量计算在计算完能带之后,可以对载流子(电子+空穴)的有效质量进行计算,考察其光电性质.以Si的能带结构为例,如图所示,其价带顶位于高对称点Gamma处,导带底大约位于X处. 由图 ...
matlab 画石墨烯能带,锯齿石墨烯能带matlab计算程序
这个程序是初步优化后的matlab版本,主要思路是先生成体系的格点坐标,再运用坐标生成体系的哈密顿量,然后进行对角化计算能带,能带的计算使用一维体系超原胞的处理方法.可以进一步优化主程序 nx=3; ...
vasp 模拟退火_VASP 计算问题小结
本文转载自博主一个人就是一个叠加态,有部分删减修改,文中对相关概念方法做了详细的总结,留坑待填... 1. 第一原理计算的一些心得 1. 第一性原理第一性原理其实是包括基于密度泛函的从头算和基于Ha ...
石墨烯能带matlab,锯齿石墨烯能带matlab计算程序
这个程序是初步优化后的matlab版本,主要思路是先生成体系的格点坐标,再运用坐标生成体系的哈密顿量,然后进行对角化计算能带,能带的计算使用一维体系超原胞的处理方法.可以进一步优化主程序 nx=3; ...
vasp 自恰计算_恰逢假期–如何最大程度地利用货币销售
vasp 自恰计算 Rewarded video ads can be an effective way to generate revenue for your mobile game. In fa ...
VASP+FDM+Phonopy 计算声子过程
以Ammonia为例,基本翻译自:https://atztogo.github.io/phonopy/procedure.html#tutorial. 1. 结构优化. 输入文件如下: INCAR ! ...
基于VASP+phonopy+shengbte计算声子相干的热学性质
声子相干的热学性质在计算材料领域有着非常重要的地位,我们可以通过第一性原理计算得到任一种晶体的热力学性质,但是很多小伙伴尤其是初学者还比较陌生,本文介绍基于DFT计算与声子相关的性质,是用VASP+p ...
vasp服务器中断,求助VASP能带计算的中断原因
我也遇到类似情况,发现服务器内存没有用到多少?是不是并行环境有问题?谢谢! DA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and K ...
vasp服务器中断,求助VASP能带计算的中断原因 - 第一原理 - 小木虫 - 学术科研互动社区...
我也遇到类似情况,发现服务器内存没有用到多少?是不是并行环境有问题?谢谢! DA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and K ...
vasp能带结构计算
Si (band structure) VASP中能带结构计算流程与态密度相同,具体不同点如下: 1. scf自洽计算与态密度dos计算中的scf自洽计算完全相同. 2. 能带结构非自洽计算能带计 ...

Vasp 石墨烯能带计算

导入用到的python模块

读入vasp的计算结果

Vasp 石墨烯能带计算相关推荐

最新文章

热门文章