当地时间11月15日下午，中国科大入围国际戈登•贝尔奖提名的超算应用成果 “2.5 Million-Atom Ab Initio Electronic-Structure Simulation of Complex Metallic Heterostructures with DGDFT”，参加了在美国得克萨斯州达拉斯市召开的全球超级计算大会（SC22）的线上答辩，报告了中国科大在新一代神威超级计算机上采用低标度高精度第一性原理计算软件DGDFT首次实现250万原子复杂金属异质结（锂钠金属合金）的超大规模材料模拟，这是中国科大第二次作为第一完成单位入围戈登•贝尔奖提名。 戈登•贝尔奖是国际上高性能计算应用领域的最高学术奖项，被誉为“超算领域的诺贝尔奖”，由ACM（美国计算机学会）每年评选和颁发，奖励在高性能计算应用领域具有突出成就的工作，具有巨大的国际影响力。此项成果由合肥微尺度物质科学国家研究中心杨金龙院士课题组胡伟团队与计算机科学与技术学院安虹教授课题组联合攻关，在崂山实验室（原青岛海洋科学与技术试点国家实验室）、中国科学院计算技术研究所、北京大学、中国科学院软件研究所、齐鲁工业大学，以及国家并行计算机工程技术研究中心相关研究人员的紧密配合下完成。
图1ACMGordonBell奖提名 先进材料是国民经济的基石，是实现制造业转型升级的重要基础。由于新型材料的研发难度较大，用于验证材料性质的物理实验设计复杂、代价昂贵。量子力学和高性能计算技术的出现从根本上改变了这种状况。通过输入材料的结构信息，利用量子力学的第一性原理计算就可以较为准确地预测已知材料的基态结构和基本的物理化学性质，并实现原子级别的精准控制。这是二十一世纪解决实验理论问题和预测新材料结构性能的最具竞争力的方法和技术途径。这种方法不需要开展真实的实验，不仅可以极大地节省实验成本，缩短新材料的开发周期，而且可以为材料的制备和改性、新材料的开发，以及对极端环境下材料性质的研究提供有效的理论指导。 1933年，薛定谔和狄拉克因发展了量子力学的薛定谔-狄拉克方程而获得诺贝尔物理学奖。尽管在量子力学建立之时狄拉克就曾预言寻求物理和化学基本规律的任务已大体完成。但是，由于描述这些基本规律的方程求解过于复杂，用基本原理来解决实际问题依然非常困难。直到上世纪60年代出现了高性能计算(HPC)，我们也只能解出有几十个原子的薛定谔方程，这离真正的复杂系统相距遥远。 1998年，Kohn和Sham发展了基于Kohn-Sham方程的密度泛函理论（DFT），将3N维波函数问题简化为三维粒子密度问题，从而将从头算电子结构模拟的计算复杂度降到O(N^3)，获得了诺贝尔化学奖。2013年，三位化学家Martin Karplus, Michael Levitt和Arieh Warshel因在现代超级计算机上开发出复杂系统的多尺度建模方法而获得诺贝尔化学奖。Kohn-Sham方程第一次将量子力学引入到包含数千个原子的真正复杂系统。但是，由于第一性原理材料模拟的计算复杂度会随着模拟材料的尺度而急剧增加，研究人员对软件的性能和计算资源的需求越来越大。 超级计算机和高性能计算技术的快速发展为第一性原理计算的发展提供了机遇，使其在凝聚态物理、化学、材料和生物等研究领域发挥着越来越重要的作用。在后摩尔时代，从头算电子结构模拟成为我们理解现代信息技术的唯一方法。全量子力学模拟对于设计10纳米以下的下一代场效应晶体管来说是必不可少的，但要求模拟的尺度至少要达到数十万个原子。近三十年来，现代高性能计算技术已经将从头算建模带到了真实的物理世界。这种能力对大多数的科学和工程领域十分重要，如能源、材料、生物医学、催化反应和强关联物理。 那么，在当代最先进的超级计算机上，利用Kohn-Sham方程我们究竟可以模拟多大尺度的系统呢?实现大规模从头算电子结构模拟主要有三种计算方法，分别是线性标度算法、人工智能算法和通用立方标度算法。采用这三种算法的研究工作都曾获得过戈登·贝尔奖或戈登·贝尔奖提名。线性标度算法和人工智能算法对分子、半导体和绝缘体的模拟过度依赖于近似原理，因而不能应用于复杂的金属体系。而目前最先进的采用通用的立方标度算法的DFT-FE软件，截止2019年，只能模拟11K个原子的固体材料，比线性标度算法和人工智能算法要小得多。 采用基组离散Kohn-Sham方程是实现第一性原理电子结构计算的基础。传统基于平面波基组的第一性原理材料模拟具有三阶的计算复杂度，很难通过大规模并行加速应用于复杂的大体系。基于局域原子轨道基组的线性标度算法可以实现大尺度模拟，但是往往精度不够并难以应用于金属体系，同时还面临非规则稀疏矩阵难以实现高度并行计算的难题。全新的自适应性局域基组（AdaptiveLocalBasis,ALB），在实空间严格截断且正交，结合了平面波（正交完备性）和数值原子轨道（局域性）各自的优点，既有媲美平面波的高精度又能利用适合于大尺度计算的线性标度算法。并且，由ALB基组构造的哈密顿矩阵具有固定的三对角块状的稀疏特性，适合于实现高度并行。 中国科大与中科院计算所、中科院软件所等单位合作，以ALB基组和间断伽辽金有限元方法为基础，发展了第一性原理计算软件DGDFT（Discontinuous Galerkin density functional theory），结合多级 并行优化设计和最先进的稀疏矩阵求解器PEXSI的低标度对角化算法，在数值算法、稀疏矩阵并行和数据通讯等方面进行一系列关键技术创新，采用在新一代神威超级计算机上自主设计的DGDFT软件，首次对复杂金属体系实现了250万原子的高精度第一性原理电子结构模拟，达到了介观尺度（>100nm）范围，可用于设计基于二维材料的新一代电子晶体管。