丙烯腈/丙烯酸甲酯/衣康酸三元共聚物序列结构的Monte Carlo模拟

以丙烯腈(AN)、丙烯酸甲酯(MA)、衣康酸(IA)为共聚单体,合成碳纤维用聚丙烯腈原丝的三元共聚物,利用Monte Carlo法模拟了AN／MA／IA三元共聚物的组成和序列结构,讨论了共聚单体对共聚物分子链序列结构的影响,探讨了增加序列结构均匀性的方案。结果表明:Monte Carlo法可以较为准确的模拟该共聚物的序列结构;AN／MA／IA以98／0.7／1.3的单体配比可以得到序列结构比较理想的共聚物;采取分批加入IA的方法,使共聚物的序列结构均匀性得到提高。     

P(CTFE-BVE-SUA)含氟无皂乳液的制备

以无皂乳液聚合法合成了聚(三氟氯乙烯-乙烯基正丁基醚-十一烯酸钠)[P(CTFE-BVE-SUA)]含氟共聚物乳液,考察了单体配比对聚合反应的影响,同时对影响乳液粒径的因素进行了研究。结果表明:以无皂乳液聚合法制得的含氟乳液稳定性好、粒径分布均匀。通过控制单体配比可以得到不同结构、不同氟含量的乳液。聚合反应中过硫酸钾(KPS)的浓度、聚合温度以及SUA用量对乳液粒径的影响较大,试验中还发现,聚合反应过程中压力的变化与反应的进程是同步的。     

高分子构象模拟研究中的Monte Carlo方法

介绍了Monte Carlo方法及其特点,进而分析了Monte Carlo用于高分子模拟的优势,并描述了两类模拟模型。论文重点综述了近年来Monte Carlo方法在高分子构象模拟中的一些研究与应用,并展望了Monte Carlo方法在高分子构象模拟中的发展趋势和前景。     

MATLAB和Monte Carlo法在连续反应动力学中的应用

基于MATLAB语设计编写了Monte Carlo法处理连续化学反应的动力学计算及其计算机模拟的应用程序,将所编程序用于模拟3-丁烯腈在盐酸催化剂溶液中的水解反应,模拟结果与按照常规化学反应动力学公式的计算结果相比较,表明Monte Carlo法模拟对于预测反应动力学过程具有较高准确性。     

丁二酰亚胺无灰分散剂作用机理的分子模拟研究

通过对30种丁二酰亚胺类无灰分散剂在水和石墨表面作用的Monte Carlo Docking模拟,得到其在水和石墨表面的低能量吸附构象及结合能.QSAR分析和Monte Carlo Docking模拟的结果表明分散剂在水表面的吸附行为,与分散剂分子的极性基团的数量、分子正负电荷的分布、分子在水中的溶解自由能有关;而分散剂在石墨表面的吸附,与分散剂分子的链长度和支链数目有关.在此基础上合成了9个丁二酰亚胺无灰分散剂,对其分散烟炱(碳黑)的性能进行了模拟评定,通过对这些分散剂作用的计算机模拟,验证了模拟方法的可行性、可操作性.     

无皂乳液聚合制备AANa-VAc-TFEMA纳米含氟乳液

以丙烯酸钠(AANa)、乙酸乙烯酯(VAc)和甲基丙烯酸三氟乙酯(TFEMA)为单体进行无皂乳液聚合,合成了P(AANa-VAc-TFEMA)纳米含氟乳液,考察了反应温度、引发剂用量、单体配比对聚合的影响,并着重对聚合物乳液的粒径进行了研究。结果表明:该聚合单体的转化率可达到97.00%以上,聚合得到的无皂乳液稳定性好、粒径分布均匀。通过控制单体配比中AANa的含量及VAc与TFEMA的比例,可以得到纳米含氟乳液。     

超临界CO_2中CTFE和BVE的共聚合研究

以超临界CO2为反应介质,制备得到三氟氯乙烯(CTFE)和乙烯基正丁基醚(BVE)共聚物,研究了单体配比、反应温度、引发剂用量、单体浓度对聚合反应及聚合物的结构与性能的影响情况。结果表明:通过控制单体配比及其他实验条件,可以得到不同结构、不同氟含量的含氟聚合物;当单体配比中CTFE的比例增大到一定程度(BVE/CTFE质量比<6/10),实验得到的聚合物是交替共聚物。同时实验发现,聚合反应过程中压力的变化与反应的进程是同步的。     

丁醇和高级烷烃气液相平衡的分子模拟

用Gibbs系综的Monte Carlo法模拟了正丁醇。正己烷和正丁醇-正庚烷体系的气液相平衡。在恒定压力0.496 6 MPa、0.294 0 MPa和0.1 MPa下,分别采用恒温恒压的Gibbs系综(NPT-Gibbs)模拟了不同温度下的气液相平衡。计算结果表明,模拟的数据比较准确,在一定温度范围内与实验结果吻合。     

催化裂解过程分子尺度反应动力学模型研究

以工业数据为基础,运用结构导向集总方法来构造烃分子,对催化裂解原料油进行了分子尺度上的Monte Carlo模拟.在深入研究催化裂解反应机理的基础上,以原料油模拟产生的分子矩阵作为反应物分子,将结构导向集总方法与Monte Carlo方法相结合,建立了催化裂解过程分子尺度的反应动力学模型.结果表明:Monte Carlo方法可以在分子尺度上实现对催化裂解原料较好的模拟,产品产率和汽油组成与实际值能较好地拟合.且随着虚拟分子数的增加,对原料油性质和反应结果的模拟计算精度提高;模型具有较好的适应性和外推性.     

乳液聚合成核阶段的模拟与分析

建立了乳液聚合成核阶段的Monte Carlo模型，并用计算机对一个体积为10^-17m^3的微型反应器中苯乙烯的乳液聚合进行了模拟。以计算机生成随机数作为自由基被胶束和乳胶粒捕获的几率，模拟了在微型反应器中每一个自由基的生成、被胶束或乳胶粒捕获的过程以及每一个乳胶粒的生成及增长过程。通过对每一个乳胶粒在增长过程中各参数的统计计算，研究了乳液聚合成核阶段诸参数(乳胶粒数目、乳胶粒直径与粒径分布、单体转化率、聚合反应速率等)与乳化剂浓度[S]及引发剂浓度[I]的关系。结果表明，苯乙烯的乳液聚合体系中乳胶粒数目与[S]^0.5996[I]^0.4016成正比：在成核阶段乳胶粒直径分布先变宽后变窄，乳液聚合过程中乳胶粒直径分布有自动变窄的趋势；成核阶段持续时间t12与[S]^0.60[I]^0.60成正比，成核阶段结束时的单体转化率X12与[S]^1.20[I]^0.20成正比。     

Monte Carlo simulation of sequence distributions of acrylonitrile copolymers

A Monte Carlo simulation examining the effect of monomer ratios on the composition and sequence distribution of acrylonitrile(AN) copolymers with N‐vinyl pyrrolidone (NVP), itaconic acid (IA), and acrylic acid (AA) as comonomers has been developed. The Kelen–Tudos method was used to estimate monomer reactivity ratios. The results of the simulation are consistent with the academic conclusion and are as foreseen by the experimental data. The average number of NVP identical monomers in a sequence length of AN/NVP copolymer chain increases continuously and the average number of AN identical monomers in a sequence length shows a prominent decrease with an increase of NVP concentration in the feed. Changes in the monomer average number of AN/IA and AN/AA copolymers in a sequence length were the same as those of AN/NVP copolymer with an increase of comonomer concentration in the feed. The optimum weight ratio of AN with comonomers for manufacturing carbon fibers is 98/2. © 2005 Wiley Periodicals, Inc. J Appl Polym Sci 96: 483–488, 2005     

Monte Carlo simulation of sequential structure control of AN-MA-IA aqueous copolymerization by different operation modes

The regulation of polyacrylonitrile (PAN) copolymer composition and sequence structure is the precondition for producing high-quality carbon fiber high quality. In this work, the sequential structure control of acrylonitrile (AN), methyl acrylate (MA) and itaconic acid (IA) aqueous copolymerization was investigated by Monte Carlo (MC) simulation. The parameters used in Monte Carlo were optimized via machine learning (ML) and genetic algorithms (GA) using the experimental data from batch copolymerization. The results reveal that it is difficult to control the aqueous copolymerization to obtain PAN copolymer with uniform sequence structure by batch polymerization with one-time feeding. By contrary, it is found that the PAN copolymer with uniform composition and sequence structure can be obtained by adjusting IA feeding quantity in each reactor of a train of five CSTRs. Hopefully, the results obtained in this work can provide valuable information for the understanding and optimization of AN copolymerization process to obtain high-quality PAN copolymer precursor.     

Herstellung von Copolymeren mit definierter Zusammensetzungsverteilung im Zulaufverfahren

The possibility of the computation of binary and ternary copolymerization for semi-batch reactions is demonstrated. For that purpose the conversion must be known as a function of time and length of the monomer feed program. The computer simulations are successfully applied to the binary system methacrylic acid (N-dimethylaminoethyl)ester/styrene and to the ternary systems methacrylic acid (N-dimethylaminoethyl)ester/styrene/acrylonitrile as well as to methacrylic acid (2-hydroxyethyl)ester/styrene/acrylic acid n-butylester. The computation of the chemical conversion heterogeneity is possible with these simulations. Furthermore by these computations it is possible to find monomer feed programs by which copolymers with defined distribution of the chemical composition are generated. This procedure allows to generate a more homogeneous terpolymer methacrylic acid (2-hydroxylethyl)ester/styrene/acrylic acid butylester by computed monomer feed programs instead of the normally very heterogeneous terpolymer system. It is shown for the ternary system methacrylic acid (N-dimethylaminoethyl)ester/styrene/acrylonitrile which initial monomer compositions must be chosen in order to generate chemically homogenous terpolymers. In the latter system an internal ternary azeotropic region supposedly exists.     

核壳型聚丙烯酸酯乳液的耐水性能研究

采用种子乳液二阶段聚合法,制备了经丙烯酸(AA)/丙烯酰胺(AM)交联、具有核壳结构的甲基丙烯酸甲酯(MMA)/丙烯酸丁酯(BA)/丙烯腈(AN)复合乳液,研究了核壳软硬单体比、AN及壳层中交联单体AM含量对乳胶膜吸水率的影响。结果表明:乳胶膜吸水率主要受壳层软硬单体比例的影响,软硬单体比值越小,乳胶膜的吸水率越低;丙烯酸/丙烯酰胺复合交联剂对乳胶膜吸水率的影响不明显;功能单体丙烯腈可明显降低乳胶膜的吸水率,乳胶膜吸水率随丙烯腈含量的增加而逐渐降低,最低达到9.40%。     

Grafting of acrylonitrile and methyl methacrylate from their binary mixtures on cellulose using ceric ions

Ceric ion‐initiated grafting on cellulose from a binary mixture of acrylonitrile and methyl methacrylate was carried out in heterogeneous and acidic conditions at 30 ± 0.1°C in a nitrogen atmosphere. To avoid the complexation of water molecules with Ce(IV) ions, the concentration of the nitric acid was taken to be more than the concentration of ceric ions. The effect of the feed concentration, reaction time, and ceric ions concentration on grafting were investigated at a fixed composition. To investigate the effect of monomer–monomer interactions on grafting, the graft copolymerization was also studied, using different feed compositions (f_AN) ranging from 0.25 to 0.80. In this range of feed composition, the synergistic effect of methyl methacrylate molecules has shown an important effect on acrylonitrile monomer and facilitate the incorporation of the acrylonitrile monomer into the grafted chains. The reactivity ratios of acrylonitrile and methyl methacrylate were calculated using the Mayo and Lewis method and were found to be 0.74 and 1.03, respectively. The average sequence lengths of the monomers (m̄M) were found to be dependent on the feed compositions and found to be arranged in alternate fashion in the grafted chains. The probability of the addition of a monomer (P_1,1) to the growing radicals on cellulose ended with its own type of monomer was found to be dependent on the feed composition. The composition of the grafted copolymers, homocopolymers, was determined by IR and elemental analysis for nitrogen. None of the grafted chain on cellulose was found to be made of a single type of monomer. The ceric ion consumption during grafting was found to be independent of the molarity of the feed but shown an appreciable change in the initial few hours of grafting. The variation in the values of the grafting parameters as a function of the reaction conditions is suitably explained. © 2000 John Wiley & Sons, Inc. J Appl Polym Sci 79: 767–778, 2001     

Simulation of reversible chain transfer catalyzed polymerization (RTCP): effect of different iodide based catalysts

Monte Carlo kinetic simulation method was performed to study the reversible chain transfer catalyzed polymerization (RTCP) of styrene in 80°C. The effect of different iodide based catalysts (AI) was investigated on RTCP systems by simulating the chain length distribution, the depletion rate of polymerization ingredients, average molecular weights and the monomer conversion. In RTCP systems a narrow distribution was obtained compared to iodide-mediated polymerization. Superior reversible chain transfer reaction constants (k_a and k_da) led to a more uniform chain length distribution and a faster PDI decrement. In each RTCP system k_da/k_a ratio, designates the concentration of A* by dictating the dominant side of exchange equilibriums which specifies the number of cross-termination reactions. Addition of cross-termination reactions to the iodide-mediated polymerization system decreases the number of combination reactions leading to lower average molecular weights. Higher k_da/k_a ratios also consequences in faster catalyst depletions and lessen monomer conversions. On the other hand, PhE-I was found to deplete rapidly in RTCP systems, changing the nonlinear increase of the number average molecular weight to a linear pattern. The Monte Carlo simulation results were in a fine agreement with experimental data which were obtained from different RTCP systems.     

New aspects of unsaturated polyester resin synthesis. Part 2. Reactant sequence distribution and its effect on cure kinetics

The distribution of unsaturations in the prepolymer of a typical unsaturated polyester (UP) resin (maleic anhydride, phthalic anhydride and 1,2‐propylene glycol) has been shown to influence the kinetics of the cure process with styrene monomer. Segments containing double bonds in close proximity appear to lower the reactivity of the resin due to steric hindrance, as indicated by the fact that the rate of cure and the final degree of cure, measured by differential scanning calorimetry (DSC), increase as the average sequence length (SL) of maleic units decreases. This implies that the reactivity of UP resins may be improved by synthesis of prepolymers with certain reactant sequence‐length distributions. The copolymer formed by the melt condensation process of maleic anhydride, phthalic anhydride and 1,2‐propylene glycol in the absence of a transesterification catalyst has a non‐random structure with a tendency towards blockiness. This was established using ¹H NMR analysis in tandem with deterministic and Monte Carlo modelling techniques. Copyright © 2003 Society of Chemical Industry     

Monte Carlo simulation of terpolymerization: Optimizing the simulation and post-processing times

Monte Carlo simulations are a useful and easy way to understand a polymerization reaction process properly. However, achieving reliable results with Monte Carlo simulations can also lead to prohibitive computational times and a considerable amount of data to be processed afterward. The present study analyses the Monte Carlo simulation of a steady-state terpolymerization process to reduce the overall computational time of the simulation and the post-processing of its results. Different sorting algorithms (Bubble, Insertion, Selection, and Tim) and Python libraries (Joblib and Numba) were used. The chain composition distribution and the micro-structures resultant of different scenarios were assessed by processing the simulated mechanism results. The simulation time results indicate the Tim sorting algorithm as the best to use in the post-processing step and the Numba library as the best suited for both the simulation and the post-processing step.