24小时热门版块排行榜    

查看: 663  |  回复: 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一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见