24小时热门版块排行榜    

查看: 2795  |  回复: 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-06-19 6/300 2024-06-26 16:38 by 大大熊呀
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[论文投稿] 我在写关于多维空间的论文,希望能与大家交流,准备在science上发表 +8 lgf519 2024-06-24 10/500 2024-06-26 16:25 by 梦渺岚烟
[育儿交流] 暑假孩子一个人在家可以干什么 +4 唠叨的映 2024-06-25 5/250 2024-06-26 15:35 by 361997417
[教师之家] 神奇的中医 +8 水冰月月野兔 2024-06-24 10/500 2024-06-26 12:10 by lyfbangong
[有机交流] 大佬们,打的核磁氢谱与chemdraw预测的有偏差 +4 来了个 2024-06-26 4/200 2024-06-26 08:45 by czyzsu
[基金申请] 这样的说辞是上会了吗 +12 学员d3zYCz 2024-06-24 15/750 2024-06-25 22:05 by 3115321
[有机交流] 对苯乙烯磺酰氯的合成机理 25+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[考博] 一名额难求啊 +14 pinbo拼搏 2024-06-19 14/700 2024-06-25 16:55 by wjykycg
[基金申请] 能看出是否上会了吗 +10 articlefan 2024-06-23 15/750 2024-06-25 16:05 by 请慎重修改昵称
[基金申请] 2024安徽省哲社立项名单公示 +3 robin_work 2024-06-24 6/300 2024-06-25 14:37 by Pickfoot
[基金申请] 焦虑没有毛线用,默默前行是王道! +4 漠上藜梭 2024-06-24 8/400 2024-06-25 14:32 by 漠上藜梭
[教师之家] 复旦夏同学退学理由说明,哪儿可以下载? +7 苏东坡二世 2024-06-21 8/400 2024-06-25 12:42 by 药嘿1233
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[考博] 申请25博士,可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 教育部基金 +5 m1393 2024-06-21 5/250 2024-06-21 21:13 by odes
[基金申请] 面上博士后 +11 jsqy 2024-06-19 13/650 2024-06-21 17:12 by sizhouyi
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
信息提示
请填处理意见