24小时热门版块排行榜    

查看: 1903  |  回复: 13
【奖励】 本帖被评价6次,作者csfn增加金币 5
当前主题已经存档。

csfn

荣誉版主 (知名作家)


[资源] CHGCAR三维图的做法(于ms中实现)

记得前面有虫子提过关于dmol电荷密度解释的问题。另外,有一些朋友也想解决对vasp的电荷三维化的问题。
这里,帖出一个关于dmol电荷密度例子的解释。
想解决vasp三维化问题的朋友,也可以考虑下。只要把chgcar用程序简单转化成grd格式,即可以在ms中进行三维化观察,非常漂亮,方便.
Example GRD file:
Each .grd file begins with a title, which is the descriptive of the volumetric field. A typical title would be "DMol3 total electron density".

The second line of each file contains a FORTRAN format descriptor of the data as written by the program that generated it. This format is typically (1p,e12.5).

The third line specifies the volume of space containing the data to be plotted: the cell lengths a, b, and c in angstroms, followed by the cell angles α, β, and γ. The unit cell is the size of the grid used for plotting. By convention, the first axis, a, is taken to be parallel to the Cartesian x axis.

The fourth line indicates the number of grid spaces in each direction, a, b, and c.

The fifth line contains seven integers indicating which index varies the fastest and the location of the Cartesian origin relative to the edges of the volumetric data.

The first integer in the fifth line must be 1, corresponding to a, or 3, corresponding to b. The fastest varying index cannot be c.

The next six integers in the fifth line indicate the number of grid spaces to the left and right of the origin in directions a, b, and c.

The subsequent lines contain the actual data. If the number of spaces indicated on line 4 are Na, Nb, and Nc, then there will a total of (Na+1) x (Nb+1) x (Nc+1) data points.

Below is a sample .grd file.

DMol3 total electron density            
(1p,e12.5)  
   9.000   8.200   9.400  90.000  90.000  90.000
   45   41   47
    1  -23   22  -20   21  -22   25
0.00000E+00
0.00000E+00
2.57683E-13
1.43774E-12
2.01988E-10
8.71334E-10
1.28986E-09
1.11599E-09
4.72026E-10
1.59762E-11
2.83878E-12
...

[ Last edited by wuchenwf on 2009-6-13 at 23:08 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)


★★★★★ 五星级,优秀推荐

不错 正要学习呢
2楼2008-05-16 20:09:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)



wuchenwf(金币+1):thank you
CHGCAR直接用vaspview画就可以
想好看一点,用v2xsf+Xcrysden处理好了。
3楼2008-05-16 23:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

q68

铁杆木虫 (著名写手)


★★★★★ 五星级,优秀推荐

就喜欢这种经验交流!
4楼2008-05-17 12:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ztwang

铜虫 (小有名气)


引用回帖:
Originally posted by y1ding at 2008-5-16 23:05:
CHGCAR直接用vaspview画就可以
想好看一点,用v2xsf+Xcrysden处理好了。

下了一个vaspview,但总是编译不成功CINCLUDE  LIBS,怎么设?谢谢!!
gcc /usr/X11R6/bin  -c -o ds3.o ds3.c
In file included from ds3.c:18:
ds3.h:26:21: error: GL/glut.h: No such file or directory
In file included from glw.h:24,
                 from ds3.h:29,
                 from ds3.c:18:
5楼2008-05-20 09:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ypf1449

银虫 (初入文坛)


把chgcar用程序简单转化成grd格式???
怎么转化啊??
能不能把步骤具体一点啊啊??
谢谢了!!!
6楼2008-06-24 16:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylin78

金虫 (小有名气)


0.5

能具体介绍下怎么用Xcrysden处理么??
7楼2008-06-24 16:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cationly

木虫 (正式写手)



wuchenwf(金币+1):thank you 1-27 16:04
引用回帖:
Originally posted by ztwang at 2008-5-20 09:42:


下了一个vaspview,但总是编译不成功CINCLUDE  LIBS,怎么设?谢谢!!
gcc /usr/X11R6/bin  -c -o ds3.o ds3.c
In file included from ds3.c:18:
ds3.h:26:21: error: GL/glut.h: No such file or direct ...

这个估计是缺少opengl库。
8楼2008-06-24 19:48:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)


★★★ 三星级,支持鼓励

是啊,怎么转换VASP的chgcar 呢?
9楼2008-06-24 19:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuxinxing1

铁杆木虫 (著名写手)


推荐


wuchenwf(金币+1):thank you 1-27 16:04
引用回帖:
Originally posted by kylin78 at 2008-6-24 16:10:
能具体介绍下怎么用Xcrysden处理么??

http://www.xcrysden.org/

此软件,不是很复杂; 可以免费下载使用。
10楼2008-07-12 03:24:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)


★★★ 三星级,支持鼓励

好的经验呀,有大人编写一个Fortran或者Matlab程序进行转换的没有呀!?
11楼2008-08-04 11:20:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 csfn 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +52 医学老男孩 2026-08-13 115/5750 2026-08-19 12:09 by 医学老男孩
[基金申请] filecode,4个jtjc了 +7 ziyangfang 2026-08-19 8/400 2026-08-19 11:40 by h初学者f
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 今天放榜吗? +11 布布和一二 2026-08-19 11/550 2026-08-19 10:06 by xter9665
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见