体系势能为正有没有什么解决办法？

appleliu

请问有没有人解决过系统弛豫稳定后总体势能为正的问题呀，我最近遇到了这个问题，但不知道解决办法，各位大神求指点！！！按理说总体势能为负值呀，我用的OPLS-AA力场。新手求求求教！！！

sobereva

发帖时候别打那么多叹号，显得极其浮躁，令人反感。本论坛是正经、专业的讨论场所，注意描述问题时的语气

先具体交代清楚体系是什么样，认真看
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）
没具体前提只能问算命的

为正本来就未必有什么必然性的错误，诸如算的是环张力特别大的分子，按照力场势函数的定义，显然势能会明显为正

appleliu

谢谢老师提醒。我的计算体系是用OPLS-AA力场模拟直链C6F14液相分子在铜基底表面的蒸发行为，采用的是参考文献https://pubs.acs.org/doi/10.1021/acs.langmuir.7b01410中描述的详细力场参数，我计算出来的密度可以和NIST数据库对应到，对应附件中的图1，就是我看许多类似长链分子的蒸发体系，计算得到的液相势能都是负值，而且我理解的是势能为负才表示系统稳定，所以我有点不确定我的体系是否正确。进一步我模拟了纯液相C6F14体系的弛豫过程，并计算了非键相互作用、成键相互作用，得到evdwl为几百的正值 ecoul为正值 ebond eangle也为正值，edihed为负值，具体数据如附件图2，示例代码如附件图3，代码按照文献中数据完全还原。请大家有空帮我看看我的计算是否合理，感谢大家百忙之中抽空查看，谢谢。