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的优化，文献里的力场如果兼容、可以直接用文献的。

sobereva

cc-pVnZ系列结合DFT使用又贵又烂（相对于合适的def/def2等合适的搭配来说），做DFT计算的情况没事甭用它。下文明确说了
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
你师兄缺乏基组选用的常识

MD前并不需要对小分子做优化。至于牵扯到Multiwfn计算1.2*CM5电荷的情况，下文的脚本里已经自带了优化步骤，里面自动用的计算级别就是非常合适的
计算适用于OPLS-AA力场做模拟的1.2*CM5原子电荷的懒人脚本
http://sobereva.com/585（http://bbs.keinsci.com/thread-21462-1-1.html）
ORCA结合Multiwfn计算RESP、RESP2和1.2*CM5原子电荷的懒人脚本
http://sobereva.com/637（http://bbs.keinsci.com/thread-28178-1-1.html）

有机分子优化撑死了用def-TZVP也足矣，def2-TZVP严重浪费

另外，不建议优化时用SMD，下文说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
建议用IEFPCM

文献里有靠谱的对某个体系专门搞的原子电荷的情况可以用文献里的。