论文部分内容阅读
水在生命体中扮演着极其重要的角色,水分子间存在大量的二体、三体等多体相互作用。通过研究这些多体作用,能够帮助我们进一步发展高精度分子模拟方法。已有研究表明,水团簇中的二体作用能占总作用能比例最大,三体、四体及更高阶相互作用贡献较小,因此,二体作用能的准确预测非常重要。研究表明,使用CCSD（T）/CBS方法可以得到非常准确的二体作用能,但对于大的水团簇,其存在的二体作用数量庞大,使用CCSD（T）方法计算所有二体作用仍很耗时。所以找到一种快速准确预测大水团簇二体作用能的方法十分必要。本文使用包含基组重叠误差（BSSE）校正的MP2/aDZ、MP2/aTZ方法以及MP2/CBS、MP2.5/CBS方法计算了水团簇（H2O）n（n=2-7,10,12,14,16,18,20）中的二体作用强度,还采用CCSD（T）/CBS方法计算了水团簇（H2O）n（n=2-7）中的二体作用强度。探讨了水团簇二体作用强度与分子间距(用两氧原子距离ROO表示)的关系,提出了一种能够快速准确预测水团簇二体作用强度的组合方法。本文对较小水团簇（H2O）n（n=2-7）的计算结果表明,MP2.5/CBS方法计算得到的水团簇的总二体作用能可以媲美CCSD（T）/CBS方法的结果,所以本文选取MP2.5/CBS方法的计算结果为基准。本文研究发现,与MP2.5/CBS方法计算得到的二体作用能结果相比,对于水团簇中远程(ROO>0.41 nm)二体作用,MP2/aDZ、MP2/aTZ、MP2/CBS三种方法均能给出很准确的结果。但是对于水团簇中近程(ROO≤0.41 nm)二体作用,MP2/CBS方法的准确性要明显优于MP2/aDZ和MP2/aTZ方法。本文计算结果还表明,对本文研究的所有32个水团簇体系,使用MP2/aDZ方法计算水团簇中ROO>0.41 nm的二体作用、使用MP2/CBS方法计算水团簇中ROO≤0.41 nm的二体作用,这样得到的总二体作用能结果与MP2.5/CBS方法的总二体作用能的相对平均偏差仅为1.1%,最大相对偏差仅为1.5%,说明可以使用组合方法来较快速地计算得到水团簇的总二体作用能:即对近程二体作用使用较费时的MP2/CBS方法计算,对远程二体作用使用较省时的MP2/aDZ方法计算。本文研究还发现,以上组合方法计算得到的总二体作用能,加上MP2/aDZ方法计算得到的总三体作用能,结果得到的近似总作用能与MP2.5/CBS方法总作用能相比较,均方根偏差仅为1.43 kJ/mol,相对平均偏差仅为0.6%。