| 查看: 1487 | 回复: 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) **** 不知该如何解决，不知是哪里出现问题导致计算停止了？此外，从图片可以看出优化后的结构和优化前有很大区别，不知道这样算出来的表面能数据是否可靠？是不是需要增加表面厚度呢？ |
回复此楼» 猜你喜欢
0703化学调剂 348分
已经有13人回复
-
327求调剂
已经有8人回复
-
材料307分求大佬组收留
已经有9人回复
-
求调剂
已经有33人回复
-
293调剂
已经有15人回复
-
求调剂求调剂
已经有11人回复
-
22408 266求调剂
已经有10人回复
-
一志愿211，化学310分，本科重点双非，求调剂
已经有8人回复
-
一志愿0703化学招61最终排名62化学求调剂
已经有10人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
doudouding
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2306.2
- 红花: 1
- 帖子: 182
- 在线: 47.2小时
- 虫号: 64507
- 注册: 2005-04-16
- 专业: 凝聚态物性I:结构、力学和
2楼2009-02-11 21:01:11
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
3楼2009-02-12 09:43:33
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
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
★ ★ ★
michelle8812(金币+2,VIP+0):非常感谢，初始接触GULP,很多地方都不太懂。。。 2-13 20:51
wuli8(金币+1,VIP+0):(*^__^*) 2-14 09:52
michelle8812(金币+2,VIP+0):非常感谢，初始接触GULP,很多地方都不太懂。。。 2-13 20:51
wuli8(金币+1,VIP+0):(*^__^*) 2-14 09:52
|
gulp的力场可以从网上下载，具体网址我也记不清， gulp官方网站下载源代码，里面带的力场应该全一些， 用google都可以搜索到。 top值影响的是表面的起始层面，试一试就可以看见一个虚线框在变动。 建议看看http://muchong.com/html/200612/364805.html表面建模 |
6楼2009-02-13 13:22:43
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
7楼2009-02-13 17:08:04
8楼2009-02-13 20:56:12
doudouding
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2306.2
- 红花: 1
- 帖子: 182
- 在线: 47.2小时
- 虫号: 64507
- 注册: 2005-04-16
- 专业: 凝聚态物性I:结构、力学和
9楼2009-02-22 20:44:24
10楼2009-02-23 16:25:44
