24小时热门版块排行榜    

查看: 3115  |  回复: 16

shengxiang

木虫 (著名写手)

[交流] 【求助】势能为正值

请教个问题,跑完动力学之后我用e_energy分析势能随时间的变化时发现势能都是正的,这怎么解释啊?

势能图如下:
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yahoohoo

铁杆木虫 (著名写手)

★ ★
shengxiang(金币+1): 2010-09-05 08:29:41
lei0736(金币+2):谢谢 2010-09-05 09:33:16
引用回帖:
Originally posted by shengxiang at 2010-09-03 21:39:31:



是这样的,我的体系中只有醇,我算的这个势能是用GROMACS中的g_energy选择势能“potentia”,我看文献中求出的势能都是负的~这个势能是以什么为标准呢?可能是正的吗?请不吝赐教~万分感谢~l

最简单的LJ势能就可以是正值
10楼2010-09-04 14:34:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

345320587

银虫 (初入文坛)


lei0736(金币+1):谢谢 2010-09-02 19:46:09
shengxiang(金币+5): 2010-09-03 21:36:01
你跑的是个什么体系?我从前跑过肽段,得到的势能也是正值,是没有问题的 不过看你的势能值貌似也太高了点了啊 是不是计算方面有什么错误呢~
2楼2010-09-02 16:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

★ ★
lei0736(金币+2):谢谢 2010-09-02 19:46:16
shengxiang(金币+5): 2010-09-03 21:36:08
引用回帖:
Originally posted by shengxiang at 2010-09-02 10:45:00:
请教个问题,跑完动力学之后我用e_energy分析势能随时间的变化时发现势能都是正的,这怎么解释啊?

你可以先考虑下这个问题:

势能为什么不能是正的?

然后是:

你模拟的体系中势能由哪些部分组成?
3楼2010-09-02 17:23:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

★ ★
lei0736(金币+2):谢谢 2010-09-02 22:10:43
shengxiang(金币+1): 2010-09-03 21:36:13
4楼2010-09-02 20:32:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)

我跑的是醇体系,那要是计算中可能有错误的话,会出现在哪个环节呢?哪些参数?

我跑的是醇体系,那要是计算中可能有错误的话,会出现在哪个环节呢?哪些参数?

[ Last edited by shengxiang on 2010-9-3 at 21:34 ]
5楼2010-09-03 21:32:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)

谢谢

引用回帖:
Originally posted by 345320587 at 2010-09-02 16:37:45:
你跑的是个什么体系?我从前跑过肽段,得到的势能也是正值,是没有问题的 不过看你的势能值貌似也太高了点了啊 是不是计算方面有什么错误呢~

我跑的是醇体系,那要是计算中可能有错误的话,会出现在哪个环节呢?哪些参数?
6楼2010-09-03 21:37:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)

谢谢

引用回帖:
Originally posted by yahoohoo at 2010-09-02 17:23:27:


你可以先考虑下这个问题:

势能为什么不能是正的?

然后是:

你模拟的体系中势能由哪些部分组成?

是这样的,我的体系中只有醇,我算的这个势能是用GROMACS中的g_energy选择势能“potentia”,我看文献中求出的势能都是负的~这个势能是以什么为标准呢?可能是正的吗?请不吝赐教~万分感谢~l
7楼2010-09-03 21:39:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)

谢谢

引用回帖:
Originally posted by 043114076 at 2010-09-02 20:32:45:
看一下这个帖子
http://muchong.com/bbs/viewthread.php?tid=2297056

感谢你的回复,那个帖子我看了,如果势能出现负值的话我认为是正确的,因为取无穷远处和真空为势能零点,但是出现正值的话我真不知道该怎么解释了,望帮助~
8楼2010-09-03 21:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

★ ★ ★ ★ ★
ghcacj(金币+5):谢谢 2010-09-04 12:10:57
shengxiang(金币+2): 2010-09-05 08:28:43
引用回帖:
Originally posted by shengxiang at 2010-09-03 21:41:15:



感谢你的回复,那个帖子我看了,如果势能出现负值的话我认为是正确的,因为取无穷远处和真空为势能零点,但是出现正值的话我真不知道该怎么解释了,望帮助~

对于分子力学来说,能量各部分的0点定义是不一样的。对于谐振势,键长键角0点都是在平衡位置,而且能量一定不会小于0. 库仑势能的0点是无穷远处,能量一定是负值。LJ VDW势除了无穷远处的0点外在r0的位置还有一个零点,当原子距离小于r0时能量大于0,大于r0时能量为负。
总的势能是这几项的加和,能量可能大于0也可能小于0,特别是当分子构型扭曲比较厉害时,成键能贡献很大。
9楼2010-09-04 11:25:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shengxiang 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 cu9Nq1xK233Z 2026-08-03 14/700 2026-08-04 14:44 by jzdxtiUmJ6PL
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +20 低垂的野花 2026-07-31 26/1300 2026-08-04 14:43 by Ldrop2023
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 5/250 2026-08-04 13:31 by jzdxtiUmJ6PL
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 6/300 2026-08-04 13:28 by jzdxtiUmJ6PL
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 jRl3mE6ddZGq 2026-08-03 11/550 2026-08-04 12:58 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见