vasp算表面吸附流程_VASP表面吸附计算实例分析

本文主要为了记录在学习Vasp计算过程中如何设置INCAR中的一些参数,并不着重讲解其含义,详见可自行到vasp官网查阅。

其次关于建模部分也不做细致讨论,一般结构可从Materials Studio里自带的晶体数据库导出,或从三大数据库网站上下载:

1、Materaials Project网址:https://materialsproject.org/(强烈推荐,邮箱注册即可使用)

2、CCDC网址:https://www.ccdc.cam.ac.uk/

3、ISCD:http://www2.fiz-karlsruhe.de/icsd_home.html

一、表面吸附计算

这里以W(100)表面吸附CO分子为例,暂不做收敛性测试,主要为了记录一些参数的设定,计算流程如下:

1.对W晶胞进行结构优化

INCAR设置:

#### initial I/O ####SYSTEM = WISTART = 0ICHARG = 2LWAVE = .FALSE. #优化晶胞时,不开波函数,节省计算时间LCHARG = .FALSE. #优化晶胞时,不开电荷密度,节省计算时间#### Ele Relax 电子步####ENCUT = 225 #截断动能根据相关文献或赝势确定,也可自己做收敛性测试ISMEAR = 1 #金属体系SIGMA = 0.2 #POTCAR值确定EDIFF = 0.1E-4NELM = 60LREAL = .F.PREC = Normal#### Geo opt 离子步####EDIFFG = -0.05IBRION = 2POTIM = 0.2NSW =

  • 0
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
VASP是一种常用的第一性原理计算软件,可用于模拟各种材料的物理和化学性质,包括表面吸附。下面是利用VASP进行表面吸附计算的一般步骤和代码示例: 1. 准备模型 首先需要准备表面模型,在VASP中通常采用超胞模型,即将一个小的表面单元重复堆叠成一个大的超胞。可以使用VASP的POSCAR文件格式来描述超胞的结构,其中包括晶格常数、原子种类和坐标等信息。 2. 初步松弛 对于表面吸附计算,需要先进行初步的结构松弛,以确定表面的稳定结构。可以使用VASP的INCAR文件设置参数,如ISIF=2表示同时对晶格和原子位置进行优化,IBRION=2表示采用自洽赝势方法进行松弛。以下是一个简单的INCAR文件示例: ``` SYSTEM = surface_adsorption ISTART = 0 ICHARG = 2 ENCUT = 400 EDIFF = 1E-6 ISMEAR = 0 ISIF = 2 IBRION = 2 NSW = 100 POTIM = 0.2 ``` 3. 吸附分子构型生成 接下来需要生成吸附分子在表面上的构型。可以使用VASP的POSCAR文件描述吸附分子在表面上的位置和朝向等信息。通常需要对吸附分子进行初步的几何优化,以确定其在表面上的最稳定位置。 4. 吸附计算 利用VASP计算结果,可以得到吸附分子在表面上的能量。吸附能可以通过下式计算得到: ``` E_ads = E_surf+mol - E_surf - E_mol ``` 其中E_surf+mol是吸附后的表面-分子复合体的能量,E_surf和E_mol分别是表面和分子的能量。通过计算吸附能可以得到吸附过程的稳定性和动力学信息。可以使用VASP的OUTCAR文件获取表面和复合体的能量信息,然后计算吸附能。 以上只是表面吸附计算的大致步骤和代码示例,具体的计算方法和参数设置需要根据具体的模型和吸附分子进行调整。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值