本人做Li+电池，使用gmx模拟固态电解质和Li金属之间的界面性质，所以需要通过gmx构建金属的界面。
我在阅读这篇文献时，Accurate Simulation of Surfaces and Interfaces of Face-Centered Cubic Metals Using 12-6 and 9-6 Lennard-Jones Potentials （ J. Phys. Chem. C 2008, 112, 17281–17290 ），
发现该文献提出金属可以仅通过平衡非键距离r0和平衡非键能eplison0进行模拟（例如opls-AA力场等）。
我想：那是不是可以直接用晶体数据库中已知的金属原子间的非键平衡距离r0和平衡非键能eplison0代入到力场中，就不需要再经过反复的调试和模拟了？
还是说，由于力场的函数形式和实际情况存在差异，r0和eplison0的真实值不适合直接用作力场的参数？
十分感谢好心人解惑！作者Author: sobereva 时间: 2022-11-4 08:43
opls-AA根本不适合描述金属单质