Stanislaw Sosnowski

Theoretical analysis and computer simulations (Monte Carlo and numerical integration of differential equations) indicate that the statistical effect of a small number of reacting molecules is affected by transfer of reagent molecules between nanoreactors (droplets in this study). As a model reaction, a simple reversible association A + B ⇆ C was chosen, which was studied previously without reagent transfer processes. For sufficiently fast exchange of reactant molecules and a sufficiently high number of nanoreactors, the studied systems virtually do not differ from large volume systems if overall kinetics and thermodynamics of the chemical process are concerned. However, if either reagent molecule exchange is not fast or the number of exchanging nanoreactors is low, then the stochastic effect is clearly visible, influencing the kinetics of reaching reaction equilibrium. In systems with a low number of nanoreactors, the apparent (average) equilibrium constant is affected as well. The distribution of reactant molecules in the nanoreactors is governed by stochastic processes, dependent on stochastic rate constants of all processes, chemical as well as physical (transfer of molecules outside and into droplets). When accumulation of reactant molecules in the continuous phase cannot be neglected, then the partition coefficients of reactants between the continuous and dispersed phases have to be taken into account. Distributions of reactant molecules described in the paper for systems composed of few nanoreactors can be especially important for some biochemical processes in living cells or devised corresponding artificial reactors. If the reactant molecules predominate in a continuous phase, then the kinetics and overall conversion of reversible association may be significantly affected by the presence of the dispersed phase and its catalytic behavior can be observed. It has been shown that Monte Carlo simulations, applying a devised simple algorithm, give reliable results, allowing theoretical studies of nanoscale-droplet systems.

A process of irreversible copolymerization with segmental exchange reactions is simulated using a Monte Carlo algorithm. The terminal (dyad) model is assumed for eight considered reactions: propagations and attacks of active centers on monomeric units in chains with subsequent chain scissions. It is established that microstructure of macromolecules at any given moment depends on the values of reaction rate constants and evolves from the beginning one, defined mainly by propagation rates, to the final one, determined by chain transfer rates. Dispersity and randomization degree of unit sequences in macromolecules can change nonmonotonically with progress of the process. Equilibration of the microstructure needs usually significantly longer time than needed for attaining the equilibrium chain lengths distribution. The proposed new Monte Carlo algorithm allows to more efficiently investigate the studied process.

The effect of configuration of an asymmetric bulky initiator 2,2′-[1,1′-binaphtyl-2,2′-diyl- bis-(nitrylomethilidyne)]diphenoxy aluminum isopropoxide (Ini) on structure of copolymer of asymmetric monomer l,l-lactide (Lac) with symmetric comonomer trimethylene carbonate (Tmc) was studied using polarimetry, dilatometry, Size Exclusion Chromatography (SEC), and Carbon Nuclear Magnetic Resonance (¹³C NMR). When the S-enantiomer of Ini was used the distribution in copolymer chains at the beginning of polymerization is statistical, with alternacy tendency, changing next through a gradient region to homoblocks of Tmc. However, when R-Ini was used, the product formed was a gradient oligoblock one, with Tmc blocks prevailing at the beginning, changing to Lac blocks dominating at the end part of chains. Initiation of copolymerization with the mixture of both initiator enantiomers (S:R = 6:94) gave a multiblock copolymer of similar features but shorter blocks. Analysis of copolymerization progress required complex analysis of dilatometric data, assuming different molar volume contraction coefficients for units located in different triads. Comonomer reactivity ratios of studied copolymerizations were determined.