计算化学公社

标题: 关于 LOCPOT 文件格式的几个疑惑 [打印本页]

作者
Author: 北大-陶豫 时间: 2023-2-13 11:18
标题: 关于 LOCPOT 文件格式的几个疑惑
本帖最后由北大-陶豫于 2023-2-13 11:22 编辑

LOCPOT 的格式是，开头是体系的结构信息（和POSCAR差不多），之后一行三个数字，分别是 NGXF NGYF NGZF，然后是所有的格点数据，我的疑问如下：
（1）这些格点数据是按什么顺序排列的呢，是x变化最慢、z变化最快，还是什么别的顺序？
x变化最慢、z变化最快，是指：

double LOCPOT[NGXF][MGYF][NGZF]; // 三维数组，三个维度的下标范围分别是 0,1, ..., NGXF-1 0,1, ..., NGYF-1 0,1, ..., NGZF-1
for (i=0; i<NGXF; ++i) // python: for i in range(NGXF)
for (j=0; j<NGYF; ++j) // python: for j in range(NGYF)
for (k=0; k<NGZF; ++k) // python: for k in range(NGZF)
scanf("%lf",&LOCPOT[i][j][k]); // 读入一个实数存进数组

复制代码

（2）文件中的 (NGXF*NGYF*NGZF) 个格点，对应晶胞中的空间位置是什么呢？是不是，NGXF方向第i个格点（i=0,1, ..., NGXF-1），NGYF方向第j个格点（j=0,1, ..., NGYF-1），NGZF方向第k个格点（k=0,1, ..., NGZF-1），分数坐标为 (i/NGXF, j/NGYF, k/NGZF)？
谢谢！

作者
Author: sobereva 时间: 2023-2-13 11:34
这是Multiwfn读取VASP的格点数据的Fortran代码的一部分：
read(10,*) (((cubmat(ix,iy,iz),ix=1,nx),iy=1,ny),iz=1,nz)
ix循环完一轮变一次iy，iy循环完一轮变一次iz

是

欢迎光临计算化学公社 (http://bbs.keinsci.com/)