24小时热门版块排行榜    

查看: 1555  |  回复: 16

lusostty

金虫 (小有名气)

[求助] 【动力学-扩散】关于模拟X时间 舍弃前Y时间的理解 已有2人参与

关于扩散系数的模拟

很多文献都是这样描述的:例如,模拟进行了1000ps,前200ps-1000ps的动力学构型计算扩散系数

我之前的做法是:直接进行1000ps的动力学模拟,得到一串数据,取其中200ps-1000ps进行线性拟合,计算扩散系数。虽然相关性还不错,但这条直线应该是不过原点的。

今天忽然想起来这种方法应该值得商榷,会不会是这样的:先进行200ps的模拟,然后在此基础上(current)进行800ps的模拟,这样得到的MSD~T线性拟合应该是过原点的

以上 两种方法哪一种正确,或者都存在问题?

谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
9楼: Originally posted by lusostty at 2014-07-12 12:26:29
那如果计算扩散系数呢?如果取中间段MSD数据  它也是不过原点的 那怎么办...

根据爱因斯坦方程,求出直线的斜率即可。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
11楼2014-07-12 12:28:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lusostty

金虫 (小有名气)

自己顶帖 :
2楼2014-07-10 16:29:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

polypro

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
对于Dc计算,前200ps是用于消除初始的人为因素干扰的。理论上讲,Dc是动态平衡值,与初始点的选取无关。
泉涸,鱼相与处于陆,相呴以湿,相濡以沫,不如相忘于江湖。
3楼2014-07-11 08:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
个人认为,两种方法计算均方位移的时候,都应该取中间时间段的数据。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2014-07-11 08:41:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 月只蓝 at 2014-07-11 08:41:11
个人认为,两种方法计算均方位移的时候,都应该取中间时间段的数据。

那拟合的直线方程 就不会过原点了  总觉得这样做似乎有点问题
5楼2014-07-11 09:49:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)

引用回帖:
3楼: Originally posted by polypro at 2014-07-11 08:03:17
对于Dc计算,前200ps是用于消除初始的人为因素干扰的。理论上讲,Dc是动态平衡值,与初始点的选取无关。

谢谢你的回复 我考察吸附量对扩散系数的影响 拟合的都是一条条不过原点的直线 这样具有可比性么?
6楼2014-07-11 09:51:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 月只蓝 at 2014-07-11 08:41:11
个人认为,两种方法计算均方位移的时候,都应该取中间时间段的数据。

谢谢你的回复 我考察吸附量对扩散系数的影响 拟合的都是一条条不过原点的直线 这样具有可比性么?
7楼2014-07-11 09:52:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
7楼: Originally posted by lusostty at 2014-07-11 09:52:04
谢谢你的回复 我考察吸附量对扩散系数的影响 拟合的都是一条条不过原点的直线 这样具有可比性么?...

文献上常直接比较MSD曲线

[ 发自小木虫客户端 ]
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2014-07-11 12:19:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)

引用回帖:
8楼: Originally posted by 月只蓝 at 2014-07-11 12:19:18
文献上常直接比较MSD曲线
...

那如果计算扩散系数呢?如果取中间段MSD数据  它也是不过原点的 那怎么办
9楼2014-07-12 12:26:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
9楼: Originally posted by lusostty at 2014-07-12 12:26:29
那如果计算扩散系数呢?如果取中间段MSD数据  它也是不过原点的 那怎么办...

为什么要求过原点呢?
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
10楼2014-07-12 12:27:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lusostty 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见