| 查看: 2356 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
baoluolove铜虫 (小有名气)
|
[求助]
请问编程计算配位数时如何处理周期性边界条件？ 已有1人参与
|
|||
| 自己按照老虎大王编程思路尝试编程计算了配位数，可是发现和文献有一定的出入，想想很有可能是没有处理周期性边界条件，请问专家们周期性边界条件如何在编程中进行处理的，十分感谢~ |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在有不均一介质的电场中电场是怎样分布的
已经有4人回复
- 请教 径向分布函数（ＲＤＦ） 问题
已经有12人回复
ChemiAndy
木虫 (正式写手)
- 模拟EPI: 40
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
5楼2014-07-16 02:10:17
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
baoluolove: 金币+5, ★★★很有帮助 2014-07-12 08:39:37
baoluolove: 金币+5, ★★★★★最佳答案, 已经搞定啦~ 2014-07-12 09:40:07
月只蓝: 金币+3, 应助指数+1, 感谢指导！ 2014-07-16 16:34:48
感谢参与，应助指数 +1
baoluolove: 金币+5, ★★★很有帮助 2014-07-12 08:39:37
baoluolove: 金币+5, ★★★★★最佳答案, 已经搞定啦~ 2014-07-12 09:40:07
月只蓝: 金币+3, 应助指数+1, 感谢指导！ 2014-07-16 16:34:48
|
Assuming the simulation box lengths along x,y,z directions are lx,ly,lz, you want to calculate the distance between particle i and j while removing the periodic boundary conditions This is what I did in Fortran xij = xi - xj - nint((xi-xj)/lx)*lx yij = yi - yj - nint((yi-yj)/ly)*ly zij = zi - zj - nint((zi-zj)/lz)*lz rij = sqrt(xij*xij + yij*yij + zij*zij) nint is a fortran function that gives the nearest integer of a real number. I'm sure you can find similar functions in C++ or other programing languages. |
» 本帖已获得的红花（最新10朵）
baoluolove3楼2014-07-12 00:58:19
baoluolove
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 105.8
- 红花: 2
- 帖子: 134
- 在线: 43.4小时
- 虫号: 2939613
- 注册: 2014-01-19
- 性别: MM
- 专业: 玻璃材料
4楼2014-07-12 08:38:44
baoluolove
铜虫 (小有名气)
 |
6楼2014-07-17 10:54:44
投票:
10