计算化学公社

标题: 如何准确地计算有大分子和水和离子的体系？ [打印本页]

作者
Author: kazusa 时间: 2019-12-18 14:24
标题: 如何准确地计算有大分子和水和离子的体系？
大家好，我想要计算的是两个六方BN板和板间放入的水合离子,两个BN板有216个原子，体系很大。
因为主要是为了计算离子-pi相互作用，所以只能用量化计算。
我先用的是ORCA的B97-3c进行的计算，但是两板会把水合离子往板的边缘挤，然后计算会不收敛或者难以收敛。我觉得可能是由于没加入隐式溶剂模型，导致没有计入背景水的效应才会有这么不合理的结构，于是我在ORCA中加入了SMDsolvent “water”的隐式溶剂进行再次计算，结果发现算不了几步就会出现SCF报错。

我现在在尝试用moleclus和MOPAC2016产生多种初始结构，再用ORCA的B97-3C进一步算。用的是70个水和两个BN板还有一个锂离子，但是我有两个担心：
1、 70个水我感觉不够表现BN板在水中的效应，但是加太多我感觉ORCA的优化又难以进行。
2、 MOPAC2016计算了快一天了才出了3个结果，我计划的是100个结果，这需要的时间太长，而且也不知道之后ORCA的计算要花多长时间。

所以想问下各位对这种体系有没有好的办法进行相对准确的计算？

作者
Author: sobereva 时间: 2019-12-19 16:12
可尝试用Gaussian，SCF容易收敛得多

是molclus不是moleclus

只需要在离子周围加上显式水，更远部分的溶剂效应通过隐式溶剂模型表达就够了

MOPAC还慢就用xtb

作者
Author: kazusa 时间: 2019-12-19 16:31

sobereva 发表于 2019-12-19 16:12
可尝试用Gaussian，SCF容易收敛得多

是molclus不是moleclus

好的，社长

但是这么大的体系Gaussian里的方法算需要太长时间了

作者
Author: sobereva 时间: 2019-12-19 18:01

kazusa 发表于 2019-12-19 16:31
好的，社长

但是这么大的体系Gaussian里的方法算需要太长时间了

两百多个原子，有个三四十核的机子，用Gaussian照样能优化得动

作者
Author: kazusa 时间: 2019-12-20 09:07

sobereva 发表于 2019-12-19 18:01
两百多个原子，有个三四十核的机子，用Gaussian照样能优化得动

嗯嗯，我已经在用Gaussian优化了，谢谢社长

欢迎光临计算化学公社 (http://bbs.keinsci.com/)