24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1759  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见