J-OCTA的溶解度模块可根据超额化学势值计算溶解度系数。
目的和方法
溶解度系数是评估聚合物渗透率的一个重要因素。J-OCTA的溶解度模块可以计算出超额化学势 μex和分子溶解在聚合物中时的溶解度系数S。在计算超额化学势时,该模块中的采样效率使用排除体积图采样(EVMS)法[1]。
s: 无量纲溶解度 即亨利常数
kB: 玻尔兹曼常数
T: 聚合物体系的温度
T0: 标准条件温度
p0: 标准条件下的压力
在本案例中,我们计算了三种聚合物(图1)中的气体分子(N2, O2, CO2和CH4)和水的溶解度系数。聚合物有50个单体单元,使用GAFF作为力场。气体分子采用Dreiding力场,水分子采用SPC-FW模型。用NPT系综对聚合物进行了3ns的MD模拟弛豫,并计算溶解度系数。
图1 用于计算的聚合物
模拟结果
我们计算了聚合物中气体分子和水的溶解度系数,结果如图2所示。在聚合物中,PS的溶解度较高。在气体分子中,非极性N2和O2分子的溶解度较低,CO2分子溶解度较高,模拟结果与实验结果基本吻合。
图2 溶解度系数计算值与实验值的比较
PE、PS、PB和水的实验值分别取自文献[2-5]。
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删