loading...
设置晶格常数
Published in:2021-12-30 |

设置结构常数信息

计算需要知道:原胞尺寸,原子坐标和晶格常数

(这些从实验**文献中获得 或者COD数据库@http://www.crystallography.net/cod/

COD数据库使用说明:

COD主页–search–搜索“1 to 8 elements”–“number of distinct elements min and max”一栏输入2和2表示搜索二元化合物–点击cif开始下载cif文件

由cif文件得到pw.x输入文件:

1)使用vesta,MS

用vesta打开,点击File–export Data—选择vasp格式

注意:导入cif文件一般是晶胞,而非原胞(对于非简单格子,如面心,体心等),为了转换为原胞,最简单的方法是MS,也可以在 vesta中输入晶胞到原胞的转换矩阵,即Edit–edit data

data–unit cell–transform,这样不太容易找对称性,通常会输出下三角格式

新版的vesta增加了生成原胞的功能

2)POSCAR格式与pw.x格式互相转换可以用QEtoolkit中的P2P工具@http://www.densityflow.com/p2p.php

3)使用QEtoolkit@http://www.densityflow.com/

4)使用cif2cell

命令:cif2cell *.cif -p pwscf --supercell=[3,3,1]

设置pw.x输入文件(ibrav,celldm,nat,ntyp)

1)法1:

设置ibrav和celldm(1-6)指定CELL,这时不写CELL_PARAMETERS

2)法2:

设置ibrav=0和CELL_PARAMETERS (angstrom),这时不允许写出celldm(1),程序会自动内部设置celldm(1) = alat = sqrt(v1*v1),这里celldm代表晶格常数,单位是Bohr

在结构驰豫(calculation=’relax’)过程中,ATOMIC_POSITIONS是根据力而变化的,如果是vc-relax,原子坐标改变的同时,CELL_PARAMETERS根据应力变化(celldm在vc-relax时是不变的)

注意:自6.4.1版本,官方不推荐celldm(1)=1.88972613(任何<2的值)的做法,>这里也修正为celldm(1)设置为晶格常数,或用ibrav≠0

Prev:
QE计算磁性
catalog
catalog