24小时热门版块排行榜    

查看: 1743  |  回复: 7

amynihao

金虫 (正式写手)


[交流] 【讨论】很奇怪的能量,温度曲线。

我用MS软件中的discover跑动力学,跑了100ps,最后出来的能量曲线和温度曲线如下图,从图可以看出温度和能量已经几乎不变了,看来体系已经平衡了,但是为什么能量却变大了,不是应该能量最低状态才是平衡态吗?而且温度下降了。好奇怪的现象,希望大家帮忙解释一下原因吧,谢了。
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

wally8962

木虫 (著名写手)


★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+2): 谢谢 2011-02-05 20:53:35
NVE?要看总能量而不是某项。
potential energy 和Nonbond都变小了。很正常。
2楼2011-01-25 11:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amynihao

金虫 (正式写手)


引用回帖:
Originally posted by wally8962 at 2011-01-25 11:20:28:
NVE?要看总能量而不是某项。
potential energy 和Nonbond都变小了。很正常。

不是NVE。是NVT.potential energy 和Nonbond都变大了啊,你怎么说变小了呢。
3楼2011-01-25 11:26:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)


★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+2): 谢谢 2011-02-05 20:53:42
amynihao(金币+5): 2011-07-24 11:45:18
严格点说,两个能量都变低了。
你别看绝对值
4楼2011-01-25 11:58:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)


★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+10): 谢谢 2011-02-05 20:53:55
amynihao(金币+5): 2011-07-24 11:45:27
引用回帖:
Originally posted by amynihao at 2011-01-25 10:37:08:
我用MS软件中的discover跑动力学,跑了100ps,最后出来的能量曲线和温度曲线如下图,从图可以看出温度和能量已经几乎不变了,看来体系已经平衡了,但是为什么能量却变大了,不是应该能量最低状态才是平衡态吗?而 ...

仔细盯着看了你的图, 超过一分钟, 也愣是没有明白你是怎么从图中看出"能量变大", "温度减小"的结论.

Anyway, 几点考虑一下:

(1) 不管你的体系怎样, 不管你是NVT还是NVE, 100ps似乎有些太短了.
(2) 看是否到达平衡, 不需要将 E vs. t 从头到尾都在图上体现出来. 否则, 你看看你的Y轴, 两个点之间的能量差是500Kcal/mol, 所以就算你看起来是一条直线了, 也还是伪直线, 说不定震荡还是在 200Kcal/mol左右. 同理, 温度的震荡在你的scale下很难判断.

所以,

(a) 尝试跑个300ps的dynamics, 用作平衡体系(我假设你的体系不是很大, 体系比较容易到达平衡);
(b) 在(a)基础上再跑5ps, 统计能量和温度(我再假设你的步长在10fs左右, 那么5ps能给你带来500组数据, 应该比较直观见效了).

[ Last edited by qphll on 2011-1-25 at 13:16 ]
5楼2011-01-25 13:15:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+2): 谢谢 2011-02-05 20:54:23
amynihao(金币+10): 2011-07-24 11:45:34
拿时间滚动平均值出来比较更容易看出趋势。
6楼2011-01-26 11:32:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倚天舰

木虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖
amynihao(金币+5): 2011-07-24 11:45:40
楼主的意思是不是说初始时的势能是接近-73000,而平衡之后是不到-71500,因此你觉着能量变大了是么?我也有这样的疑问啊。不知道楼主解决了这个问题没
7楼2011-07-23 11:10:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

supermicro

新虫 (初入文坛)


amynihao(金币+25): 2011-07-24 11:46:00
??
8楼2011-07-24 11:43:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 amynihao 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 321求调剂 +4 大米饭! 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
信息提示
请填处理意见