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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何显示3D spin density
131/2返回列表12下一页
查看: 959  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +8 不愿透露姓名的� 2026-04-02 8/400 2026-04-08 20:12 by luoyongfeng
[考研] 269求调剂 +6 啊啊我我 2026-04-07 6/300 2026-04-08 20:04 by 我减肥1
[考研] 328求调剂 +16 lftmya 2026-04-07 17/850 2026-04-08 18:45 by 环化材-小生
[考研] 284求调剂 +17 梵@@ 2026-04-06 17/850 2026-04-08 11:35 by 1shin_ichi
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 材料调剂 +12 一样YWY 2026-04-02 14/700 2026-04-07 12:52 by Sammy2
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 第一志愿东南大学物理313，有科研竞赛获奖经历，希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 302分 085601求调剂推荐 +11 zyx上岸！ 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 341求调剂 +3 学无止境，冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 324求调剂 +14 想上学求调 2026-04-02 15/750 2026-04-04 20:31 by 无际的草原
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 能源动力 调剂 +3 不破不立0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[考研] 377求调剂 +3 RASKIN 2026-04-02 3/150 2026-04-02 09:45 by zzchen2000
信息提示
关闭
请填处理意见
关闭