xo扩展oiom方法用于大分子体系的精确高效计算.pdf

摘 要 量子化学模拟的一个发展方向是用越来越精确的方法处理越来越大的体系．但精确 的量子化学方法限于严峻的计算标度，在大分子体系的应用中面临巨大障碍．组合方法 (hy嘶d 提出的ONIOM方法由于兼具技术上的简单性和应用上的灵活性，是量化程序包中集成 度最高和应用领域最广泛的组合方法之一．但由于链接原子边界效应的影响太过显著， oNIOM方法计算结果极其依赖于用户的设置．另一方面，ONIOM方法虽允许系统性地 改进计算精度，但往往需以明显增加计算量为代价．这些缺陷增加了用户在oNIOM计算 中的盲目性，应用于大分子筛体系计算也存在诸多限制．与ONIOM方法不同，分子片方 法较少受边界误差的影响．鉴于分子片方法与ONIOM方法实质的相通性，我们扩展了 oNIOM方法的灵活性以突破其固有的缺陷，而没有牺牲其自身的简单性．本博士论文工 作主要包括如下三个方面： 1．首次通过容斥原理将分子片方法与ONIOM方法结合起来，并据此提出扩展的 既能以oNl0M“洋葱皮”一样的形式组合，也能以分子片方法中部分重叠的方式划分． 同时XO的任一子模型，既可使用同一水平的方法描述，也可用完全不同的方法来处理． XO方法给予用户在体系划分和计算方法组合上最大的自由度，能够针对性地设计出比 单纯使用分子片方法或oNIoM方法精度更高效率更好的计算方案．从可模拟的性质而 言，XO方法不仅适合处理传统的具有活性区域的体系，也可有效计算体系整体的性质， 包括整个分子的几何构型，光谱频率以及热力学性质等．为了避免用户在组合方法计算 上的随意性和减少对计算方案可靠程度的盲目性，我们提出多种定性规则用于指导XO 计算方案的设计，并建议基于Force误差分析来定量地测试特定X0计算方案在几何构 型优化和频率计算上的表现，同时基于能量误差分析来定量测试特定计算方案在能量计 算和生成焓计算的表现． 代优化算法，用于减少X02构型优化所需要的步数，但不损失任何计算精度．特别地，对 于X02(H：H)计算，配合IFC算法可以有效地得到全簇模型方法相一致的几何构型，从而 达到XO方法所允许的精度上限．我们以分子筛为测试体系，对IFC算法进行详细的测 I 度比标准优化方式提高约十倍，这极大地提升了X02结构优化任务所能处理体系的大 小． 水平下，使用167T簇模型(包含557个原子)研究了钛硅分子筛TS．1中Ti原子的落位， T和126 在B3LYP／6．311G事牛水平下，使用了124T簇模型(分别包含416个原子和419 个原子)研究了MOR分子筛的酸性位结构．对于TS．1分子筛体系，我们构建单一的大簇 模型直接比较了各个T位的相对稳定性，避免使用多个小模型间接比较所受人为选择的 影响．我们计算得出各Ti位热力学分布，首次合理地解释了最新的实验结果．对于MOR 分子筛体系，我们得出的灿在四个T位取代的热力学相对稳定性与实验分析得出的分 布完全一致，同时计算出的桥羟基振动频率无论质子的具体落位还是数值大小均与实验 结果吻合一致．XO方法在分子筛体系的成功应用展示了XO计算策略的高效性和理论 预言的可靠性，也同时预示了XO方法在实践应用中的巨大前景． 关键词：ONIOM，XO，钛硅分子筛，MOR分子筛，酸性位 II Abstract isa demandformeoreticaJ of T0day，meregrowing inter】)retation System诵tll ever-la唱er costofaccurate scale precision．The eVer-greatcr qu砒咖【mmechallical(QM)metllodsuSllally thesizeofthe of en． ul