本论文针对目前传统试验手段用在研究分子扩散过程以及扩散系数方面的局限性,通过分子动力学模拟(Molecular Dynamic Simulation)的手段,深入详细地研究了小分子在反渗透复合膜以及水溶液内的扩散过程。主要进行了四方面的研究工作:分子在体系内扩散系数的分子模拟方法研究;甘氨酸等12种氨基酸分子在水溶液内扩散系数的模拟;氧气等9种气体分子在水溶液内扩散系数的模拟;水和盐分子(NaCl, MgCl2,CaSO4,K2SO4)在8种反渗透复合膜中扩散状态及扩散系数的研究,膜的结构单体与水和盐分子在膜内扩散系数关系的分析与讨论。
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.10×10-9m2s-1,1.80×10-9m2s-1,2.00×10-9m2s-1,聚类分析方法的误差要小于微分-区限变分法,充分证明聚类方法是一个有效的分析方法。
③模拟甲烷和氨气在水中的扩散过程时,对0.1fs的步长值延长了2-3倍的模拟时间,结果表明模拟时间的延长对最后的模拟结果没有明显的影响,从而排除了模拟时间对扩散系数模拟结果的影响。另外,论文对甲烷,氨气和二氧化碳气体在水中的扩散过程采用了0.1fs-1.5fs范围内几个不同的步长值,模拟结果表明,步长值对扩散系数的模拟结果有较大影响。随着步长的增加,扩散系数有一个极值,与该极值相邻的两个步长对应的扩散系数与实验值最为接近。这两个步长值为极值对应的步长值的70%和150%。结合步长在模拟计算中的作用,本论文认为对于一个特定体系,应在0.1fs到1.5fs之间取多个值进行模拟,取模拟值最大的步长的70%即为模拟的合理值。
2)本论文首次通过分子模拟的手段研究了氨基酸在水溶液内的扩散过程,并利用微分-区限变分法计算了氨基酸在水中的扩散系数。目前氨基酸在水中扩散系数主要通过实验手段来获得,对于氨基酸在水中的扩散系数的模拟,目前尚未见报道。本文通过分子动力学模拟的手段模拟计算了甘氨酸、丙氨酸等12种氨基酸在水中的扩散系数。模拟值与实验值的误差小于7%,证明了分子模拟方法应用在扩散系数计算方面的有效性。
3)本论文通过分子模拟的手段模拟了氧气等9种气体在水溶液内扩散过程,并利用聚类分析方法计算了气体在水中的扩散系数。对于O2,CO2,CH4,NH3四种气体在水中的扩散系数的模拟,已有文献的模拟值与实验值的误差最大超过50%,而本论文的模拟数据与实验值的误差小于30%。充分证明了本文所采用的分子模拟途径的可靠性。
4)本文以8种反渗透复合膜为主要研究对象,模拟了海水中主要成分(H2O,NaCl, MgCl2,CaSO4和K2SO4)在膜内的扩散过程并计算了各粒子在膜内的扩散系数。
①模拟了TMC/MPD膜在体系中的扩散系数,模拟结果为0.0115A2/ps,与已有的对于该膜的模拟结果0.015 A2/ps基本吻合。证明了本文模拟过程中对于模型建立,模拟方法以及数据分析方法的有效性。
②水分子在TMC/MPD膜内的扩散状态表明水分子与高分子链上的功能基团有明显的分子间相互作用。水分子在靠近苯环的区域内与苯环有排斥运动;在靠近酰胺基的区域内有吸引运动。
③在所有膜中,水分子在膜HT/MMPD中的扩散系数最大,远远高于在其他膜内的扩散系数。对比水分子在膜TMC/MMPD和TMC-HT/MMPD的扩散系数发现水分子在加入了HT的膜中扩散系数要高一些了;对比水分子在膜TMC/MPD和膜TMC/MMPD内的扩散系数,水分子在膜TMC/MMPD的扩散系数要小很多。由于MMPD单体在苯环上含有一个甲基,在同样数量和比例生成聚合物膜的情况下,MMPD比MPD的刚性要强一些,从而不利于水分子的通通过。对比水分子在膜CFIC/MMPD和膜CFIC/MPD的扩散系数也可得出同样的结论。
④NaCl分子在膜HT/MMPD中的扩散系数最大,而在膜ICIC-HT/MMPD中的扩散系数则最小,其次分别是膜ICIC-IPC/MPD ,膜TMC-HT/MMPD,膜CFIC/MPD以及膜TMC/MMPD。模拟结果中只有ICIC-IPC/MPD与实验结果有较大差距。
⑤在所模拟的8种膜内,Na+和Cl-在膜内的扩散系数值相差较大,NaCl分子的运动受限于运动较慢的离子。膜的种类不同,制约NaCl在膜内扩散的离子便不同:对于膜TMC/MPD以及膜CFIC/MMPD,NaCl分子的扩散受Cl-的制约;而其他6种类型的膜则受Na+的制约。
⑥对于膜ICIC-HT/MPD,水分子在MgCl2,CaSO4和K2SO4体系中的扩散系数与在NaCl体系中的结果基本一致,该结果表明在同种膜中,水分子的扩散过程不受体系中盐分子类型的影响。
⑦在ICIC-HT/MPD反渗透膜中,阴离子(Cl-,SO42-)的扩散系数均远大于阳离子(Mg2+,Ca2+,K+)的扩散系数。这与NaCl水体系中,Cl-的扩散系数远大于Na+的扩散系数的运动状况一致,该结果说明当盐分子单独存在时,其在膜内的扩散过程受哪个离子制约只与膜种类有关,与盐分子本身无关。
