分子统计热力学方法（MST）：高效计算纳米材料等温准静.pdfVIP

中国科学院力学研究所 非线性力学国家重点实验室《简讯》20 12第二期(总第203 期) 科研进展方法论 分子统计热力学方法 (MST)：高效计算纳米材料 等温准静态力学行为 国际复合材料期刊（Composites: Part B）, 2012 年 B 卷上发表了中国科学院 力学研究所非线性力学国家重点实验室肖攀博士等介绍分子统计热力学方法 (MST) 的论文，该方法可在几类纳米材料力学中应用，突显出它的可靠性和高效性。 纳米固体力学中常用的计算方法是分子动力学(MD)。MD 能有效表征原子的动力 学过程，但在时间尺度上受原子振动特征时间尺度(约 10-12 s)制约，难以模拟准静 态的实际过程；在空间尺度上，难以应用到百纳米以上体系。分子/集团统计热力学 (简称 MST/CST)多尺度计算方法，适用于有限温度下微/纳米材料的准静态行为分析。 这一计算框架包含：采用分子表象的 MST 方法、采用准连续表象的 CST 方法以及采 用分子/准连续联合表象的 MST/CST 耦合方法。发表在 Composites: Part B 上的论 第 1 页 中国科学院力学研究所 非线性力学国家重点实验室《简讯》20 12第二期(总第203 期) 文主要讨论该框架的核心：MST。 论文综述性地介绍了 MST 方法在几类纳米材料的力学性能研究中的应用：面心 立方晶体纳米材料的单轴加载及纳米压入过程；半导体纳米材料的结构相变及力学 性能尺寸效应。计算结果表明，MST 方法能有效刻画材料在力学加载中的微结构变化 过程，如位错、相变等行为，并具有相比 MD 方法更高的计算效率(如图 1(a))。在 ZnO 纳米线力学性能计算中，MST 方法成功刻画了拉伸过程中纳米线由纤锌矿向四角 结构的转变，其结果与 MD 及第一原理结果一致。更重要的是，MST 表现出比 MD 高约 60 倍的计算效率，因而能计算更大尺寸的体系。例如，采用 MST 计算到了目前实验 所能测量的 ZnO 纳米线尺寸，其计算结果与实验结果一致性地揭示了弹性常数的尺 寸效应 (如图 1(b)所示，实验结果由共同合作的清华大学朱静院士课题组完成)。 图 1. (a) MST计算纳米杆压缩过程的耗时随原子总数的变化以及与 MD 方法的对比； (b) ZnO纳米杨氏模量尺寸效应计算值与实验结果的对比。 该研究成果将为微/纳米材料力学性能研究提供高效可靠的计算工具，便于模拟 第 2 页 中国科学院力学研究所 非线性力学国家重点实验室《简讯》20 12第二期(总第203 期) 与实验结果的直接对比。 论文中所提及的 MST 计算软件已进行软件著作权登记（登记号：2011SR003778）。 MST 软件并行版本采用 C++面向对象程序开发，近 2 万行代码，具有良好的通用性和 扩展性。软件可以免费下载使用。 该研究得到科技部 973 计划 (2007CB814803)、国家自然科学基金 10732090 和和中国科学院创新基金 (KJCX2-YW-M04)的资助。 论文信息： Pan Xiao, Jun Wang, Fujiu Ke and Yilong Bai, “Molecular statistical thermodynamics –A distinct and efficient numerical approach to quasi-static analysis of nanomaterials at finite temperature”, Composites: Part B, 43 (2012), 57-63. 实验研究进展 力学所在 “电弹性毛细”实验研究中取得进展 2012 年 4 月 25 日，美国《科学》杂志网站主页以 “液滴容器将使血液检测不再 痛苦”为题，报道了中科院力学所赵亚溥研究员 “纳微系统力学与物理力学”课题 组关于 “电弹性毛细”的实验研究结果（图 2）。该报道指出：电弹性毛细 (Electro-Elasto-Capillarity, EEC)将有效地避免待检测