问题：gmx应用LJ12-6势模拟金属界面与表面

王纪峰

本人做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的真实值不适合直接用作力场的参数？
十分感谢好心人解惑！

sobereva

opls-AA根本不适合描述金属单质

这跟晶体数据库里的原子间距也没关系

你实际模拟时，直接用金属的cif文件扩胞成你需要的slab或其它什么模型，并在gromacs的拓扑文件里直接用文献里描述相应纯金属适合的LJ参数就完了（可以是C6、C12形式，也可以是r0、epsilon形式，也可以是sigma、epsilon形式，具体看文章）