| 查看: 963 | 回复: 12 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xiaowenzhi16铁杆木虫 (正式写手)
|
[交流]
【求助】如何显示3D spin density
|
||
|
请问VASP里的自旋密度如何显示为3D spin density 的电子云状的图像?用什么软件可以做到呢？？ 谢谢 |
回复此楼» 猜你喜欢
生物学308求调剂（一志愿华东师大）
已经有6人回复
-
调剂求收留
已经有7人回复
-
272分材料子求调剂
已经有36人回复
-
275求调剂
已经有8人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有20人回复
-
070300化学学硕311分求调剂
已经有19人回复
-
材料与化工调剂
已经有13人回复
-
材料与化工调剂
已经有33人回复
-
复试调剂
已经有7人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有17人回复
xiaowenzhi16
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 7096.7
- 散金: 50
- 红花: 1
- 帖子: 564
- 在线: 311.7小时
- 虫号: 523628
- 注册: 2008-03-12
- 性别: GG
- 专业: 半导体物理
3楼2009-05-13 22:19:11
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
4楼2009-05-14 09:21:53
xiaowenzhi16
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 7096.7
- 散金: 50
- 红花: 1
- 帖子: 564
- 在线: 311.7小时
- 虫号: 523628
- 注册: 2008-03-12
- 性别: GG
- 专业: 半导体物理
5楼2009-05-14 10:14:22
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xiaowenzhi16(金币+6,VIP+0):实在是太感谢啦。以后还要麻烦你的。 5-14 12:00
wuchenwf(金币+10,VIP+0):谢谢 5-14 23:01
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)"  scaledo 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,zend 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 ] |
6楼2009-05-14 11:25:56
5