24小时热门版块排行榜

>论坛更新日志 (2258)
>虫友互识 (128)
>导师招生 (96)
>招聘信息布告栏 (36)
>硕博家园 (25)
>考博 (23)
>公派出国 (22)
>博后之家 (21)
>文献求助 (16)
>考研 (15)
>休闲灌水 (13)
>论文道贺祈福 (12)
>找工作 (12)
>药学 (11)
>论文投稿 (11)
>计算模拟 (10)

返回列表

当前主题已经存档。

wfbgc

木虫 (小有名气)

应助: 3 (幼儿园)
金币: 1802.1
散金: 10
红花: 2
帖子: 220
在线: 172小时
虫号: 448198
注册: 2007-11-01
性别: GG
专业: 原子和分子物理

[交流] 【求助】请问哪位达人能介绍一下在MS下建模，在abinit下计算？

首先，我是做实验的，最近特想做些理论计算。
abinit的建模相对比较麻烦，我有个师弟给我说过，可以在MS下建模，然后在abinit下计算，可是，一直不太懂，故请教论坛的各位达人。

PS：能否给一个比较详细的教程？

回复此楼

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

1楼2009-03-15 16:22:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

BigBlackBear

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 563.5
散金: 10
帖子: 124
在线: 106.2小时
虫号: 112883
注册: 2005-11-22
专业: 半导体材料

★ ★ ★ ★ ★ ★
wuchenwf(金币+6):谢谢 O(∩_∩)O

在MS下建模，得到原胞或超胞的参数（a,b,c大小，各原子坐标等信息），将这些信息拷贝到abinit输入文件相应的位置。

比如想计算一个2×2×2的AB超胞能带结构：

1.在MS中import AB，如果没有现成的，就找到AB的晶体学参数，手动建立。如果是import的，将AB转变成primitive cell，如果是手动建立的，估计你找到的晶体学参数就是按primitive cell给的。
2.将primitive cell转变成2×2×2 supercell，菜单上有这个功能。
3.export这个supercell的晶体结构参数，可以输出.car或.cif文件。

下例为一个.cif文件
－－－－－－－－
data_AB
_audit_creation_date             2008-11-21
_audit_creation_method          'Materials Studio'
_symmetry_space_group_name_H-M 'P1'
_symmetry_Int_Tables_number    1
_symmetry_cell_setting          triclinic
loop_
_symmetry_equiv_pos_as_xyz
  x,y,z
_cell_length_a                   6.1490
_cell_length_b                   6.1490
_cell_length_c                   6.1490
_cell_angle_alpha                60.0000
_cell_angle_beta                60.0000
_cell_angle_gamma                60.0000
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_adp_type
_atom_site_occupancy
A1    A    0.12500 0.12500 0.12500 0.00000  Uiso 1.00
B2 B -0.00000 0.00000  -0.00000 0.00000  Uiso 1.00
A3    A    0.62500 0.12500 0.12500 0.00000  Uiso 1.00
B4 B 0.50000 0.00000  -0.00000 0.00000  Uiso 1.00
A5    A    0.12500 0.62500 0.12500 0.00000  Uiso 1.00
B6 B -0.00000 0.50000  -0.00000 0.00000  Uiso 1.00
A7    A    0.62500 0.62500 0.12500 0.00000  Uiso 1.00
B8 B 0.50000 0.50000 0.00000 0.00000  Uiso 1.00
A9    A    0.12500 0.12500 0.62500 0.00000  Uiso 1.00
B10 B -0.00000 0.00000 0.50000 0.00000  Uiso 1.00
A11 A    0.62500 0.12500 0.62500 0.00000  Uiso 1.00
B12 B 0.50000  -0.00000 0.50000 0.00000  Uiso 1.00
A13 A    0.12500 0.62500 0.62500 0.00000  Uiso 1.00
B14 B 0.00000 0.50000 0.50000 0.00000  Uiso 1.00
A15 A    0.62500 0.62500 0.62500 0.00000  Uiso 1.00
B16 B 0.50000 0.50000 0.50000 0.00000  Uiso 1.00
loop_
－－－－－－－－－
4.从.cif文件中找出supercell的a,b,c以及各原子坐标
即
_cell_length_a                   6.1490
_cell_length_b                   6.1490
_cell_length_c                   6.1490

_cell_angle_alpha                60.0000
_cell_angle_beta                60.0000
_cell_angle_gamma                60.0000

A1    A    0.12500 0.12500 0.12500 0.00000  Uiso 1.00
B2 B -0.00000 0.00000  -0.00000 0.00000  Uiso 1.00
A3    A    0.62500 0.12500 0.12500 0.00000  Uiso 1.00
B4 B 0.50000 0.00000  -0.00000 0.00000  Uiso 1.00
A5    A    0.12500 0.62500 0.12500 0.00000  Uiso 1.00
B6 B -0.00000 0.50000  -0.00000 0.00000  Uiso 1.00
A7    A    0.62500 0.62500 0.12500 0.00000  Uiso 1.00
B8 B 0.50000 0.50000 0.00000 0.00000  Uiso 1.00
A9    A    0.12500 0.12500 0.62500 0.00000  Uiso 1.00
B10 B -0.00000 0.00000 0.50000 0.00000  Uiso 1.00
A11 A    0.62500 0.12500 0.62500 0.00000  Uiso 1.00
B12 B 0.50000  -0.00000 0.50000 0.00000  Uiso 1.00
A13 A    0.12500 0.62500 0.62500 0.00000  Uiso 1.00
B14 B 0.00000 0.50000 0.50000 0.00000  Uiso 1.00
A15 A    0.62500 0.62500 0.62500 0.00000  Uiso 1.00
B16 B 0.50000 0.50000 0.50000 0.00000  Uiso 1.00
注意这里的a、b、c单位为埃，原子坐标为fractal 格式(见_atom_site_fract_x)

5 abinit输入文件

# computation of the total energy using a conventional cell(2 2 2)
# Computation of the band structure.
# First, a SCF density computation, then a non-SCF band structure calculation.

ndtset 2

#Dataset 1 : usual self-consistent calculation
kptopt1 1       # Option for the automatic generation of k points,
               # taking into account the symmetry
nshiftk1 4
shiftk1  0.5 0.5 0.5  # These shifts will be the same for all grids
      0.5 0.0 0.0
      0.0 0.5 0.0
      0.0 0.0 0.5
ngkpt1  4 4 4
prtden1  1       # Print the density, for use by dataset 2
toldfe1  1.0d-6

#Dataset 2 : the band structure
iscf2 -2
getden2  -1
kptopt2  -3
nband2 8
ndivk2 10 12 17    # 10, 12 and 17 divisions of the 3 segments, delimited
                     # by 4 points.
kptbounds2  0.5  0.0  0.0 # L point
         0.0  0.0  0.0 # Gamma point
         0.0  0.5  0.5 # X point
         1.0  1.0  1.0 # Gamma point in another cell.
tolwfr2  1.0d-12
enunit2  1          # Will output the eigenenergies in eV

#Definition of the unit cell
chkprim 0

acell 3*11.6199 # the primitive AB supercell(2 2 2) a= 11.61992783 Bohr
angdeg 60 60 60

#Definition of the atom types
ntypat 2       # There is two types of atoms
znucl A  B       # The keyword "znucl" refers to the atomic number of the
               # possible type(s) of atom. The pseudopotential(s)
               # mentioned in the "files" file must correspond



#Definition of the atoms
natom 16          # There are 16 atoms in the 2*2*2 supercell
typat 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 #
xred             # This keyword indicate that the location of the atoms
               # will follow, one triplet of number for each atom
0.12500 0.12500 0.12500
0.00000 0.00000 0.00000
0.62500 0.12500 0.12500
0.50000 0.00000 0.00000
0.12500 0.62500 0.12500
0.00000 0.50000 0.00000
0.62500 0.62500 0.12500
0.50000 0.50000 0.00000
0.12500 0.12500 0.62500
0.00000 0.00000 0.50000
0.62500 0.12500 0.62500
0.50000 0.00000 0.50000
0.12500 0.62500 0.62500
0.00000 0.50000 0.50000
0.62500 0.62500 0.62500
0.50000 0.50000 0.50000

#Definition of the planewave basis set
ecut  10       # Maximal kinetic energy cut-off, in Hartree

#Definition of the SCF procedure
nstep 100       # Maximal number of SCF cycles
toldfe 1.0d-6    # Will stop when, twice in a row, the difference
               # between two consecutive evaluations of total energy
               # differ by less than toldfe (in Hartree)
diemac 12.0    # Although this is not mandatory, it is worth to
               # precondition the SCF cycle. The model dielectric
               # function used as the standard preconditioner
               # is described in the "dielng" input variable section.
               # Here, we follow the prescription for bulk silicon.
6 计算
7 分析结果

[ Last edited by BigBlackBear on 2009-3-19 at 16:35 ]

赞一下(2人)

回复此楼

9年

2楼2009-03-19 16:31:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)

1ST强帖: 1
应助: 61 (初中生)
贵宾: 0.33
金币: 5959.3
散金: 1
红花: 21
帖子: 1884
在线: 491.1小时
虫号: 142265
注册: 2005-12-21
专业: 凝聚态物性 II ：电子结构

★
wuchenwf(金币+1):谢谢 O(∩_∩)O

引用回帖:

Originally posted by wfbgc at 2009-3-15 16:22:
首先，我是做实验的，最近特想做些理论计算。
abinit的建模相对比较麻烦，我有个师弟给我说过，可以在MS下建模，然后在abinit下计算，可是，一直不太懂，故请教论坛的各位达人。

PS：能否给一个比较详细的教程？

还是找个理论组合作吧。

赞一下(1人)

回复此楼

3楼2009-03-19 20:13:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖