|
|
[交流]
[求助]GULP计算YAG晶体表面能
YAG,分子式为Y3Al5O12,属于立方晶系,空间群为Ia3d
单位晶胞中包含8个YAG分子,共160个原子,其中阳离子都处在特殊位置,O2-离子处在一般位置:
24个Al 离子处于O离子构成的正四面体的中心
16个Al 离子位于O离子构成的正八面体中心
这些正八面体和正四面体所占据的空间,其间形成一些十二面体的空隙,呈畸变立方,角上为O离子,中心位置为Y离子,单位晶胞中有24个Y
原子坐标:
Al 16 a (0, 0, 0),
Al 24 d (0.375, 0, 0.25),
Y 24 C (0.125, 0, 0.25),
O 96 h (-0.02985, 0.05056, 0.14878).
采用MS的GULP模块计算(001),(-110)(-112)三个晶面的表面能,采用“Bush”力场,创建的表面模型和优化后的结构如图所示



我建的表面模型厚度都是最小值,就是模型所含原子数与单位晶胞的相同,计算出来(-112)和(-110)的表面能分别为-0.06352和-0.1792 J/m2, 但是(001)的计算过程很快停止,无法继续进行,原文如下
Start of surface optimisation :
Cycle: 0 Energy: -2439.461143 Gnorm: 3.619983 CPU: 2.672
** Hessian calculated **
Cycle: 1 Energy: -2789.807186 Gnorm: 2.999379 CPU: 10.391
** Hessian calculated **
Cycle: 2 Energy: -3268.548968 Gnorm: 2.141373 CPU: 19.500
** Hessian calculated **
Cycle: 3 Energy: -3496.883530 Gnorm: 2.090267 CPU: 28.875
** Hessian calculated **
Cycle: 4 Energy: -3962.867941 Gnorm: 1.027086 CPU: 38.156
** Hessian calculated **
Cycle: 5 Energy: -4309.378260 Gnorm: 0.708665 CPU: 47.344
** Hessian calculated **
Cycle: 6 Energy: -4310.303948 Gnorm: 0.710238 CPU: 56.625
** Hessian calculated **
Cycle: 7 Energy: -4310.359825 Gnorm: 0.710336 CPU: 66.016
** Hessian calculated **
Cycle: 8 Energy: -4311.057138 Gnorm: 0.711578 CPU: 75.125
** Hessian calculated **
Cycle: 9 Energy: -4311.057138 Gnorm: 0.711578 CPU: 84.266
** Hessian calculated **
**** Conditions for a minimum have not been satisfied. However ****
**** no lower point can be found - treat results with caution ****
**** unless gradient norm is small (less than 0.1) ****
不知该如何解决,不知是哪里出现问题导致计算停止了?此外,从图片可以看出优化后的结构和优化前有很大区别,不知道这样算出来的表面能数据是否可靠?是不是需要增加表面厚度呢? |
|