24小时热门版块排行榜    

查看: 1904  |  回复: 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的回帖

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的回帖
查看全部 14 个回答

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的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 12/600 2026-08-19 13:45 by gloomy6159
[基金申请] filecode,4个jtjc了 +8 ziyangfang 2026-08-19 9/450 2026-08-19 13:08 by blacksandy
[基金申请] 今天基金会出结果吗?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 医学老男孩
[基金申请] 朋友圈看到的 +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 启萌科技
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +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
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +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
信息提示
请填处理意见