CN109473144A - 一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 - Google Patents
一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 Download PDFInfo
- Publication number
- CN109473144A CN109473144A CN201811276321.9A CN201811276321A CN109473144A CN 109473144 A CN109473144 A CN 109473144A CN 201811276321 A CN201811276321 A CN 201811276321A CN 109473144 A CN109473144 A CN 109473144A
- Authority
- CN
- China
- Prior art keywords
- macromolecule
- carbon dioxide
- supercritical carbon
- monomer
- atom
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 title claims abstract description 141
- 229920002521 macromolecule Polymers 0.000 title claims abstract description 83
- 229910002092 carbon dioxide Inorganic materials 0.000 title claims abstract description 70
- 239000001569 carbon dioxide Substances 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000000137 annealing Methods 0.000 claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 239000000178 monomer Substances 0.000 claims description 71
- 125000004429 atom Chemical group 0.000 claims description 63
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Natural products C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 claims description 15
- 229920000642 polymer Polymers 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 7
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 6
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 5
- 229910052799 carbon Inorganic materials 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 239000004408 titanium dioxide Substances 0.000 claims description 3
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 claims description 2
- 230000009881 electrostatic interaction Effects 0.000 claims description 2
- 229910052731 fluorine Inorganic materials 0.000 claims description 2
- 239000011737 fluorine Substances 0.000 claims description 2
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 2
- 239000002245 particle Substances 0.000 claims description 2
- QQONPFPTGQHPMA-UHFFFAOYSA-N Propene Chemical compound CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 claims 2
- DHJMZSUXHHPWGK-UHFFFAOYSA-N 2-fluoroprop-2-enoic acid;prop-2-enoic acid Chemical class OC(=O)C=C.OC(=O)C(F)=C DHJMZSUXHHPWGK-UHFFFAOYSA-N 0.000 claims 1
- HIBWGGKDGCBPTA-UHFFFAOYSA-N C=CC1=CC=CC=C1.C=CC1=CC=CC=C1 Chemical compound C=CC1=CC=CC=C1.C=CC1=CC=CC=C1 HIBWGGKDGCBPTA-UHFFFAOYSA-N 0.000 claims 1
- 239000002253 acid Substances 0.000 claims 1
- 150000002148 esters Chemical class 0.000 claims 1
- 125000001153 fluoro group Chemical group F* 0.000 claims 1
- 125000003011 styrenyl group Chemical group [H]\C(*)=C(/[H])C1=C([H])C([H])=C([H])C([H])=C1[H] 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 2
- 229960004424 carbon dioxide Drugs 0.000 description 56
- 150000001252 acrylic acid derivatives Chemical class 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 102100033176 Epithelial membrane protein 2 Human genes 0.000 description 2
- 101000851002 Homo sapiens Epithelial membrane protein 2 Proteins 0.000 description 2
- GCNLQHANGFOQKY-UHFFFAOYSA-N [C+4].[O-2].[O-2].[Ti+4] Chemical compound [C+4].[O-2].[O-2].[Ti+4] GCNLQHANGFOQKY-UHFFFAOYSA-N 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002904 solvent Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000790917 Dioxys <bee> Species 0.000 description 1
- YCKRFDGAMUMZLT-UHFFFAOYSA-N Fluorine atom Chemical compound [F] YCKRFDGAMUMZLT-UHFFFAOYSA-N 0.000 description 1
- 229910002090 carbon oxide Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000000329 molecular dynamics simulation Methods 0.000 description 1
- 238000000302 molecular modelling Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 125000001997 phenyl group Chemical group [H]C1=C([H])C([H])=C(*)C([H])=C1[H] 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000751 protein extraction Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Manufacture Of Macromolecular Shaped Articles (AREA)
- Treatments Of Macromolecular Shaped Articles (AREA)
Abstract
本发明涉及一种构建高分子与超临界二氧化碳溶液混合体系全原子模型的方法,属于计算机技术领域。所述方法包括如下步骤:1.构建超临界二氧化碳溶液体系的全原子模型;2.构建高分子的初始空间坐标;3.构建高分子的全原子模型对应的拓扑参数;4.对高分子进行能量优化;5.将高分子加入超临界溶液体系;6.对体系进行高温退火模拟以得到稳定的高分子与超临界二氧化碳混合体系的全原子模型。通过该方法,可以构建各种不同种类的高分子与超临界二氧化碳溶液组成的混合体系所对应的全原子模型,为进一步研究不同高分子对超临界二氧化碳溶液微观动力学的影响规律以及超临界二氧化碳的应用开发奠定有效的技术基础。
Description
技术领域
本发明涉及计算机技术领域,具体涉及一种构建高分子与超临界二氧化碳溶液混合体系全原子模型的方法。
背景技术
超临界二氧化碳是指二氧化碳在温度和压力大于临界点时的独特状态。它具有接近于空气的低粘度和高扩散系数,同时具有接近于液体的高密度。它制备方法简单且易于循环利用,属于环境友好型绿色溶剂。它不仅可以作为绿色溶剂参与多种化学反应,还可以用于蛋白质提取和其他生物分子的萃取。更重要的是,超临界二氧化碳可以用于连续油管钻井,同时还可以提高原油采收率,对储集层没有任何污染。这对于油气藏的开发有着重要的意义。
然而由于二氧化碳分子直线式的分子构型与较高的分子极性,导致其超临界状态下的粘滞系数可调区间非常窄,导致其无法满足生物分子萃取和石油开采等领域对高粘滞溶液的需求。在对超临界二氧化碳粘性的研究过程中,大量的实验研究表明由带有苯环非极性侧链和具有氟烃极性侧链构成的特定比例的高分子能够有效地增加超临界二氧化碳的粘滞系数。因此,研究人员希望从分子角度设计不同种类的高分子来调控超临界二氧化碳的粘性,从而拓展其在不同领域的应用。然而由于高分子单体种类繁多,且高分子中单体的配比难以控制,导致实验上无法实现对高分子与超临界二氧化碳之间相互作用规律定量研究。
随着计算机技术以及分子建模方法的成熟,人们开始尝试利用精确的全原子分子动力学方法来研究各种体系中分子之间的微观动力学以及相互作用规律。该方法具有环境参数精确可调和分辨率高等优点。通过全原子分子动力学模拟,我们能够在原子尺度观察二氧化碳分子在高分子不同官能团附近的微观动力学过程,同时能够定量地研究它们与二氧化碳分子的微观相互作用机制,最终阐明高分子调控超临界二氧化碳粘性的定量规律。因此,我们提供了基于全原子分子动力学构建高分子与超临界二氧化碳混合体系模型的方法,为实现计算机辅助设计不同粘性的高分子与超临界二氧化碳的混合体系奠定基础。
发明内容
本发明要解决的问题是提供一种构建高分子与超临界二氧化碳溶液混合体系全原子模型的方法。可以利用成功构建的全原子模型进一步研究高分子与超临界二氧化碳液体在原子尺度的动力学和相互作用问题,如高分子调控超临界二氧化碳粘性规律等等。
为实现上述目的,一种基于全原子分子动力学构建高分子与超临界二氧化碳混合体系模型的方法,包含如下步骤:
(1)利用Gromacs软件,对二氧化碳的EPM2全原子在NPT系综下进行高温退火模拟,构建稳定的超临界二氧化碳溶液体系,其中二氧化碳分子约为11000个;
(2)利用Gauss View构建由苯乙烯和含氟丙烯酸酯构成的特定比例的嵌段高分子,并保存为pdb文件格式;
(3)基于成熟的OPLSAA全原子力场以及高分子原子的化学成键类型,确定高分子中所有原子对应的原子类型,同时确定高分子包含的单体内原子之间的拓扑关系,并产生两个单体对应的top文件格式;
(4)基于链接单体之间原子成键作用的强弱,确定高分子单体之间的连接拓扑关系,然后根据高分子包含的单体之间的排列方式,将单体对应的拓扑和单体之间的拓扑关系按顺序关系进行组合,生成该高分子对应完整拓扑关系,并产生相应的top文件格式;
(5)利用Gromacs以及选定的全原子力场对初始的pdb文件中的高分子构型进行能量最小化,从而得到高分子的能量最优构型
(6)将能量优化后的高分子加入超临界二氧化碳溶液体系,并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟以消除二氧化碳和高分子的不合理接触,最终得到稳定的高分子与超临界二氧化碳混合体系的全原子模型;
其中,步骤(1)和(6)进行NPT系综下的分子动力学模拟时,动力学参数设置为:采用Leap-Frog算法进行积分求解,对应的积分步长为0.5fs;采用Particle Mesh Ewald算法对长程静电作用进行求解,其截断半径为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。该值与实验密度值869.2Kg/m3接近。
(2)利用Gauss View构建由极性含氟丙烯酸酯和非极性苯乙烯构成的29:1的嵌段高分子,并保存为pdb文件格式,其中高分子单体的全原子构型及其标注如附图2所示,而高分子的三维原子构型示意图如附图3所示;
(3)基于成熟的OPLSAA全原子力场以及高分子中原子的化学成键类型,确定高分子中所有原子对应的原子类型,同时确定高分子中单体内原子之间的拓扑关系,并产生两个单体对应的top文件;
(4)基于链接单体之间原子成键作用的强弱,确定高分子单体之间的连接拓扑关系。由于该高分子的单体比例为29:1,因此其包含两类连接关系:第一类是极性单体与极性单体之间的连接拓扑,第二类是极性单体与非极性单体之间的连接拓扑。针对第一类连接拓扑,我们选定含氟丙烯酸酯中如附图2a中所示C1,C2,C3,O1,O2原子与含氟丙烯酸酯中如附图2a中所示C1,C2,C3,O1,O2原子之间的连接关系作为连接拓扑,针对第二类连接拓扑,我们选定含氟丙烯酸酯中如附图2a中所示C1,C2,C3,O1,O2原子与苯乙烯中如附图2b中所示C1,C2,C3,C4,C8原子之间的连接作为拓扑。然后根据高分子包含的单体之间的嵌段排列方式,将单体对应的拓扑和两类连接拓扑按顺序关系进行组合,生成该高分子对应完整拓扑关系,并产生相应的top文件;
(5)基于选定的全原子力场模型,利用Gromacs对初始pdb文件中的高分子构型进行能量最小化;
(6)将结构优化后的高分子加入超临界二氧化碳溶液体系,并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟,消除二氧化碳和高分子的不合理接触,最终得到稳定的高分子与超临界二氧化碳混合体系的全原子构型如附图4所示;
至此,从步骤(1)到(6),完成了对极性单体和非极性单体比例为29:1的高分子与超临界二氧化碳混合体系对全原子模型的构建。
实施例2
本实施例详细阐述了本发明“一种构建高分子与超临界二氧化碳混合体系全原子模型的方法”在具体实施时针对极性单体和非极性单体比例为15:15的高分子与超临界二氧化碳混合体系对全原子模型的构建,具体步骤如下:
(1)采用EMP2全原子模型对二氧化碳进行全原子建模,利用Gromacs软件对二氧化碳体系在NPT系综下进行高温退火模拟,其退火模拟过程依次从(1000bar,800K,)至(800bar,600K,),至(600bar,400K,)至(400bar,320K,)最后至(280bar,320K)。得到稳定的超临界二氧化碳体系,其密度随时间变化曲线如附图1所示,平均值为830.2Kg/m3。该值与实验密度值869.2Kg/m3接近。
(2)利用Gauss View构建由极性含氟丙烯酸酯非极性苯乙烯构成的15:15的嵌段高分子,并保存为pdb文件格式,其中高分子单体的全原子构型及其标注如附图2所示,而高分子的三维原子构型示意图如附图5所示;
(3)基于成熟的OPLSAA全原子力场以及高分子中原子的化学成键类型,确定高分子中所有原子对应的原子类型,同时确定高分子中单体内原子之间的拓扑关系,并产生两个单体对应的top文件;
(4)基于链接单体之间原子成键作用的强弱,确定高分子单体之间的连接拓扑关系。由于该高分子的单体比例为15:15,因此其包含三类连接关系:第一类是极性单体与极性单体之间的连接拓扑,第二类是极性单体与非极性单体之间的连接拓扑,第三类是非极性单体与非极性单体之间的连接拓扑。针对第一类连接拓扑,我们选定含氟丙烯酸酯中如附图2a中所示C1,C2,C3,O1,O2原子与含氟丙烯酸酯中如附图2a中所示C1,C2,C3,O1,O2原子之间的连接关系作为连接拓扑;针对第二类连接拓扑,我们选定含氟丙烯酸酯中如附图2a中所示C1,C2,C3,O1,O2原子与苯乙烯中如附图2b中所示C1,C2,C3,C4,C8原子之间的连接作为拓扑;针对第三类连接拓扑,我们选定苯乙烯中如附图2b中所示C1,C2,C3,C4,C8原子与苯乙烯中如附图2b中所示C1,C2,C3,C4,C8原子之间的连接关系作为连接拓扑。然后根据高分子包含的单体之间的嵌段排列方式,利用自主开发的工具程序自动将单体对应的拓扑和三类连接拓扑按顺序关系进行组合,生成该高分子对应完整拓扑关系,并产生相应的top文件
(5)基于选定的全原子力场模型,利用利用Gromacs对初始pdb文件中的高分子构型进行能量最小化;
(6)将能量优化后的高分子加入超临界二氧化碳溶液体系,并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟以消除二氧化碳和高分子的不合理接触。最终得到稳定的高分子与超临界二氧化碳混合体系的全原子构型如附图6所示;
至此,从步骤(1)到(6),完成了对极性单体和非极性单体比例为15:15的高分子与超临界二氧化碳混合体系对全原子模型的构建。
在已构建的两种高分子与超临界二氧化碳混合体系全原子模型的基础之上,利用Gromacs软件对这两种混合体系进行4ns的全原子分子动力学模拟,并分析高分子对二氧化碳分子的微观动力学的影响,如二氧化碳分子的微观MSD随时间变化曲线如附图7所示。可以看出,高分子的加入使得超临界二氧化碳的粘滞性发生了变化,并且该变化与加入的高分子的种类有关。这说明,依据本发明“一种构建高分子与超临界二氧化碳混合体系全原子模型的方法”构建的高分子与超临界二氧化碳混合体系的全原子模型具备合理性。
以上所述为本发明的两个典型实施例而已,本发明不应该局限于该实施例和附图所公开的内容。任何属于本领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (8)
1.一种构建高分子与超临界二氧化碳混合体系全原子模型的方法,其特征在于:
所述方法具体步骤为:
(1)利用Gromacs软件,对二氧化碳的EPM2全原子在NPT系综下进行高温退火模拟,构建稳定的超临界二氧化碳溶液体系;
(2)利用Gauss View构建由苯乙烯(Styrene)和含氟丙烯酸酯(Fluoroacrylate)构成的特定比例的嵌段高分子;
(3)基于成熟的OPLSAA全原子力场以及高分子原子的化学成键类型,确定高分子中所有原子对应的原子类型,同时确定高分子包含的单体内原子之间的拓扑关系;
(4)基于链接单体之间原子成键作用的强弱,确定高分子单体之间的连接拓扑关系,然后根据高分子包含的单体之间的排列方式,将单体对应的拓扑和单体之间的拓扑关系按顺序关系进行组合,生成该高分子对应完整拓扑关系;
(5)利用Gromacs以及选定的全原子力场对高分子构型进行能量最小化,从而得到高分子的能量最优构型;
(6)将能量优化后的高分子加入超临界二氧化碳溶液体系,并利用Gromacs软件对高分子与超临界二氧化碳混合体系进行NPT系综下的高温退火模拟,最终得到稳定的高分子与超临界二氧化碳混合体系的全原子模型。
2.根据权利要求1所述的方法,其特征在于:步骤(1)和(6)中,采用Leap-Frog算法进行积分求解,对应的积分步长为0.5fs;采用Particle Mesh Ewald算法对长程静电作用进行求解,其截断半径为0.9nm;Lenard-Jones相互作用的截断半径为1.3nm;采用Berendsen热浴算法来恒定模拟的温度,其耦合常数为0.5ps;采用Berendsen压浴算法来恒定模拟的压强,其耦合常数为1.0ps。
3.根据权利要求1所述的方法,其特征在于:在步骤(1)和(6)中,进行NPT系综下的高温退火模拟,具体退火过程依次为(1000bar,800K)-(800bar,600K)-(600bar,400K)-(400bar,320K)-(280bar,320K)。
4.根据权利要求1所述的方法,其特征在于步骤(2)中所述拟定的高分子极性单体与非极性单体之间以嵌段模式进行组合,且极性单体与非极性单体之间的比例为29:1~15:15。
5.根据权利要求1所述的方法,其特征在于:步骤(2)中高分子的极性单体为含氟丙烯酸酯,而非极性单体为苯乙烯,两个单体对应的全原子模型为OPLSAA力场。
6.根据权利要求1所述的方法,其特征在于:步骤(3)中高分子单体之间的连接拓扑包括所有与主链原子之间拓扑距离小于等于3的非氢原子之间的成键关系构成,同时还包括最近邻原子间的成键关系,次近邻原子间的键角关系,以及次次近邻原子间的二面角关系。
7.根据权利要求1所述的方法,其特征在于:步骤(6)中,将高分子加入到超临界二氧化碳时,高分子表面的范德瓦尔半径设置为0.35nm,且体系的高分子与二氧化碳的混合体系对应的高分子占体系的质量分数设定为0.19wt%左右。
8.根据权利要求1所述的方法,其特征在于步骤(6)中消除高分子和二氧化碳之间的不合理接触的过程中,采用了NPT系综下的高温退火模拟实现初始构型的弛豫。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811276321.9A CN109473144B (zh) | 2018-10-30 | 2018-10-30 | 一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811276321.9A CN109473144B (zh) | 2018-10-30 | 2018-10-30 | 一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109473144A true CN109473144A (zh) | 2019-03-15 |
CN109473144B CN109473144B (zh) | 2021-07-09 |
Family
ID=65666432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811276321.9A Expired - Fee Related CN109473144B (zh) | 2018-10-30 | 2018-10-30 | 一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109473144B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110136782A (zh) * | 2019-05-28 | 2019-08-16 | 北京科技大学 | 一种金属氧化物团簇结构模型的构造方法 |
CN111681702A (zh) * | 2020-06-09 | 2020-09-18 | 青岛科技大学 | 一种基于分子动力学的生物燃料提纯方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1649040A (zh) * | 2004-12-04 | 2005-08-03 | 西南师范大学 | 可调变粘滞度磁性液体的制备方法 |
CN1662571A (zh) * | 2002-06-14 | 2005-08-31 | Ami-阿格罗林茨三聚氰胺国际有限公司 | 用于制备氨基树脂产品的组合物和用这些组合物制备产品的方法 |
CN101033269A (zh) * | 2007-02-12 | 2007-09-12 | 杭州师范学院 | 一种原子转移自由基聚合引发体系及其在制备高分子化合物的应用 |
US20110061742A1 (en) * | 2009-09-11 | 2011-03-17 | Tyco Valves & Controls Lp | Cryogenic liquid cylinder manifold |
CN102350229A (zh) * | 2011-08-24 | 2012-02-15 | 北京中电加美环境工程技术有限责任公司 | 一种金属元素改性zsm-5分子筛复合膜及其制备方法和应用 |
CN104657564A (zh) * | 2015-03-16 | 2015-05-27 | 长春理工大学 | 基于分子动力学磨粒流加工数值模拟研究方法 |
CN104778329A (zh) * | 2015-04-24 | 2015-07-15 | 中国石油大学(华东) | 一种分析co2/ch4在褐煤中吸附机理的方法 |
WO2015166275A1 (en) * | 2014-05-02 | 2015-11-05 | Oxford Nanopore Technologies Limited | Mutant pores |
WO2016033368A1 (en) * | 2014-08-27 | 2016-03-03 | Ohio State Innovation Foundation | Improved peptidyl calcineurin inhibitors |
CN105694358A (zh) * | 2016-04-08 | 2016-06-22 | 天津工业大学 | 磺化聚醚醚酮-磺化氧化石墨烯杂化膜及制备和应用 |
CN106432615A (zh) * | 2016-09-20 | 2017-02-22 | 新辉(中国)新材料有限公司 | 一种通co2气体能够絮凝的聚合物乳液及制备、絮凝的方法 |
CN107330235A (zh) * | 2017-05-08 | 2017-11-07 | 华南理工大学 | 一种多肽链自组装过程的分子动力学研究方法 |
CN108108540A (zh) * | 2017-12-12 | 2018-06-01 | 华南理工大学 | 一种模拟tbpb半笼形水合物膜分离混合气体的方法 |
CN108226080A (zh) * | 2018-01-17 | 2018-06-29 | 中国石油大学(北京) | 一种高压原位红外样品池及其在co2-有机液体体系中的应用 |
-
2018
- 2018-10-30 CN CN201811276321.9A patent/CN109473144B/zh not_active Expired - Fee Related
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1662571A (zh) * | 2002-06-14 | 2005-08-31 | Ami-阿格罗林茨三聚氰胺国际有限公司 | 用于制备氨基树脂产品的组合物和用这些组合物制备产品的方法 |
CN1649040A (zh) * | 2004-12-04 | 2005-08-03 | 西南师范大学 | 可调变粘滞度磁性液体的制备方法 |
CN101033269A (zh) * | 2007-02-12 | 2007-09-12 | 杭州师范学院 | 一种原子转移自由基聚合引发体系及其在制备高分子化合物的应用 |
US20110061742A1 (en) * | 2009-09-11 | 2011-03-17 | Tyco Valves & Controls Lp | Cryogenic liquid cylinder manifold |
CN102350229A (zh) * | 2011-08-24 | 2012-02-15 | 北京中电加美环境工程技术有限责任公司 | 一种金属元素改性zsm-5分子筛复合膜及其制备方法和应用 |
WO2015166275A1 (en) * | 2014-05-02 | 2015-11-05 | Oxford Nanopore Technologies Limited | Mutant pores |
WO2016033368A1 (en) * | 2014-08-27 | 2016-03-03 | Ohio State Innovation Foundation | Improved peptidyl calcineurin inhibitors |
CN104657564A (zh) * | 2015-03-16 | 2015-05-27 | 长春理工大学 | 基于分子动力学磨粒流加工数值模拟研究方法 |
CN104778329A (zh) * | 2015-04-24 | 2015-07-15 | 中国石油大学(华东) | 一种分析co2/ch4在褐煤中吸附机理的方法 |
CN105694358A (zh) * | 2016-04-08 | 2016-06-22 | 天津工业大学 | 磺化聚醚醚酮-磺化氧化石墨烯杂化膜及制备和应用 |
CN106432615A (zh) * | 2016-09-20 | 2017-02-22 | 新辉(中国)新材料有限公司 | 一种通co2气体能够絮凝的聚合物乳液及制备、絮凝的方法 |
CN107330235A (zh) * | 2017-05-08 | 2017-11-07 | 华南理工大学 | 一种多肽链自组装过程的分子动力学研究方法 |
CN108108540A (zh) * | 2017-12-12 | 2018-06-01 | 华南理工大学 | 一种模拟tbpb半笼形水合物膜分离混合气体的方法 |
CN108226080A (zh) * | 2018-01-17 | 2018-06-29 | 中国石油大学(北京) | 一种高压原位红外样品池及其在co2-有机液体体系中的应用 |
Non-Patent Citations (3)
Title |
---|
GURINA D L等: "Structural Features of Binary Mixtures of Supercritical CO2 with Polar Entrainers by Molecular Dynamics Simulation", 《RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY》 * |
史俊勤: "超临界二氧化碳溶胀原油及降低油水界面张力的分子模拟研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
孙鑫等: "二氧化碳压裂技术研究综述", 《油田化学》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110136782A (zh) * | 2019-05-28 | 2019-08-16 | 北京科技大学 | 一种金属氧化物团簇结构模型的构造方法 |
CN111681702A (zh) * | 2020-06-09 | 2020-09-18 | 青岛科技大学 | 一种基于分子动力学的生物燃料提纯方法 |
CN111681702B (zh) * | 2020-06-09 | 2022-09-13 | 青岛科技大学 | 一种基于分子动力学的生物燃料提纯方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109473144B (zh) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | A FTIR and 2D-IR spectroscopic study on the microdynamics phase separation mechanism of the poly (N-isopropylacrylamide) aqueous solution | |
CN109473144A (zh) | 一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 | |
Thickett et al. | Emulsion polymerization: State of the art in kinetics and mechanisms | |
Charleux | Theoretical aspects of controlled radical polymerization in a dispersed medium | |
CN107417861B (zh) | 一种耐盐性堵水剂的制备方法 | |
CN106947019A (zh) | 一种水凝胶微球的制备方法及装置 | |
Thomson et al. | Nitroxide-mediated surfactant-free emulsion polymerization of n-butyl methacrylate with a small amount of styrene | |
Brunel et al. | New insight into cluster aggregation mechanism during polymerization-induced self-assembly by molecular dynamics simulation | |
JP4665614B2 (ja) | ポリマー材料の設計システム | |
Chang et al. | Simultaneously construct IRT-based parallel tests based on an adapted CLONALG algorithm | |
Gromer et al. | Simulation of latex film formation using a cell model in real space: Vertical drying | |
Lazzari et al. | Interplay between aggregation and coalescence of polymeric particles: experimental and modeling insights | |
CN110167993A (zh) | 乳液颗粒、包含它的乳液及乳液的制备方法 | |
Saldívar et al. | Mathematical modeling of emulsion copolymerization reactors: Experimental validation and application to complex systems | |
CN113012763B (zh) | 一种基于四族组分的原油氧化反应动力学模型建立方法 | |
Gurr et al. | Synthesis, characterization, and direct observation of star microgels | |
CN104525061B (zh) | 基于离子液体和两亲性Fe3O4纳米粒子的磁响应性皮克林乳液 | |
Sood | Particle size distribution control in emulsion polymerization | |
CN102786608B (zh) | 一种界面自由基乳液聚合方法 | |
CN105482029A (zh) | 一种节能环保生物基水性乳液聚合物及其制备方法 | |
CN114783537A (zh) | 一种降阻剂的优选及性能评价方法及其应用 | |
CN105160208A (zh) | 用于疾病亚型问题的基于网络的聚类方法 | |
CN106055887A (zh) | 基于相场模型下控制Ag2Ga纳米针长径比的研究方法 | |
CN108334737A (zh) | 一种基于分子模拟的破乳剂筛选方法 | |
Wilding et al. | Development and Experimental Validation of a Dispersity Model for In Silico RAFT Polymerization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210709 |