| 查看: 930 | 回复: 2 | |||
[交流]
【求助】紧急求教读取CHGCAR的matlab脚本，谢谢
|
| 现在需要深入研究vasp的CHGCAR，希望用matlab读取它，但对矩阵的读取一直没查到相关资料，求教哪位有读能提供取CHGCAR的matlab小脚本，谢谢。 |
回复此楼» 猜你喜欢
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有12人回复
-
359求调剂
已经有7人回复
-
070300化学专业279调剂
已经有10人回复
-
环境工程297分求调剂一志愿杭高院
已经有9人回复
-
一志愿武理材料工程302调剂环化或化工
已经有9人回复
-
309分085801求调剂
已经有5人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有5人回复
-
材料调剂
已经有4人回复
-
材料调剂
已经有3人回复
-
311（085601）求调剂
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 继续求助Matlab绘图！
已经有8人回复
- 弱弱问一下：stereochemical activity 是什么意思？ 谢谢
已经有5人回复
- 紧急求助翻译一篇德文文献~翻译其中合成的一部分~急求~谢谢~
已经有4人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 帮翻译几段文献，请不要在线翻译，谢谢
已经有2人回复
- 用matlab统计图片中某种颜色成分占总图片中所有颜色成分的比率
已经有4人回复
- 【求助】弱弱的问一下Matlab数据提取的问题
已经有5人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】CHGCAR用什么来读取
已经有10人回复
- 【求助】matlab如何只读取图像的一部分？
已经有11人回复
- 【求助】有哪位大侠做过次氯酸叔丁酯啊？帮帮我啊！谢谢
已经有6人回复
- 【求助】matlab读取fortran输出的txt文件【已解决】
已经有3人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
武汉纺织大学全国重点实验室招收材料、化学、纺织相关研究生
+1/67
-
山东青岛-导师招生
+1/52
-
留学--博士招生
+1/50
-
重庆交通大学 2026硕士调剂，材料化工缺额59+8人！能源动力缺额11人！！
+1/48
-
中国石油大学（华东）化学化工学院学硕、专硕接收调剂
+1/41
-
【#上海调剂急录#】985联合培养！优质课题组，带你发一区文章！
+1/40
-
天津科技大学-生物质资源利用与制浆技术团队2026年硕士招生
+1/10
-
烟台大学先进陶瓷方向 招收材料专业调剂研究生
+1/9
-
【博士招生】广东工业大学国家优青课题组招收2026年环境方向博士生
+1/6
-
哈工大赫晓东院士团队隋超教授招收材料学博士1名
+1/5
-
招生材料与化工、环境科学调剂研究生
+1/5
-
澳门理工大学2026人工智能药物发现博士（奖学金+申请考核）
+1/3
-
广东医科大学 李雪萌教授团队 招收申请考核制博士（临床检验诊断）
+1/3
-
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生）
+1/3
-
武汉纺织大学-材料科学与工程学院-杰青团队招收调剂研究生
+3/3
-
【上海调剂捡漏】985联合培养！别错过能送你去交大/同济/中科院的神仙导师！
+1/3
-
2026年上海工程技术大学材料学院智能材料与器件课题组招收硕士调剂生。
+1/2
-
英国全额资助博士招生｜射频与无线系统方向（全球招生）
+1/1
-
南京林业大学化工院有机合成方向 硕士一名
+1/1
-
2026年青岛农业大学作物应用基因组团队招收考研调剂生
+1/1
★ ★
sunyang1988(金币+2): 谢谢帮助。这样直接粘贴可能有格式错误，如果方便请上传到网盘 2011-01-16 14:10:22
雁儿霏霏(金币+20): 谢谢大虾，能给提供个程序文件吗？粘贴有一些格式错误。邮箱：shangyan2009@gmail.com ，谢谢！ 2011-01-16 19:37:24
sunyang1988(金币+2): 谢谢帮助。这样直接粘贴可能有格式错误，如果方便请上传到网盘 2011-01-16 14:10:22
雁儿霏霏(金币+20): 谢谢大虾，能给提供个程序文件吗？粘贴有一些格式错误。邮箱：shangyan2009@gmail.com ，谢谢！ 2011-01-16 19:37:24
|
% Read in CHG file to 3D matrix. fid = fopen('CHG','r'); % Read POSCAR part of CHG file: sysname = fgetl(fid); lconst = fscanf(fid,'%f',1); basisvec = fscanf(fid,'%f',[3,3])'; temp1 = fgetl(fid); temp2 = fgetl(fid); species = sscanf(temp2,'%i'); Natoms = sum(species); temp3 = fgetl(fid); atompos_f3 = fscanf(fid,'%f',[3,Natoms])'; volume = abs(dot(basisvec(1,  ,cross(basisvec(2,,basisvec(3,)))*lconst^3;% Read charge density data. (i,j,k) value in i+(j+k*Ny)*Nx, with coordinate % system defined by basis vectors!: gridsize = fscanf(fid,'%i',3); Nx = gridsize(1); Ny = gridsize(2); Nz = gridsize(3); temp4 = fscanf(fid,'%f',[Nx*Ny*Nz,1]); density_sub_h1_surf=reshape(temp4,Nx,Ny,Nz)/volume; % % Magnitization (if present): % temp5 = fscanf(fid,'%i',3); % temp6 = fscanf(fid,'%f',[Nx*Ny*Nz,1]); % magn_lowacc_t2_N�蕆eshape(temp6,Nx,Ny,Nz)/volume; % minval = min(min(min(density))); % maxval = max(max(max(density))); % For rectangular unit cell, unit axes: % Plot isosurface: % prect = patch(isosurface(density,7.5), 'FaceColor', 'red', 'EdgeColor', 'none'); % isonormals(density,prect); % %axis tight; % camlight; lighting phong; % Non-rectangular (general) unit cell: % Need to created position 3D matrices xtemp = zeros(Nx*Ny*Nz,1); ytemp = zeros(Nx*Ny*Nz,1); ztemp = zeros(Nx*Ny*Nz,1); for k = 0:Nz-1 for j = 0:Ny-1 for i = 0:Nx-1 xtemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,1)+j/Ny*basisvec(2,1)+k/Nz*basisvec(3,1); ytemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,2)+j/Ny*basisvec(2,2)+k/Nz*basisvec(3,2); ztemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,3)+j/Ny*basisvec(2,3)+k/Nz*basisvec(3,3); end end end X = reshape(xtemp,Nx,Ny,Nz)*lconst; Y = reshape(ytemp,Nx,Ny,Nz)*lconst; Z = reshape(ztemp,Nx,Ny,Nz)*lconst; z = squeeze(Z(1,1, );ave_z_h1 = squeeze(sum(sum(density_sub_h1_surf,1),2))/(Nx*Ny); %(sqrt(2)*Nx*Ny)*lconst^2; figure %subplot(211) plot(z,ave_z_h1,'k') % figure % p = patch(isosurface(X,Y,Z,density,7.5), 'FaceColor', 'red', 'EdgeColor', 'none'); % isonormals(density,p); % % maxx = max(basisvec(:,1))*lconst; % % minx = min(basisvec(:,1))*lconst; % % maxy = max(basisvec(:,2))*lconst; % % miny = min(basisvec(:,2))*lconst; % % maxz = max(basisvec(:,3))*lconst; % % minz = min(basisvec(:,3))*lconst; % % axis([minx maxx miny maxy minz maxz]) % camlight; lighting phong; % Extract charge density along certain lines: ------------------------------ |
2楼2011-01-16 12:51:04
3楼2011-12-30 22:59:32
