DL_POLY2.16使用速成

分子模拟软件DL_POLY速成指导

DL_POLY2.16使用速成

一 软件的编译

1 安装visual c 6.0 对安装目录无要求，必须注意安装最后阶段，当问及是否在路径中设置环境变量时，一定选择是。

2 安装Compaq Visual Fortran 6.5 对安装目录无要求，必须注意安装最后阶段，当问及是否在路径中设置环境变量时，一定选择是。

3 在windowsxp下安装linux 模拟环境 Cygwin，只能安装在C盘跟目录下 。具体安装办法参看installcygwin.doc。很多人选择默认全部安装，因为cygwin的默认安装很多组建是不安装的。所以一定请参看我给的安装方法安装cygwin。

以上步骤1，2可以相互交换次序，但是cygwin的安装必须在前两步完成后才可以进行。不过请先看installcygwin.doc下载cygwin，很慢的。

4 进入C:/cygwin/bin/找到link.exe将其改名为link.exe.bak

5 在以上各步骤完成的基础上，把dlpoly拷贝到cygwin目录下，然后把dlpoly/build目录下的makewin拷贝到srcf90目录下改名为makefile。

6 打开cygwin的模拟终端（即双击winxp下的那个可执行程序cygwin.exe）转到srcf90下，命令如下

cd c:/cygwin/dlpoly216/srcf90 ！该命令假设所有dlpoly的文件夹都放在dlpoly216文件夹内。 然后在命令提示符后输入如下命令

make Compaq-f90

回车即开始编译，在前述操作无误，且C, fotran编译器设置正确的情况下通过后会在execute目录下产生DLPLOY.X可执行程序。以后运行时候只需要把DLPOLY.X拷贝到你提交运行任务的目录下在终端窗口输入如下指令

./DLPOLY.X

回车，即开始运算，如果有错误程序会终止，在OUTPUT文件里写有错误代号，可以从使用手册里面查找以确定出错原因。

注意：Linux操作系统区分大小写，所以字母大写小写不代表同样的文件或者操作命令，使用者可以在网上随便找到简单的linux操作指令。

二 软件包的运行

1 请在dlpoly216/下建立相关的运算目录比如Si表示该目录计算Si,且在该目录下准备好常用的输入文件，共三个（所有输入输出文件均可用记事本打开修改，我用的UltraEditer）

(1) CONFIG文件

注意，文件名一定要大写。文件格式如下

Si sample Config file ！文件头，字符形变量，长度100

0 2 ！整形变量，0表示输入文件仅包含坐标2表示长方体正交边界条件

12.882225867378 0.000000000000 0.000000000000 ！模拟材料长的三个分量

0.000000000000 12.882225867378 0.000000000000 ！模拟材料宽的三个分量

0.000000000000 0.000000000000 12.882225867378 ！模拟材料高的三个分量

Si 1 ！原子名称与序数，注意名称要与FIELD文件中的名称一致

-5.997878668000 -5.951354114000 -5.974887095000 ！原子坐标

Si 2

-2.744948035000 -4.885863816000 -4.846913142000

。

。

分子模拟软件DL_POLY速成指导

。

Si 215

4.878605683000 2.717825432000 4.875233647000

Si 216

5.980280337000 5.956520906000 6.010595727000

具体的文件格式参看手册，注意：原点必须取在模拟材料的几何中心

（2）CONTROL文件

以NPT系综为例，格式如下

Si nanoware with Tersoff potential ！文件头100长度字符形变量

integrate leapfrog ！积分采用蛙跳法

pressure 0.0 ！压强设置为0

temperature 200.0 ！温度设置为200K

ensemble npt ber 0.5 0.5 ！采用npt系综调压调温采用beresend调压调温器

steps 8010 ！总步数为8010步，这样设置是为了输出RDF

equilibration 8000 ！其中平衡用8000步多余10步不调节用于输出RDF

timestep 0.001 ！每步一飞秒

scale 1 ！每一步进行调温调压

cut 3.1 ！截断半径

delr 0.5 ！径向分布函数每次搜索步长0.5Å

no elec ！无电荷

no vdw ！无范得瓦而斯作用

print every 10 ！每10步输出

stats 10 ！每10步统计

rdf 1 ！每1步统计rdf

print rdf ！单另输出RDF，其实OUTPUT文件最后也有的

collect ！收集

traj 1 500 0 ！HISTORY输出指令从第一步开始每500步输出，只输出坐标

job time 999999.9 ！总时间

close time 99.9 ！关闭时间

belt ！禁止样品转动

finish ！到此控制文件结束

具体格式参看手册

（3）FIELD文件（以硅为例）

Si nanoware WITH TERSOFF POTENTIAL ！文件头

UNITS ev ！能量单位

MOLECULES 1 ！一种原子

Si ！名字叫硅

NUMMOL 1 ！一种分子

ATOMS 216 ！原子总数216

Si 28.0855 0.00 216 ！原子名称，质量，电荷，数目

FINISH ！该原子数据完

TERSOFF 1 ！tersoff势一种，以下为参数

Si ters 1830.8 2.4799 471.18 1.7322 2.7

3.0 1.1000e-6 0.78734 100390.0 16.217 -0.59825

分子模拟软件DL_POLY速成指导

Si Si 1.0 1.0

CLOSE

具体势参数需要根据选择的经验势的不同仔细参看手册

以上三个文件为基本输 …… 此处隐藏：2394字，全部文档内容请下载后查看。喜欢就下载吧 ……