应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【讨论】电子结构图的制作
北京石油化工学院2026年研究生招生接收调剂公告
182/2返回列表上一页12
查看: 1637  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

laohaoren2802

银虫 (正式写手)
期待高手指点。
回复此楼
11楼2009-06-28 18:42:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meimeidyh

铁杆木虫 (著名写手)

进城务工人员
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
spur(金币+1,VIP+0): 3Q！！谢谢参与！欢迎常来！ 7-1 12:24
是用castep做的。
先做好彩色的电荷密度切面，然后选择color maps，在color mapping上选取合适的具有代表性的值，就会出现如图所示的等高线了，再在spectrum栏选择black-white就可以了。不过还是建议用matlab做电荷密度等高线，更直观。
赞一下(11人) 回复此楼
12楼2009-07-01 09:18:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sss3012

银虫 (正式写手)
zhfuchun  朋友，能否介绍一下你是如何用MATLAB 画出电荷密度的，如何从CASTEP电荷密度计算中获得数据，谢谢
赞一下 回复此楼
13楼2009-07-01 11:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meimeidyh

铁杆木虫 (著名写手)

进城务工人员
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+5,VIP+0):谢谢 7-2 22:28
转一个好帖给你看看（用matlab做电荷密度等高线）：
1.找到电荷密度图的存储文件，一般是在自己所建的计算project中，隐藏的，castep模块是.charg_frm格式的，可以用txt打开后另存为txt格式。
2.读取数据，在matlab命令里输入一下命令：
[path,fn]=uigetfile('*.txt','Open');
fp=fopen([fn,path],'r');
head=fscanf(fp,'%s',4);
data=fscanf(fp,'%f',[4,18*18*144]);
fclose(fp);
这是读取刚才保存的txt数据的，其中18*18*144是可以改的，分别对应a,b,c的重复单元。
3.作矩阵。
a=reshape(data(1,,[18,18,144]);
b=reshape(data(2,,[18,18,144]);
c=reshape(data(3,,[18,18,144]);
d=reshape(data(4,,[18,18,144]);
其中18*18*144是可以改的，分别对应a,b,c的重复单元，即长宽高的值。
4.作等高线图
layer=1;
contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8);
layer=1（2，3，4....） 可以修改作图的层数，分别是从前面到后面。最后的8表示显示线条数目，可以自己修改成9，10或其它。
hold on 表示累积作图，hold off表示清空前面的 作图
举例；做表面层和第九层的等高线图
>> layer=9;
contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8);
>> hold on
>> layer=1;
contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8);

这是沿侧面做等高线（即xz平面），其它方向作图类似。只要修改其中的参数就可以了，如把a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8);
改为a(layer,:,,[18,144]),reshape(c(layer,l:,,[18,144]),reshape(d(layer,:,,[18,144]),8);

[ Last edited by meimeidyh on 2009-7-2 at 13:17 ]
赞一下(11人) 回复此楼
14楼2009-07-02 13:14:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
很好啊，初学者路过学习啦
赞一下(1人) 回复此楼
15楼2009-07-04 13:10:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zi_shu

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
还没见过这样的图，学习了。
赞一下(1人) 回复此楼
16楼2009-07-07 14:35:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mink

至尊木虫 (著名写手)
呵呵，只要会编程，处理数据不难，只是细心点。不过要多摸索，我用matla做过
赞一下 回复此楼
╱�?█◣要厚道!!
17楼2009-07-08 08:37:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qli0218

木虫 (小有名气)

问问楼上的

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我从dmol算的结果里也画出了电子图，不过先在也只是把数据导出来了，没有进行切面，想问一下用matlab如何进行切面阿？怎么改变画出的电荷密度图线的疏密程度？
谢谢
赞一下(1人) 回复此楼
18楼2009-07-08 08:56:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sss3012 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 347材料专硕求调剂 +15 zj8215216 2026-04-06 15/750 2026-04-08 00:22 by JourneyLucky
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 土木水利专硕276分求调剂 +4 我想上学！！6 2026-04-05 7/350 2026-04-07 22:10 by Delta2012
[考研] 085404 293求调剂 +8 勇远库爱314 2026-04-06 9/450 2026-04-07 13:05 by flydream1314
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 求助071001调剂！！！ +4 黄守松 2026-04-05 5/250 2026-04-06 10:55 by 1028907439
[考研] 298求调剂 +3 manman511 2026-04-05 3/150 2026-04-05 18:09 by kk112233
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 324求调剂 +9 想上学求调 2026-04-03 9/450 2026-04-04 23:57 by 果冻大王
[考研] 一志愿重庆大学085404，总分314分，求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 一志愿陕西师范大学生物学317分 +5 1563日。 2026-04-02 5/250 2026-04-03 06:58 by ilovexiaobin
[考研] 一志愿大工学硕，求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +4 爱读书的小鳄鱼 2026-04-02 4/200 2026-04-02 18:35 by 求调剂zz
[考研] 318求调剂 +8 七忆77 2026-04-01 8/400 2026-04-01 10:37 by Jaylen.
信息提示
关闭
请填处理意见
关闭