24小时热门版块排行榜    

查看: 3085  |  回复: 16

yaneehan2010

银虫 (正式写手)


[交流] 扩散系数模拟的结果比资料上的小一个数量级

如题:我用MS进行分子扩散系数模拟时,结果比在手册上查的小一个数量级,这是为啥呀?请各位给予解答呀、、
另外,想请教一下,最终出来MSD曲线,那个直线段的斜率数值单位是埃^2/ps,除以6之后,是不是要进行单位换算呀,手册上是cm^2/s。
先行谢过呀、、
回复此楼

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

先进材料与分子模拟 感悟人生

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水晶蚂蚁

金虫 (小有名气)


★ ★
yaneehan2010(金币+1): 谢谢参与
chaizhm: 金币+1, 谢谢~ 2012-07-11 08:47:11
你的扩散系数模拟MSD 曲线是取得哪一段做为斜率?不同位置,数据不一样,单位换算肯定是要换算的,过程不麻烦。而且你的扩散系数与手册有差别也是正常的,你的计算设置也影响结果的,你做的是什么系综?
2楼2012-07-11 08:33:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaneehan2010(金币+1): 谢谢参与
保证单位统一以后再计算
3楼2012-07-11 08:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaneehan2010

银虫 (正式写手)


引用回帖:
2楼: Originally posted by 水晶蚂蚁 at 2012-07-11 08:33:58
你的扩散系数模拟MSD 曲线是取得哪一段做为斜率?不同位置,数据不一样,单位换算肯定是要换算的,过程不麻烦。而且你的扩散系数与手册有差别也是正常的,你的计算设置也影响结果的,你做的是什么系综?

感谢回复呢,我先用了NPT,后来是NVE各跑了100000,MSD曲线我是取了拖尾前的那段近似直线的部分计算斜率。另外还用matlab将数据拟合了一条,这个单位换算是10^(-4)对吧。
我模拟的是298K,1atm时的情况,我建晶胞以及以后的系综时,温度都写得是298K。密度我写的是其中一种分子在298K时的密度,这样设置有什么问题么?该怎么改呢?谢谢啦
4楼2012-07-11 10:27:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaneehan2010

银虫 (正式写手)


引用回帖:
3楼: Originally posted by chaizhm at 2012-07-11 08:47:36
保证单位统一以后再计算

谢谢您的回复,我单位换算了呀,10^(-4)对吧?
5楼2012-07-11 10:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水晶蚂蚁

金虫 (小有名气)


引用回帖:
4楼: Originally posted by yaneehan2010 at 2012-07-11 10:27:55
感谢回复呢,我先用了NPT,后来是NVE各跑了100000,MSD曲线我是取了拖尾前的那段近似直线的部分计算斜率。另外还用matlab将数据拟合了一条,这个单位换算是10^(-4)对吧。
我模拟的是298K,1atm时的情况,我建晶胞 ...

单位换算关系是对的,你的set up 设置中T和B 是如何设置的,运算了100000步么?那么你步长是多少,运算时间除以步长再乘以1000是步数,而且你跑完动力学,你选取的近似直线部分是多长时间内,建议中间段不同时间多取几次,然后对得到的斜率取平均值。希望对楼主有所帮助。
6楼2012-07-11 15:59:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaneehan2010

银虫 (正式写手)


引用回帖:
6楼: Originally posted by 水晶蚂蚁 at 2012-07-11 15:59:09
单位换算关系是对的,你的set up 设置中T和B 是如何设置的,运算了100000步么?那么你步长是多少,运算时间除以步长再乘以1000是步数,而且你跑完动力学,你选取的近似直线部分是多长时间内,建议中间段不同时间多 ...

哪里的set up?
我的time step:1.0fs,Dynamics time:100ps,number of steps:100000.
我是直接把那段时间的都算了、、
7楼2012-07-11 17:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gxx2012

新虫 (初入文坛)



yaneehan2010(金币+1): 谢谢参与
引用回帖:
7楼: Originally posted by yaneehan2010 at 2012-07-11 17:22:48
哪里的set up?
我的time step:1.0fs,Dynamics time:100ps,number of steps:100000.
我是直接把那段时间的都算了、、...

麻烦问一下楼主,你说的扩散系数的手册名字是什么啊?我一直在寻找。。。还请指点。。。万分感谢!
8楼2012-10-19 15:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zx2456

木虫 (正式写手)



yaneehan2010(金币+1): 谢谢参与
不知道楼主做的是什么体系的扩散系数?我做过高温熔盐的最后一段的MSD都平了,意味着不流动了,但是实际上那时应该是流动的状态,不知道是怎么回事,请教高手回答。
9楼2012-12-12 20:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytudou

银虫 (小有名气)



yaneehan2010(金币+1): 谢谢参与
送红花一朵
请问查得是什么手册啊?
10楼2015-02-07 21:00:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶味果冻

新虫 (小有名气)



yaneehan2010(金币+1): 谢谢参与
你查的手册单位cm^2/s,是说的MS计算出来的MSD图 用Tools中plot graph 出来的斜率值的单位么,
11楼2015-03-17 15:51:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaneehan2010

银虫 (正式写手)


不好意思,我这块工作只做了一点儿,就换了课题,很长时间不做了
12楼2015-03-18 15:29:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)



yaneehan2010(金币+1): 谢谢参与
引用回帖:
10楼: Originally posted by ytudou at 2015-02-07 21:00:46
请问查得是什么手册啊?

请问手册查到了吗?我也在查手册,一直没找到。
13楼2015-09-08 13:07:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)


引用回帖:
12楼: Originally posted by yaneehan2010 at 2015-03-18 15:29:49
不好意思,我这块工作只做了一点儿,就换了课题,很长时间不做了

那请问你当时查的手册还有吗?我一直在找,没有找到。有的话能不能分享一下,谢谢你。
14楼2015-09-08 13:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxfysdx2009

金虫 (小有名气)



yaneehan2010(金币+1): 谢谢参与
引用回帖:
14楼: Originally posted by 说好de at 2015-09-08 13:08:12
那请问你当时查的手册还有吗?我一直在找,没有找到。有的话能不能分享一下,谢谢你。...

你找到手册了么?我为查一个数据找了好几天了。
15楼2015-11-15 11:05:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

说好de

新虫 (小有名气)


引用回帖:
15楼: Originally posted by wxfysdx2009 at 2015-11-15 11:05:43
你找到手册了么?我为查一个数据找了好几天了。...

没有找到……

发自小木虫Android客户端
16楼2015-11-16 08:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mailuoshan

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by yaneehan2010 at 2012-07-11 10:27:55
感谢回复呢,我先用了NPT,后来是NVE各跑了100000,MSD曲线我是取了拖尾前的那段近似直线的部分计算斜率。另外还用matlab将数据拟合了一条,这个单位换算是10^(-4)对吧。
我模拟的是298K,1atm时的情况,我建晶胞 ...

我想问一下用DL_POLY软件仿真得到的MSD数据是单个原子不同时刻的MSD,怎么把构成阳离子的不同原子的MSD等效成阳离子的MSD
17楼2016-09-23 20:25:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yaneehan2010 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 请问大家,我有希望进浙大普博么 +5 高艺文 2024-09-28 8/400 2024-09-29 00:08 by 小红豆
[论文投稿] applied catalysis B 投稿模版 5+4 午后夏日 2024-09-27 5/250 2024-09-28 22:59 by 莱茵润色
[硕博家园] 希望有师兄师姐帮忙引荐申博 +6 zhanyaqian 2024-09-25 11/550 2024-09-28 22:24 by 自强2019
[硕博家园] 毕业论文的数据能否再发小论文 +6 20081002 2024-09-28 7/350 2024-09-28 22:20 by 鱼翔浅底1
[教师之家] 东南大学电气学院电力电子方向招收硕士博士 +4 蜡笔小鑫1989 2024-09-26 7/350 2024-09-28 21:16 by 伍文童
[论文投稿] 求推荐医学类杂志 35+4 wshxtim1 2024-09-25 4/200 2024-09-28 13:56 by Wormaciae
[考博] 25有机合成申博 +3 再等你啊 2024-09-24 4/200 2024-09-28 13:09 by majun521
[论文投稿] 想跟大家聊聊对国产SCI期刊的看法 +14 wleizl 2024-09-26 17/850 2024-09-28 11:59 by wleizl
[电化学] 刚装好的电池电压掉很快怎么解决?? +3 花花你最棒 2024-09-26 6/300 2024-09-28 08:20 by cnofc
[考博] 材料/电信/生物-2025普博生自荐-985本双非硕一区一作 +4 enowei0127 2024-09-23 8/400 2024-09-27 22:40 by enowei0127
[论文投稿] Journal of colloid and interface science期刊 10+4 1821588220 2024-09-24 10/500 2024-09-27 20:37 by 1821588220
[硕博家园] 当前读博士还有性价比吗? +8 苏东坡二世 2024-09-23 8/400 2024-09-27 11:39 by 半生梦君
[考博] 数学博导 +4 学术霸王 2024-09-25 6/300 2024-09-27 11:14 by 青古
[有机交流] 二氯甲烷的去除 +3 cgsa吧 2024-09-24 8/400 2024-09-27 10:55 by bear2007
[考博] 有机转材料申博难度 30+3 昕散 2024-09-22 10/500 2024-09-26 22:42 by 824282658
[基金申请] 请教2024后期资助到哪个阶段了 +3 拾光者5566 2024-09-23 8/400 2024-09-26 09:27 by 拾光者5566
[论文投稿] NC投稿多久有消息呀? +4 Jordanblood 2024-09-24 6/300 2024-09-25 16:58 by FZX_2024
[有机交流] chemdraw结构式复制到word变形 +5 笑看人生1993 2024-09-22 11/550 2024-09-25 14:44 by 冰蓝夜游神
[论文投稿] 爱斯维尔旗下Postharvest Biology and Technology,采用什么系统查查 20+4 sdsdsxh 2024-09-22 9/450 2024-09-25 07:33 by sdsdsxh
[有机交流] 做聚酰亚胺一步法合成,机械搅拌怎么密封不漏气 +3 哼哈EN 2024-09-24 5/250 2024-09-24 14:16 by mrzhl1986
信息提示
请填处理意见