| 查看: 1426 | 回复: 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) **** 不知该如何解决，不知是哪里出现问题导致计算停止了？此外，从图片可以看出优化后的结构和优化前有很大区别，不知道这样算出来的表面能数据是否可靠？是不是需要增加表面厚度呢？ |
回复此楼» 猜你喜欢
心脉受损
已经有5人回复
-
博士读完未来一定会好吗
已经有15人回复
-
Springer期刊投稿求助
已经有4人回复
-
读博
已经有3人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
谢谢大家的讨论，不胜感激 我下了GULP[1]3.1,压缩文件2.19M,解压后有个exe文件夹，点gulp.exe后出现了好像是要输命令的窗口，不懂。。。看了readme: Distribution for GULP 3.0.1 : ============================+ The following is a guide to the directories that should be present in this distribution : Docs : documentation and information relating to the program Examples : example input files for the program Libraries : library files for the program Src : source code and files relating to the installation of the program To install GULP on most Unix/Linux machines : (1) Go to Src/ (2) Edit local.F to set the pointers to where the files will be kept (3) Type "make" (4) If the above fails - edit "getmachine" and correct the settings for your machine NOTE for Linux users : In order to compile the program you need to have an f90 compiler installed. By default, the Intel ifort compiler is assumed. 请问是什么意思呢？ [ Last edited by michelle8812 on 2009-2-24 at 21:38 ] |
14楼2009-02-24 12:18:04
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