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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何显示3D spin density
131/2返回列表12下一页
查看: 958  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

xiaowenzhi16

铁杆木虫 (正式写手)

[交流] 【求助】如何显示3D spin density

请问VASP里的自旋密度如何显示为3D spin density 的电子云状的图像?用什么软件可以做到呢？？
谢谢
回复此楼

» 猜你喜欢
1楼 2009-05-13 20:43:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaofenglili

银虫 (小有名气)
★ ★ ★ ★ ★ ★
qasd(金币+2,VIP+0):谢谢~ 5-13 21:07
xiaowenzhi16(金币+3,VIP+0):谢谢啦，过一下再给哦。 5-13 22:20
xiaowenzhi16(金币+1,VIP+0):谢谢 5-14 12:00
VESTA可以！
回复此楼
2楼2009-05-13 20:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowenzhi16

铁杆木虫 (正式写手)
具体怎么操作呢，我试了很久，能做错电荷密度三维图，但是自旋三维图不行呀，能够详细点不？先谢谢啦
赞一下 回复此楼
3楼2009-05-13 22:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)
★ ★ ★ ★
xiaowenzhi16(金币+3,VIP+0)::谢谢啦，过一下再给哦 5-14 10:15
xiaowenzhi16(金币+1,VIP+0):谢谢 5-14 12:00
引用回帖:
Originally posted by xiaowenzhi16 at 2009-5-13 22:19:
具体怎么操作呢，我试了很久，能做错电荷密度三维图，但是自旋三维图不行呀，能够详细点不？先谢谢啦

你要自己先处理下CHGCAR文件,得到自旋电荷密度,在画图.
赞一下(2人) 回复此楼
4楼2009-05-14 09:21:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowenzhi16

铁杆木虫 (正式写手)
你要自己先处理下CHGCAR文件,得到自旋电荷密度,在画图.

楼上的专家，谢谢你的提示，我也看出来了，自旋电荷密度是在CHGCAR文件当中,就是不会把它的数据提取出来。请不吝赐教。
谢谢。
赞一下 回复此楼
5楼2009-05-14 10:14:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

铁杆木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xiaowenzhi16(金币+6,VIP+0):实在是太感谢啦。以后还要麻烦你的。 5-14 12:00
wuchenwf(金币+10,VIP+0):谢谢 5-14 23:01
给你一段我自己编的程序，提取电荷或自旋密度图

program readcharge4vasp
!..............................................................
!    read density from name=CHG file of VASP
!..............................................................
!    spin =0 - total charge \____ only for CHGCAR
!         =1 - spin density /
!..............................................................

character title*20,type_coodinates*20
integer :: natoms=0,NGX,NGY,NGZ,n=1,totalatoms=0,spin
integer,allocatable :: ntypat(
real*8 :: scale, basic_vectorx,basic_vectory,basic_vectorz
real*8, allocatable :: grid(:,:,
write(*,*)'Number of atom types in the unit cell'
read(*,*) n
write(*,*)'Charge density or spin density? (0 or 1)'
read(*,*) spin
allocate(ntypat(n))
open(unit=11,file='CHG')
if (spin==0) then
   open(unit=12,file='Charge_Density.vasp')
else if (spin==1) then
   open(unit=12,file='Spin_Density.vasp')
end if

read(11,'(A15)') title
write(12,'(a15)') title
read(11,*) scale
write(12,"(F6.3)" scale

do i=1,3
read(11,*) basic_vectorx,basic_vectory,basic_vectorz
write(12,'(3F13.7)')basic_vectorx,basic_vectory,basic_vectorz
end do
read(11,*) (ntypat(i),i=1,n)
write(12,*) (ntypat(i),i=1,n)

read(11,*) type_coodinates
write(12,*) type_coodinates

do i=1,n
totalatoms=ntypat(i)+totalatoms
end do

do i=1,totalatoms
read(11,*) x,y,z
write(12,"(3F11.7)" x,y,z
end do
read(11,*)
write(12,*)

read(11,*) NGX,NGY,NGZ
write(12,'(3I5)') NGX,NGY,NGZ
allocate(grid(NGX,NGY,NGZ))

if (spin==1) then
   read(11,"(10E12.5E2)" (((grid(i,j,k),i=1,NGX),j=1,NGY),k=1,NGZ)
   read(11,*)
else if (spin==0) then
end if

read(11,"(10E12.5E2)" (((grid(i,j,k),i=1,NGX),j=1,NGY),k=1,NGZ)
write(12,"(10E13.5E2)" (((grid(i,j,k),i=1,NGX),j=1,NGY),k=1,NGZ)
close(11)
close(12)
end


之后用vesta打开。

[ Last edited by 灯塔守望者 on 2009-5-14 at 12:18 ]
赞一下(12人) 回复此楼
6楼2009-05-14 11:25:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)
★
xiaowenzhi16(金币+1,VIP+0):谢谢跟帖 5-14 12:01
跟着学习：我想问我是在主机上的linux下计算的，我如和用在windows下装的vesta软件处理，因为我没有权限在主机上装软件，谢谢
赞一下 回复此楼
7楼2009-05-14 11:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
谢谢！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
回复此楼
8楼2009-05-14 12:16:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

子虚乌有5388

木虫 (正式写手)
★
wuchenwf(金币+1,VIP+0):谢谢 5-14 23:02
把数据拷贝出来  然后用gunplot可以作三维图吧
赞一下(10人) 回复此楼
房无一间，地无一亩，身无分文，胸无大志。
9楼2009-05-14 15:16:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)
★
wuchenwf(金币+1,VIP+0):谢谢 5-14 23:02
引用回帖:
Originally posted by ddang100 at 2009-5-14 11:30:
跟着学习：我想问我是在主机上的linux下计算的，我如和用在windows下装的vesta软件处理，因为我没有权限在主机上装软件，谢谢

VESTA 可以在自己的windows使用的.
赞一下(10人) 回复此楼
10楼2009-05-14 16:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowenzhi16 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334求调剂 +13 Riot2025 2026-04-08 14/700 2026-04-09 00:12 by lature00
[考研] 311求调剂 +5 surte 2026-04-08 11/550 2026-04-08 22:10 by 土木硕士招生
[考研] 265求调剂 +3 风说她早忘了 2026-04-07 3/150 2026-04-08 20:20 by 逆水乘风
[考研] 320分人工智能调剂 +9 振—TZ 2026-04-03 10/500 2026-04-08 19:56 by 振—TZ
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 284求调剂 +17 梵@@ 2026-04-06 17/850 2026-04-08 11:35 by 1shin_ichi
[考研] 本科郑州大学，一志愿华东师范大学282求调剂 +21 熊哥xtk 2026-04-07 24/1200 2026-04-08 09:48 by 123考研123
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 0854求调剂 +4 assdll 2026-04-05 4/200 2026-04-06 12:29 by 中飞院空管学院�
[考研] 工科08-机械专硕-求调剂 +3 雷欧飞踢 2026-04-02 3/150 2026-04-05 18:49 by 蓝云思雨
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
[考研] 278求调剂 +3 依旧！ 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 一志愿0817化学工程与技术，求调剂 +24 我不是只因 2026-04-02 28/1400 2026-04-04 15:15 by dongzh2009
[考研] 357求调剂 +13 1050389037 2026-04-03 13/650 2026-04-03 22:27 by 无际的草原
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭