| 查看: 4463 | 回复: 26 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
[交流]
（文献学习）一篇计算HfO2表面 结构与原子性质的文章
|
||||||
|
Structural and electronic properties of cubic HfO2 surfaces G.H. Chen, Z.F. Hou, X.G. Gong Surface Physics Laboratory and Department of Physics, Fudan University, Shanghai 200433, China 文章不错，很有学习价值 现在想跟大家交流一下此文中的各部分计算的方法，向大家请教一下Projected bulk band structures 的计算方法。此贴打算一直写下去，直到把这篇文章搞懂  HfO2 表面计算 1. (偷偷懒，不做测试了,参数取文献中的值) 先以（100）-Hf 为例计算 2. 在ISCD中找到cubic-HfO2,优化晶体结构 晶格常数a=4.7 4.8 4.9 5.0 5.1 5.2 5.3 5.4 5.5 结果(E-V) 103.8200 -115.935673 110.5900 -119.148061 117.6500 -121.097380 125.0000 -122.006153 132.6500 -122.060947 140.6100 -121.419122 148.8800 -120.206692 157.4600 -118.539643 166.3800 -116.507645  3. 取a=5.1, 设ISIF=3，做结构优化，收敛精度要高 优化后 length of vectors 5.054564137 5.054564137 5.054564137 4. 建立(1 0 0)-Hf表面  5. 优化-relaxation ISIF=2 中间三层固定住 Hf O 1.000000000000000 3.574116577215030 0.000000000000000 0.000000000000000 0.000000000000000 3.574116577215030 0.000000000000000 0.000000000000000 0.000000000000000 25.163692411000000 7 12 Selective dynamics Direct 0.000000000000000 0.000000000000000 0.198698979399951 T T T 0.000000000000000 0.000000000000000 0.399566326466650 T T T 0.000000000000000 0.000000000000000 0.801301020600049 T T T 0.000000000000000 0.000000000000000 0.600433673533350 T T T 0.500000000000000 0.500000000000000 0.299132652933301 T T T 0.500000000000000 0.500000000000000 0.700867347066700 T T T 0.500000000000000 0.500000000000000 0.500000000000000 F F F 0.000000000000000 0.500000000000000 0.248915816166626 T T T 0.000000000000000 0.500000000000000 0.449783163233325 F F F 0.000000000000000 0.500000000000000 0.650650510300024 T T T 0.500000000000000 0.000000000000000 0.248915816166626 T T T 0.500000000000000 0.000000000000000 0.449783163233325 F F F 0.500000000000000 0.000000000000000 0.650650510300024 T T T 0.000000000000000 0.500000000000000 0.751084183833374 T T T 0.000000000000000 0.500000000000000 0.550216836766675 F F F 0.000000000000000 0.500000000000000 0.349349489699976 T T T 0.500000000000000 0.000000000000000 0.751084183833374 T T T 0.500000000000000 0.000000000000000 0.550216836766675 F F F 0.500000000000000 0.000000000000000 0.349349489699976 T T T 优化后结构  6.不算表面能了，直接算 projected band structures（可是怎么算呢，不会啊， ）未完待续，另提供本人蹂躏版文献下载，欢迎大家前来拍砖，拍砖有奖，打酱油没有  [ Last edited by huangll99 on 2012-5-2 at 17:32 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Structural and electronic properties of cubic HfO2 surfaces.pdf
2012-04-30 17:28:26, 453.86 K
» 收录本帖的淘帖专辑推荐
 精华网帖收集 | VASP算例 | 计算-vasp |
» 本帖已获得的红花（最新10朵）
legodgz» 猜你喜欢
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有5人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
招博士
已经有6人回复
-
限项规定
已经有8人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp中怎样可视化原子轨道
已经有24人回复
- vasp计算界面结构时原子弛豫的问题
已经有12人回复
- 如何处理在体结构中掺杂小原子（如B，N，C等）并做结构弛豫？
已经有39人回复
- 两个表面之间作用能或原子吸附能的计算疑惑
已经有9人回复
- 求助以下原子（分子）最稳定晶体结构的空间群
已经有4人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 【求助】多原子结构supercell的k点选择
已经有3人回复
- 【求助】为什么纳米粒子较小时，增加了表面原子的无序度和表面缺陷？
已经有11人回复
- 【求助】Ms计算磁性要不要设置每个原子的磁矩？？
已经有10人回复
- 【求助】请问各位虫友，从哪可以查物质晶体结构及原子间距啊
已经有5人回复
- 【求助】castep结构优化与性质计算
已经有5人回复
- 【求助】计算表面能要不要固定中间的原子
已经有12人回复
» 抢金币啦！回帖就可以得到:
-
哈尔滨工程大学青岛创新发展基地招聘青年教师
+1/475
-
“超分子材料交叉研究团队”联合诚聘博士后 [清华/吉大/复旦/北大]
+1/79
-
西北工业大学民航学院复合材料领域招聘两名博后
+1/79
-
【2026/2027 哈工大计算机类博士招生】
+1/78
-
娃娃们今儿考试喽。。。。
+1/77
-
西北工业大学民航学院招博士与硕士复合材料方向
+1/74
-
南京林业大学特聘教授团队招聘博后和2026博士研究生
+1/73
-
北京理工大学 珠海校区全职院士招数名博士生--申请考核制-半导体、光学、微电子
+2/64
-
华南理工大学(985、广州) 氢能源技术课题组诚招博士后（不限专业）
+1/49
-
东北大学-招收2026年硕士研究生2-3名（金属材料3D打印方向）
+5/40
-
中山大学农业与生物技术学院周潇峰课题组诚聘微生物/植物病理学方向科研助理
+1/28
-
诚征女友，非诚勿扰
+3/26
-
西班牙巴塞罗那访学、博后、留学互动
+1/15
-
青岛大学化学院泰山学者冯佳课题组申请考核制博士招生
+1/13
-
2026年天津大学杰青团队招收化学合成、计算机和器件的方面博士
+1/11
-
深圳大学信息功能聚合物电介质方向“申请-考核制”博士生招生
+2/10
-
【陕西师范大学】催化化学课题组2026年招收博士后/讲师/副高
+1/8
-
华南理工大学宋波教授招聘材料和化学方向博士后（长期有效）
+1/6
-
福建师范大学柔性电子学院 院士团队招2026级博士 光电器件、发光传感忆阻器
+1/3
-
New Jersey Institute of Technology 陈浩老师团队诚招26级秋季博士生
+1/2
huangll99: 回帖置顶 2012-05-01 10:31:30
|
7.接下来算bulk 的表面布里渊区投影能带结构，疑问是：处理ENGEIVAL是个大问题 初步计划，建五个目录，做五个KPOINTS,稍显笨拙，但我没有办法啊  K1----- kpoints 20 Line-mode rec 0 0 0 0 0.5 0 0 0.5 0 0.5 0.5 0 0.5 0.5 0 0 0 0 K2---- KPOINTS 20 Line-mode rec 0 0 0.125 0 0.5 0.125 0 0.5 0.125 0.5 0.5 0.125 0.5 0.5 0.125 0 0 0.125 K3---- KPOINTS 20 Line-mode rec 0 0 0.25 0 0.5 0.25 0 0.5 0.25 0.5 0.5 0.25 0.5 0.5 0.25 0 0 0.25 K4---- KPOINTS 20 Line-mode rec 0 0 0.375 0 0.5 0.375 0 0.5 0.375 0.5 0.5 0.375 0.5 0.5 0.375 0 0 0.375 K5--- KPOINTS 20 Line-mode rec 0 0 0.5 0 0.5 0.5 0 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0 0 0.5 |
8楼2012-05-01 10:30:47
huangll99: 回帖置顶 2012-05-01 09:55:33
|
参考了一下李明宪的教程，projected band structure 计算分了两步，一个是沿着表面布里渊区的高对称K点计算slab的能带结构；一个是沿着表面布里渊区的高对称K点计算体材料的能带，这时，对于体材料的计算，K点的选取很特别，K1和K2和计算slab时一样，K3依次从0取到0.5，中间取三个点应该够了 哈哈，试一下 6.算slab的能带结构，k点 kpoints 20 Line-mode rec 0 0 0 0 0.5 0 0 0.5 0 0.5 0.5 0 0.5 0.5 0 0 0 0 然后处理EIGENVAL, bbscon-plotband.f90是个好程序，处理的.dat每个带都分好列了。嘿嘿 能带图  [ Last edited by huangll99 on 2012-5-1 at 10:06 ] |
3楼2012-04-30 21:14:26
楼主继续 4楼2012-05-01 01:47:02
6楼2012-05-01 10:00:52