水分子通过形变碳纳米管的运输行为分析论文.pdfVIP

水分子通过形变碳纳米管的运输行为分析论文 Hummer 等人[5]在 2022 年采纳分子动力学模拟的方法，发觉 水分子不仅可以在单壁碳纳米管中运输，而且会加快其运输速度。一 部分讨论人员发觉假如把碳纳米管口作部分扭曲后，能很大地加快水 分子通过碳纳米管的速率。 2022 年李敬源等人[6]将电荷放置于碳 纳米管之外，讨论了外加电荷对管内水分子的影响，观看到了电荷对 于水分子运输的开关方面特别灵敏。 此后，大量的学者[79]在这些 重大成果启发之下，对水在碳纳米管中动力学性质作了很多讨论。如 周毅[7]在碳纳米管中水的行为的分子动力学模拟中讲解并描述：在 纳米尺度下的受限空间中，水分子所表现的行为与宏观中截然不同。 何俊霞等人[8]通过转变形变位置探究纳米通道中的窄结构对水输运 产生的开关特性。 由于碳纳米管生长和合成试验中[10],总有一些拉伸或压缩的 形变碳纳米管存在，但形变碳纳米管内水分子的运输特性未见到相关 报道。 故本文运用分子动力学的方法讨论碳纳米管拉伸或压缩时其 形变程度对其内部水的运输行为影响。 即选择手性矢量 〔7,7〕和 〔8,8〕2 种碳纳米管，拉伸或压缩量 k 分别为：0.5,0.6,0.8,1. 0,1. 2,1. 5,1. 8 倍的碳纳米管进行模拟计算，进而分析形变碳纳 米管中水分子的流量、偶极矩概率分布、轴向分布以及径向分布。 1 讨论模型与方法 本文采纳分子动力学软件 NAMD 和 CHARMM力场文件。 讨论 模型如图1 所示，它主要由 1 个形变的单壁碳纳米管和 2 个石墨 第 1 页 烯层构成，并且 2 个石墨烯成处于碳纳米管的 2 个端口处，在石墨 烯层外同时添加 2 个 TIP3P 水体，其大小为 1. 5 nm ×1. 5nm × 1. 5 nm,其中形变碳纳米管由 NanotubeModeler 软件生成，其原长 为 6. 05 nm. 模拟时，设定温度为 300 K,并采纳 Berenden方法进行控温， 截断半径为 1. 2 nm,利用 Particle Mesh Ewald 的方法计算系统 的静电互相作用，网格划分精度设置为 0. 1 nm,当超过 1 nm 时采 纳过渡函数修正计算，采纳周期性边界条件其大小为1. 5 nm × 1. 5 nm ×9. 2 nm,模拟步长 2 fs,模拟时长均为 10 ns. 为解决错误 结构的影响本文执行10 000步的能量最小化处理，同时为了削减运 算量和简化模拟体系，将碳纳米管和石墨烯层进行固定。 为使水分 子通过碳纳米管，本文沿 z 轴正方向施加 0. 5 v/nm 的电场。 全 部的计算任务均是在高性能服务器上完成的。 2 结果与分析 2. 1 水分子通过碳纳米管的流量 水分子通过形变碳纳米管的流量如图 2 所示，可以发觉：当 碳纳米管未发生形变时，即 k 1 时，管径大的流量大，该结论与 王俊等人[11]全都； 当碳纳米管发生较小的形变时，即 k 0.8 或 1. 2 时，管径大的水分子流量反而减小，但均大于未发生形变时的 流量； 当形变程度 k≥1. 5 时，2 种管径的水分子流量几乎为零， 这说明 〔7,7〕 和 〔8,8〕 型碳纳米管被拉伸到 1. 5 倍以上时， 几乎没有水分子能通过碳纳米管。 分析认为碳纳米管的轴向形变程 第 2 页 度对管内水分子的运输有重要的影响作用，同时可以通过转变碳纳米 管的形变程度使其起到分子开关的作用。 2. 2 水分子通过碳纳米管的径向分布 碳纳米管中水分子的径向分布函数如图 3 所示，其中图 3 〔a〕 、〔b〕 分别表示 〔7,7〕 和 〔8,8〕 型碳纳米管中水分 子的径向分布图。 由图 3 〔a〕可以发觉不同形变程度下，碳纳米 管内水分子的径向分布函数均只消失 1 个峰值并位于 r 0. 16 nm 处，说明水分子在该处消失的概率最大，这说明水分子在管内是以水 柱面的形式通过碳纳米管； 同时还可以发觉，当 k 越小时其峰越尖 锐，这说明水分子仅仅从柱面处通过，当 k 增加时其水柱内也有水 分子通过。 由图 3 〔b〕 同样可以发觉径向分布函数仅存在 1 个 位于 r 0. 23 nm 处的峰，同样也说明水分子是以水柱面的形式通 过碳纳米管。