论文部分内容阅读
针对中大体系的能量和电子结构计算,我们提出了一种将高精度的WFT方法和低开销的DFT方法相结合的局域杂化分块计算(LHDC)的DFT方法,该方法首先将整个体系划分为小的子体系,然后利用高精度的WFT方法能量对每个子体系的DFT方法杂化参数进行优化。为进一步提高杂化参数的优化精度,我们还引入了子体系之间的相互作用校正。最终,这些优化后的局部杂化参数被应用在整个分子的DFT方法计算中。在本工作中,我们对一系列的长链烷烃长链、烯烃以及苯二聚体进行了测试。测试结果表明,在能量计算中,LHDC方法可以获得与所使用的WFT方法相近的结果,即使是针对离域化很强的全共轭体系,势能曲线与WFT方法的相对误差也令人满意(小于1kcal/mol)。LHDC方法除了受到所使用的WFT方法和DFT方法影响外,还受到DC分区策略和基组的影响。在DC分区策略中,选择较大的子体系区域或者引入更高阶的多体校正都可以提高精度。其原因在于这两种方案都可以降低被忽略掉的非相邻区域间的相互作用。在基组的影响测试中,我们发现杂化参数的优化值随着基组逐步增大表现出与能量相似的收敛性。此外含不同元素的分子具有差异较大的杂化参数优化值。LHDC方法的计算量包含整个体系的DFT计算、子体系WFT能量和相互作用的计算以及杂化参数的自洽场优化计算。我们的计算表明,当使用二体校正方法LHDC2时,计算开销仅仅比标准DFT方法略高,但可得到精确的WFT方法的精度。从WFT方法的角度看,我们的LHDC方法使用DFT作为不同区域间的胶合方案,并通过二体或者三体的校正来修正区域之间的胶合强度,这样可以确保大体系WFT能量值的精确重现。另外最终的整体DFT计算还可以有效解决DC方法难以解决的电子密度在区域间的重新分配问题,并且可以在整体上给出分子的各种性质。从DFT方法的角度看,每个子体系的杂化参数都被重新优化过,这可以进一步提高密度泛函在分子局部的精度,因此我们的方法可以被看做是一种新的局域的杂化参数DFT方法。