24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 2448  |  回复: 8

说好de

新虫 (小有名气)

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

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

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

MS相关

» 猜你喜欢

加油好姑娘!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪

【答案】应助回帖

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

» 本帖已获得的红花(最新10朵)

#有问题发论坛,不接受私聊#爱晶体,爱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朵)

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朵)

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朵)

#有问题发论坛,不接受私聊#爱晶体,爱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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学 +29 妮妮ninicgb 2026-04-04 32/1600 2026-04-08 22:50 by 猪会飞
[基金申请] 有爆料,一个青年教师卖房得400万,然后换了一个四青帽子 +5 babu2015 2026-04-08 5/250 2026-04-08 22:23 by god_tian
[考研] 0703化学求调剂 +16 不知名的小卅 2026-04-08 16/800 2026-04-08 21:53 by 旺旺豆包
[考研] 材料调剂 +13 一样YWY 2026-04-06 13/650 2026-04-08 20:44 by LiuXingYu星星
[考研] 085600材料与化工专硕329 求调剂 +21 额cc 2026-04-06 22/1100 2026-04-08 18:33 by 环化材-小生
[考研] 347材料专硕求调剂 +18 zj8215216 2026-04-06 18/900 2026-04-08 16:27 by luoyongfeng
[考研] 材料科学与工程320求调剂,080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 372求调剂 +4 jj涌77 2026-04-02 4/200 2026-04-07 09:31 by 白云朵朵飞
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +4 江淮北月 2026-04-05 4/200 2026-04-06 22:44 by chenzhimin
[考研] 287求调剂 +3 通信学硕081000 2026-04-03 4/200 2026-04-06 21:03 by going home
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 生物与医药086000调剂一志愿西北农林320分 +3 美美女士 2026-04-03 3/150 2026-04-05 21:55 by 学员8dgXkO
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 一志愿北京科技大学材料工程085601,求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
信息提示
请填处理意见