反渗透膜及其水溶液内扩散过程的分子模拟的研究.pdf

反渗透膜及水溶液内扩散过程的分子模拟研究 摘要 本论文针对目前传统试验手段用在研究分子扩散过程以及扩散系数方面的 Simulation)的手段，深入 局限性，通过分子动力学模拟(MolecularDynamic 详细地研究了小分子在反渗透复合膜以及水溶液内的扩散过程。主要进行了四方 面的研究工作：分子在体系内扩散系数的分子模拟方法研究；甘氨酸等12种氨 基酸分子在水溶液内扩散系数的模拟；氧气等9种气体分子在水溶液内扩散系数 散状态及扩散系数的研究，膜的结构单体与水和盐分子在膜内扩散系数关系的分 析与讨论。 1)论文首次从模型的分子数样本以及模拟步长出发，对影响分子在体系内 扩散系数的主要因素进行了研究。并针对动力学模拟数据提出了新的分 析方法。 ①对浓度分别为0．7mol／L和0．9mol／L的甘氨酸水体系进行了动力学模 拟。建立模型时分别加入了不同数目的分子：对0．7mol／L的甘氨酸模 型的构建，分别在立方晶格中加入了3、5、7、10、12个甘氨酸分子； 对于0．9mol／L的甘氨酸则分别在立方晶格中加入了4、7、9、13、15 个甘氨酸分子。扩散系数模拟值与实验值的对比表明分子数样本对于扩 散系数的模拟结果有较大的影响。分子数越多，则模拟结果与实验值越 接近，当分子数达到一定值时，扩散系数趋于一个稳定的值。 论文还针对不同温度下氧气在水中的扩散过程作了关于样本数的模拟。 分别考察了温度为303K和313K条件下的气体扩散。模拟体系中氧气分 子数分别采用1，2，3个。结果表明随着模型中分子数的增多，得到的 扩散系数的值也逐渐增大，并逐渐恒定于正确值。 ②通过对比传统分析方法以及微分一区限变分法用在均方位移(MSD)曲线 上的优缺点，本论文提出了新的数据分析方法一聚类分析方法。该方法 克服了传统方法用在MSD线性不好时给分析结果带来误差的缺点，对微 分一区限变分法进行了改进。将聚类方法和微分一区限变分法应用在氧气 在水中的扩散系数的分析中，实验值，微分一区限变分法值和聚类方法 值分别为2．fOX 101酽s一，1．80×lO’n2s一，2．OOX10％2s～，聚类分析方法 的误差要小于微分一区限变分法，充分证明聚类方法是一个有效的分析 方法。 ③模拟甲烷和氨气在水中的扩散过程时。对0．1fs的步长值延长了2—3 倍的模拟时间，结果表明模拟时间的延长对最后的模拟结果没有明显的 影响，从而排除了模拟时间对扩散系数模拟结果的影响。另外，论文对 甲烷，氨气和二氧化碳气体在水中的扩散过程采用了0．1fs-I．5fs范围 内几个不同的步长值，模拟结果表明。步长值对扩散系数的模拟结果有 较大影响。随着步长的增加，扩散系数有一个极值，与该极值相邻的两 个步长对应的扩散系数与实验值最为接近。这两个步长值为极值对应的 步长值的70％和150％。结合步长在模拟计算中的作用，本论文认为对于 一个特定体系，应在0．Ifs到1．5fs之间取多个值进行模拟，取模拟值 最大的步长的70％即为模拟的合理值。 2本论文首次通过分子模拟的手段研究了氨基酸在水溶液内的扩散过程。 并利用微分一区限变分法计算了氨基酸在水中的扩散系数。目前氨基酸 在水中扩散系数主要通过实验手段来获得，对于氨基酸在水中的扩散系 数的模拟，目前尚未见报道。本文通过分子动力学模拟的手段模拟计算 了甘氨酸、丙氨酸等12种氨基酸在水中的扩散系数。模拟值与实验值 的误差小于7％，证明了分子模拟方法应用在扩散系数计算方面的有效 性。 3)本论文通过分子模拟的手段模拟了氧气等9种气体在水溶液内扩散过 程，并利用聚类分析方法计算了气体在水中的扩散系数。对于02，c02， CH,。NH。四种气体在水中的扩散系数的模拟，已有文献的模拟值与实验 值的误差最大超过50％，而本论文的模拟数据与实验值的误差小于3096。 充分证明了本文所采用的分子