| 查看: 4533 | 回复: 98 | |||
| 本帖产生 1 个 模拟EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】DL_POLY的config文件如何生成?? 已有9人参与
|
|||
|
我想模拟聚合物在金属团簇表面的吸附,想请教大家这样的体系的config文件如何生成? 另外我把手册已经看了一遍,但是希望运行一下.有么有人根据手册的算例给写一个模拟过程,给传一个谢谢!! |
回复此楼» 猜你喜欢
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有9人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
|
您好,又来麻烦了.. 因为HISTORY的坐标不能直接用来统计分析(从前面的几个截图可以看出来,有的分子是断开的). 所以需要一定的转换吧,下面这段程序是在官方论坛上面有人这样写的,虽然讲的很清楚,可是我不明白他的算法是什么意思? 大家能给解释一下吗?谢谢!! here is a way to generate image flags (or unwrapped coordinates) from the HISTORY file which contains the wrapped coordinates: Assumptions: 1. All the atoms listed in the first configuration listed in the HISTORY file are assumed to be inside the box 2. The time between successive configurations in the HISTORY file is small enough for the atoms to not cross the simulation box boundaries more than once 3. Simulation box is cuboid, should be easily modifiable for other geometries Procedure: 1. Initialize the image flags ix, iy and iz in the x, y and z directions to 0, i.e. Code:ix = 0; iy = 0; iz = 0; where i = index of atoms going from 1 to nAtoms 2. Read the coordinates of all the atoms in the first configuration (config_old) and label them x_old, y_old and z_old 3. Read the coordinates of all the atoms from the next configuration (config_new) and label them x_new, y_new and z_new 4. The image flags for the configuration denoted as 'config_new' are then given as follows Code: ix = ix+int(2*(x_old-x_new)/L_x); iy = iy+int(2*(y_old-y_new)/L_y); iz = iz+int(2*(z_old-z_new)/L_z); where, int(k) = the integer part of the floating point number k L_x = Simulation box length in the x direction L_y = Simulation box length in the y direction L_z = Simulation box length in the z direction Note 1. If the box were cubic, L_x = L_y = L_z = L 2. Do not reset the image flags to zero, just keep on incrementing them 5. The unwrapped coordinated for the configuration denoted as 'config_new' are then given as Code: x_uw = x_new+L_x*ix y_uw = y_new+L_y*iy z_uw = z_new+L_z*iz 6. You can then use the unwrapped coordinates to fine the diffusion, etc. |
78楼2010-04-08 15:17:20
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-01-24 16:53:03
3楼2010-01-24 20:21:22
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2010-01-25 09:45:24