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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何显示3D spin density
北京石油化工学院2026年研究生招生接收调剂公告
131/2返回列表12下一页
查看: 957  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 353求调剂 +6 晴空万里air 2026-04-07 6/300 2026-04-08 10:06 by 不吃魚的貓
[考研] 326分，一志愿沪9，求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 328求调剂 +13 lftmya 2026-04-07 14/700 2026-04-07 22:45 by JourneyLucky
[考研] 调剂 +18 不逢春 2026-04-05 19/950 2026-04-07 22:04 by lijunpoly
[考研] 326求调剂 +5 9ahye 2026-04-02 6/300 2026-04-07 21:37 by lijunpoly
[考研] 085602调剂 初试总分335 +10 19123253302 2026-04-05 10/500 2026-04-07 15:23 by 小乔同学ya
[考研] 085600材料与化工301分求调剂院校 +25 刺痛jk 2026-04-06 26/1300 2026-04-07 15:07 by -迷了路啊路
[考研] 081700，311，求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 生物学求调剂 一志愿沪9，326分 +6 刘墨墨 2026-04-06 6/300 2026-04-06 19:36 by lijunpoly
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 312求调剂 +6 小小墨123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
信息提示
关闭
请填处理意见
关闭