24小时热门版块排行榜    

查看: 1829  |  回复: 15

jianjunjiang

铁杆木虫 (正式写手)

[求助] 请问如何计算某个组合(比如说水分子)的msd 已有2人参与

各位大侠:
     请问如何利用分子动力学计算某个组合的msd,比如说对于水分子,是不是先利用水分子中氢原子和氧原子的坐标和质量首先算出水分子的质心,然后计算该质心的msd即可,期待您的答复!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
jianjunjiang: 金币+3, ★★★很有帮助, 谢谢您的答复!! 2016-08-07 10:30:32
我感觉差别不会很大吧。你把两条线画在一起比比看,差别很大吗?你说的这两个定义在数学上应该是不等价的,但我觉得区别应该很小。既然gromacs的结果和你的其中一种定义的结果相同,就说明它很可能用了你说的定义。也许手册中会有介绍?

[ 发自手机版 http://muchong.com/3g ]
7楼2016-08-07 03:36:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MS_Z

金虫 (初入文坛)

同意brucefan 我个人觉得有一点差别不一定需要太纠结 如果真的想知道差别 建议直接看看gromacs内部msd的算法 有可能他给出的msd是用了诸如多时间起点等平滑算法

发自小木虫IOS客户端
8楼2016-08-07 10:12:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

brucefan

专家顾问 (著名写手)

2楼2016-08-04 14:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by brucefan at 2016-08-04 14:59:53
yes.

谢谢您的答复,就是用“自由坐标”算的,计算单个原子的msd和gromacs的结果一样,我计算的质心位置也和gromacs一样,但计算水分子的msd却和gromacs不一样,郁闷呀!
3楼2016-08-04 18:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by brucefan at 2016-08-04 14:59:53
yes.

今天又重新算了一下msd,对于水分子而言,我分别算了单个氧原子的msd,假设为msd_o,然后又单独算了另外两个氢原子的msd,假设分别为msd_H1和msd_H2,最后计算msd=(M_O*msd_o+M_H*msd_H1++M_H*msd_H2)/(M_O+M_H+M_H),其中M_O和M_H分别代表氧原子和氢原子的质量,发现算出来的msd和gromacs计算出来的该水分子的msd结果一样。是不是这意味着计算组合体msd时,是先分别计算出各个原子的msd,然后在依照他们的质量进行加权平均???
4楼2016-08-06 21:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
4楼: Originally posted by jianjunjiang at 2016-08-06 21:02:04
今天又重新算了一下msd,对于水分子而言,我分别算了单个氧原子的msd,假设为msd_o,然后又单独算了另外两个氢原子的msd,假设分别为msd_H1和msd_H2,最后计算msd=(M_O*msd_o+M_H*msd_H1++M_H*msd_H2)/(M_O+M_H+M_H ...

而如果先算出水分子的质心,在计算msd,计算得到的结果确和gromacs的结果不一样
5楼2016-08-06 21:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
5楼: Originally posted by jianjunjiang at 2016-08-06 21:04:51
而如果先算出水分子的质心,在计算msd,计算得到的结果确和gromacs的结果不一样...

敬请各位高手给予指点!
6楼2016-08-06 21:06:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
7楼: Originally posted by brucefan at 2016-08-07 03:36:38
我感觉差别不会很大吧。你把两条线画在一起比比看,差别很大吗?你说的这两个定义在数学上应该是不等价的,但我觉得区别应该很小。既然gromacs的结果和你的其中一种定义的结果相同,就说明它很可能用了你说的定义。 ...

谢谢您的答复。这两种算法数学上确实不等价,我算的系统是膜系统。为了检验算法是否正确,只算了系统中的一个水分子,差别还比较明显。看来只有看看手册了
9楼2016-08-07 10:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
8楼: Originally posted by MS_Z at 2016-08-07 10:12:36
同意brucefan 我个人觉得有一点差别不一定需要太纠结 如果真的想知道差别 建议直接看看gromacs内部msd的算法 有可能他给出的msd是用了诸如多时间起点等平滑算法

谢谢您的指点!!
10楼2016-08-07 10:36:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jianjunjiang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 277材料科学与工程080500求调剂 +6 自由煎饼果子 2026-03-16 6/300 2026-03-21 17:21 by 学员8dgXkO
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见