编者按:中科院之声与中国科学院上海硅酸盐研究所联合开设“科普硅立方”专栏,为大家介绍先进无机非金属材料的前世今生。我们将带你——认识晶格,挑战势垒,寻觅暗物质,今古论陶瓷;弥补缺陷,能级跃迁,嫦娥织外衣,溢彩话琉璃。相信大家对于炒菜一定不陌生,就算自己没亲自炒过,也一定见看过父母、朋友或者是饭店里的厨师炒过。
百度百科里面给了“炒菜”相对普适性的的定义:炒菜,是中国菜的常用制作方法,将一种或几种菜在特制的凹型锅内,以火传导到铁锅中的热度为载体,加入油、佐料和菜品后用特制工具“锅铲”翻动将菜炒熟的烹饪过程。在传统的科学研究中采用的方法往往与“炒菜”模式类似。试验中处理的模型可类比炒菜用的锅,而物理思路和理论分析算法看作是主料,其他各种参数的调节和细节的处理大致就是调料了。
在试验中,改变“调料”,改用其他锅或主料,得到的结果也经常不尽相同。但是,这种“炒菜式”的科研方式存在着一个致命的短板,那就是需要消耗大量的人力、资源、时间和金钱,且极大地制约了新材料的发展速度。此时,计算材料学应运而生。
所谓计算材料学(ComputationalMaterials Science),是一门综合了材料科学与计算机科学且正处于快速发展的新兴学科,是关于材料组成、结构、性能、服役性能的计算机模拟与设计的学科,是材料科学研究里的“计算机实验”。它涉及材料、物理、计算机、数学、化学等多门学科。
通俗一点解释为,计算材料学是通过计算机将理论知识编程各种软件程序,然后利用该软件进行一定的科学研究,结合一定的实验数据,快速找出其本质的物理意义,从而完善理论,进一步改善和加速实验研究。计算材料学可以利用不同的模型和计算方法,从不同尺度与时域研究材料的结构演化及其与性能的相互关系,从而更好地理解和设计新材料。具体的计算方法很多,从空间尺度的不同,分为微观尺度、介观尺度和宏观尺度。
微观尺度(原子尺度及小于晶粒尺寸),利用密度泛函理论(Density Functional Theory, DFT),可以通过电子态密度(Density of Electronic States,DOS)、巴德电荷(BaderCharge)以及局域电荷密度等很好地将电子信息与材料性能和结构联系在一起。
另外,也可以选择分子动力学(Molecular Dynamics,MD)和蒙特卡罗方法(MonteCarlo, MC)来寻找原子的迁移路径和势垒。针对具有十分复杂的界面结构的问题时,用经典尖锐界面模型去跟踪界面演化,会遭遇到严重的数值困难。
并且真实材料中的相界或晶界实际上并不是严格的零厚度界面,而是具有一定厚度(纳米尺度)的边界层,这层厚度控制材料相变动力学,由此引入一个序参量场区分两相(如固相和液相),该方法为相场法(Phase-field method, PFM),相场法被广泛应用于模拟液/固相变(枝晶生长、多元多相凝固、多晶凝固等)、固态相变(沉淀相析出、马氏体相变、铁电相变等)、应力相变(薄膜生长,定向粗化)、结构缺陷相变(裂纹扩展、位错动力学)等。
在微观尺度内,还有微观动力学模拟等方法。介观尺度(晶粒尺寸大小),主要有元胞自动机法,以晶粒生长为例,元胞自动机的模拟过程主要有确定形核、选择时间步长、计算温度场、计算溶质场、计算界面元胞的生长动力学、确定胞状态、捕捉邻胞、计算曲率等过程组成。
宏观尺度(宏观试样大小),主要是有限元分析(Finite Element Analysis, FEA)方法,该方法是利用数学近似的方法对真实物理系统进行模拟,还利用简单而又互相作用的元素就可以用有限数量的未知量去逼近无限未知量的真实系统。计算材料学在化学、材料等学科中运用非常广泛,常见量子化学模拟软件包括:VASP、Materials Studio、Gaussian等。
这些软件主要针对分子层面的微观问题的模拟计算,在给定分子结构的前提下可以计算材料的光谱性质、电子结构、分子不同构象的能级等等。计算材料学的一个重要发展方向是材料基因组计划,其研究的关键是实现材料研发的“高通量”,即并发式完成“一批”而非“一个”材料样品的计算模拟、制备和表征,实现系统的筛选和优化材料,从而加快材料从发现到应用的过程。
这些原先难以通过实验完成的工作,现在可以很方便地使用超级计算机来完成。