【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 2461  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

[求助] 求助Dmol3的格点文件grd格式说明已有1人参与

鄙人不才,看了一下MS的help文档中对grd文件的介绍和例子:
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
前五行我大致知道什么意思了,不过第五行所谓的"The fastest varying index "是什么意思啊?
另外下面的第六行开始的列数据是按照怎样的顺序排列?到底是怎样与三维空间对应?
非立方格点有没有办法或是程序转换成cube或CHGCAR这样的格式?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

不一定挂在论坛,计算问题问题欢迎发帖或留言。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
2楼: Originally posted by sobereva at 2014-06-16 11:57:46
The fastest varying index通常为1,代表x方向的指数变化是最快的,也就是说,读取下面格点数据值时最内层循环的是x序号,例如这样读取
do k=1,nz
        do j=1,ny
                do i=1,nx
                        read(10,*) cubmat(i,j,k)
                end  ...

前辈,现在multiwfn处理grd文件的AIM是不是
还不能支持PBC啊?

[ 发自手机版 http://muchong.com/3g ]
不一定挂在论坛,计算问题问题欢迎发帖或留言。
3楼2014-06-16 12:38:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+5 2014-06-16 14:18:30
卡开发发: 金币+10, ★★★★★最佳答案 2014-06-16 15:01:38
The fastest varying index通常为1,代表x方向的指数变化是最快的,也就是说,读取下面格点数据值时最内层循环的是x序号,例如这样读取
do k=1,nz
        do j=1,ny
                do i=1,nx
                        read(10,*) cubmat(i,j,k)
                end do
        end do
end do

Multiwfn程序的fileIO.f90里subroutine readgrd就是读取.grd文件的代码,可供参考

非立方格点情况的.grd文件原则上可以转换为.cub,只要适当定义平移矢量即可
2楼2014-06-16 11:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


引用回帖:
3楼: Originally posted by 卡开发发 at 2014-06-16 12:38:50
前辈,现在multiwfn处理grd文件的AIM是不是
还不能支持PBC啊?
...

Multiwfn载入grd后可以做AIM盆分析(也包括寻找(3,-3)或(3,+3)临界点),但是不能做拓扑分析。
目前并没有尝试过对周期性体系的grd文件做盆分析,虽然可能能做
4楼2014-06-16 14:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
4楼: Originally posted by sobereva at 2014-06-16 14:27:22
Multiwfn载入grd后可以做AIM盆分析(也包括寻找(3,-3)或(3,+3)临界点),但是不能做拓扑分析。
目前并没有尝试过对周期性体系的grd文件做盆分析,虽然可能能做...

前辈,之前我试了一下grd做盆分析,非pbc倒是没啥问题,pbc盆积分得到的总电荷总数不太对。我先看看我自己有没有办法处理。

多谢前辈指导。

[ 发自手机版 http://muchong.com/3g ]
不一定挂在论坛,计算问题问题欢迎发帖或留言。
5楼2014-06-16 15:01:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +9 孙大大@ 2024-04-17 9/450 2024-04-19 21:50 by 刘国宁
[硕博家园] +5 s150535912 2024-04-18 6/300 2024-04-19 20:33 by s150535912
[论文投稿] 最近遇到这样一个问题 2+4 asd123gfa689 2024-04-18 10/500 2024-04-19 18:57 by asd123gfa689
[基金申请] 特别资助审核状态 +7 pantray 2024-04-17 7/350 2024-04-19 13:51 by Breezy恋
[论文投稿] 求助 100+3 guhangyu 2024-04-15 3/150 2024-04-19 08:48 by bobvan
[基金申请] 国自然青年基金只有一篇能第一标注,其他都是第二标注能结题吗? 10+4 lice_1987 2024-04-18 6/300 2024-04-19 08:29 by kudofaye
[有机交流] 兄弟们帮我看看这两个结构怎么合成 +5 xl2088131 2024-04-17 5/250 2024-04-19 08:09 by 西瓜霜华
[基金申请] 国家资助博士后BC档出校后资助的概率多大? +3 卡卡罗特哦 2024-04-16 3/150 2024-04-18 12:58 by wolfgangHugh
[考研] 一志愿北化085600 328求线上调剂 +10 Mrshark 2024-04-13 10/500 2024-04-18 08:30 by 刘国宁
[有机交流] 怎么清洗烧瓶 20+5 ww34523 2024-04-16 6/300 2024-04-17 15:20 by 591950582
[考研] 322求调剂 +7 本己上岸 2024-04-16 7/350 2024-04-17 11:49 by duanxz
[论文投稿] MDPI初审7天,是不是被拒 20+6 applegirl5762 2024-04-15 8/400 2024-04-17 08:49 by bobvan
[考博] 25年申博求助,射频方向 +4 SherlockAH 2024-04-15 10/500 2024-04-16 18:06 by ljysudaee
[考研] 求调剂 +4 鹿萌月yy 2024-04-15 5/250 2024-04-16 16:16 by chenweiwade
[考研] 320求调剂 +5 陆志伟 2024-04-15 5/250 2024-04-16 11:11 by 19862091
[考研] 广州大学光电信息工程专业调剂,招收物理学专业学生 +5 txhx4010 2024-04-14 7/350 2024-04-16 10:52 by domax
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[考研] 求调剂 +4 桃岸雪 2024-04-15 5/250 2024-04-15 18:49 by mthwyj
[论文投稿] with efitor 越久是不是越容易拒稿。我的已经一个多月了 +5 lizhengke06 2024-04-14 5/250 2024-04-15 18:33 by jonewore
[考研] 273求调剂 +5 Late婉安 2024-04-15 7/350 2024-04-15 13:01 by Late婉安
信息提示
请填处理意见