GULP 案例 2 代码解释：如何计算相变压力？

对前一个内容进行输入文件的解释

完整的 GULP 输入文件：

以 ZnO-WZ 为例进行解释：（此外，还需要计算RS和HC结构的焓值）

opti conp prop phonon

# opti— perform optimisation run（声明是结构优化）
# conp — constant pressure calculation（声明是等压计算）
# prop — calculate properties for final geometry（声明是最终的计算性质）
# phonon — calculate phonons for final geometry（声明是计算声子相关特性）

name ZnO-WZ
#声明计算物质的名字

Temperature 0.0 K   # 温度为零
Pressure 0.0 GPa   # 给定压力值

cell
3.242  3.242  5.188  90.0  90.0  120.0
#声明ZnO-WZ结构的晶格参数（晶格常数和角度）

frac
Zn core 0.333333   0.666667   0.000000  
O  core 0.333333   0.666667   0.381900
#声明Zn和O原子的归一化坐标

space
186
#声明物质的空间群

species
Zn   1.14
O   -1.14
# 指定电荷种类和电荷量

buck 
Zn  Zn  78.91   0.5177  0.00     0.0 12.0
Zn  O   257600  0.1396  0.00     0.0 12.0

Morse
O   O   0.1567  1.164   3.405    0.0 12.0 
#指定势函数的形式和势参数

dump every 1 ZnO-WZ.dump
#指定输出焓值等数据信息

计算所用代码：

opti conp prop phonon
name ZnO-WZ

Temperature 0.0 K
Pressure 15.0 GPa

cell
3.242  3.242  5.188  90.0  90.0  120.0
frac
Zn core 0.333333   0.666667   0.000000  
O  core 0.333333   0.666667   0.381900
space
186

species
Zn   1.14
O   -1.14

buck 
Zn  Zn  78.91   0.5177  0.00     0.0 12.0
Zn  O   257600  0.1396  0.00     0.0 12.0

Morse
O   O   0.1567  1.164   3.405    0.0 12.0 
dump every 1 ZnO-WZ.dump

提取数据的过程：

从输出文件中找到Final enthalpy所对应的值即是焓

  **** Optimisation achieved ****

  Final enthalpy =     -26.48176238 eV

  Final Gnorm  =       0.00021207

最后，通过Origin作图，GULP计算结果与文献一致！

欢迎大家交流！扫码即可关注微信公众号：原子与分子模拟