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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】势能为正值
北京石油化工学院2026年研究生招生接收调剂公告
172/2返回列表上一页12
查看: 2888  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★
shengxiang(金币+5): 2010-09-05 08:31:50
lei0736(金币+8):谢谢 2010-09-05 09:33:36
通常来说，许多被模拟体系平衡时其势能为负值。比如对原子体系来讲，一般来说凝聚系的势能零点为无限远原子，即能量是与孤立原子体系相比。孤立原子形成凝聚态，这个过程如果是自发的，能量应是减小，所以凝聚系平衡势能常是负值。尤其是金属体系，基本上应该都是负值。但如果是分子体系，由于势能毕竟只是自由能的一部分，而且对于通过拟合体系性质得到的经验势来说，它的势能零点有时候好像也不是特别清楚，不一定是孤立原子，比如我们拟合时常常只关心势函数的一阶导和二阶导数，或者我们只关心各原子混合性质等等，所以出现势能为正的情况，倒也不见得违背什么原则。事实上，假如处理一个分子体系，它的势函数模型中包含了分子内的部分，因为常见的分子内的势函数，三体势和四体势（二面角）都经常选择为正的形式（如K*(theta-theta0)^2这种），而分子间斥力、色散力等自然也都是正的，那么总势能也将是有可能大于零的。

总之，对于分子体系，如果计算中不出现什么异常情况，我觉得势能大于零也可以被接受。一句话，如果是负值，更让人放心一些，如果是正值，假如仔细核对无误，而且计算中无异常，可以生成合理的结构，得到合理的性质，那就不要追究势能正负的问题了。

[ Last edited by 老虎大王 on 2010-9-4 at 18:58 ]
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

xwwjd
11楼2010-09-04 18:54:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)
引用回帖:
Originally posted by 043114076 at 2010-09-04 11:25:49:


对于分子力学来说，能量各部分的0点定义是不一样的。对于谐振势，键长键角0点都是在平衡位置，而且能量一定不会小于0. 库仑势能的0点是无穷远处，能量一定是负值。LJ VDW势除了无穷远处的0点外在r0的位置还有 ...

明白了，你用GROMACS吗？GROMACS中g_energy求出的势能是哪一类呢？
赞一下 回复此楼
12楼2010-09-05 08:29:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)
引用回帖:
Originally posted by yahoohoo at 2010-09-04 14:34:16:


最简单的LJ势能就可以是正值

明白了一些了，就像楼上说的，得看是哪种势能了
赞一下 回复此楼
13楼2010-09-05 08:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)
引用回帖:
Originally posted by 老虎大王 at 2010-09-04 18:54:27:
通常来说，许多被模拟体系平衡时其势能为负值。比如对原子体系来讲，一般来说凝聚系的势能零点为无限远原子，即能量是与孤立原子体系相比。孤立原子形成凝聚态，这个过程如果是自发的，能量应是减小，所以凝聚系平 ...

老虎大王，久仰久仰，小弟新人，还多照顾，多谢你了~你对知识了解得都很透彻啊，向你学习~
赞一下 回复此楼
14楼2010-09-05 08:33:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwwjd

铜虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
引用回帖:
11楼: Originally posted by 老虎大王 at 2010-09-04 18:54:27
通常来说，许多被模拟体系平衡时其势能为负值。比如对原子体系来讲，一般来说凝聚系的势能零点为无限远原子，即能量是与孤立原子体系相比。孤立原子形成凝聚态，这个过程如果是自发的，能量应是减小，所以凝聚系平衡 ...

你好，我也遇到这种势能为正的问题，我算的就是简单的乙醇体系，液态的时候势能是负值，但是我算气态乙醇的时候就是一个正值，其中它的远程库伦静电作用能是一个正值，这是为什么呢？不能理解，求大神指教。。
赞一下 回复此楼
你一会看云，一会看我，我觉得你看云时很近，看我时很远。。。
15楼2014-05-06 10:10:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 感谢指导！ 2014-05-16 18:51:16
引用回帖:
15楼: Originally posted by xwwjd at 2014-05-06 10:10:00
你好，我也遇到这种势能为正的问题，我算的就是简单的乙醇体系，液态的时候势能是负值，但是我算气态乙醇的时候就是一个正值，其中它的远程库伦静电作用能是一个正值，这是为什么呢？不能理解，求大神指教。。...

你可以比较一下气态和和液态的各项势能之间的差别。
赞一下 回复此楼
16楼2014-05-07 08:21:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwwjd

铜虫 (初入文坛)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流 2014-05-16 18:51:21
引用回帖:
6楼: Originally posted by shengxiang at 2010-09-03 21:37:07
我跑的是醇体系，那要是计算中可能有错误的话，会出现在哪个环节呢？哪些参数？...

请问楼主，这个问题你后来是怎么解决的，我也遇到类似问题，模拟的是环己烷的一个体系，势能总是正值~
赞一下 回复此楼
你一会看云，一会看我，我觉得你看云时很近，看我时很远。。。
17楼2014-05-07 09:19:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shengxiang 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 0703化学 +13 妮妮ninicgb 2026-04-04 14/700 2026-04-06 07:12 by houyaoxu
[考研] 085600，320分求调剂 +15 大馋小子 2026-04-04 16/800 2026-04-06 06:34 by houyaoxu
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 270分求调剂 +4 maxjxbsk 2026-04-01 4/200 2026-04-05 17:04 by yulian1987
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 求调剂，一志愿郑州大学材料与化工专硕，英二数二342分，求老师收留 +18 v12abo 2026-04-02 20/1000 2026-04-05 11:37 by a8144223
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 363求调剂 +7 zh096 2026-04-04 7/350 2026-04-04 17:13 by dongzh2009
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 材料调剂 +11 吴棂颖！ 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +5 redcom 2026-04-01 6/300 2026-04-03 15:15 by Howard28
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
信息提示
关闭
请填处理意见
关闭