| 查看: 2321 | 回复: 10 | ||||
baoluolove铜虫 (小有名气)
|
[求助]
请问编程计算配位数时如何处理周期性边界条件？ 已有1人参与
|
|||
| 自己按照老虎大王编程思路尝试编程计算了配位数，可是发现和文献有一定的出入，想想很有可能是没有处理周期性边界条件，请问专家们周期性边界条件如何在编程中进行处理的，十分感谢~ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有16人回复
-
我不理解！
已经有14人回复
-
国自然面上复盘~欢迎讨论
已经有4人回复
-
系统查不到
已经有10人回复
-
基金不中，共勉
已经有3人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有9人回复
-
2026国自然函评费到账
已经有25人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有28人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在有不均一介质的电场中电场是怎样分布的
已经有4人回复
- 请教 径向分布函数（ＲＤＦ） 问题
已经有12人回复
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
- 专业: 玻璃材料
2楼2014-07-11 22:13:59
baoluolove
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 105.8
- 红花: 2
- 帖子: 134
- 在线: 43.4小时
- 虫号: 2939613
- 注册: 2014-01-19
- 性别: MM
- 专业: 玻璃材料
4楼2014-07-12 08:38:44
ChemiAndy
木虫 (正式写手)
- 模拟EPI: 40
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
5楼2014-07-16 02:10:17
baoluolove
铜虫 (小有名气)
 |
6楼2014-07-17 10:54:44
baoluolove
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 105.8
- 红花: 2
- 帖子: 134
- 在线: 43.4小时
- 虫号: 2939613
- 注册: 2014-01-19
- 性别: MM
- 专业: 玻璃材料
|
你好，我之前用你的这种方法计算配位数，实现了我的目标，我现在同样用你的这个转换方法算键角分布函数，却与没用周期性边界条件的结果是一样的，你知道这是这怎么一回事吗 |
» 本帖已获得的红花（最新10朵）
7楼2014-10-08 21:20:13
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
8楼2014-10-09 06:20:25
15216105346
金虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2062.4
- 散金: 2693
- 红花: 7
- 帖子: 1205
- 在线: 298小时
- 虫号: 3441372
- 注册: 2014-09-25
- 专业: 固体无机化学
9楼2015-01-31 22:22:43
10楼2017-03-21 15:00:25
投票:
