设置结构常数信息
计算需要知道:原胞尺寸,原子坐标和晶格常数
(这些从实验**文献中获得 或者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时是不变的)