24小时热门版块排行榜    

查看: 646  |  回复: 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文章,我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:50 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 14:14 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见