当前位置：文档之家 › VASP上机练习讲解几何结构优化分子

VASP上机练习讲解几何结构优化分子

  • 格式：pdf
  • 大小：1.10 MB
  • 文档页数：25

下载文档原格式

  / 25

合集下载

2015源资培训班-VASP上机练习讲解(MedeA-VASP相关模块MT-Phonon-Electronics-TSS)

2015源资培训班-VASP上机练习讲解(MedeA-VASP相关模块MT-Phonon-Electronics-TSS)

How to analyze?
Frequency: 6.904 THz IR
Phonon Dispersion
Right-click
How to analyze?
Frequency: 6.904 THz
Z
Move
up
Oxygen atoms move along Z direction
How to analyze?
VASP 量子力学
计算模块
LAMMPS 分子动力学
GIBBS 蒙特卡洛 MOPAC 半经验量化
Transition State Search 过渡态搜索 Mechanical Thermal 力学、热力学性质 Phonon 声子谱、光谱预测 Electronics 费米面、电子输运性质
Thermal Conductivity 热导
粗略搜索过渡态 的image个数
对插入image 结构的设置
鞍点处进行精细搜索的次数（同粗搜 索采用相同image个数）
结合CI-NEB方法，找到在设定反应 路径上更为准确的过渡态
化学反应
Step 2 和 step 3 是可选项： 在Step 2中，推荐采用Dimer方法 去更精确的找准反应路径中鞍点的位置以及准 确能量，并对过渡态结构做优化处理。 Step 3 是采用RMM-DIIS方法，再优化前两步搜索的鞍点结构。此方法的收敛半 径非常小，所以要求之前的鞍点结构必须十分接近真正过渡态的结构。
化学反应
过渡金属表面上的CO解离
确定反应初始态结构 (initial state)
创建表面吸附模型 MedeA-VASP模块 MedeA-TSS模块
确定反应末态结构 (final state)

2015源资培训班-VASP上机练习讲解(几何结构优化-分子)

2015源资培训班-VASP上机练习讲解(几何结构优化-分子)
VASP 上机练习讲解
几何结构优化（分子）
1. H2O 弛豫结构
2
1. H2O分子
MedeA创建模型：打开昨天创建好的结构
File Open File from Disk 选择保存结构目录中的H2O.sci 打开结构
3
1. H2O分子
创建输入文件：INCAR、POSCAR、KPOINTS、POTCAR
1）POSCAR：
H2O
！体系名称
1.0
！缩放系数
10.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 10.0000000000000000 0.0000000000000000 ！晶胞矢量
0.0000000000000000 0.0000000000000000 10.0000000000000000
！截断能大小 ！电子弛豫最大步数 ！能量收敛标准 ！不保留CHGCAR ！不保留WAVECAR ！积分路径 ！积分展宽 ！电子优化的算法
注：由于此次培训班输入文件以H2O结构优化为模板，因此此算例中INCAR不需改 动。在之后的算例中，其他参数也应随计算体系修改。
8
1. H2O分子
3）KPOINTS： 切换至Xshell窗口，编辑KPOINTS vi KPOINTS 不需修改，直接退出 :q
13
1. H2O分子
计算结束后，查看主要输出文件：OUTCAR、OSZICAR、 CONTCAR、XDATCAR
vi XXXXX
2）OSZICAR vi OSZICAR Ctrl+b 向前翻页 Ctrl+f 向后翻页 如右图显示每步 电子弛豫过后的 体系总能量

2015源资培训班-VASP上机练习讲解(杂化泛函)

2015源资培训班-VASP上机练习讲解(杂化泛函)

TIME = 0.4
*杂化计算时，利于收敛
8
杂化泛函计算步骤
态密度、能带结构建议：
1. 结构优化（PBE） 2. 静态自洽scf计算（PBE）
*保留波函数和电荷密度
3. 静态自洽scf计算（HSE06）
*采用第2步得到的波函数，加快收敛速度 *保留波函数和电荷密度
4. 非自洽计算（HSE06）
*采用第3步得到的波函数和电荷密度
ALGO = Damped LHFCALC = .TRUE. HFSCREEN = 0.2 TIME = 0.4 LMAXFOCK = 4
创建HSE-scf文件夹： cp -r Si-HSE-scf Si-HSE-dos
修改输入文件： INCAR（如右修改） KPOINTS（不变） POTCAR（不变） POSCAR（不变）
NBANDS = 10 LORBIT = 10 EMIN = -10 EMAX = 25 NEDOS = 1000
12
结果分析
PBE： 0.50 eV
HSE06： 1.20eV
13
3
常用杂化泛函
B3LYP
适用于分子体系 The standard hybrid functional for molecules after Becke, Lee, Yang & Parr47, Vosko, Wilk & Nusair, Stephens, Devlin, Chabalowski & Frisch
5
半导体材料GGA-PBE与 HSE06泛函计算对比研究
6
半导体
• 禁带宽度是半导体十分重要的电子性质。 • 众所周知，普通的LDA和GGA方法计算带隙存在问题，一 般低估1-2个电子伏特。解决的办法则是采用杂化泛函HSE 给出更为准确的禁带宽度。 PBE Si ZnS TiO2(anatase) 0.63eV 1.90eV 1.80eV HSE06 1.20eV 3.28eV 3.20eV 实验值 1.12eV 3.68eV 3.20eV

VASP上机练习讲解-杂化泛函

VASP上机练习讲解-杂化泛函
8
杂化泛函计算步骤
态密度、能带结构建议：
1. 结构优化（PBE） 2. 静态自洽scf计算（PBE）
*保留波函数和电荷密度
3. 静态自洽scf计算（HSE06）
*采用第2步得到的波函数，加快收敛速度 *保留波函数和电荷密度
4. 非自洽计算（HSE06）
*采用第3步得到的波函数和电荷密度
9
算例——Si
Si 应用领域：
太阳能电池 电子元器件 半导体 。。。。。。
禁带宽度：
实验值为1.20eV
10

算例——Si
3. 静态自洽scf（HSE06）
创建HSE-scf文件夹： cp -r Si-scf Si-HSE-scf 修改输入文件： INCAR（如右修改） KPOINTS（不变） POTCAR（不变） POSCAR（不变）
SYSTEM = Si ISTART = 1 ICHARG = 0 PREC=Normal LREAL = .F. IBRION = -1 ISIF=2 NSW = 0 POTIM = 0.5 EDIFFG =-0.05 ENCUT = 320 eV NELM = 60 EDIFF = 0.1E-04 LCHARG = .T. LWAVE = .T. ISMEAR = -5 SIGMA = 0.2 ALGO = Damped LHFCALC = .TRUE. HFSCREEN = 0.2 TIME = 0.4 LMAXFOCK = 4 PRECFOCK = Normal
3
常用杂化泛函
B3LYP
适用于分子体系 The standard hybrid functional for molecules after Becke, Lee, Yang & Parr47, Vosko, Wilk & Nusair, Stephens, Devlin, Chabalowski & Frisch

2015源资培训班-VASP上机练习讲解(MedeA-LAMMPS-GIBBS-MOPAC模块的使用)

2015源资培训班-VASP上机练习讲解(MedeA-LAMMPS-GIBBS-MOPAC模块的使用)

GIBBS
MedeA®——蒙特卡洛计算方法
超临界CO2流体中声速和密度计算
创建分子模型 MedeA-GIBBS
通过Monte Carlo计算能够获得的性质
Calculation of first order derivative of the thermodynamic potential (i.e. pressure, molar volume, enthalpy)
•OPLS-AA •PCFF •PCFF+ •COMPASS •CFF91 & CFF93 •CVFF •AUA-4 （仅用于Gibbs）
MedeA® LAMMPS-EAM
嵌入原子法（Embedded Atom Method, EAM）力场模拟是一种 描述金属体系结构、力学性质、 热性质的有效方法。 支持计算多种金属体系的性质： 结构、缺陷的结构和能量、力学 性质及动力学性质（熔点等）。
t-mer
-357.90 1.45 8.29
MedeA®——半经验量化计算方法
UV-vis光谱预测
MedeA®——创建热固性材料模型
关键优势： •完全集成到MedeA® Flowcharts中，并且能够直 接与以下模型结合： 与MedeA® Amorphous Builder结合创建未 经处理的原料结构 与MedeA® LAMMPS结合生成动力学轨迹 与其他性质计算模块结合，使用精确的力 场来预测密度（shrinkage），力学性质和热 导率 •在交联过程中，对材料特性进行严格检验： 网状结构和凝胶点的演变 粘结应力 自动缺陷检测（避免非物质的环连锁） •研究有/没有外加溶剂的情况下，多种树脂和固 化剂的结合
MedeA Molecular Buiห้องสมุดไป่ตู้der MedeA Polymer Builder

VASP经典学习教程-有用

VASP经典学习教程-有用

VASP 学习教程太原理工大学量子化学课题组2012/5/25目录第一章Linux命令 (1)1.1 常用命令 (1)1.1.1 浏览目录 (1)1.1.2 浏览文件 (1)1.1.3 目录操作 (1)1.1.4 文件操作 (1)1.1.5 系统信息 (1)第二章SSH软件使用 (2)2.1 软件界面 (2)2.2 SSH transfer的应用 (3)2.2.1 文件传输 (3)2.2.2 简单应用 (3)第三章VASP的四个输入文件 (3)3.1 INCAR (3)3.2 KPOINTS (4)3.3 POSCAR (4)3.4 POTCAR (5)第四章实例 (5)4.1 模型的构建 (5)4.2 VASP计算 (8)4.2.1 参数测试 (8)4.2.2 晶胞优化(Cu) (13)4.2.3 Cu(100)表面的能量 (2)4.2.4 吸附分子CO、H、CHO的结构优化 (2)4.2.5 CO吸附于Cu100表面H位 (4)4.2.6 H吸附于Cu100表面H位 (5)4.2.7 CHO吸附于Cu100表面B位 (6)4.2.8 CO和H共吸附于Cu100表面 (7)4.2.9 过渡态计算 (8)第一章Linux命令1.1 常用命令1.1.1 浏览目录cd: 进入某个目录。

如：cd /home/songluzhi/vasp/CH4 cd .. 上一层目录；cd / 根目录；ls: 显示目录下的文件。

注：输入目录名时，可只输入前3个字母，按Tab键补全。

1.1.2 浏览文件cat：显示文件内容。

如：cat INCAR如果文件较大，可用：cat INCAR | more (可以按上下键查看) 合并文件：cat A B > C (A和B的内容合并，A在前，B在后) 1.1.3 目录操作mkdir：建立目录；rmdir：删除目录。

如：mkdir T-CH3-Rh1111.1.4 文件操作rm：删除文件；vi：编辑文件；cp：拷贝文件mv：移动文件；pwd：显示当前路径。

vasp做分子动力学

vasp做分子动力学的好处，由于vasp是近些年开发的比较成熟的软件，在做电子scf速度方面有较好的优势。

缺点：可选系综太少。

尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。

主要使用的系综是NVT和NVE。

下面我将对主要参数进行介绍！一般做分子动力学的时候都需要较多原子，一般都超过100个。

当原子数多的时候，k点实际就需要较少了。

有的时候用一个k点就行，不过这都需要严格的测试。

通常超过200个原子的时候，用一个k点，即Gamma点就可以了。

INCAR：EDIFF 一般来说，用1E-4或者1E-5都可以，这个参数只是对第一个离子步的自洽影响大一些，对于长时间的分子动力学的模拟，精度小一点也无所谓，但不能太小。

IBRION=0 分子动力学模拟IALGO=48 一般用48，对于原子数较多，这个优化方式较好。

NSW=1000 多少个时间步长。

POTIM=3 时间步长,单位fs,通常1到3.ISIF=2 计算外界的压力.NBLOCK= 1 多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT. KBLOCK=50 NBLOCK*KBLOCK个步长写一次XDATCAR.ISMEAR=-1 费米迪拉克分布.SIGMA =0.05 单位:电子伏NELMIN=8 一般用6到8,最小的电子scf数.太少的话,收敛的不好.LREAL=AAPACO=10 径向分布函数距离,单位是埃.NPACO=200 径向分布函数插的点数.LCHARG=F 尽量不写电荷密度,否则CHG文件太大.TEBEG=300 初始温度.TEEND=300 终态温度。

不设的话，等于TEBEG.SMASS -3 NVE ensemble；-1 用来做模拟退火；大于0 NVT 系综。

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// 1)收敛判据的选择结构弛豫的判据一般有两种选择：能量和力。

2015源资培训班-VASP上机练习讲解(几何结构优化-表面、吸附)

Selective dynamics Direct
0.00000000 0.00000000 0.09914000 F F F 0.00000000 0.50000000 0.09914000 F F F 0.50000000 0.00000000 0.09914000 F F F 0.50000000 0.50000000 0.09914000 F F F 0.16667000 0.33333000 0.19828000 F F F 0.16667000 0.83333000 0.19828000 F F F 0.66667000 0.33333000 0.19828000 F F F 0.66667000 0.83333000 0.19828000 F F F -0.00020533 -0.00004234 0.38489174 T T T 0.00001920 0.50001702 0.38661440 T T T 0.50011017 -0.00005160 0.38497381 T T T 0.50008950 0.50022203 0.38491575 T T T 0.33338701 0.16667413 0.28884565 T T T 0.33079186 0.66543798 0.28974725 T T T 0.83459743 0.16928665 0.28973673 T T T 0.83468151 0.66531165 0.28971833 T T T
输入文件INCAR
！优化表面
！同Pd晶胞优化采用相同设置
提示： ISIF参数：晶胞优化取3，表面优化取2 ENCUT参数：由体系决定，可参考文献及 POTCAR中的值确定。
12
（3）优化Pd(111)表面

2015源资培训班-VASP上机练习讲解(几何结构优化-表面、吸附)


输ห้องสมุดไป่ตู้文件INCAR
！优化表面
！同Pd晶胞优化采用相同设置
提示： ISIF参数：晶胞优化取3，表面优化取2 ENCUT参数：由体系决定，可参考文献及 POTCAR中的值确定。
12
（3）优化Pd(111)表面
Pd111 1.0 5.63410000 0.00000000 0.00000000 -2.81704988 4.87927380 0.00000000 0.00000043 0.00000074 23.20050000
17
（6）Pd(111)吸附CO分子计算 Pd(111)-CO
SYSTEM = Pd111-CO-fcc ISTART = 0 ICHARG = 2 PREC=Normal LREAL = .F. IBRION = 2 ISIF=2 NSW = 100 POTIM = 0.5 EDIFFG =-0.05 对于在表面吸附分子，INCAR参数不变。
3
（1）对Pd晶胞进行结构优化
MedeA平台上建模：
晶格常数： a = b = c =3.86700 α=β=γ=90
（1）对Pd晶胞进行结构优化
SYSTEM = Pd ISTART = 0 ICHARG = 2 PREC=Normal LREAL = .F. IBRION = 2 ISIF=3 NSW = 100 POTIM = 0.5 EDIFFG = -0.05 ENCUT = 250 eV NELM = 60 EDIFF = 0.1E-04 LCHARG = .F. LWAVE = .F. ISMEAR = 1 SIGMA = 0.2
优化晶胞常数
！需要改变Pd的晶胞体积和晶格参数
！截断能根据相关文献或贋势确定，也可自己测试

2015源资培训班-VASP上机练习讲解 (MedeA-VASP模块)

– Clean all results – Start from fresh Input
• Restart
– Use existing results – Get files from completed tasks – Redo missing tasks
6
TaskServer 计算任务
7
JobServer vs TaskServer
• File Preferences
9
MedeA的基本设置
10
MedeA的基本设置
11
MedeA的基本设置
12
MedeA-VASP参数设置
13
Load VASP module
• Tools → VASP 5.3 → Run
14
Calculation——重要参数
计算类型
计算性质
15
Calculation——重要参数
– Step in calculating some property – Limited logic:
• STOP when X
– Generates input set for one or more code – Combines results
– input set for a code
8
MedeA的基本设置
通过VASP convergence功 能够自动确定适用于目标体 系计算的Cutoff和Kmeshing两个参数。
对两个参数进行不同取值， 对目标体系反复计算，分析 体系总能量变化趋势，当趋 于恒定不变时，即为合适的 参数取值。
36
Automated Convergence
37
Automated Convergence
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

0.44329 0.58837 0.50000
注：由于此次培训班输入文件以H2O结构优化为模板，因此此算例中只对原子坐标进 行修改，其他不需改动。 在之后的算例中，其他参数也应随计算体系修改。
7
1. H2O分子
2）INCAR： 切换至Xshell窗口，查看INCAR cat INCAR 不需修改
注：由于此次培训班输入文件以H2O结构优化为模板，因此此算例中INCAR不需改 动。在之后的算例中，其他参数也应随计算体系修改。
8
1. H2O分子
3）KPOINTS： 切换至Xshell窗口，查看KPOINTS cat KPOINTS 不需修改
mesh auto 0 G 111 000
提示：对于10x10x10Å大小的晶胞，使用Gamma点计算即可。
将模板文件夹bak直接拷贝成H2O文件夹
cp -r bak H2O
进入H2O文件夹 cd H2O
编辑POSCAR vi POSCAR
将其中模型坐标部分删掉 dd/d]]
进入插入模式 i
右键粘贴，将UltraEdit中复制的
坐标信息粘贴替换模板中的坐标部分
退出插入模式 Esc 保存POSCAR :wq
ENCUT = 400 eV NELM = 60 EDIFF = 0.1E-04 LCHARG = .F. LWAVE = .F. ISMEAR = 0 SIGMA = 0.2 ALGO = Fast
！截断能大小 ！电子弛豫最大步数 ！能量收敛标准 ！不保留CHGCAR ！不保留WAVECAR ！积分路径 ！积分展宽 ！电子优化的算法
5）投vasp作业： qsub vasp.pbs
提示：POTCAR中贋势排列顺序必须与POSCAR中原子种类及位 置顺序一致。
10
1. H2O分子
实时查看计算过程：OUTCAR、OSZICAR tail -f OUTCAR
退出此模式：ctrl + c
11
1. H2O分子
计算结束后，查看主要输出文件：OUTCAR、OSZICAR、 CONTCAR、XDATCAR
14
1. H2O分子
计算结束后，查看主要输出文件：OUTCAR、OSZICAR、 CONTCAR、XDATCAR
vi XXXXX
3）CONTCAR vi CONTCAR 如右图显示最后一步 离子弛豫过后的体系 晶格参数及原子坐标， 即最终优化好的结构。
同POSCAR的格式完全一致。 如果在计算时，结构没有收敛或者结构优化之后要接着进行scf计算，则可以将作 业结束后生成的CONTCAR替代新作业中的POSCAR（cp CONTCAR POSCAR）。 在第二次投作业时，由上一步的更新的结构继续计算。
16
1. H2O分子
记录计算结果，保存计算文件
能量：建议采用本子记下来，注明体系及重要参数。 文件备份：采用Xmanager软件的xftp传输工具将H2O文件夹由
服务器传回至本地。 查看结构：采用MedeA直接打开CONTCAR文件 其他输出文件：对于专门的性质计算，可能还需要保留其他
CHG、WAVECAR、CHGCAR等输出文件。如果单纯的结构优 化，可以不保留。
此时可以直接将H2O文件夹复制并重命名为O2文件夹，然后在O2 文件夹中将相关输入文件修改，投作业后输出文件将会覆盖之前 H2O的输出文件（cp –r H2O O2）
3. 修改POTCAR（cp PAW/O/POTCAR O2 ，只含有O的贋势） 4. 修改POSCAR、INCAR、KPOINTS（vi）
17
1. H2O分子
分析初始结构和优化后结构
初始结构： H-O：1.05Å
优化后的结构： H-O：0.9725Å
对于优化前后的结构对比，可以观察键长、键角、 晶格参数、吸附位等
18
练习1 O2 分子
19
练习1 O2 分子
MedeA中打开本地O2分子，并导出为cif文件
20
练习1 O2 分子
1. 创建O2文件夹（mkdir） 2. 准备四个输入文件和脚本（cp）
坐标部分删除并替换
注：由于此次培训班输入文件以H2O结构优化为模板，因此此算例中其他不需改动。 6
1. H2O分子
1）POSCAR：
H2O
！体系名称
1.0
！缩放系数
10.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 10.0000000000000000 0.0000000000000000 ！晶胞矢量
1）POSCAR： File Export 文件名：H2O 保存类型：cif 保存
4
1. H2O分子
1）POSCAR： 用UltraEdit打开H2O.cif文件 只复制坐标部分（采用列选择模式）
可以实现列选 择编辑模式
5
1. H2O分子
1）POSCAR：
切换至Xshell界面，登陆 ，进入vasp2015文件夹 cd vasp2015
15
1. H2O分子
计算结束后，查看主要输出文件：OUTCAR、OSZICAR、 CONTCAR、XDATCAR
vi XXXXX
4）XDATCAR vi XDATCAR 如右图显示每一步离 子弛豫过后的原子坐 标，直至最后一步。 可用于当结构优化不 收敛或硬件原因中途 终止计算时，取中间 某一步结构继续计算。
0.0000000000000000 0.0000000000000000 10.0000000000000000
12
！各种原子的个数 （一个O原子，两个H原子）
Direct
0.50000 0.50000 源自.500000.45183 0.40670 0.50000 ！原子坐标 ， 第一行为O原子，后两行为H原子
SYSTEM = H2O ISTART = 0 ICHARG = 2 PREC=Normal LREAL = .F. IBRION = 2 ISIF=2 NSW = 100 POTIM = 0.5 EDIFFG =-0.05
！体系的命名 ！ 新作业 ！保留原子的电荷密度 ！计算精度 ！采用倒易空间 ！结构优化方法:CG ! 只有原子relax ！离子弛豫最大步数 ！离子弛豫步长 ！力的收敛标准
注：由于此次培训班输入文件以H2O结构优化为模板，因此此算例中KPOINTS不需 改动。在之后的算例中，其他参数也应随计算体系修改。
9
1. H2O分子
4）POTCAR： 切换至vasp2015文件夹 cd .. 将O和H的贋势文件连接成一个文件，命名为POTCAR cat PAW/O/POTCAR PAW/H/POTCAR > H2O/POTCAR 进入H2O文件夹 cd H2O 查看文件夹中的文件是否完全 ll
13
1. H2O分子
计算结束后，查看主要输出文件：OUTCAR、OSZICAR、 CONTCAR、XDATCAR
vi XXXXX
2）OSZICAR vi OSZICAR Ctrl+b 向前翻页 Ctrl+f 向后翻页 如右图显示每步 电子弛豫过后的 体系总能量
在作业进行的过程中，可通过此文件查看每步电子弛豫后的能量变化（dE）是 否达到收敛精度（即EDIFF）。
24
练习1 O2 分子
25
vi XXXXX
1）OUTCAR vi OUTCAR ]] 到OUTCAR最末端 Ctrl+f 向后翻页 Ctrl+b 向前翻页 若作业正常结束，则 显示右图。 [[ 到最前端
12
1. H2O分子
计算结束后，查看主要输出文件：OUTCAR、OSZICAR、 CONTCAR、XDATCAR vi XXXXX
cd H2O POSCAR根据cif文件修改 INCAR文件需要填加二行
ISPIN=2
VOSKOWN=1
KPOINTS不需修改，晶胞尺寸相同，k点不变
4. 投作业（qsub）
21
练习1 O2 分子
INCAR
SYSTEM = O2 ISTART = 0 ICHARG = 2 PREC=Normal LREAL = .F. IBRION = 2 ISIF=2 NSW = 100 POTIM = 0.5 EDIFFG =-0.05
VASP 上机练习讲解
几何结构优化（分子）
1. H2O 弛豫结构
2
1. H2O分子
MedeA创建模型：打开昨天创建好的结构
File Open File from Disk 选择保存结构目录中的H2O.sci 打开结构
3
1. H2O分子
创建输入文件：INCAR、POSCAR、KPOINTS、POTCAR
23
练习1 O2 分子
5. 实时跟踪输出文件 tail –f OSZICAR 6. 保存重要输出文件 xftp回传INCAR、KPOINTS、POSCAR、CONTCAR、 OUTCAR至本地 7. 记录计算结果 OUTCAR中记录O2分子总能量 E= -9.847 eV 8. 对比优化前后结构 MedeA打开初始结构和CONTCAR对比 初始结构：O-O键长 1.46Å 优化之后：O-O键长 1.2343Å
ENCUT = 400 eV ISPIN = 2 VOSKOWN =1 NELM = 60 EDIFF = 0.1E-04 LCHARG = .F. LWAVE = .F. ISMEAR = 0 SIGMA = 0.2 ALGO = Fast
！打开自旋极化 ！计算磁性需要
22
练习1 O2 分子
POSCAR
1）OUTCAR
vi OUTCAR
Ctrl+b 向上翻页 如右图显示最后 一步离子弛豫过 后的体系总能量、 原子坐标及原子 受力。