[交流] 【求助】DLPOLY计算出来的结果，是否要编程才能分析？

» 猜你喜欢

为什么资助数各大高校都创新高，自己申请怎么就这么难已经有15人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
基金系统什么内容也没有已经有10人回复
学科评审组评审是指会评吗？已经有5人回复
梦想已经有8人回复
面上合作单位盖章已经有9人回复
面上函评意见出来了，像什么等级？已经有17人回复
小白求助投论文要求的highlights应该如何写已经有4人回复
麻烦专家们看看评委们的意见（F口面上）已经有12人回复

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

光学涡旋计算全息图的matlab程序已经有9人回复
求助循环伏安谱图分析已经有11人回复
Fluent计算时的迭代问题已经有14人回复
调和分析---著名数学家沃尔夫奖获得者 Stein 名著已经有92人回复
在XRD中如何确定哪些是超晶格峰，超晶格的峰位能通过计算算出来吗？已经有12人回复
根据某种蛋白质的氨基酸分析结果计算htot（水解度要用的常数值）已经有11人回复
求助关于用vasp计算二维石墨烯纳米带所得到的能带结果与报道对不上已经有7人回复
[原创]一种基于能带计算结果的载流子有效质量计算方法已经有174人回复
晚辈诚惶诚恐请教：制甲醇，在线色谱产物分析计算方法？已经有7人回复
计算结果显示导带底不在特殊K点上，这种情况如何解释已经有11人回复
ATK分子器件电子输运理论计算，结果分析的经典文献推荐已经有126人回复
用c语言编程怎来表达时滞！谢谢！已经有6人回复
matlab 积分计算算不出结果，求帮助！！！已经有10人回复
有三个计算物理的题，大家谁能帮帮忙啊，要求编程序。我弄了好几天了，弄不出来，恳请已经有18人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-01-12 15:57:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

043114076

木虫 (正式写手)

应助: 49 (小学生)
金币: 3364.9
帖子: 788
在线: 369.5小时
虫号: 570592

★ ★
zyj8119(金币+1):谢谢参与
zyj8119(金币+10):应该怎么做呢？ 2011-01-12 19:56:11
ghcacj(金币+1):谢谢 2011-01-13 11:43:50

vmd 自身带的分析功能可以实现

2楼2011-01-12 16:47:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

winterz

木虫 (小有名气)

模拟EPI: 1
应助: 0 (幼儿园)
金币: 3774.3
帖子: 81
在线: 918.1小时
虫号: 935931

★ ★ ★ ★
zyj8119(金币+1):谢谢参与
zyj8119(金币+18):谢谢兄弟。。 2011-01-12 19:36:41
ghcacj(金币+3):谢谢 2011-01-13 11:43:57

引用回帖:

Originally posted by zyj8119 at 2011-01-12 15:57:55:
计算出来的HISTORY文件，要怎么才能得到RDF,MSD等信息？

MSD使用自带的Java GUI和msdprg.f都可以计算。
RDF只要在CONTROL文件中写入以下两行，就会生成一个RDFDAT文件，里面是各原子RDF信息。
rdf 100
print rdf

3楼2011-01-12 19:27:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by 043114076 at 2011-01-12 16:47:18:
vmd 自身带的分析功能可以实现

VMD应该怎么做？RDF和MSD都能分析吗？

4楼2011-01-12 19:37:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

043114076

木虫 (正式写手)

应助: 49 (小学生)
金币: 3364.9
帖子: 788
在线: 369.5小时
虫号: 570592

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+2):谢谢 2011-01-13 11:44:05

引用回帖:

Originally posted by zyj8119 at 2011-01-12 19:37:32:

VMD应该怎么做？RDF和MSD都能分析吗？

rdf可以做，vmd可以做rmsd。自带的有工具

5楼2011-01-12 20:26:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by 043114076 at 2011-01-12 20:26:52:

rdf可以做，vmd可以做rmsd。自带的有工具

但是RDF选项里面的selection 1和selection 2怎么办？

6楼2011-01-12 20:33:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by winterz at 2011-01-12 19:27:24:

MSD使用自带的Java GUI和msdprg.f都可以计算。
RDF只要在CONTROL文件中写入以下两行，就会生成一个RDFDAT文件，里面是各原子RDF信息。
rdf 100
print rdf

madprg.f是否是DLPOLY的源代码？还是在某处下到的源代码文件？

7楼2011-01-12 20:39:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

winterz

木虫 (小有名气)

模拟EPI: 1
应助: 0 (幼儿园)
金币: 3774.3
帖子: 81
在线: 918.1小时
虫号: 935931

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+3):谢谢 2011-01-13 11:44:16

引用回帖:

Originally posted by zyj8119 at 2011-01-12 20:39:55:

madprg.f是否是DLPOLY的源代码？还是在某处下到的源代码文件？

在源码包中的utility目录里有很多与DLPOLY相关的工具程序，其中msprg.f、msdblk.f、msdblk2.f都是计算MSD的程序，但是DLPOLY的手册中对这些工具的介绍很粗略，一般都需要读懂它们输入部分的源代码才能知道怎么使用。

8楼2011-01-12 20:49:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

043114076

木虫 (正式写手)

应助: 49 (小学生)
金币: 3364.9
帖子: 788
在线: 369.5小时
虫号: 570592

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+2):谢谢 2011-01-13 11:44:24

引用回帖:

Originally posted by zyj8119 at 2011-01-12 20:33:54:

但是RDF选项里面的selection 1和selection 2怎么办？

记得好像是写原子名，你可以看一下vmd的说明书，或者自己试试

9楼2011-01-12 21:46:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by winterz at 2011-01-12 20:49:34:

在源码包中的utility目录里有很多与DLPOLY相关的工具程序，其中msprg.f、msdblk.f、msdblk2.f都是计算MSD的程序，但是DLPOLY的手册中对这些工具的介绍很粗略，一般都需要读懂它们输入部分的源代码才能知道怎么 ...

貌似这些也是需要编译的，使用ifort编译，然后把得到的那个编译好的文件，放到history文件夹里面。运行那个可执行文件？

10楼2011-01-12 22:16:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

貌似ifort编译不了这个文件啊，我输入ifort msdprg.f -o msd.exe ，貌似什么也没有。ifort是我路径都正确，用gfortran也不行，这个文件应该怎么用呢？

11楼2011-01-12 22:28:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by 043114076 at 2011-01-12 21:46:12:

记得好像是写原子名，你可以看一下vmd的说明书，或者自己试试

说明书写的不清楚，我先mouse->label->atoms,然后graphics->labels->貌似什么用也没有。

12楼2011-01-12 22:30:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by zyj8119 at 2011-01-12 22:28:48:

貌似ifort编译不了这个文件啊，我输入ifort msdprg.f -o msd.exe ，貌似什么也没有。ifort是我路径都正确，用gfortran也不行，这个文件应该怎么用呢？

终于把msdprg.f编译成了msdprg.exe,然后把HISTORY拷贝到一起，运行
./msdprg.exe

13楼2011-01-12 23:10:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

winterz

木虫 (小有名气)

模拟EPI: 1
应助: 0 (幼儿园)
金币: 3774.3
帖子: 81
在线: 918.1小时
虫号: 935931

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+10):谢谢 2011-01-13 11:43:23

引用回帖:

Originally posted by zyj8119 at 2011-01-12 23:10:11:

终于把msdprg.f编译成了msdprg.exe,然后把HISTORY拷贝到一起，运行
./msdprg.exe

不要把HISTORY改名，要重新编辑一个msd_input文件，按照msdprg.f里要求的格式和你模拟所设置的参数输入变量：

=======================================================================
c    read the control variables

c    name of atomic species for autocorrelation
   read(5,*)atname

c    name of selected HISTORY file
   read(5,*)fname

c    HISTORY file formatted or unformatted
   read(5,*)fmt

c    total number of atoms in a HISTORY file configuration
   read(5,*)natms

c    time interval between successive configurations
   read(5,*)tstep

c    max number of configurations to sample
   read(5,*)nconf

c    working length of msd array
   read(5,*)nmsd

c    sampling interval
   read(5,*)isampl

c    interval between origins
   read(5,*)iomsd

c    average cell vectors (taken from dl_poly OUTPUT file)
   read(5,*)avcell
=========================================================================

然后直接执行./msdprg.exe，会生成一个msd_output

以下是一个msd_input样例

=========================================================================
C
HISTORY
formatted
720
1.0
5000
500
10
1
25.345 0.00 0.00
  0.00  25.345 0.00
  0.00 0.00  25.345
=========================================================================

14楼2011-01-13 08:57:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:

Originally posted by winterz at 2011-01-13 08:57:50:

不要把HISTORY改名，要重新编辑一个msd_input文件，按照msdprg.f里要求的格式和你模拟所设置的参数输入变量：

=======================================================================
c read the ...

牛人啊，膜拜！！！

15楼2011-01-13 10:08:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

meixinqing

金虫 (著名写手)

应助: 46 (小学生)
贵宾: 0.122
金币: 5691.4
帖子: 1923
在线: 369.4小时
虫号: 280674

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+5): 谢谢 2011-02-05 20:56:26

可以用vmd分析

load的时候选择DLPOLY HISTORY
然后在graphical representations 里的selection 里看name 知道你要的原子名

然后在extension -- analysis --rdf 里，selection写name 原子名

eg

name O
name H

不过不太喜欢用vmd的，精度不够高，设置的delta r小，出来的图就怪怪的。。。。DL-POLY可以设置输出RDF，精度可以设置的很高