药学专业高分子化学的虚拟仿真教学探索——基于Python的自由基聚合反应模拟

通过Python语言开发了一款教学用自由基聚合反应数值模拟器FRPython。该程序通过对自由基聚合反应体系的常微分方程组进行数值模拟,得到引发剂自由基浓度、链自由基浓度、动力学链长,以及聚合度分布等信息。以聚甲基丙烯酸甲酯聚合体系的数值模拟为例,FRPython可直观有效地演示自由基聚合反应过程,同时避免了稳态假设、方程推演等较为深奥的内容,为药学专业高分子化学相关知识的教学提供有益参考。     

我国高分子合成化学的发展（二）

二、聚合反应与共聚合聚合反应与共聚合反应一直是高分子合成与改性的重要内容与手段,近年仍有重要发展。下面就自由基聚合、共聚合、正负离子聚合等方面的研究概况进行简要地阐述。1.自由基聚合近年的研究主要在开发含功能基的烯类单体、新型自由基引发体系和引发机理、聚合动力学尤其高转化的聚合动力学等。李笃信等研究了 N-乙烯基咔唑在 znO 为分散     

基于Monte Carlo模拟的化学反应动力学参数估算

提出并采用基于MonteCarlo模拟与动力学实验相结合的化学反应动力学参数估算方法,由基元反应确定MonteCarlo模拟具体做法,将MonteCarlo模拟结果与动力学实验结果相比较,根据比较结果自动调整和优化动力学参数,从而无需事先确定动力学方程即可有效估算各种化学反应的动力学参数值.采用该方法估算了丙烯氨氧化反应动力学参数,并对估算结果进行了分析与讨论.     

红外光谱法研究苯乙烯原子转移自由基聚合反应动力学

利用红外光谱法测定反应转化率研究原子转移自由基聚合反应动力学。在环己酮溶液中,以卤代烃为引发剂,过渡金属卤化物与配位剂２,２‘－联吡啶为催化体系,进行了苯乙烯聚合。分别就反应温度,反应物浓度和引发体系对苯乙烯聚合速率的影响进行了动力学测定,证实了原子转移自由基聚合具有活性聚合的特征,同时计算并讨论了苯乙烯原子转移自由基聚合反应的动力学数据。     

网格化教学设计课程“自由基聚合中的引发剂”

自由基聚合生产工艺基础是《聚合物制备工程》课程中非常重要的一章,主要涉及本体聚合、乳液聚合、悬浮聚合、溶液聚合四种聚合方法,而引发剂的选择和利用将极大地影响这四种聚合的效果。正确掌握引发剂的种类、特性、动力学过程以及不同种类引发剂在不同聚合反应中的选择规则,对于控制自由基聚合非常重要。但是由于引发剂种类繁多,且特性各不相同,学生在理解和记忆时往往会出现偏差,最终在实际生产中导致聚合失败。本文将利用网格化的教学方法,帮助学生建立引发剂的记忆网格。本着循序渐进与"高难度高速度"相结合的原则,利用网格化的课程设计,通过精心的"导课",深入浅出的"正课"以及探索求知的"结课",将繁杂的内容通过网格呈现给学生,并将复杂的动力学公式,分解成简单容易理解的基本方程呈现给学生,克服学生在学习过程中的畏难情绪,促进学生发挥主观能动性,实现"生本课堂"的教学理念。     

甲基丙烯酸甲酯本体自由基聚合的引发效率f

基于ESR技术测定的自由基浓度和终止速率常数,用严格的非稳态动力学公式第一次获得了AIBME引发MMA本体聚合体系各个转化率下的引发效率f数据,提出了扩散控制引发反应模型。在反应全程中理论结果与实验结果相符。     

极性单体阴离子型聚合反应产物的分子量分布

关于任意官能度多官能团引发剂瞬时引发并有单体链终止的阴离子型聚合反应体系,本工作通过非稳态动力学分析,求得了分子量分布函数和平均聚合度的一般表示式,讨论了单体的最大消耗量和聚合物的官能度分布问题,建立了从反应的初始条件和单体转化率计算产物的各种分子参数的方法。上述理论结果适用于甲基丙烯酸甲酯等极性单体在极性溶剂中的阴离子型聚合反应。数值计算的结果表明:当引发剂的官能度为2时,除了少数例外,所得聚合物的分子量分布一般具有双峰。     

过氧化物复合引发剂分解动力学模型及应用

过氧化物分解动力学是控制自由基聚合反应的重要因素 .单一的过氧化物在某一温度下只有固定的分解速率常数 ,而聚合反应往往需要在较宽的温度范围调节聚合反应速度和聚合物质量 .单一引发剂无法满足要求 ,工业上常常使用多种 (高、中、低活性 )引发剂复合 [1,2 ] .过氧化物引发剂分解速率常数测定方法主要有碘量法 [3] 和示差扫描量热法 (DSC) [4 ] 两种 .单一引发剂分解动力学基本符合Arrhenius模型 ,有时也会出现诱导分解现象 [5] . VC悬浮聚合单一引发剂用量已有理论估算 [6] ,但多元复合引发剂用量未见研究 ,史悠彰 [7] 曾提出用“折…     

亚硝酸钠引发丙烯腈聚合的研究

研究了亚硝酸钠引发硝酸溶液中丙烯腈的聚合反应。测得表现聚合速度 Rp=Ae~(-10,800/RT)[AN]~2.2[NaNO_2]~(0.17-1.0)[HNO_3]~(1.0-0.67 丙烯腈-丙烯酸甲酯共聚合反应中竞聚率分别是γ_An=0.96,γ_MA=1.17,表明聚合反应是按自由基机理进行。 根据聚合动力学和红外光谱分析,认为以亚硝酸钠引发硝酸溶液中丙烯腈的聚合反应与电解或金属溶蚀过程中的次级引发相同。     

碳纳米管的聚合物功能化与结构控制:聚合物接枝碳纳米管

碳纳米管高分子化是发展高性能的聚合物基纳米功能材料的重要研究方向,本文从"grafting-to"和"grafting-from"两种方式对聚合物接枝碳纳米管的最新进展进行了系统综述。"Grafting-to"方法主要包括羧基衍生反应(酰化、酯化)、加成反应(大分子自由基加成、叠氮环加成)和硫醇偶联反应。"Grafting-from"方法包括普通自由基聚合、可控/活性自由基聚合、离子聚合、开环聚合和逐步聚合反应,其中碳纳米管表面引发活性自由基聚合进一步分为原子转移自由基聚合、氮氧稳定自由基聚合和可逆加成-断链转移聚合。此外,本文还简述了碳纳米管自身的聚合反应,并探讨了目前聚合物修饰碳纳米管所面临的问题和今后的发展方向。     

苯乙烯在顺1,4-聚丁二烯上的接枝共聚反应动力学的Monte Carlo模拟

用ＭｏｎｔｅＣａｒｌｏ模拟方法对苯乙烯在顺１，４－聚丁二烯上的接枝反应动力学进行了研究。模拟结果表明，对于不同的橡胶浓度，接枝的和均聚的聚苯乙烯的数均聚合度、接枝效率与苯乙烯在顺１，４－聚丁二烯上的接枝共聚反应的动力学的解析解十分吻合，证明本方法能够有效地应用于自由基型接枝共聚合反应体系。     

An equation-free probabilistic steady-state approximation: dynamic application to the stochastic simulation of biochemical reaction networks

Stochastic chemical kinetics more accurately describes the dynamics of "small" chemical systems, such as biological cells. Many real systems contain dynamical stiffness, which causes the exact stochastic simulation algorithm or other kinetic Monte Carlo methods to spend the majority of their time executing frequently occurring reaction events. Previous methods have successfully applied a type of probabilistic steady-state approximation by deriving an evolution equation, such as the chemical master equation, for the relaxed fast dynamics and using the solution of that equation to determine the slow dynamics. However, because the solution of the chemical master equation is limited to small, carefully selected, or linear reaction networks, an alternate equation-free method would be highly useful. We present a probabilistic steady-state approximation that separates the time scales of an arbitrary reaction network, detects the convergence of a marginal distribution to a quasi-steady-state, directly samples the underlying distribution, and uses those samples to accurately predict the state of the system, including the effects of the slow dynamics, at future times. The numerical method produces an accurate solution of both the fast and slow reaction dynamics while, for stiff systems, reducing the computational time by orders of magnitude. The developed theory makes no approximations on the shape or form of the underlying steady-state distribution and only assumes that it is ergodic. We demonstrate the accuracy and efficiency of the method using multiple interesting examples, including a highly nonlinear protein-protein interaction network. The developed theory may be applied to any type of kinetic Monte Carlo simulation to more efficiently simulate dynamically stiff systems, including existing exact, approximate, or hybrid stochastic simulation techniques.     

Long-chain branching in free-radical polymerization due to chain transfer to polymer

The Monte Carlo sampling technique is used to investigate the branched structure formation during free-radical polymerization that involves chain transfer to polymer. This method accounts for the history of the generated branched structure and can provide virtually any structural information, because one can observe each polymer molecule directly. In this paper, we investigate the whole molecular weight distribution (MWD) for both pre- and postgelation periods, the MWDs for polymer molecules containing 0, 1, 2, 3, … branch points, the branching density of polymer molecules as functions of both size and the number of branch points, the spatial distribution of the branched chains at the theta state, etc. Contrary to the term ‘long-chain’ branching, many branch chains are relatively small, and the branched structures formed are significantly different from those usually depicted to introduce ‘branched polymers’ in many introductory textbooks. The radii of gyration at the theta state can be approximated by the Zimm-Stockmayer equation for random branching, in spite of various violations against the assumptions used in deriving the equation © 1995 John Wiley & Sons, Inc.     

D_4/N-β-氨乙基-γ-氨丙基甲基二甲氧基硅烷本体共聚反应聚合度变化的Monte Carlo模拟

用MonteCarlo方法模拟八甲基环四硅氧烷 (D4 )与N β 氨乙基 γ 氨丙基甲基二甲氧基硅烷(APAEDMS)的本体开环共聚动力学 .模拟过程采用自由体积理论简化处理扩散效应并与本征反应动力学耦合 .本征动力学常数通过模拟主要共聚基元反应得到 ,基于优化的动力学常数通过模拟从分子水平揭示D4 APAEDMS本体开环共聚反应过程主要表现为活性阴离子聚合行为 ,同时又伴有部分逐步聚合特征 .     

Path ensembles and path sampling in nonequilibrium stochastic systems

Markovian models based on the stochastic master equation are often encountered in single molecule dynamics, reaction networks, and nonequilibrium problems in chemistry, physics, and biology. An efficient and convenient method to simulate these systems is the kinetic Monte Carlo algorithm which generates continuous-time stochastic trajectories. We discuss an alternative simulation method based on sampling of stochastic paths. Utilizing known probabilities of stochastic paths, it is possible to apply Metropolis Monte Carlo in path space to generate a desired ensemble of stochastic paths. The method is a generalization of the path sampling idea to stochastic dynamics, and is especially suited for the analysis of rare paths which are not often produced in the standard kinetic Monte Carlo procedure. Two generic examples are presented to illustrate the methodology.     

Initialization of RAFT Agents with Different Leaving Groups – Determination of the Transfer Coefficients

The leaving group of the initial RAFT agent is one key factor during the initialization period of the RAFT polymerization. Benzyl-, Phenylethyl- and Cumyl Dithiobenzoate RAFT mediated bulk polymerization of styrene at 120 °C was investigated by direct measurement of the concentration of the RAFT agent and determination of the molecular weight in the initial phase of the RAFT polymerization using Size Exclusion Chromatography (SEC). The initialization period can be described by means of two transfer coefficients that are determined by the stability of the free-radical leaving group. By means of a Monte Carlo simulation the RAFT initialization process was characterized and the transfer coefficients from the experimental data were determined.     

O—酰基—α—酮肟光分解反应的Monte Carlo处理

本文首次对光化学反应体系用Monte Carlo方法进行模拟处理。通过5个O-酰基-α-酮肟光分解反应的Monte Carlo模拟,可避免解析解中由于对吸收光强须采用一级近似求解动力学微分方程组,而造成拟合反应在后期产生与实验结果的偏差。     

Simulation of styrene free radical polymerization over bi-functional initiators using Monte Carlo simulation method and comparison with mono-functional initiators

In the present study, based on a complete mechanism, a Monte Carlo simulation method is employed to investigate the kinetics of styrene free radical polymerization over bi-functional initiators in a bulk medium. The effects of the concentration of initiator and the monomer, of the temperature on monomer conversion, average molecular weights, polydispersity index, and molecular weight distribution are inspected and compared with mono-functional initiators. According to the simulation results, an increase in either the concentration of initiator or the temperature leads to the rise of the monomer conversion and to the reduction of the average molecular weights, while the increase of the monomer concentration results in the rise of both monomer conversion and molecular weights, which is in accord with predictions of the theory of free-radical polymerization. In addition, application of bi-functional initiators increases both monomer conversion and average molecular weight and results in narrower chain length distributions.     

Derivation of distributions of the degree of polymerization by probability theory

A method of determining distributions of the degree of polymerization by means of probability theory is derived. This method, based on Kolmogorov's forward differential equation, is illustrated by application to two familiar kinetic schemes and is used to derive the distribution of degree of polymerization for free-radical polymerization with diffusion-controlled termination.     

A simple approach to treat anisotropic elastic collisions in Monte Carlo calculations of the electron energy distribution function in cold plasmas

A simplified method allowing one to treat anisotropic electron heavy species elastic scattering in MonteCarlo models of gas discharges with the proper value for collision frequency is proposed The method is applied to an electric discharge in a Ne · Xe/HCl mixture, and the results are compared with the solution of the two-term expansion of the Boltzmann equation under the same conditions. Methods for reduction of computational time in Monte Carlo codes and the use of the Monte Carlo flux method are also discussed.