关键词:
第一性原理计算
结构预测
机器学习势场模型
高斯过程回归算法
硼团簇
CALYPSO
摘要:
团簇作为连接微观世界与宏观世界的桥梁,其由数量不等的分子或是原子组成。团簇由于表面效应和量子限制效应,其表现出来的物理与化学性质与晶体迥然不同,而且其物理化学性质以及基态结构的构型与团簇的尺寸是紧密相关的。团簇领域的重要问题之一是团簇的生长模式的探索,然而对于小尺寸体系而言,我们可以通过第一性原理计算方法进行理论模拟,但是对于大尺寸体系,我们就显得束手无策了,这是一直困扰着科学家们的问题。硼元素由于其“缺电子”特性,导致原子之间的成键模式非常复杂,不同尺寸下的硼团簇的基态结构构型吸引着诸多科学家对其进行研究。目前,硼团簇已经发现了许多的基态候选结构的构型,且构型之间存在着竞争的关系,这使大尺寸的硼团簇的基态结构探索变得极具挑战性。但是,受限于计算开销,我们很难对大尺寸的硼团簇体系展开大规模的基于第一性原理计算的理论模拟。基于第一性原理计算的理论模拟已经成为结构预测领域的一种主要手段,其能够准确地预测结构的许多物理性质,能够帮助我们克服苛刻的实验条件的桎梏。如今,随着理论模拟的需求不断增加,即大尺寸体系的大规模模拟,我们在采用第一性原理计算方法进行理论模拟的过程中不得不面对其昂贵的计算开销,很难对大尺寸的体系展开计算,而采用经验势模型,其可移植性以及精度都大大限制了其应用范围。至此,对于大尺寸体系而言,发展一个精度高且计算开销小的原子间势能面模型是必不可少的。近年来,机器学习方法的蓬勃发展为科学家们解决此问题提供了一种全新的思路,即发展机器学习势场模型,不但可大大减小计算开销,提高计算效率,同时可以保证很高的精度。在本文的研究工作中,我们主要实现了以下两项工作:第一:我们搭建了基于高斯过程回归算法的机器学习势场模型,即高斯回归势场模型。通过验证集的表现,以及与成熟的神经网络模型之间的对比,证明了其有效性。第二:我们针对大尺寸硼团簇体系,展开了大规模的基态结构搜索。我们采取课题组自主研发的CALYPSO结构预测方法与第一性原理计算方法相结合的方案,产生了17616个硼团簇B_n(n=10–84)的数据集以此来进行机器学习势场模型的拟合。首先,我们在硼团簇Β101体系上,进行了测试,证明方法的可信度。随后,对于硼团簇B_n(n=58-108),我们计算得到了不同尺寸下的基态候选结构,并以此谈论了大尺寸硼团簇B_n(n=82-108)的生长模式。本文的方法为结构预测领域的理论模拟提供了一种新的思路,本文的数据结果则为大尺寸硼团簇体系的生长模式的探索提供了参考价值。