一种构建高分子与超临界二氧化碳混合体系全原子模型的方法技术

本发明专利技术涉及一种构建高分子与超临界二氧化碳溶液混合体系全原子模型的方法，属于计算机技术领域。所述方法包括如下步骤：1.构建超临界二氧化碳溶液体系的全原子模型；2.构建高分子的初始空间坐标；3.构建高分子的全原子模型对应的拓扑参数；4.对高分子进行能量优化；5.将高分子加入超临界溶液体系；6.对体系进行高温退火模拟以得到稳定的高分子与超临界二氧化碳混合体系的全原子模型。通过该方法，可以构建各种不同种类的高分子与超临界二氧化碳溶液组成的混合体系所对应的全原子模型，为进一步研究不同高分子对超临界二氧化碳溶液微观动力学的影响规律以及超临界二氧化碳的应用开发奠定有效的技术基础。

【技术实现步骤摘要】
一种构建高分子与超临界二氧化碳混合体系全原子模型的方法
本专利技术涉及计算机
，具体涉及一种构建高分子与超临界二氧化碳溶液混合体系全原子模型的方法。
技术介绍
超临界二氧化碳是指二氧化碳在温度和压力大于临界点时的独特状态。它具有接近于空气的低粘度和高扩散系数，同时具有接近于液体的高密度。它制备方法简单且易于循环利用，属于环境友好型绿色溶剂。它不仅可以作为绿色溶剂参与多种化学反应，还可以用于蛋白质提取和其他生物分子的萃取。更重要的是，超临界二氧化碳可以用于连续油管钻井，同时还可以提高原油采收率，对储集层没有任何污染。这对于油气藏的开发有着重要的意义。然而由于二氧化碳分子直线式的分子构型与较高的分子极性，导致其超临界状态下的粘滞系数可调区间非常窄，导致其无法满足生物分子萃取和石油开采等领域对高粘滞溶液的需求。在对超临界二氧化碳粘性的研究过程中，大量的实验研究表明由带有苯环非极性侧链和具有氟烃极性侧链构成的特定比例的高分子能够有效地增加超临界二氧化碳的粘滞系数。因此，研究人员希望从分子角度设计不同种类的高分子来调控超临界二氧化碳的粘性，从而拓展其在不同领域的应用。然而由于高分子单体种类繁多，且高分子中单体的配比难以控制，导致实验上无法实现对高分子与超临界二氧化碳之间相互作用规律定量研究。随着计算机技术以及分子建模方法的成熟，人们开始尝试利用精确的全原子分子动力学方法来研究各种体系中分子之间的微观动力学以及相互作用规律。该方法具有环境参数精确可调和分辨率高等优点。通过全原子分子动力学模拟，我们能够在原子尺度观察二氧化碳分子在高分子不同官能团附近的微观动力学过程，同时能够定量地研究它们与二氧化碳分子的微观相互作用机制，最终阐明高分子调控超临界二氧化碳粘性的定量规律。因此，我们提供了基于全原子分子动力学构建高分子与超临界二氧化碳混合体系模型的方法，为实现计算机辅助设计不同粘性的高分子与超临界二氧化碳的混合体系奠定基础。
技术实现思路
本专利技术要解决的问题是提供一种构建高分子与超临界二氧化碳溶液混合体系全原子模型的方法。可以利用成功构建的全原子模型进一步研究高分子与超临界二氧化碳液体在原子尺度的动力学和相互作用问题，如高分子调控超临界二氧化碳粘性规律等等。为实现上述目的，一种基于全原子分子动力学构建高分子与超临界二氧化碳混合体系模型的方法，包含如下步骤：(1)利用Gromacs软件，对二氧化碳的EPM2全原子在NPT系综下进行高温退火模拟，构建稳定的超临界二氧化碳溶液体系，其中二氧化碳分子约为11000个；(2)利用GaussView构建由苯乙烯和含氟丙烯酸酯构成的特定比例的嵌段高分子，并保存为pdb文件格式；(3)基于成熟的OPLSAA全原子力场以及高分子原子的化学成键类型，确定高分子中所有原子对应的原子类型,同时确定高分子包含的单体内原子之间的拓扑关系，并产生两个单体对应的top文件格式；(4)基于链接单体之间原子成键作用的强弱，确定高分子单体之间的连接拓扑关系，然后根据高分子包含的单体之间的排列方式，将单体对应的拓扑和单体之间的拓扑关系按顺序关系进行组合，生成该高分子对应完整拓扑关系，并产生相应的top文件格式；(5)利用Gromacs以及选定的全原子力场对初始的pdb文件中的高分子构型进行能量最小化,从而得到高分子的能量最优构型(6)将能量优化后的高分子加入超临界二氧化碳溶液体系，并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟以消除二氧化碳和高分子的不合理接触，最终得到稳定的高分子与超临界二氧化碳混合体系的全原子模型；其中，步骤(1)和(6)进行NPT系综下的分子动力学模拟时，动力学参数设置为：采用Leap-Frog算法进行积分求解，对应的积分步长为0.5fs；采用ParticleMeshEwald算法对长程静电作用进行求解，其截断半径为0.9nm；Lenard-Jones相互作用的截断半径为1.3nm；采用Berendsen热浴算法来恒定模拟的温度，其耦合常数为0.5ps；采用Berendsen压浴算法来恒定模拟的压强，其耦合常数为1.0ps。其中，步骤(1)和(6)进行NPT系综下的高温退火模拟时，具体退火过程依次为(1000bar，800K)-(800bar，600K)-(600bar，400K)-(400bar，320K)-(280bar，320K)。其中，步骤(2)中所述拟定的高分子极性单体与非极性单体之间以嵌段模式进行组合，且极性单体与非极性单体之间的比例为29：1～15：15。其中，步骤(3)中高分子的极性单体为含氟丙烯酸，而非极性单体为苯乙烯，两个单体对应的全原子模型为OPLSAA力场。其中，步骤(4)中，高分子单体之间的连接拓扑包括所有与主链原子之间拓扑距离小于等于3的非氢原子之间的成键关系构成，包括最近邻原子间的成键关系，次近邻原子间的键角关系，以及次次近邻原子间的二面角关系。其中，步骤(6)中将高分子加入到超临界二氧化碳时，高分子表面的范德瓦尔半径设置为0.35nm，且体系的高分子与二氧化碳的混合体系对应的高分子占体系的质量分数设定为0.19wt％左右。其中，步骤(6)中采用了NPT系综下对混合体系实施高温退火模拟，消除了高分子和二氧化碳之间的不合理接触的过程，实现了对初始构型进行弛豫。本专利技术的技术效果在于：利用本方法能够更简洁且高效地构建的各种高分子与超临界二氧化碳混合体系的全原子模型；该方法构建的模型具有更加精确和更稳定的特点；该方法为我们从原子尺度研究高分子与二氧化碳分子在极端条件下的相互作用以及利用高分子调节超临界二氧化碳的粘滞性质奠定技术基础。附图说明下面结合附图对本专利技术做进一步地说明：图1为实施例1中超临界二氧化碳体系密度随时间的变化图，其中体系压强为280bar，温度为320K；图2为实施例1和实施例2中极性单体含氟丙烯酸酯与非极性单体苯乙烯的原子结构以及原子标注示意图；图3为实施例1中极性单体和非极性单体比例为29：1的高分子对应的原子构型示意图；图4为实施例1中高分子与超临界二氧化碳混合体系在能量优化之后的原子构型示意图；图5为实施例2中极性单体和非极性单体比例为15：15的高分子对应的原子构型示意图；图6为实施例2中高分子与超临界二氧化碳混合体系在能量优化之后的原子构型示意图；图7为实施例1和实施例2中超临界二氧化碳分子的微观均方位移(MSD)随时间变化规律图；具体实施方式下面结合实施例和附图对本专利技术进一步加以说明。实施例1本实施例详细阐述了本专利技术“一种构建高分子与超临界二氧化碳混合体系全原子模型的方法”在具体实施时针对极性单体和非极性单体比例为29：1的高分子与超临界二氧化碳混合体系对全原子模型的构建，具体步骤如下：(1)采用EMP2全原子模型对二氧化碳进行全原子建模，利用Gromacs软件对二氧化碳体系在NPT系综下进行高温退火模拟，其退火模拟过程依次从(1000bar，800K，)至(800bar，600K，)，至(600bar，400K，)至(400bar，320K，)最后至(280bar，320K)。得到稳定的超临界二氧化碳体系，其密度随时间变化曲线如附图1所示，平均值为830.2Kg/m3。该值与实验密度本文档来自技高网...

【技术保护点】
1.一种构建高分子与超临界二氧化碳混合体系全原子模型的方法，其特征在于：所述方法具体步骤为：(1)利用Gromacs软件，对二氧化碳的EPM2全原子在NPT系综下进行高温退火模拟，构建稳定的超临界二氧化碳溶液体系；(2)利用Gauss View构建由苯乙烯(Styrene)和含氟丙烯酸酯(Fluoroacrylate)构成的特定比例的嵌段高分子；(3)基于成熟的OPLSAA全原子力场以及高分子原子的化学成键类型，确定高分子中所有原子对应的原子类型,同时确定高分子包含的单体内原子之间的拓扑关系；(4)基于链接单体之间原子成键作用的强弱，确定高分子单体之间的连接拓扑关系，然后根据高分子包含的单体之间的排列方式，将单体对应的拓扑和单体之间的拓扑关系按顺序关系进行组合，生成该高分子对应完整拓扑关系；(5)利用Gromacs以及选定的全原子力场对高分子构型进行能量最小化,从而得到高分子的能量最优构型；(6)将能量优化后的高分子加入超临界二氧化碳溶液体系，并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟，最终得到稳定的高分子与超临界二氧化碳混合体系的全原子模型。...

【技术特征摘要】
1.一种构建高分子与超临界二氧化碳混合体系全原子模型的方法，其特征在于：所述方法具体步骤为：(1)利用Gromacs软件，对二氧化碳的EPM2全原子在NPT系综下进行高温退火模拟，构建稳定的超临界二氧化碳溶液体系；(2)利用GaussView构建由苯乙烯(Styrene)和含氟丙烯酸酯(Fluoroacrylate)构成的特定比例的嵌段高分子；(3)基于成熟的OPLSAA全原子力场以及高分子原子的化学成键类型，确定高分子中所有原子对应的原子类型,同时确定高分子包含的单体内原子之间的拓扑关系；(4)基于链接单体之间原子成键作用的强弱，确定高分子单体之间的连接拓扑关系，然后根据高分子包含的单体之间的排列方式，将单体对应的拓扑和单体之间的拓扑关系按顺序关系进行组合，生成该高分子对应完整拓扑关系；(5)利用Gromacs以及选定的全原子力场对高分子构型进行能量最小化,从而得到高分子的能量最优构型；(6)将能量优化后的高分子加入超临界二氧化碳溶液体系，并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟，最终得到稳定的高分子与超临界二氧化碳混合体系的全原子模型。2.根据权利要求1所述的方法，其特征在于：步骤(1)和(6)中，采用Leap-Frog算法进行积分求解，对应的积分步长为0.5fs；采用ParticleMeshEwald算法对长程静电作用进行求解，其截断半径为0.9nm；Lenard-Jones相互作用的截断半径为1.3nm；采用Berendsen热浴...

【专利技术属性】
技术研发人员：闫永凤，张军，赵明伟，戴彩丽，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东,37