材料大规模计算的需求与挑战-北京应用物理与计算数学研究所.PDF

专题 第 10 卷  第 1 期  2014 年 1 月 材料大规模计算的需求与挑战 宋海峰　高兴誉　邵建立　王帅创 关键词 ：材料计算　第一原理　分子动力学 中国工程物理研究院 引言 的时间复杂度是体系原子数的立 两部分，主要包括大规模计算的 方。正因如此，大体系的第一原 需求与挑战、软件研制与应用。 材料作为社会进步和人类文 理计算一直是高性能计算的重点 明的物质基础和先导，与能源、 研究课题。如今，大规模第一原 大规模第一原理计算 信息并称为现代社会的三大支 理计算不仅面向大体系，长时间、 柱。材料的发现和变革影响了人 多尺度、高通量的计算需求以及 类文明的进程。随着信息科学和 强关联材料理论方法的发展都对 需求与挑战 技术的发展，计算技术在科学领 高性能计算提出了新的挑战。 密度泛函理论计算的核心与 域的广泛应用正改变着科学研究 分子动力学是基于材料晶体 热点是求解科恩- 沈(Kohn-Sham) 的方式和面貌。计算科学与科学 结构建模的确定性模拟方法。它 [1] 方程 。Kohn-Sham 方程是一个 实验、科学理论共同推动着现代 忽略电子运动的自由度，电子效 非线性特征值问题，通过自洽迭 科学技术的进步。材料的大规模 应体现在原子间的相互作用势函 代求解。每个自洽步计算一个线 计算基于第一原理方法，结合热 数中，原子的运动遵从牛顿定律。 性化的特征值问题，时间复杂度 力学、动力学、多尺度等方法， 因此，分子动力学可以直接追踪 2 为O (N M ) ，其中N 是需要求解 充分发挥超级计算和信息技术的 原子的运动轨迹，从而获得现象 的特征对（特征值及其特征函数） 优势，将会极大地加速材料的研 的各种微观变形机理，尤其在极 个数，M 是离散的自由度。N 和 究和新材料的开发。 端加载条件 （高温、高压）下显 M 均与体系的原子数或价电子数 第一原理方法通过求解量子 示出其不可替代的作用。已有研 成正比，M 一般远大于N 。对于 力学方程，从微观电子结构研究 究表明，对于大多数金属，基于 狭义的大规模第一原理计算，N 材料物性。密度泛函理论以电子 原子间多体势的经典分子动力学 3 6 通常大于10 ，M 达到10 ，单自 密度为基本参量，针对多电子体 方法可以给出与实验结果契合的 洽步的计算量在 1012 以上。大体 系建立了一种可计算的模型。在 热力学和动力学性质。20 世纪 系的第一原理计算主要源于三方 过去的20 多年时间里，基于密 90 年代以来，分子动力学模拟在 面需求 ：(1) 热力学模拟中精确 度泛函理论的第一原理方法得到 现代科学研究中一直扮演着重要 描述原子的受力运动，周期单元 了迅猛发展与广泛应用，成为材 角色，已成为沟通理论与实践、 需要包含足够多的不等价原子 ； 料科学、凝聚态物理、纳米科 宏观与微观的主要桥梁之一。 (2) 缺陷、掺杂、非电中性破坏 学、化学等领域不可或缺的计算 本文介绍的材料大规模计算 了晶体的周期性，需要选取足够 方法。标准的密度泛函理论方法 分为第一原理和经典分子动力学 大的周期单元 ；(3) 纳米科学中 26