md计算步骤以及电荷选取提问

GHR1111

最近在计算电解液的md，使用的是lammps opls力场，在进行md前统一对有机小分子进行优化，在Gaussian中使用的def2tzvp进行的结构优化，并添加了smd，设置了溶剂的eps。然后导入ligpargen获取力场参数，multiwfn获取1.2cm5的电荷。想请问以上步骤是否合理。def2tzvp精度是否太高，因为组内师兄也有在md前用ccpvdz进行结构优化的，精度太高是不是也不太好。第二个问题是，阴阳离子是否也需要在溶剂的smd下进行结构优化，对于阴阳离子的力场参数和电荷我直接采用的文献中opls力场下的值，想请问阴阳离子电荷是否可以直接用文献中的，还是需要结构优化完后重新计算cm5电荷呢

Stardust0831

def2tzvp够用的，计算级别统一就行、不存在“精度太高是不是也不太好”的问题，你能算动结构优化那说明分子也都不大(我个人倾向于B3LYP-D3/def2-SVP结构优化、B3LYP-D3/def2-TZVP单点，省一些算力开销，此处基组可以不统一的原因见：浅谈为什么优化和振动分析不需要用大基组）。ccpvdz是给post-HF等理论方法用的，给普通杂化泛函有些浪费。如果只是用于动力学模拟初始结构的建模，只需要结构合理即可，不需要细致的dft的优化，文献里的力场如果兼容、可以直接用文献的。