24小时热门版块排行榜    

查看: 623  |  回复: 1

xiaowu787

木虫 (正式写手)

[交流] 【求助】怎么计算Al2O3的状态方程

现在做分子动力学,希望对比一下VASP优化计算Al2O3的状态方程。前面也看过几个计算状态方程的较为简单的例子,现在对VASP刚接触,不知道这种晶格Al2O3的POSCAR该怎么写?设置KPOINTS中对于Al2O3需要注意什么?INCAR都需要那些参数?
有的说对于非立方晶系还要考虑压强,该怎么弄呢?请高手指点一下?


Al2O3的cif文件如下:

#### CIF created by Crystallographica 2 ####

data_AluminiumOxideAlpha

_audit_creation_method                    'Crystallographica 2'
_cell_angle_alpha                         90
_cell_angle_beta                          90
_cell_angle_gamma                         120
_cell_formula_units_Z                     6
_cell_length_a                            4.7602
_cell_length_b                            4.7602
_cell_length_c                            12.993
_cell_volume                              254.971
_cgraph_comments                         'Recovering experimental and theoretical electron densities       in corundum using the multipolar model : IUCr multipole       refinement project '
_cgraph_title                            'Aluminium Oxide - Alpha'
_chemical_formula_sum                    'Al2 O3'
_symmetry_space_group_name_H-M           'R -3 c'
_symmetry_space_group_name_Hall          '-R 3 2"c'

loop_
        _symmetry_equiv_pos_as_xyz
        'x, y, z'
        'x-1/3, y+1/3, z+1/3'
        'x+1/3, y-1/3, z-1/3'
        '-x, -y, -z'
        '-x-1/3, -y+1/3, -z+1/3'
        '-x+1/3, -y-1/3, -z-1/3'
        '-y, x-y, z'
        '-y-1/3, x-y+1/3, z+1/3'
        '-y+1/3, x-y-1/3, z-1/3'
        'y, -x+y, -z'
        'y-1/3, -x+y+1/3, -z+1/3'
        'y+1/3, -x+y-1/3, -z-1/3'
        '-x+y, -x, z'
        '-x+y-1/3, -x+1/3, z+1/3'
        '-x+y+1/3, -x-1/3, z-1/3'
        'x-y, x, -z'
        'x-y-1/3, x+1/3, -z+1/3'
        'x-y+1/3, x-1/3, -z-1/3'
        'x-y, -y, -z+1/2'
        'x-y-1/3, -y+1/3, -z-1/6'
        'x-y+1/3, -y-1/3, -z+1/6'
        '-x+y, y, z+1/2'
        '-x+y-1/3, y+1/3, z-1/6'
        '-x+y+1/3, y-1/3, z+1/6'
        '-x, -x+y, -z+1/2'
        '-x-1/3, -x+y+1/3, -z-1/6'
        '-x+1/3, -x+y-1/3, -z+1/6'
        'x, x-y, z+1/2'
        'x-1/3, x-y+1/3, z-1/6'
        'x+1/3, x-y-1/3, z+1/6'
        'y, x, -z+1/2'
        'y-1/3, x+1/3, -z-1/6'
        'y+1/3, x-1/3, -z+1/6'
        '-y, -x, z+1/2'
        '-y-1/3, -x+1/3, z-1/6'
        '-y+1/3, -x-1/3, z+1/6'

loop_
        _atom_site_label
        _atom_site_type_symbol
        _atom_site_fract_x
        _atom_site_fract_y
        _atom_site_fract_z
        _atom_site_U_iso_or_equiv
        _atom_site_thermal_displace_type
        _atom_site_occupancy
        Al1 Al3+ 0 0 0.35217 0 Uani 1
        O1 O2- 0.30634 0 0.25 0 Uani 1

loop_
        _atom_site_aniso_label
        _atom_site_aniso_U_11
        _atom_site_aniso_U_22
        _atom_site_aniso_U_33
        _atom_site_aniso_U_12
        _atom_site_aniso_U_13
        _atom_site_aniso_U_23
        Al1 0.00255 0.00255 0.00266 0.00128 0 0
        O1 0.00305 0.00305 0.0035 0.00153 0.00034 0.00068

_eof

#### End of Crystallographic Information File ####
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)


zzy870720z(金币+1):谢谢建议 2010-10-29 08:33:14
xiaowu787(金币+2): 谢谢 2011-04-25 18:57:53
问题太多,不好回答。建议多看VASP资料。
2楼2010-10-28 22:14:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowu787 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 6/300 2026-08-04 13:28 by jzdxtiUmJ6PL
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 13:22 by jzdxtiUmJ6PL
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 jRl3mE6ddZGq 2026-08-03 11/550 2026-08-04 12:58 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐,不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见