石墨烯掺入对PVDF膜模型纳米结构与力学响应的影响：分子动力学研究

《Membranes》：Effects of Graphene Incorporation on the Nanoscale Structure and Mechanical Response of PVDF Membrane Models: A Molecular Dynamics Study

　　

聚偏二氟乙烯（PVDF）膜需要兼具合适的形貌与力学完整性，但纳米尺度几何结构的变化并不必然改善所有载荷响应。研究人员采用原子尺度分子动力学模拟，比较了由五重复单元低聚物构建的原始PVDF膜模型与石墨烯/PVDF膜模型。采用力驱动沉降程序考察侧向元胞尺寸与聚合物链数对沉积形貌的影响。结构指标为二维投影孔隙面积分数。增加聚合物链数降低了投影孔隙分数，而增加侧向元胞尺寸则产生更薄的沉积层及更大的未覆盖面积。在相应构建条件下，含石墨烯结构表现出更大的投影覆盖率。另行制备的力学试样在力学比较中采用相同的PVDF参数进行拉伸与剪切加载。所报道的平均表观拉伸割线模量在复合材料中约高14.6%，属描述性有限速率比较。其三次独立运行的平均峰值应力为122.46 ± 32.57 MPa（均值 ± 样本标准差），而其平均曲线的峰值为109.38 MPa。这些结果将所报道的初始斜率差异与峰值响应变异性区分开来，并阐明了峰值承载阻力的统计意义。它们为讨论膜骨架几何与载荷依赖响应提供了原子尺度基础，但不构成对水渗透性、选择性或抗污染性能的预测。

研究背景方面，膜分离技术在水净化、废水回用及复杂水流深度处理中具有重要地位，膜材料除分离性能外还需在操作过程中保持足够的结构完整性。污染与反复清洗带来额外操作需求，而制备路线决定了膜性能所依赖的孔形貌与聚合物结构，因此膜构建、纳米尺度结构与力学响应之间的关系对耐久聚合物膜设计具有重要意义。聚偏二氟乙烯（PVDF）因耐化学性、热稳定性及成膜能力广泛用于微滤与超滤，其性能取决于聚合物化学与成膜过程中形成的形貌。较大的孔隙分数可能为传输提供更多空间，但膜厚度、孔隙几何与抗变形能力应分别表征而非由单一结构描述符推断。石墨烯基添加剂已被研究作为PVDF膜改性剂，但原始石墨烯与氧化石墨烯（GO）必须区分：GO含氧官能团改变其表面化学，而本研究关注嵌入PVDF的纯碳石墨烯片。已有GO/PVDF实验研究考察了共混与填料浓度、化学功能化添加剂、膜生物反应器运行及天然有机物去除、混凝与静电纺丝制备等，这些研究展示了含石墨烯膜的制备依赖性变化，但其传输与抗污染结果不应直接转移至未功能化石墨烯/PVDF模型。添加剂即使不假设润湿性改善也可影响聚合物堆叠与载荷传递，因此表面化学与过滤行为的变化应与承载聚合物组装体的变化区分。对本研究未氧化石墨烯/PVDF体系，直接的力学问题是初始拉伸斜率增加是否伴随峰值阻力增加。石墨烯复合材料力学取决于填料排列与聚合物–填料相互作用，拉曼测量表明应力可传递至石墨烯且界面具有有限承载能力，含单壁碳纳米管的PVDF模拟也表明刚度变化依赖加载方向与纳米管间距，这些研究促使将初始刚化与组装体最大拉伸阻力分开考察。分子动力学（MD）可获取补充膜尺度测量的局部构型与相互作用，原子模型允许明确指定构建历史、分子排列与加载协议。本研究问题为：侧向限域与聚合物载量如何调控沉积几何，以及添加石墨烯是否在初始刚度与峰值拉伸阻力上产生相同变化。结构扫描采用共同力驱动构建程序比较膜形貌、投影覆盖率与几何孔隙尺寸描述符；另行制备的力学试样在共同力场框架下比较拉伸与剪切响应。其贡献在于模型特定的构建几何与不同力学可观测量的比较，而非重建实验相转化或完整结构–力学–传输关系。研究方法上，主要关键技术方法包括：采用LAMMPS进行分子动力学模拟，力学模型使用固定拓扑OPLS-AA力场框架，长程静电采用PPPM处理；PVDF链以五重复单元低聚物表示（C10H11F11），结构复合体系含15个纯碳C54片，力学试样含150条PVDF链（4800个聚合物原子），复合材料含810个石墨烯原子、共5610个原子；结构模型通过力驱动沉降程序构建，采用二维投影孔隙面积分数与周期性Voronoi几何孔隙尺寸描述符表征；拉伸加载以1×10?? fs?1工程速率沿x方向伸长，剪切加载以0.002 ? fs?1驱动顶部区域；使用OVITO可视化。研究结果方面，3.1 PVDF与石墨烯/PVDF膜模型的构建：膜形貌受侧向几何比Λ与聚合物链数n调控，固定n=150时增大L使沉积层变薄，固定L=6.8 nm时增加n产生更厚结构；φ2D量化显示增大侧向元胞尺寸降低聚合物覆盖、增加未覆盖投影区域，增加链数提高投影原子覆盖并降低φ2D；石墨烯/PVDF模型在相应构建条件下投影孔隙面积分数更低，这反映碳片直接投影及沉积过程中周围PVDF组织变化；sV分布显示增大Λ使分布向更大几何尺度移动并变宽，增加n抑制大尺寸尾部。3.2 拉伸与剪切变形：拉伸应力–应变曲线显示应力增加至最大值后下降，区分初始加载区与峰后软化；所报道平均运行级表观拉伸割线模量为原始PVDF 2.9270 GPa、石墨烯/PVDF 3.3530 GPa，描述性差异约14.6%；石墨烯/PVDF个体峰值真应力为111.4527、96.8164、159.1100 MPa，均值与样本标准差为122.4597 ± 32.5728 MPa（n=3），平均曲线峰值为109.3834 MPa；剪切采用顶部刚性驱动、底部弹性锚定，非线性应力–应变曲线表征载荷依赖变形而非单一线性剪切模量；拉伸加载起始时19.78 ± 1.81%的PVDF原子接触石墨烯，150条链中116.3 ± 7.4条含至少一个接触原子，105个可能片对中11.0 ± 3.6对最小原子间距≤4.0 ?。3.3 与已发表研究及模型范围的关系：Moradi等研究静电纺PVDF/石墨烯纤维膜，Raji等比较熔融加工PVDF中石墨烯纳米片与GO，Wu等研究纳米石墨烯/PVDF复合膜力学浓度依赖性，Lee等考察拉伸PVDF/GO纤维，这些基准因孔隙率、拉伸、填料化学与浓度差异而不可互换；Yang等400单体每链PCFF模型在300 K给出密度1.59 g cm?3与杨氏模量5.6 ± 1.4 GPa；模型未预设PVDF晶格，此构建选择不同于结晶度定量测定；原子加载时间尺度相对实验室测试可限制链重排。讨论部分总结：结果将初始刚度与峰值响应分别评价，不赋予统计显著性、统计等价或整体增强效应；接触统计量化加载起始时聚合物–石墨烯邻近性但未确立唯一载荷传递机制；结构沉积与力学模型家族在电荷分配、边界条件与制备协议上不同，结构投影覆盖率未指定为个别力学响应之原因。研究结论翻译：采用原子模拟比较了原始PVDF与石墨烯/PVDF短链膜模型的沉积几何与力学响应。增加聚合物链数产生更厚沉积层与更低投影孔隙面积分数，而在固定聚合物量下增大侧向元胞尺寸产生更薄沉积层与更大未覆盖面积。含石墨烯结构在相应构建条件下表现出更大投影覆盖率。在另行制备的DP=5力学模型中，所报道平均运行级1–2%拉伸割线模量在复合材料中约高14.6%，此描述性比较仍受曲线分辨率与采样限制。其三次运行平均峰值应力为122.46 ± 32.57 MPa（均值 ± 样本标准差），平均曲线峰值为109.38 MPa。因此初始刚度与峰值响应需要分别的估计量与不确定性陈述。接触统计量化加载起始时聚合物–石墨烯邻近性但未确立唯一载荷传递机制。相同PVDF参数支持力学家族内比较。两家族均使用C54片，而结构沉积在电荷分配、边界条件与制备协议上不同。这些结果未确立链长或速率无关性、因果性试样级结构–力学关系，或定量水传输、选择性与抗污染性能。