金属银熔体快冷过程的计算机模拟

利用分子动力学模拟技术研究了由周期性边界条件控制的500个银(Ag)原子构成的金属Ag熔体快速冷凝过程.冷却速率为8×1013K/s.模拟在FS(Finnis-Sinclair)相互作用势的基础上,通过双体分布函数、键对分析技术、键取向序等多种方法,对液银快速冷凝过程的微观结构转变特性作了分析,给出了连续快速冷凝过程中液银原子间依靠相互作用力形成的独特的微观结构图像.并考察了冷却过程中体系能量和元胞体积随温度的变化.模拟结果表明在快速冷凝过程中液Ag没有形成bcc结构的倾向.     

〈111〉晶向冲击加载下单晶铜中纳米孔洞增长的早期动力学行为

利用分子动力学方法模拟计算了单晶铜中纳米孔洞在沿〈111〉晶向冲击加载下增长的早期过程.测量发现不同加载强度下等效孔洞半径随时间近似成线性变化.观测到单孔洞增长的两种位错生长机理：加载强度较低时，只在沿着冲击加载方向的孔洞顶点附近区域有位错的成核和运动；而随着加载强度超过一定阈值，在沿冲击加载和其垂直方向的孔洞顶点区域都观察到位错的成核和运动.在前一种机理作用下，孔洞只沿加载方向增长；在后一种机理作用下，孔洞同时沿加载和垂直于加载方向增长.分析孔洞表面原子的位移历史，发现沿加载及与其垂直方向的孔洞顶点沿径向的速度基本恒定，由此提出了一个孔洞生长模型，可以解释孔洞增长的线性生长规律. 关键词： 纳米孔洞 分子动力学 冲击加载 位错     

燃料电池:从电催化到关键材料

介绍了在燃料电池关键材料与电催化方面的一些研究进展,特别是理论计算与实验研究相结合对燃料电池催化剂构效关系和聚合物电解质的研究.主要内容包括1 nm Pt颗粒的构效关系;Pd催化活性的电子调控与形貌剪裁;以及使用非贵金属催化剂的碱性聚合物电解质燃料电池研究等.     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构     

分子动力学和丙氨酸变异研究MDM2与抑制剂P4的结合模式

抑制p53-MDM2相互作用已经成为癌症治疗的新途径.采用分子动力学模拟和MM-PBSA(molecular mechanics-Possion-Boltzmann surface area)方法研究了肽类抑制剂P4与MDM2的结合模式.研究表明P4与MDMD2疏水性裂缝的范德华作用是抑制剂结合的主导力量。采用丙氨酸变异计算研究了P4与MDM2的作用热区.结果表明残基Lys51,Leu54,Leu57,Ile61,Met62,Tyr67,Gln72,His73,Val93,His96和Ile99的丙氨酸变异导致了范德华作用的降低,而对极性相互作用几乎没有产生影响,同时也证明这些残基处于P4与MDM2作用表面的热区,对抑制剂的结合有重要贡献,这能为抗癌药物的设计提供理论上的指导.     

非晶氧化钛薄膜形成过程中钛离子能量对表面结构影响的机理

研究了非晶氧化钛薄膜沉积过程中入射钛离子能量对表面结构形成机理以及薄膜特性的影响.模拟结果表明,通过提高入射钛离子能量,可以有效降低成膜表面粗糙度,从而减小薄膜表面的光学散射损耗.研究发现,当入射离子能量提高后,薄膜生长模式从"岛"状生长过渡到了"层"状生长,且离子入射点附近的平均扩散系数也有显著增加,这有利于形成更加平整的高质量薄膜表面.     

(n,n)-(2n, 0)碳纳米管异质结的扭转力学特性

相近直径的锯齿型和扶手椅型碳纳米管可以共轴组合形成5-7碳环交替出现的柱形对称异质结. 本文利用分子动力学方法研究了直径相近且等长锯齿型和扶手椅型碳纳米管形成的(n, n)-(2n, 0)结在扭转过程中的扭矩和轴向应力随扭转角度的变化规律以及应力传递过程. 研究发现, (n, n)-(2n, 0)结扭转应变在达弹性限度内不会产生轴向应力, 该效应对基于碳纳米管扭转特性的纳米振荡器件的设计具有重要意义.     

α-Fe裂纹的分子动力学研究

通过分子动力学方法，模拟了α-Fe裂纹的单轴拉伸实验中的形变过程.研究了不同晶体取向裂纹的形变特点和断裂机理，观察到各种形变现象，如位错形核和发射，位错运动，堆垛层错或孪晶的形成，纳米空洞的形成与连接等.计算结果表明，裂纹扩展是塑性过程和弹性过程相结合的过程，其中塑性过程表现为由裂尖发射的位错导致的原子切变行为，而弹性过程的发生则是由无位错区中的原子断键所导致.同时还研究了α-Fe裂纹的形变特点和断裂机理与温度场和应力场的依赖关系.     

冷速对液态合金Ca50Zn50快速凝固过程中微观结构演变的影响

采用分子动力学方法对六种不同冷速对原子尺寸相差较大的液态合金Ca50Znso凝固过程中微观结构演变的影响进行了模拟研究,并采用双体分布函数、Honeycutt—Andersen（HA）键型指数法、原子团类型指数法（CTIM一2）、可视化等方法进行了深入分析,结果表明：系统存在一个临界冷速,介于和5×10^11K／s与1×10^11K／s之间,在临界冷速以上（如1×10^11K／s,1×10^11K／s,1×10^11K／s和5×10^11K／s）时,系统形成以1551,1541,1431键型或二十面体基本原子团（12012000）等为主体的非晶态结构;在临界冷速以下时,系统形成以1441和1661键型或bcc基本原子团（1460800）为主体（含有少量的hcp（1200066）和fcc（12000120）基本原子团）的部分晶态结构．在非晶形成的冷速范围内,其总双体分布函数的第一峰明显分裂成与近邻分别为Zn—Zn,Ca—Zn,Ca—Ca相对应的三个次峰;且随着冷速的下降,同类原子近邻的次峰峰值升高、异类原子近邻的次峰峰值下降;Zn原子容易偏聚,随着冷速降低,二十面体的数量增多,非晶态结构也越稳定．在晶态形成的冷速范围内,Zn原子己大量偏聚形成大块bcc晶态结构,Ca原子也部分形成hcp和fcc晶态结构．     

碳纳米管阵列水渗透性质的研究

碳纳米管阵列组成的碳纳米管分子膜在生物学分子器件等方面有重要应用. 本文利用分子动力学方法计算研究水分子对(11, 11)碳纳米管阵列的渗透过程. 结果发现, 只有当阵列间隙面积大于57.91 Ų时, 水分子才能进入阵列间隙中, 并揭示了碳管内部、阵列间隙内水分子结构随相邻碳管间距变化的演化趋势以及管内外水分子电偶极矩的分布特性.