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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » 扩散系数的计算及体系平衡的判断
51/1返回列表
查看: 2254  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yida

木虫 (小有名气)

[求助] 扩散系数的计算及体系平衡的判断

我想模拟糖分子在水体系中与金表面的相互作用，采用Discover跑的动力学，扩散系数是许多分子的统计平均值的话，我所构建的模型中只含一个糖分子，想问下这个的模型计算扩散系数合适吗?如何判断体系有没有达到平衡态？我看到许多判断标准，有看动力学时间的，有看MSD对t曲线的，还有看lnMSD对lnt作图看是否成线性的，还有看RMSD是否平行于x轴呈一条直线，想问下如何得到RMSD随时间变化的曲线？我用的是MS 6.0，用的是Forcite模块分析，还有就是如果看是否呈线性的话r值是多少？比如我运行了5ns的动力学，3-4.5ns的呈线性这样是否达到平衡？问题比较多，希望高手指导
回复此楼

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

晶体模型计算分析 扩散系数

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-07-01 15:48:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yida

木虫 (小有名气)
引用回帖:
4楼: Originally posted by sciencejoy at 2013-07-02 08:56:09
你最开始跑的500ps，为什么是500ps？这个时间应该选取至少2倍弛豫时间，不是随便设个500ps，有时候，跑500ps未必体系平衡。很少有模拟选NVE系综，你为什么选NVE系综？NVT、NPT无法满足你的体系？

所需要弛豫时间怎么计算？我看MS教程上举的例子是先用NVT系综平衡，后用NVE系综跑动力学，解释是分析热力学性质最好采用NVE系综，我的体系是下边一层金，上边是水的体系中有糖单元，最上边有一层真空层，采用NPT的话盒子会越来越小，真空层直接消失，水分子也会到金中去，感觉很不合理，我很不理解，就咨询了个师兄，师兄说在密度合理的情况下尽量少用NPT系综，我看了关于模拟的书，包括小木虫上的讨论，对这几个系综还不是很了解，请高手解惑
赞一下 回复此楼
5楼2013-07-02 09:18:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

sciencejoy

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2013-07-02 08:18:13
yida: 金币+10, ★★★很有帮助, 非常感谢 2013-07-02 08:34:21
计算弛豫时间，至少2倍弛豫时间，体系一般是平衡态了。同时算算能量、温度等，是不是也不随时间变化。
赞一下 回复此楼
2楼2013-07-01 22:53:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yida

木虫 (小有名气)
引用回帖:
2楼: Originally posted by sciencejoy at 2013-07-01 22:53:32
计算弛豫时间，至少2倍弛豫时间，体系一般是平衡态了。同时算算能量、温度等，是不是也不随时间变化。

请教下弛豫时间如何计算？
我一般建好模型后会把模型的结构优化，然后跑500ps NVT系综的动力学，最后跑NVE系综分析热力学性质，但是NVE系综的问题就是能量虽然是守恒的，动力学过程中温度一直在上升，比如我跑NVE系综设置的初始温度是298K,跑过NVE后的温度可能就上升到了500K左右，这时的性质是不是就不能代表298K时的性质了？看到有的帖子上提到可以采用NVT来对温度进行校正，想问下具体怎么操作？
问题比较多，请高手指导！
赞一下 回复此楼
3楼2013-07-02 08:39:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
yida: 金币+2, ★有帮助 2013-07-02 09:19:00
chaizhm: 金币+2, 谢谢~ 2013-07-02 18:38:05
你最开始跑的500ps，为什么是500ps？这个时间应该选取至少2倍弛豫时间，不是随便设个500ps，有时候，跑500ps未必体系平衡。很少有模拟选NVE系综，你为什么选NVE系综？NVT、NPT无法满足你的体系？
赞一下 回复此楼
4楼2013-07-02 08:56:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +6 来了个 2024-06-26 7/350 2024-06-27 13:09 by 张小生子
[基金申请] 博后面上今天下午会公布吗？大家有无消息？ +10 地理学1995 2024-06-24 14/700 2024-06-27 11:24 by jackiddd
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +7 来了个 2024-06-26 7/350 2024-06-27 11:11 by wubiansiya
[基金申请] 省基金太难了 +8 aasahr 2024-06-26 8/400 2024-06-27 10:35 by 小狮子2010
[硕博家园] 联培博士文章第一单位署名问题交流~ +10 橙成成c 2024-06-23 23/1150 2024-06-27 10:11 by 橙成成c
[基金申请] 专科院校申请国自然 +17 孤单12站 2024-06-24 26/1300 2024-06-27 08:45 by 3001160025
[论文投稿] 第一篇论文投稿进程记录 +7 慎独的小花卷 2024-06-20 12/600 2024-06-27 00:49 by 玩耍绝缘体
[考研] 刚当完兵回来想考研 +7 五官在线 2024-06-23 18/900 2024-06-26 17:05 by 1158057902
[基金申请] 国产期刊影响因子大于12的有20多个个了 +19 babu2015 2024-06-20 24/1200 2024-06-26 16:57 by deliciou
[基金申请] 省基金对国自然申请是有利还是有弊 +6 zycumt04 2024-06-26 7/350 2024-06-26 15:02 by 潇湘之迷
[硕博家园] 难道我真的要放弃吗？ +42 133456 2024-06-20 44/2200 2024-06-26 12:32 by lyfbangong
[基金申请] 评审专家会不会很在意申请人的单位啊 +11 lancet0903 2024-06-24 13/650 2024-06-26 11:42 by 漠上藜梭
[有机交流] 过柱子，产品，杂质在是 pe:ea=100:1 也一起出来？ +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[基金申请] 这样的说辞是上会了吗 +12 学员d3zYCz 2024-06-24 15/750 2024-06-25 22:05 by 3115321
[基金申请] 博后面上到底什么时候出结果？？？ +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 青年和面上，哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[公派出国] 博士csc联培会看重第一学历学校层次吗 +4 也就这样 2024-06-23 4/200 2024-06-24 08:18 by 晓目崇
[博后之家] 在国内某高校做全职博士后2年，现在找到新的单位，出站或退站对新工作有什么影响？ +10 nxplfcc 2024-06-20 10/500 2024-06-22 07:52 by 徐长安
信息提示
关闭
请填处理意见
关闭