| 查看: 1492 | 回复: 15 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
[求助]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) **** 不知该如何解决，不知是哪里出现问题导致计算停止了？此外，从图片可以看出优化后的结构和优化前有很大区别，不知道这样算出来的表面能数据是否可靠？是不是需要增加表面厚度呢？ |
回复此楼» 猜你喜欢
291 求调剂
已经有22人回复
-
296求调剂
已经有3人回复
-
一志愿北理工298英一数二求调剂
已经有11人回复
-
材料与化工调剂
已经有18人回复
-
080500求调剂
已经有14人回复
-
284求调剂
已经有7人回复
-
085800 能源动力求调剂
已经有5人回复
-
287求调剂
已经有20人回复
-
298求调剂
已经有10人回复
-
求调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
3楼2009-02-12 09:43:33
doudouding
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2306.2
- 红花: 1
- 帖子: 182
- 在线: 47.2小时
- 虫号: 64507
- 注册: 2005-04-16
- 专业: 凝聚态物性I:结构、力学和
2楼2009-02-11 21:01:11
4楼2009-02-13 09:55:33
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
5楼2009-02-13 13:12:19
