| 查看: 5450 | 回复: 98 | |||
| 本帖产生 1 个 模拟EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】DL_POLY的config文件如何生成?? 已有9人参与
|
|||
|
我想模拟聚合物在金属团簇表面的吸附,想请教大家这样的体系的config文件如何生成? 另外我把手册已经看了一遍,但是希望运行一下.有么有人根据手册的算例给写一个模拟过程,给传一个谢谢!! |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有16人回复
-
40+中年人失业后，还能找到好工作吗？
已经有16人回复
-
人工智能虚构史料、伪造文献频现，多家期刊声明抵制
已经有4人回复
-
CsPbI2Br全无机钙钛矿SEM图像异常分析
已经有8人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有9人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
|
您好,又来麻烦了.. 因为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