Molecular Weight Distribution Simulation in Equilibrium Ring‐Opening Polymerization: A New Macroscopic Model

Compared with other types of polymerization, the molecular weight and its distribution (MWD) of equilibrium ring‐opening polymerization (EROP) are complicated and have not been extensively studied. By using statistic method, a series of equations based on polymerization mechanism is established to describe numbers of rings, chains, and active centers. Using this new model, the predicted results of polydimethylsiloxane synthesized by EROP agree well with experimental results. This model has advantages in molecule number, calculation speed, and stability when compared with the Monte Carlo simulation. It also has the potential to replace Monte Carlo simulation in MWD prediction.     

Analytical calculus and Monte Carlo simulation of crosslinked polymer formation in the copolymerization of tetraethoxysilane and poly(dimethylsiloxane)

To understand the fundamental aspects of the polycondensation reaction of hydrolyzed tetraethoxysilane (TEOS) and silanol‐terminated poly(dimethylsiloxane) (PDMS), we modeled the reaction system as a step‐growth polymerization of A₄ and polydisperse A₂, assuming the reactivities of all functional groups are equal. The analytical solution for the weight‐average molecular weight is developed, and in addition, a Monte Carlo simulation is conducted to investigate the detailed structural development. It was found that as long as the molecular weight of PDMS is much larger than TEOS, the apparent behavior is significantly different from usual gelling systems. The gel point is relatively insensitive to the weight fraction of crosslinker (TEOS), the polydispersity index may decrease during polymerization before the rapid increase to infinity, and the molecular weight distribution profile may not show a significant broadening toward gelation. Even though the present model assumes a complete random reaction process among functional groups, formation of a heterogeneous structure in which a tight core consisting of TEOS‐based molecules is surrounded by soft PDMS chains was observed in the Monte Carlo simulation.     

Off‐lattice Monte Carlo simulation of hyperbranched polymers, 1 Polycondensation of AB2 type monomers

The polycondensation of hyperbranched polymers, based on AB₂ type monomers, was simulated using an off‐lattice Monte Carlo method in order to investigate the polymerization kinetics and microstructures of hyperbranched polymers. The effects of temperature and activation energy of reaction on the conversion rate seem qualitatively acceptable, indicating that our simulation model properly describes the polymerization reaction of hyperbranched polymers. Number average degree of polymerization and polydispersity index were calculated as a function of conversion and compared with Flory's expectation. Fractions of dendritic, linear and terminal unit were also determined from simulation and compared with theoretical predictions. As the hyperbranched polymer grows, it is observed that the molecular shape changes from a regular fan‐shape structure to an edge‐curled up structure.     

自由基聚合反应的Monte Carlo模拟方法

本文将化学反应动力学的MonteCarlo模拟方法运用到引发剂引发的自由基聚合反应的非稳态动力学,针对自由基聚合反应动力学数值模拟所特有的"无伸缩问题",采用"偏倚抽样法"解决了MonteCarlo模拟中的"无伸缩问题",模拟结果与非稳态动力学解的结果完全一致,此算法易推广到研究更复杂的自由基聚合反应体系。     

Numerical simulations of the polymerization reaction of deuterated diacetylene 2,4‐hexadiynylene bis(p‐toluenesulfonate) crystal

In this paper, we developed two types of programs in order to simulate the polymerization reaction of a fully deuterated crystal of diacetylene 2,4‐hexadiynylene bis(p‐toluenesulfonate) (pTS‐D). The first simulation is based on a modification of Baughman's model, a classical model for simulating the polymerization of diacetylene crystals. The agreement between the simulated and experimental results concerning the reaction kinetics is satisfactory. With this simulation algorithm, we take into account the experimental observation that the polymerization of pTS‐H and pTS‐D crystals is really a random process of formation of polymer chains along the crystallographic axis b . The second simulation is based on the Monte Carlo method, which permits not only to simulate the kinetics of the reaction, but also the chain‐length distribution in the hydrogenated and deuterated compounds. These two types of simulations were already developed for the hydrogenated crystal of diacetylene, named pTS‐H. Two main modifications are applied in the case of pTS‐D for taking into account experimental results: in the first the rate constants of chain‐terminating microscopic processes are different in pTS‐H and pTS‐D which must be considered. The second modification concerns the evolution of the lattice deformation during the course of polymerization. The experimental variation of the b parameter as a function of polymer content X in pTS‐D is different from that in pTS‐H; this result is important to consider when calculating the activation energy of the initiation and propagation microscopic processes.     

Bayesian Modeling and Markov Chain Monte Carlo Simulations for a Kinetic Study of Homo‐ and Co‐ Polymerization Systems

A Bayesian modeling and Markov Chain Monte Carlo simulation was developed for a kinetic study of homopolymerization and copolymerization systems at the molecular scale. Two copolymerization models – the terminal unit model and the penultimate unit model – were considered. Prior estimates of the kinetic parameters were obtained by L₁‐norm robust statistics. Using the structure of experimental data through a likelihood function, Bayesian modeling was employed to update the prior estimates. The joint posterior probability regions and shimmer bands were calculated for updated reactivity ratios. A method for assessing the power of experimental data in discrimination between copolymerization models is presented. This method was validated for free radical polymerization in binary systems. The evolution of species and radical populations during the course of polymerization were determined. The computational time was considerably decreased by calculating the propagation step from lifetime of the polymer chain and local monomer concentration. To avoid inaccuracies in the results caused by poor choice or false computation of the time step, the time step between successive Monte Carlo events was adapted to the time scale of the fastest reaction. The simulation algorithm is exact, in the sense that it takes full account of the fluctuations and correlations.

     

Monte‐Carlo simulation of secondary electron emission by x‐ray irradiation—an application of x‐ray absorption near‐edge structure (XANES)

A Monte‐Carlo simulation program has been developed for describing x‐ray absorption near‐edge structure (XANES) observed by synchrotron radiation. The Monte‐Carlo simulation was applied for interpreting XANES spectroscopy on a polycrystalline Ag specimen under synchrotron irradiation with photon energy 3340–3390 eV around the absorption edge of the Ag L_α line at 3352 eV. The results clearly indicate that Monte‐Carlo simulation describes the experimental results with considerable success. Dependence of secondary electron yield on the incident angle of synchrotron radiation was also studied. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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

Development of an Integrated Framework for Multiscale,Multiphase Modeling of Industrial Slurry‐Phase Reactors for Polyethylene Production

A framework based on the Monte Carlo/random‐pore polymeric flow model is proposed to simulate both single‐particle and continuous slurry reactor industrial polymerizations. The Sanchez–Lacombe equation of state describes the distributions of components in the different phases of these systems. The developed process model is applied to describe heterogeneously catalyzed polymerizations of ethylene in n‐hexane diluent with or without 1‐hexene as a comonomer, but the proposed methodology is applicable to any ethylene/1‐olefin copolymerization in slurry reactors. In addition to the effects of catalyst particle size and reactor residence time distributions, the proposed hybrid model is used to investigate the impact of several catalyst characteristics under different process conditions on polymer yield and microstructure. Particular attention is paid to the catalyst fragmentation process and active center distribution through the particle. These simulations demonstrate the versatility and thoroughness of combining Monte Carlo simulation with single‐particle models to analyze and predict the behavior of commercial polyolefin reactors.     

Flexible docking simulations: Scaled collective variable Monte Carlo minimization approach using Bezier splines,and comparison with a standard Monte Carlo algorithm

An algorithm for docking a flexible ligand onto a flexible or rigid receptor, using the scaled‐collective‐variables Monte Carlo with energy minimization approach, is presented. Energy minimization is shown to be one of the best techniques for distinguishing between native‐ and nonnative‐generated conformations. Incorporation of this technique into a Monte Carlo procedure enables one to distinguish the native conformation directly during the conformational search. It avoids the generation of a large number of ligand conformers for which more sophisticated energy evaluation tools would have had to be applied to identify the nativelike conformations. The efficiency of the Monte Carlo minimization was greatly improved by incorporating a new grid‐based energy evaluation technique using Bezier splines for which the energy function, as well as all of its derivatives, can be deduced from the values at the grid points. Comparison between our ECEPP/3‐based algorithm and the Monte Carlo algorithm presented elsewhere (Hart, T. N.; Read, R. J. Prot Struct Funct Genet 1992, 13, 206–222) has been made for docking NH₂ D Phe Pro Arg COOH, the noncovalent analog of NH₂ D Phe Pro Arg chloromethylketone (PPACK), onto the active site of human α‐thrombin. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 244–252, 1999     

Mathematical Modeling of the Microstructure of Poly(propylene) Made with Ziegler‐Natta Catalysts in the Presence of Electron Donors

Two different modeling techniques, the method of moments and Monte Carlo simulation, were compared for propylene polymerization with coordination catalysts including a new mechanistic step, site transformation by electron donors. We used the models to show how the molecular weight and tacticity distributions of several poly(propylene) chain populations were affected by changing the concentration of hydrogen, electron donor, and propylene in the reactor, under steady‐state or dynamic operating conditions. The Monte Carlo simulation describes the molecular weight and tacticity distributions for the whole polymer and chain populations with distinct microstructural characteristics. We have also applied the Monte Carlo model to simulate the pentad sequence distributions and its equivalent ¹³C NMR spectra.

     

Simulation of degradation processes. II. Testing some relations for random crosslinking without cyclization

By use of a Monte Carlo method the course of crosslinking without cyclization was simulated. For the simulated case the following was determined: number-average molecular weight and degree of polymerization of the system; weight-average molecular weight and degree of polymerization of the system; number of molecules of a given degree of polymerization; quantity of gel; number-average degree of polymerization of sol and branching. The values obtained by simulation were compared with some equations published to date. It was shown that the dependence of weight-average polymerization degree and, consequently, of the weight-average molecular weight on time and the dependence of the number of molecules of a given polymerization degree on time is probably more complicated than hitherto supposed.     

Polyethylene Made with Combinations of Single‐Site‐Type Catalysts: Monte Carlo Simulation of Long‐Chain Branch Formation

The formation of long‐chain branches (LCBs) during ethylene polymerization with a combination of catalysts was studied by Monte Carlo simulation. The model describes polymerization with a non‐branching catalyst that produces linear macromonomers, and a branching catalyst that produces linear and branched macromonomers. The LCBs are formed when the branching catalyst incorporates a macromonomer. The discussion is based on the three types of chain topology obtained during the synthesis: linear, comb‐branched, or hyperbranched. Simulation results show how the chain length distribution and the number of LCBs change according to the ratio between the two catalysts present in the reactor. The ratio hyperbranched/comb‐branched is defined to evaluate the system composition and the contribution of each catalyst.     

A Parallelised High Performance Monte Carlo Simulation Approach for Complex Polymerisation Kinetics

A novel, parallelised approach to Monte Carlo simulations for the computation of full molecular weight distributions (MWDs) arising from complex polymerisation reactions is presented. The parallel Monte Carlo method constitutes perhaps the most comprehensive route to the simulation of full MWDs of multiple chain length polymer entities and can also provide detailed microstructural information. New fundamental insights have been developed with regard to the Monte Carlo process in at least three key areas: (i) an insufficient system size is demonstrated to create inaccuracies via poor representation of the most improbable events and least numerous species; (ii) advanced algorithmic principles and compiler technology known to computer science have been used to provide speed improvements and (iii) the parallelisability of the algorithm has been explored and excellent scalability demonstrated. At present, the parallel Monte Carlo method presented herein compares very favourably in speed with the latest developments in the h‐p Galerkin method‐based PREDICI software package while providing significantly more detailed microstructural information. It seems viable to fuse parallel Monte Carlo methods with those based on the h‐p Galerkin methods to achieve an optimum of information depths for the modelling of complex macromolecular kinetics and the resulting microstructural information.

     

Chain Length Distributions of Polyolefins Made with Coordination Catalysts at Very Short Polymerization Times – Analytical Solution and Monte Carlo Simulation

We developed an analytical solution to describe how the CLD of polymers made with coordination polymerization catalysts vary as a function of time for very short polymerization times before the CLD becomes completely developed. We compared the analytical solution with a dynamic Monte Carlo model for validation, obtaining excellent agreement. Our analytical solution can be used to determine when the steady‐state hypothesis, commonly used in polymerization models, becomes valid as a function of polymer chain length. We also extended our model to describe polymerization with multiple‐site‐type catalysts. Depending on the polymerization kinetic parameters of the different site types on the catalyst, the fully developed CLD is reached through very different intermediate CLDs. This modeling approach, although rather simplified, can be used to interpret results from short polymerization time experiments such as the ones done in stopped‐flow reactors.

     

Multicanonical Monte Carlo calculation of the free‐energy map of the base–amino acid interaction

The specific interactions between base pairs and amino acids were studied by the multicanonical Monte Carlo method. We sampled numerous interaction configurations and side‐chain conformations of the amino acid by the multicanonical algorithm, and calculated the free energies of the interactions between an amino acid at given C_α positions and a fixed base pair. The contour maps of free energy derived from this calculation represent the preferred C_α position of the amino acid around the base, and these maps of various combinations of bases and amino acids can be used to quantify the specificity of intrinsic base–amino acid interactions. Similarly, enthalpy and entropy maps will provide further details of the specific interactions. We have also calculated the free‐energy map of the orientations of the C_α C_β bond vector, which indicates the preferential orientation of the amino acid against the base. We compared the results obtained by the multicanonical method with those of the exhaustive sampling and canonical Monte Carlo methods. The free‐energy map of the base–amino acid interaction obtained by the multicanonical simulation method was nearly identical to the accurate result derived from the exhaustive sampling method. This indicates that a single multicanonical Monte Carlo simulation can produce an accurate free‐energy map. Multicanonical Monte Carlo sampling produced free‐energy maps that were more accurate than those produced by canonical Monte Carlo sampling. Thus, the multicanonical Monte Carlo method can serve as a powerful tool for estimating the free‐energy landscape of base–amino acid interactions and for elucidating the mechanism by which amino acids of proteins recognize particular DNA base pairs. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 954–962, 2000     

GPGPU for orbital function evaluation with a new updating scheme

We have accelerated an ab initio quantum Monte Carlo electronic structure calculation using general purpose computing on graphical processing units (GPGPU). The part of the code causing the bottleneck for extended systems is replaced by Compute Unified Device Architecture‐GPGPU subroutine kernels which build up spline basis set expansions of electronic orbital functions at each Monte Carlo step. We have achieved a speedup of a factor of 30 for the bottleneck for a simulation of solid TiO₂ with 1536 electrons. To improve the performance with GPGPU we propose a new updating scheme for Monte Carlo sampling, quasi‐simultaneous updating, which is intermediate between configuration‐by‐configuration updating and the widely used particle‐by‐particle updating. The error in the energy due to by the single precision treatment and the new updating scheme is found to be within the required accuracy of ～10^?3 hartree per primitive cell. © 2012 Wiley Periodicals, Inc.     

Monte‐Carlo simulation of x‐ray photoelectron emission from silver

Silver 3d x‐ray photoelectron spectroscopy (XPS) spectra were simulated with the Monte‐Carlo method using an effective energy‐loss function that was derived from a reflected electron energy‐loss spectroscopy (REELS) analysis based on an extended Landau approach. After confirming that Monte‐Carlo simulation based on the use of the effective energy‐loss function can successfully describe the experimental REELS spectrum and Ag 3d XPS spectrum, we applied Monte‐Carlo simulation to predict the angular distribution of Ag 3d x‐ray photoelectrons for different x‐ray incidence angles and different photoelectron take‐off angles. The experimental photoelectron emission microscope that we are constructing was confirmed as being close to the optimum configuration, in which the x‐ray incident angle as measured from the surface normal direction is 74° and the photoelectron take‐off angle is set normal to the surface. The depth distribution functions of the Ag 3d X‐ray photoelectrons for different energy windows suggest that the photoelectron emission microscope will exhibit greater surface sensitivity for narrower photoelectron energy windows. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of the Intramolecular Chain Transfer to Polymer on PLP/SEC Experiments of Alkyl Acrylates

Pulsed‐laser polymerization (PLP) has been adopted by IUPAC as the method of choice for the determination of propagation rate constants (k_p). However, the method has failed in the polymerization of alkyl acrylates at temperatures above 30 °C. In this work, the PLP experiments were analyzed by simulation using a Monte Carlo algorithm. It was found that the experimental difficulties encountered to accurately determine k_p at temperatures above 30 °C were caused by extensive intramolecular chain transfer. This mechanism is not operative at lower temperatures because of its high activation energy.

Pulsed‐laser polymerization of BA in bulk at temperatures between −41 and +40 °C: Simulated MWD trace.     

Molecular weight distribution in free radical polymerization with long-chain branching

A new theory to predict the molecular weight distribution in free radical polymerization that includes chain transfer to polymer is proposed. This theory is based on the branching density distribution of the primary polymer molecules. The branching density distribution provides the information on how each chain is connected to other chains, and therefore, a full molecular weight distribution can be calculated by application of the Monte Carlo simulation. The present theory accounts for the history of the generated branched structure and can be applied to various reaction systems that involve branching and crosslinking regardless of the reactor types used. The present simulation confirmed the validity of the method of moments in a batch polymerization proposed earlier. It was shown clearly why gelation never occurs by chain transfer to polymer without the assistance of other interlinking reaction such as bimolecular termination by combination. © 1993 John Wiley & Sons, Inc.