嵌段共聚物微相分离的计算机模拟.pptVIP

元胞动力学方法（CDS） 嵌段共聚物的微相分离过程可用不含随机热力学噪声的Ginzburg-Landau方程描述： 为系统的局部序参量。 为系统粗粒化的自由能泛函。 元胞动力学方法（CDS） 对称两嵌段共聚高分子熔体薄膜微区形态随时间的演化 两嵌段共聚物和均聚物的共混物AB/C 薄膜所形成的微相结构 AB链中A嵌段的体积分数fA =0.3， ?C 从0.3 ? 0.7变化 Feng J., Liu H. L., Hu Y., Macrom. Theory Simul., 11(5), 2002: 549-565 耗散粒子动力学方法（DPD） 每个粗粒化链节代表了很多个分子。 耗散粒子动力学方法（DPD） 耗散粒子动力学方法（DPD）由Hoogerbrugge和Koelman首次提出，后经他人多次完善，现已成功应用于高分子熔体和表面活性剂溶液等具有介观结构的系统的模拟。 在DPD方法中，假设所有粒子在保守力FC、耗散力FD与随机力FR的共同作用下做牛顿运动： 耗散粒子动力学方法（DPD） 通过数值求解牛顿运动方程，可获得系统的微观构象和分子聚集状态的变化，从而得到系统的各项宏观性质和微相结构的演变。 牛顿运动方程的数值求解格式有很多种，例如velocity-Verlet算法、蛙跳算法、预报-校正算法等。 一种修正的velocity-Verlet算法由以下几步构成： 耗散粒子动力学方法（DPD） 受限于纳米圆孔的A5B5嵌段共聚高分子微相形态 压力增加时，受限于纳米圆孔的A5B5微相形态的转化 Feng J., Liu H. L., Hu Y., Macrom. Theory Simul., 15(9), 2006: 674-685 耗散粒子动力学方法（DPD） (a) (b) (c) (A5)4(B5)4星形嵌段共聚高分子的几种介观结构 (a): ?N=20; (b): ?N=30; (c): ?N=40 Xu Y. Feng J., Liu H. L., Hu Y., Mol. Simul., 32(5), 2006: 375-383 Block copolymer, Neutral wall, f=0.5, 0.4,0.3,0.2,0.1 动态密度泛函方法（DDFT） 动态密度泛函理论（DDFT）又称为介观动力学（Meso-Dyn），是荷兰学者 Fraaije 于上世纪90年代提出的一种建立在粗粒化高斯链模型基础上的介观模拟方法。 动态密度泛函方法（EOS-DDFT） 为解决DDFT的计算效率低下以及不能反映压力对嵌段共聚物微观分相的影响这两个问题，基于本课题组多年来在流体的分子热力学方面的研究积累，我们提出了基于状态方程的动态密度泛函理论。在该理论中，我们首先通过引入键势来获得链的连接性对体系自由能的贡献，进而得到理想链的自由能泛函；然后提出了一个重要假设：体系局部密度分布满足流体状态方程。在此基础上，我们通过本课题组发展的方阱链流体状态方程来推导得到体系的非理想自由能泛函，从而避免了传统DDFT所涉及的烦琐的自洽场迭代过程，并首次将压力参数引入到自由能泛函的表达式中。 动态密度泛函方法（EOS-DDFT） 动力学方程可以用与时间相关的Ginzburg-Landau方程表示： Xu H., Liu H. L., Hu Y., Chem. Eng. Sci., 62, 2007: 3494-3501 动态密度泛函方法（EOS-DDFT） 两嵌段共聚物在不同对比压力下微相分离的相图 图中符号“■”为无序结构，“●”为球状结构，“▲”为柱状结构，“◆”为双连续状结构 Xu H., Liu H. L., Hu Y., Macrom. Theory Simu., 16(3), 2007: 262-268 谢 谢 ! * * （3）微相分离的平衡结构受系统松弛过程的影响很大，控制系统的演化条件可以得到形态和有序尺度不同的结构。 嵌段共聚物相分离的特点 例如：控制嵌段共聚物的膜厚可得到不同的微相结构 三嵌段共聚高分子的自组装结构S.I. Stupp et al.,Science,276, 384(1997) 嵌段共聚物相分离的特点 嵌段共聚物微相分离的应用 AB两嵌段共聚物和针状纳米颗粒复合物的自组装G.W. Peng et al., Science, 288, 1802(2000) 嵌段共聚物微相分离的应用 Monte Carlo 模拟（MC） Monte Carlo