精华内容
下载资源
问答
  • 在分析态密度(DOS)时,常常用p4vasp这个软件来读取vasprun.xml,进而获得需要的各元素或轨道对应的图形。然后,利用p4vasp的export功能将数据输出画图。具体操作参见 https://www.bigbrosci.com/2018/02/08/ex39/ ...

    背景


    在分析态密度(DOS)时,常常用p4vasp这个软件来读取vasprun.xml,进而获得需要的各元素或轨道对应的图形。然后,利用p4vasp的export功能将数据输出画图。具体操作参见

    https://www.bigbrosci.com/2018/02/08/ex39/
    注:本例中的数据来源也是该博文,感谢原作者的分享

    然而,在将数据输出为dat格式后,发现一个问题,即:p4vasp把所有的数据都由上到下排在两列中,画图时面临着只有一个x, 一个y的情况,画出来的图如下:
    在这里插入图片描述

    在这里插入图片描述
    而软件中的样子(我们想画成的样子)起码有颜色区别:
    在这里插入图片描述

    观察了数据后发现,每组数据是301行,之后有一个空行,然后再接301行,如此反复。当然数据量少可以手动剪切粘贴,但是数据很多的时候就会很繁琐。因此,想着利用python中的pandas功能来将这些数据每302行(含最后一个空行)为一组,取出,放在新的两列,然后再进行画图。 整理好思路,我们先把dat文件放到excel转为csv格式(记得分列),然后运行Python程序(py文件和csv在一个文件夹即可)

    代码如下:

    # -*- coding: utf-8 -*-
    """
    Created on Sat Dec 12 20:15:08 2020
    
    @author: fya
    """
    
    #观察数据,301行后有一个空行。
    import pandas as pd
    
    df0=pd.read_csv('testRu.csv',header=None) #记得更改为自己的csv文件名
    l=len(df0)  #获取长度
    l=l+1 #把最后一行空行也算上,要不然后面除以302会有错误
    print(l)
    print(df0.head())
    
    df=df0.iloc[0:301,:].reset_index(drop=True) #前301行数据为第一组,0:301是取从第0行到第300行,第301不取。reset_index(drop=True) 是为了去除原有的index,重新从0生成index
    
    print(df.tail())
    
    
    df_values=locals()
    l=int(l/302) #用除法后要变为int,不然默认float,到后面range部分会出错
    print(l)
    for i in range(2,l+1):     #range括号左含右不含
        l1=302*i-302
        l2=302*i-1
        df_values['df'+str(i)]=df0.iloc[l1:l2,:].reset_index(drop=True)   #随着i的不同,每301行都取出,index不受原始数据影响,都是从0开始
        df=pd.concat([df,df_values['df'+str(i)]],axis=1,join='outer') #将取出的数据横着排列,axis=1,形成新的df,然后与下一个i产生的新的列继续合并直到结束
    
    print(df)
    df.to_csv('dos.csv',header=None,index=None) #不要column和index,以免画图时出错
    print('Done!')
    

    用Python分好组是如下图这样的:
    在这里插入图片描述

    放入origin,设置对应的x,y,画图如下:

    在这里插入图片描述

    方法2


    此方法无需再转存为csv文件,可以直接用dat文件,需要注意读入数据后没有空行

    # -*- coding: utf-8 -*-
    """
    Created on Sat Dec 12 20:15:08 2020
    
    @author: fya
    """
    
    
    import pandas as pd
    
    df0=pd.read_table('testRu.dat',sep='\s+',header=None,engine='python') #直接读取dat文件,读取后无空行
    l=len(df0)  #获取长度
    #l=l+1 #把最后一行空行也算上
    print(l)   # 得到值:903
    print(df0.head())
    
    df=df0.iloc[0:301,:].reset_index(drop=True) #前301行数据为第一组(0-300行,第301不取),reset_index(drop=True) 是为了去除原有的index,重新从0生成index
    
    print(df.tail())
    
    
    df_values=locals()
    l=int(l/301) #用除法后要变为int,不然默认float,到后面range部分会出错
    print(l)
    for i in range(2,l+1):     #range括号左含右不含
        l1=301*(i-1)
        l2=301*i
        df_values['df'+str(i)]=df0.iloc[l1:l2,:].reset_index(drop=True)   #随着i的不同,每301行都取出,index不受原始数据影响,都是从0开始
        df=pd.concat([df,df_values['df'+str(i)]],axis=1,join='outer') #将取出的数据横着排列,axis=1,形成新的df,然后与下一个i产生的新的列继续合并直到结束
    
    print(df)
    df.to_csv('dos.dat',header=None,index=None)
    print('Done!')
    

    此法输出dat文件,拖到origin直接画图,得到和方法一同样的结果

    展开全文
  • 二:计算态密度的参数准备 对于态密度计算INCAR参数的一些思考 ISTAT=1 这里我们选择1,但可以选择的还有0,在INCAR中不填,系统应该会默认为1 。他的目的是继承上一步的静态计算当中得到的WAVECAR,以减少计算的...

    一:计算前的准备工作

    1,得到一个稳定的弛豫过后的POSCR结构

    2,对结构进行一次静态计算

    二:计算态密度的参数准备

    对于态密度计算INCAR参数的一些思考

    1. ISTAT=1 这里我们选择1,但可以选择的还有0,在INCAR中不填,系统应该会默认为1 。他的目的是继承上一步的静态计算当中得到的WAVECAR,以减少计算的耗时。
    2. ICHARG=11 这一项还可以选择 0,1,2和12 选择用11进行计算的话就会读取我们在上一步静态计算过程中得到的CHGCAR并进行态密度计算这样,可以节约时间,以便平衡加大K点等其他参数所所带来的耗时增加。
    3. IBRION=-1 这里选择 -1 我们在结构优化的时候得到的体系已经是最稳定的这里选择-1就是把所有原子固定住。(思考: 为什么要固定?如果不固定是不是也能算呢?固定是为了节约计算量以节约时间吗?如果不固定也能算,那不固定对态密度峰值精度的影响大吗?)
    4. NSW=0 因为不需要调整结构所以不需要弛豫步数
    5. LCHCAR=.F. 数据有点大,可以输出,但一般不选择输出节约空间
    6. LWAVE=.F. 数据有点大,可以输出,但一般不选择输出节约空间
    7. ISMEAR=0 (疑问:我都是用了0,但我也不知道为什么)
    8. SIGMA=0.1 这个参数对能带或态密度中间的间隙影响很大,越接近0能带间隙就会越大,但相应计算的耗时也会增加(所以有时也不经怀疑自己,这弄计算的纯粹就是在凑)
    9. LORBIT=11 这个参数可以将每一个原子对应的轨道的贡献输出到OUTCAR里
    10. ISPIN=2 这个参数用于计算磁性,慎用,开启后会态密度图会成上下分布,对于磁性体系会出现上下分裂的情况
    11. NEDOS=3000 加大态密度计算的撒点数,这样可以让态密度的峰值更平顺
    12. EMIN=-20 决定态密度中x轴(能量)的最小取值范围
    13. EMAX=20 决定态密度中x轴(能量)的最大取值范围
      其余的值没有写在INCAR里了,系统会自动使用默认值

    对于态密度计算KPOINTS取值的一些思考

    取值尽量大一点

    三:计算态密度需要的文件

    INCAR
    KPOINTS
    POSCAR (上一步静态得到的CONTCAR)
    POTCAR
    WAVCAR (上一步静态得到的)
    CHGCAR (上一步静态得到的)

    四:计算态密度

    将任务投入计算列队

    五:处理数据

    将计算得到的vasprun.xml使用p4vasp.exe处理,导出数据
    后使用matlab或origin绘图

    最后:留疑

    静态计算的数据也可以绘制态密度图,但经过态密度加强过的峰值会更平缓,

    展开全文
  • 态密度图OriginPro导出数据 p4vasp 能带&态密度图 画能带图:使用文件vasprun.xml(input-bands中的) 画态密度图:使用文件vasprun.xml(input-dos文件夹下) 打开vasprun.xml文件,请注意文件路径中不能包含...

    p4vasp

    能带&态密度图

    1. 画能带图:使用文件vasprun.xml(input-bands中的)
    2. 画态密度图:使用文件vasprun.xml(input-dos文件夹下)
    3. 打开vasprun.xml文件,请注意文件路径中不能包含英文
      打开vasprun.xml文件
      4.选择画图类型
      在这里插入图片描述
      在这里插入图片描述
      dos:态密度;band:能带
      用input-dos文件夹下的vasprun.xm没有高对称点。

    安装及更多操作参考:
    P4VASP的安装 & 使用P4VASP画DOS(态密度)图.

    OriginPro

    导出数据

    在这里插入图片描述
    在这里插入图片描述
    然后点export,文件将储存在p4vasp的文件夹下

    数据导出:左边能带的横坐标 k矢量或者说是k点距离 右边的是能量
    能量一般取-3~3 ev

    展开全文
  • 这样,得到一隐藏文件.cell, 将它用编辑器打开,从中得到vasp所需的POSCAR信息,修改得到POSCAR。 AgGaS2 bulk 1.000000000000000 -2.7934999465942410 2.7934999465942370 5.2045001983642580 2.7934

    第一步:构型优化

    1.准备四个输入文件

    POSCAR INCAR POTCAR KPOINTS

    POSCAR: 从ms中导入AgGaS2结构,选择CASTEP,file,save,并保存成原包。这样,得到一隐藏文件.cell, 将它用编辑器打开,从中得到vasp所需的POSCAR信息,修改得到POSCAR。

    AgGaS2 bulk

    1.000000000000000

    -2.7934999465942410 2.7934999465942370 5.2045001983642580

    2.7934999465942370 -2.7934999465942390 5.2045001983642590

    2.7934999465942380 2.7934999465942380 -5.2045001983642580

    4 2 2

    Direct

    0.3750000000000000 0.4072000086307526 0.5322000086307526

    0.8750000000000000 0.8427999913692474 0.4677999913692474

    . . . . . . . . . . . .

    2.POTCAR:用PBE-GGA的赝势,提取,Ag Ga O的赝势合并成一个赝势。(一般VASP有自带)

    3.下面是INCAR

    SYSTEM = optimization of AgGaS2

    LPLANE=.TRUE.

    NPAR= 8

    Elecronic minimisation

    ISTART = 0

    LREAL = .FALSE.

    PREC = Low

    EDIFF = 1e-4

    EDIFFG = -0.03

    IALGO = 48

    NELMIN = 4

    ISYM = 0

    GGA = PBE

    ISPIN = 1

    NBANDS = 120

    OUTPUT CONTROL

    LCHARG = .TRUE.

    LWAVE = .TRUE.

    LVTOT = .FALSE.

    IONIC RELAXATION

    NBLOCK = 1

    NSW = 1

    IBRION = -1

    DOS RELATED (disabled)

    ISMEAR = 0 (tetrahedron/gaussian/m-p)

    SIGMA = 0.05

    4.再下面KPOINTS

    A

    0

    M

    4 4 4

    0 0 0 至此,四文件已准备好,进行计算,

    mpiexec -np 8 ~/bin/vasp.4.5-mk-mp-pgi < /dev/null > vasp.out &(运行命令)

    (问题:计算前对某些参数的测试,比如截断能,晶格参数等,标准是什么?vasp运行后,不知道结果是否满足要求呢?)

    第二步:将构型优化后产生的CONTCAR文件拷贝为POSCAR文件

    mpiexec -np 8 ~/bin/vasp.4.5-mk-mp-pgi < /dev/null > vasp.out &(运行命令)

    第三步:能带的计算.建立新的计算目录

    mkdir band 创建目录

    cd band 进入目录

    cp ../INCAR . 复制

    cp ../PO* .

    cp ../CHG* .

    对INCAR文件作如下修改:

    ICHARG =11

    NBANDS=120 与上面构型优化时一致

    建立新的KPOINTS 手动定义K点

    AgGaS2

    51

    Rec(关键字要对齐第一行)

    0.5000 0.5000 0.5000 1.0000

    0.4500 0.4500 0.4500 1.0000

    . . . . . .

    K点坐标获得途径

    布里渊区各特殊点的坐标可利用CASTEP获得

    然后编辑

    vi inp

    lin1 特殊点

    lin2 特殊点之间的取点数

    lin3 第一个K点的坐标

    lin4 第二个K点的坐标

    . . . . (同上..)

    . . . .

    保存

    运行 vasp_kpoints<inp

    就可以得到坐标

    手动复制到KPOINTS

    运行VASP

    mpiexec -np 8 ~/bin/vasp.4.5-mk-mp-pgi < /dev/null > vasp.out &(运行命令)

    在运行

    Vaspband

    Draw the band structure of VASP

    Two files, EIGENVAL and OUTCAR are needed

    ----------------------------------------------------------------

    Spin-polarized calculation(T/F)? (是否自旋极化计算)

    f

    Number of lines along the BZ: (布里渊区特殊线数目)

    5

    K point sequence used to separate lines: (各线段始终点序号)

    1 11 21 31 41 51

    Set the Fermi level to zero(T/F)? (是否费米设置为零点)

    t

    可以得到相应的 band.dat 用origin绘图

    然后标注各个特殊点

    Fermi energy=1.2182v exp=2.76ev (可能是参数设计太粗糙…) 不过放大还算清楚..

    (用MS作图 K网格相同)

    第四步:DOS的计算及其绘图

    其POSCAR POTCAR KPOINTS可以相同于第一步优化

    设置INCAR

    SYSTEM = optimization of AgGaS2

    NPAR= 8

    LPLANE = .TRUE.

    Elecronic minimisation

    ISTART = 3 (从WAVECAR读入波函数)

    LREAL= .FALSE.

    PREC = Medium

    EDIFF = 1e-4

    EDIFFG = -0.02

    IALGO = 48

    NELMIN = 3

    ISYM = 1

    ISIF = 3

    ISPIN = 1

    ISMEAR = 0

    SIGMA = 0.1

    NBANDS=120

    LORBIT=11 (投影到 s.p.d)

    NEDOS=501 (DOS的点数)

    ICHARG=11 (从CHG读入电荷密度)

    EMIN=-20

    EMAX=20

    OUTPUT CONTROL

    LCHARG = .TRUE.

    LWAVE = .TRUE.

    LVTOT = .FALSE.

    IONIC RELAXATION

    NBLOCK = 1

    NSW = 1

    IBRION = -1

    POTIM = 0.05

    运行vasp 命令同上

    然后编写DOS

    F (是否自旋极化计算 )

    9 (轨道数)

    6 (要计算的DOS)

    -1 (为负数时,表示所考察的为原子的DOS)

    1 (原子编号)

    -1 (为正正数时,表示为原子轨道的DOS)

    2 (下面都同上)

    -1

    3

    1

    1

    1

    5

    1

    8

    保存 运行vaspdos (如果没有设定路径的话 。要自己指定路径)

    可以得到fort.25 文件

    运行cdos98 不过需要加一条DOS数

    就可以得到 DOS1 DOS2 DOS3 . . . . . .

    用origin分别绘图…..

    (可能是DOS编号出现错误了的原因 导致Ga(p)和S(p) 没有画出..)

    第五步:光学性质的计算

    (一)线性的计算

    1.也是进行构型优化(同能带计算时结构优化)..只是计算中INCAR文件需包含关键词:ISIF = 3

    所有计算均需采用PAW型赝势

    2.构型优化后将CONTCAR拷贝成POSCAR

    在INCAR里面

    NPAR=1 (计算光学性质时,必须按能带依次处理)

    LOPTICS = .TRUE. (计算光学性质)

    NSW=1

    IBRION= -1(无需进行构型优化)

    计算结束后 (正常), 会得到OPTIC(用于计算线性光学性质)和momentum_matrix

    (用于计算NLO即非线性光学性质)

    3.编辑OPTCRT文件 常用关键词及其含义

    ISYMM = 2

    OMMIN = 0

    OMMAX = 20

    NEDOS = 4000

    NBCON = 200

    LJDOS = .TRUE.

    LDOS = .TRUE.

    LKRAMERS = .TURE.

    GAMMA = 0.002

    LSEARCH =.TRUE.

    EMINSEARCH = 0

    EMAXSEARCH = 20

    SCISSOR=0.5

    4.创建一个新目录,将OPTIC OPTCTR,KPOINTS,POSCAR 复制到其下

    命令 cp ../ OPTIC . OPTCTR,KPOINTS,POSCAR .

    修改KPOINTS

    创建个临时目录复制vasp计算所需要的四个输入文件.运行单机版vasp。

    (由于56机子没有装单机版,所以这一步COPY到其他的机子上做的)

    单机运行命令

    ~/bin/vasp.4.5-mk-sp-pgi</dev/null>vasp.out&

    将其中的IBZKPT复制为KPOINTS

    运行vasp

    命令:mpiexec –np 8 ~/bin/vasp.4.5-mk-mp-pgi</dev/null>vasp.out&

    正常运行将产生2个文件

    EPS----存放介电函数数据

    JDOS----存放JDOS数据

    在运行vasp_lo

    用origin软件根据EPS数据绘图

    在此基础上才可以得到其他线性光学性质:

    (1)折射率(refractive index)

    (2)吸收系数(adsorption constant)

    (3)能量损失系数(energy-loss coefficient) :

    (4)消光系数(extinction coefficient)

    (5)双折射率曲线

    (6)光导电变化曲线

    (二)非线性光学性质

    1. 非线的运行需要以下三个文件

    KPOINTS、EIGENVAL和momentum_matrix三个文件,

    其中KPOINTS文件内容与计算线性光学性质相同.

    EIGENVAL文件为vasp计算输出文件

    运行vasp_nlo, 需输入以下数据:

    momentum_matrix 存放跃迁矩阵元的文件名

    177.32 单胞体积

    2.76 实验带隙(若为0则无需对带隙进行校正)

    1 起始能带(一般均为1)

    38 终止能带

    F 为False时表示只计算静态倍频系数

    2.运行后会得到静态倍频系数的计算结果.(此数据好像没有自己保存)

    d( 1, 1) SHG (real part)= 1.60120 and (imag part) 0.00000(pm/V)

    d( 1, 2) SHG (real part)= -2.37106 and (imag part) 0.00000(pm/V)

    d( 1, 3) SHG (real part)= -2.72264 and (imag part) 0.00000(pm/V)

    d( 1, 4) SHG (real part)= 10.98424 and (imag part) 0.00000(pm/V)

    d( 1, 5) SHG (real part)= -0.00510 and (imag part) 0.00000(pm/V)

    d( 1, 6) SHG (real part)= -2.46438 and (imag part) 0.00000(pm/V)

    d( 2, 1) SHG (real part)= -2.46438 and (imag part) 0.00000(pm/V)

    d( 2, 2) SHG (real part)= 0.34388 and (imag part) 0.00000(pm/V)

    d( 2, 3) SHG (real part)= -1.40097 and (imag part) 0.00000(pm/V)

    d( 2, 4) SHG (real part)= -0.42303 and (imag part) 0.00000(pm/V)

    d( 2, 5) SHG (real part)= 10.98424 and (imag part) 0.00000(pm/V)

    d( 2, 6) SHG (real part)= -2.37106 and (imag part) 0.00000(pm/V)

    d( 3, 1) SHG (real part)= -0.00510 and (imag part) 0.00000(pm/V)

    d( 3, 2) SHG (real part)= -0.42303 and (imag part) 0.00000(pm/V)

    d( 3, 3) SHG (real part)= 1.21340 and (imag part) 0.00000(pm/V)

    d( 3, 4) SHG (real part)= -1.40097 and (imag part) 0.00000(pm/V)

    d( 3, 5) SHG (real part)= -2.72264 and (imag part) 0.00000(pm/V)

    d( 3, 6) SHG (real part)= 10.98424 and (imag part) 0.00000(pm/V)

    由上可得最大倍频系数为d(1.4)

    动态倍频系数计算 运行vasp

    需要输入以下数据

    momentum_matrix 存放跃迁矩阵元的文件名

    42.34 单胞体积

    2.26 实验带隙

    1 起始能带

    26 终止能带

    T 为True时表示要计算动态倍频系数

    10 考察动态倍频系数时的能量范围

    500 格点数

    1 4 倍频系数分量d14

    最后,有相关需求,欢迎通过微信公众号联系我们。

    微信公众号:320科技工作室。

    展开全文
  • 今天主要记录强关联电子体系+U、DOS和能带计算过程中参数的设置,相信大家在读一些计算类文献时会看到文章中计算部分有The GGA + U calculations were performed using the model proposed by Dudarev et al这样一段...
  • 本文介绍用VASP如何求硅的电子态密度和能带,共分为如下5个部分: (1). 生成4个输入文件: POSCAR POTCAR INCAR KPOINTS (2). 优化晶格参数,求出能量最低所对应的晶格参数 (3). 固定晶格参数, 求出能态密度(DOSCAR...
  • P4VASP的安装 & 使用P4VASP画DOS(态密度)图

    千次阅读 多人点赞 2020-12-29 16:40:47
    P4vasp 下载: ...这里以vasp计算为例 首先导出DOS计算的文件夹中的vasprun.xml文件 从p4vasp导入该文件进行作图(直接将vasprun.xml拖入p4vasp一般不会有反应) 具体步骤: 点击open: 在弹出的弹窗
  • 固体与表面-从零学习vasp计算(1)

    千次阅读 2018-10-12 20:52:36
    固体与表面-从零学习vasp计算1固体与表面-从零学习vasp计算(1)理论催化计算实战教学目的和理念课程用到的程序课程准备工作: 固体与表面-从零学习vasp计算(1) 今天开一个坑,做一系列的vasp实战教学课程,同时在b站...
  • vasp计算能带轨道成分

    千次阅读 2020-01-19 01:13:36
    vasp计算需要使用课题组服务器做计算,涉及到SSH连接登陆问题。 1. 对于Windows用户,使用Xshell + WinSCP 解决命令行操作,以及文件传输 2. Mac用户直接使用Item终端连接操作。 3. 使用vscode 的remote加上...
  • 固体与表面-从零学习vasp计算4固体与表面-从零学习vasp计算(4)Setting up a VASP calculationINCARINCAR的注意事项:关键词讲解: SYSTEM关键词讲解: ISTART = 0 | 1 | 2 | 3关键词讲解: ICHARG = 0 | 1 | 2 | 4 |...
  • pymatgen生成能带和态密度数据

    千次阅读 2020-08-11 20:20:56
    VASP计算好能带和态密度后,往往需要拖到windows中利用p4vasp来画图,或者用其导出能带和态密度图,这样做实在太麻烦了。这里介绍用Pymatgen程序生成能带和态密度的矩阵,直接放到origin中。只需要一个python程序...
  • 利用VASP计算能带结构,这是英文版的,是有兴趣可以翻译成中文的。
  • 利用vasp工具计算石墨烯的能带结构,态密度,以及自选密度分布和光学性质。这篇笔记是自己在写毕业论文时跟着师哥师姐在学校集群上做计算时整理的笔记,对vasp计算的原理不太了解,只会这个简单的过程,里面有错误之...
  • VASP计算AgGaS2能带及态密度及光学性质全过程
  • VASP表面计算步骤小结

    千次阅读 2019-10-02 05:18:35
    一、概述 vasp用“slab”模型来模拟表面...vasp计算表面的大概步骤是:材料体性质的计算;表面模型的构造;表面结构的优化;表面性质的计算。二、分步介绍1、材料体性质计算:本步是为了确定表面计算时所需的一些...
  • 态密度图的基本分析方法.xps
  • pymatgen绘制能带态密度图(band+DOS)

    千次阅读 热门讨论 2020-05-24 21:26:41
    pymatgen是一个API包,可以连接materials project进行高通量计算。现在我用pymatgen画一下能带+DOS图,炒鸡好看。 安装请见https://pypi.org/project/pymatgen/ 前提是需要有vasprun.xml,这是VASP运行时就会生成...
  • vasp可以做哪些计算及命令

    千次阅读 2020-05-06 14:58:14
    电子性质:电荷密度分布、电子态密度、能带结构 光学性质:折射率、消光系数、反射率、吸收系数 力学性质:弹性常数、体积模量、泊松比、切变模量 磁学性质:顺磁、铁磁、反铁磁 密度泛函理论、Linux...
  • Vasp二维材料单胞基本计算范例

    千次阅读 2019-12-10 17:23:52
    3.能带和态密度计算 上面几个过程,实际需要改动的参数只有几个,比如cut energy,k点大小,还有费米展宽参数,所有严格来说需要做批量测试。 对于简单二维体系来说,结构优化不是特别重要,可以使用实验结果,...
  • vasp-电荷局域密度图(ELF)

    千次阅读 2020-07-29 15:19:09
    运行vasp得到ELFCAR文件。 将ELFCAR拖入Vesta软件中,Vesta软件下载地址为:http://www.jp-minerals.org/vesta/en/download.html 点击 Utiltties - 2D Data Display. 点击Slice,确定切面和fractional coordinates...
  • pymatgen画能带投影的态密度

    千次阅读 2019-11-20 02:02:55
    import matplotlib.pyplot as plt from pymatgen.io.vasp.outputs import Vasprun from pymatgen.electronic_structure.plotter import BSDOSPlotter,BSPlotter,BSPlotterProjected,DosPlotter bs_vasprun = Va...
  • Vasp 石墨烯能带计算

    千次阅读 2019-10-16 14:18:30
    首先做自洽计算产生正确的基态电子密度 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.00.....
  • 1.VASP简介 VASP是一个基于量子力学的第一性原理软件包,内容包含VASP的代码包和赝势,需要在Linux下安装和运行...VASP主要是做结构优化,电子性质(态密度,能带结构,电荷密度,光学性质等),第一性原理分子动力学。
  • 背景介绍:(Spanish Initiative for Electronic Simulations with Thousands ofAtoms)是一种实现电子结构计算和第一性原理分子动力学模拟的程序,同时也是一种实现的方法。软件主要由剑桥大学地球科学系Emilio教授、...
  • vasp-入门帖.rar

    2009-07-18 11:45:05
    1、vasp软件包使用入门指南 2、电子态密度的各种计算方式 3、Calculate Band Structure Using VASP 4、vasp中文说明书.pdf 5、处理vasp结果的程序 6、vasp应用实例文件(run就ok)
  • VASP 结构优化、静态自洽、非自洽计算

    万次阅读 多人点赞 2017-03-01 17:29:02
    结构优化是指对输入体系进行原子弛豫,得到一个稳定的体系结构(通常以 前后两次总自由能之差 或 原子所受最大的力作为原子的迟豫收敛标准),当然结构优化包含了电子自洽计算过程(电子自洽迭代计算以总能作为收敛...
  • VASP赝势简介

    2021-09-13 16:24:26
    vasp计算中用到的三种赝势:模守恒赝势,超软赝势,PAW赝势(按产生顺序)。 按方法不同分为USPP(ultrasoft pesudopotential,超软赝势)和PAW(projector augmented wave,投影缀加平面波),两种方法都可以相当程度...
  • VASP ----

    热门讨论 2008-05-09 23:00:00
    VASP是使用赝势和平面波基组,进行第一定律分子动力学计算的... 设计定位,旋转和局域波的态密度; 费米水平的部分态密度; 电荷密度旋转密度; 轨道和轨道密度。 操作平台:Linux/Unix

空空如也

空空如也

1 2 3 4 5
收藏数 91
精华内容 36
关键字:

vasp计算态密度