楼主有没有联系方式留一个，想请教一下molclus+xtb的问题

楼主还有个问题，你这个模拟就是直接GAFF力场，有没有做其他修改，我使用GAFF力场30ns MD之后，阴阳离子分不开，Li周围都是PF6，几乎没有PC和FEC，这是咋回事呢？能看下你的itp以及mdp么，非常感谢，可以加我qq863618483.

nianbin 发表于 2021-6-23 23:06
这篇文献能发来看下么

Molecular Dynamics Simulation of Li+BF4- in Ethylene Carbonate, Propylene Carbonate, and Dimethyl Carbonate Solvents

h840473807 发表于 2019-4-21 11:50
看到一篇文献，电解液用的OPLS力场，我按照文章里Li原子的参数，将非键项改成sigma=1.46e-1 epsilon=7.94e- ...

这篇文献能发来看下么

看到一篇文献，电解液用的OPLS力场，我按照文章里Li原子的参数，将非键项改成sigma=1.46e-1 epsilon=7.94e-1（原来为sigma=2.03e-1 epsilon=7.66e-2），Li+就主要是四配位了的，就不知道这样做合适不合适

h840473807 发表于 2019-4-5 15:58
确实之前没有考虑过这个体积摩尔浓度的问题，只是随便跑了一个类似的浓度。我感觉和模拟的实际体系的质量 ...

好的，谢谢。

wbn 发表于 2019-4-5 04:29
MD金属离子的点电荷模型不太符合实际情况，所以它和配体的nonbonded interaction 和实际不符是很常见的...  ...

谢谢，我先看看文献有什么适合这一体系的力场吧

Frank 发表于 2019-4-5 06:57
谢谢回复。

我现在也在尝试用packmol构建一个0.1M的电解质溶液体系，我按照文献中的方法用的是11.8A c ...

确实之前没有考虑过这个体积摩尔浓度的问题，只是随便跑了一个类似的浓度。我感觉和模拟的实际体系的质量摩尔浓度保持一致就差不多吧。如果有什么好方法可以讨论

h840473807 发表于 2019-4-4 17:33
http://bbs.keinsci.com/thread-12549-1-1.html
照此博文，用packmol构建，我加了500个EC分子,53个LiPF6 ...

谢谢回复。

我现在也在尝试用packmol构建一个0.1M的电解质溶液体系，我按照文献中的方法用的是11.8A cubic的盒子里面放100个电解质分子加20K个溶剂分子，但是我在想构建完模型之后跑NPT的时候V就变了，电解质分子的数目不变，就不能维持0.1M的电解质浓度了，请问您有没有考虑过这个问题？

MD金属离子的点电荷模型不太符合实际情况，所以它和配体的nonbonded interaction 和实际不符是很常见的... 你最好查查文献看有没有更符合你体系的Li原子力场。实在没办法的话只能scale原子的charge或者调节LJ参数直到和quantum的结果对上。要不可以试试在Li原子或EC氧原子周围加几个mass和LJ参数为0的ghost原子模拟一下空SP3轨道或者孤电子对，就是不知道这样会不会被审稿人怼

Frank 发表于 2019-4-4 17:14
借楼请教一下，楼主是跑MD之前怎么设置电解质和溶剂分子的数目的？

http://bbs.keinsci.com/thread-12549-1-1.html
照此博文，用packmol构建，我加了500个EC分子,53个LiPF6,盒子46 A^3

借楼请教一下，楼主是跑MD之前怎么设置电解质和溶剂分子的数目的？