应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » 扩散系数模拟值与实验值、公式计算值有误差
91/1返回列表
查看: 2808  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

说好de

新虫 (小有名气)

[求助] 扩散系数模拟值与实验值、公式计算值有误差 已有2人参与

我在用MS模拟有机蒸气（例如正己烷）在空气中的扩散系数。
我是这样做的：建立模型，能力最小化，常温下跑NPT，之后再动力学跑NVT，最后根据MSD-TIME斜率求出扩散系数，但是扩散系数模拟值与实验值、公式计算值之间有较大误差，实验值、公式计算值是模拟值的3~5倍，请问还要怎么样进一步再改进，使得模拟值跟实验值、公式计算值更接近一些？请求大家指导，谢谢！
还有就是请问我想要用现在的模拟结果发表文章的话，会好发吗？
回复此楼

» 收录本帖的淘帖专辑推荐

MS相关

» 猜你喜欢
加油好姑娘！
1楼 2015-09-09 16:21:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪

【答案】应助回帖

感谢参与，应助指数 +1
考虑这几个方面的影响吧
力场的挑选
NVT平衡的时间
MSD计算的取值区间
前两个因素可以通过查文献参考，msd计算不是取整个时间段，是要扣除前面不稳定的区间的
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

说好de
#有问题发论坛，不接受私聊#爱晶体，爱MaterialsStudio#MS专辑点击右边#http://emuch.net/taotie/41673.html#物性数据#http://webbook.nist.gov/chemistry/#
2楼2015-09-10 10:45:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujc1983

木虫 (著名写手)

苦逼青椒一枚

【答案】应助回帖

感谢参与，应助指数 +1
我个人认为NVE系综下的结果会更好一些，这样没有温度的外界影响
MSD计算属于相关函数计算，需要有大量的数据点：长的模拟时间，小的输出间隔
MSD取值范围是个比较麻烦的事情，但是最好符合一个原则：log（MSD）- log（t）作图，其斜率尽可能接近1。这是扩散系数公式（MSD/t /6）成立的前提，如果偏离很远的话，除以6算出来的扩散系数是不对的
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

说好de
3楼2015-09-11 16:31:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by pigrass at 2015-09-10 10:45:26
考虑这几个方面的影响吧
力场的挑选
NVT平衡的时间
MSD计算的取值区间
前两个因素可以通过查文献参考，msd计算不是取整个时间段，是要扣除前面不稳定的区间的

恩，谢谢您的回复。
我计算的时候是选择的后一段区间进行计算的，并且多取了几个进行了平均值计算，这样的话，应该也会减少些误差吧。
那请问这样的话，3~5倍的模拟结果能好发文章吗？
赞一下 回复此楼
加油好姑娘！
4楼2015-09-11 19:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by xujc1983 at 2015-09-11 16:31:21
我个人认为NVE系综下的结果会更好一些，这样没有温度的外界影响
MSD计算属于相关函数计算，需要有大量的数据点：长的模拟时间，小的输出间隔
MSD取值范围是个比较麻烦的事情，但是最好符合一个原则：log（MSD）-  ...

谢谢您的回复。
NVE系综下能量是基本不变的，并且不会受系统热力学过程的干扰，是这样吗？
那我就要作log（MSD）- log（t）图，看看其斜率是不是尽可能接近1，才能进行下一步计算。
那我用现在的模拟结果发表文章能好发吗？
赞一下 回复此楼
加油好姑娘！
5楼2015-09-11 19:17:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujc1983

木虫 (著名写手)

苦逼青椒一枚
引用回帖:
5楼: Originally posted by 说好de at 2015-09-11 19:17:40
谢谢您的回复。
NVE系综下能量是基本不变的，并且不会受系统热力学过程的干扰，是这样吗？
那我就要作log（MSD）- log（t）图，看看其斜率是不是尽可能接近1，才能进行下一步计算。
那我用现在的模拟结果发表文 ...

是的，NVE系综下体系总能应该平衡的
是计算完后选取MSD区间的时候要做双log验证

发自小木虫Android客户端
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

说好de
6楼2015-09-11 19:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪
引用回帖:
4楼: Originally posted by 说好de at 2015-09-11 19:06:53
恩，谢谢您的回复。
我计算的时候是选择的后一段区间进行计算的，并且多取了几个进行了平均值计算，这样的话，应该也会减少些误差吧。
那请问这样的话，3~5倍的模拟结果能好发文章吗？...

这倒不清楚 你可以查查其他计算扩散的文献 看看他们大概的差异是多少 如果直接发的话肯定会有审稿的问的
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

说好de
#有问题发论坛，不接受私聊#爱晶体，爱MaterialsStudio#MS专辑点击右边#http://emuch.net/taotie/41673.html#物性数据#http://webbook.nist.gov/chemistry/#
7楼2015-09-14 08:15:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by xujc1983 at 2015-09-11 19:44:36
是的，NVE系综下体系总能应该平衡的
是计算完后选取MSD区间的时候要做双log验证
...

恩，好的，谢谢您的回复。
赞一下 回复此楼
加油好姑娘！
8楼2015-09-14 08:55:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)
送红花一朵
引用回帖:
7楼: Originally posted by pigrass at 2015-09-14 08:15:58
这倒不清楚 你可以查查其他计算扩散的文献 看看他们大概的差异是多少 如果直接发的话肯定会有审稿的问的...

恩，好的，谢谢您的回复。
赞一下 回复此楼
加油好姑娘！
9楼2015-09-14 08:55:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 说好de 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭