24小时热门版块排行榜    

查看: 2992  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

说好de

新虫 (小有名气)

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

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 BCdvBAEpt3zm 2026-11-01 4/200 2026-11-02 10:23 by dcc2yOIRGGTT
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 10:02 by dcc2yOIRGGTT
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 09:31 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +3 iaeyuan 2026-11-01 3/150 2026-11-02 07:49 by kudofaye
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 06:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:05 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
请填处理意见