Coarse‐grained molecular dynamics simulations of stereoregular poly(methyl methacrylate)/poly(vinyl chloride) blends

The stereoregular poly(methyl methacrylate)/poly(vinyl chloride) blends with a wide formulation range are extensively simulated using the coarse‐grained (CG) molecular dynamics (MD) method. To improve the representability, the bonded CG potentials are re‐parameterized against the atomistic simulated melt systems whereas the nonbonded CG potentials are adopted as developed in our previous work. Based on the CG potentials, the MD simulations reproduce all the local distributions of pure systems and the miscibility of mixed systems. Moreover, the global conformational properties are also closer to the target ones than those obtained using the previous CG potentials. The changes in density and volume upon mixing are computed together with the energies of mixing. They are all negative over the entire composition range and indicate stronger intermolecular interactions between distinct components than those between identical components. In particular, it is found that upon mixing the changes in density are insensible to chain tacticity but the changes in volume and the energies of mixing do, which quantitatively confirms that both inter‐molecular interactions and free‐volumes mainly contribute to the observed phase behaviors. Such models and methods reported herein can be used to quickly optimize formulations of polymer blends. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 203–212     

Molecular dynamics simulations of atomistic hydration structures of poly(vinyl methyl ether)

Molecular dynamics simulations have been performed on the aqueous solutions of poly(vinyl methyl ether) (PVME) at various concentrations. Both radial and spatial distribution functions are used to investigate the detailed hydration structures. The structures of water are found to get increasingly concentrated when polymers are introduced and the water motions are severely hindered by the polymer matrix. At low concentrations, larger populations of tt conformers in meso dyads than those at higher concentrationsare found and this phenomenon is believed to be due to the increasing in bonding of water molecule to two ether oxygens in meso dyad. At higher concentrations, the size and conformations of polymers are quite similar to those in bulk. A transition of hydrogen bond fractions between PVME and water at around the concentration of 0.3 is observed and this value is perfectly in agreement with the results of conformational analysis and Raman spectra. Second neighbor hydrogen bond statistics revealed the domination of complicated hydrogen bond networks at low concentrations, but single hydrogen bonds as well as isolated clusters composed of 2-4 water molecules are usual around each polymer repeat unit.     

Comparison of selected polarizable and nonpolarizable water models in molecular dynamics simulations of ice I(h)

We present a molecular dynamics simulation study in which we determined the melting point of ice I(h) for the polarizable SWM4-NDP water model (Lamoureux et al., Chem. Phys. Lett., 2006, 418, 245-249) and compared the performance of several popular water force fields, both polarizable and nonpolarizable, in terms of melting temperature, stability and orientational structuring of ice. The simulations yield the melting temperature of SWM4-NDP ice as low as T(m) = 185 ± 10 K, despite the quadrupole moment of a SWM4-NDP water molecule being close to the experimental gas phase value. The results thus show that the dependence of T(m) on the molecular quadrupole, observed for the three- and four-site water models, is generally lost if polarization is explicitly included. The study also shows that adding polarizability to a planar three-charge water model increases orientational disorder in hexagonal ice. In addition, analysis of the tetrahedral order in bulk ice reveals a correlation between the pre-existing degree of orientational disorder in ice simulated using different polarizable and nonpolarizable models and the melting temperature of the models. Our findings thus suggest some new considerations regarding the role of polarization forces in a crystalline solid that may guide future development of reliable polarizable water models for ice.     

Free energy of a trans-membrane pore calculated from atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations of a lipid bilayer were performed to calculate the free energy of a trans-membrane pore as a function of its radius. The free energy was calculated as a function of a reaction coordinate using a potential of mean constraint force. The pore radius was then calculated from the reaction coordinate using Monte Carlo particle insertions. The main characteristics of the free energy that comes out of the simulations are a quadratic shape for a radius less than about 0.3 nm, a linear shape for larger radii than this, and a rather abrupt change without local minima or maxima between the two regions. In the outer region, a line tension can be calculated, which is consistent with the experimentally measured values. Further, this line tension can be rationalized and understood in terms of the energetic cost for deforming a part of the lipid bilayer into a hydrophilic pore. The region with small radii can be described and understood in terms of statistical mechanics of density fluctuations. In the region of crossover between a quadratic and linear free energy there was some hysteresis associated with filling and evacuation of the pore with water. The metastable prepore state hypothesized to interpret the experiments was not observed in this region.     

Electric field poled polymeric nonlinear optical systems: molecular dynamics simulations of poly(methyl methacrylate) doped with disperse red chromophores

We demonstrate a complete procedure for simulations of electric field poled polymeric nonlinear optical systems with the purpose to evaluate the macroscopic electro-optic coefficients. The simulations cover the electric field poling effects on the chromophore order at the liquid state, the cooling procedure from the liquid to the solid state in the presence of the poling field, and the back-relaxation of the system after the removal of the field. We use Disperse Red chromophore molecules doped in a poly(methyl methacrylate) matrix for a numerical demonstration of the total procedure. On the basis of the simulation results, the polymer mobility and the static properties of the dopant chromophores are derived. In the liquid state, the chromophore molecules are closer to the side chains than to the backbones of the polymer matrix, and after the simulated annealing, the polymer matrix tends to be closely packed, leading to a significant change in the polymer structure around the chromophore molecules. Besides predicting the absolute macroscopic electro-optic coefficient values, the results are used to derive the microscopic origin of these values in terms of geometric and electronic structure, loading, poling, and back-relaxation effects, thereby aiding to establish design principles for optimum guest-host configurations.     

Polymerization of monolayers. V. Tacticity of the insertion poly(methyl methacrylate)

Insertion poly(methyl methacrylate) (PMMA) formed within a monolayer of monomer adsorbed on montmorillonite was studied by means of NMR spectroscopy and shown to be composed of short stereosequences with a predominant isotactic component. The stereospecificity of the insertion PMMA can be understood in terms of the organization of the monomer adsorbed on the surface of the mineral. The monolayer of monomer is assumed to be composed of monads and isotactic diads due to dipole-ion interactions between the molecules of methyl methacrylate and exchangeable cations. Results of experiments relating factors such as temperature, density of ion population of the surface, and the nature of ions are discussed. It is shown that the populations of isotactic and syndiotactic triads P_i and P_s determined experimentally in the insertion PMMA are in a reasonable agreement with values calculated from the model. It is shown, furthermore, that P_i increases with increase in the density of the ion population, in agreement with the equations derived from the model.     

Mechanical properties of surfactant bilayer membranes from atomistic and coarse-grained molecular dynamics simulations

We use simulations to predict the stability and mechanical properties of two amphiphilic bilayer membranes. We carry out atomistic MD simulations and investigate whether it is possible to use an existing coarse-grained (CG) surfactant model to map the membrane properties. We find that certain membranes can be represented well by the CG model, whereas others cannot. Atomistic MD simulations of the erucate membrane yield a headgroup area per surfactant a(0) of 0.26 nm(2), an elastic modulus K(A) of 1.7 N/m, and a bending rigidity kappa of 5 k(B)T. We find that the CG model, with the right choice for the size and potential well depth of the head, correctly reproduces a(0), kappa, as well as the fluctuation spectrum over the whole range of q values. Atomistic MD simulations of EHAC, on the other hand, suggest that this membrane is unstable. This is indicated by the fact that kappa is of the order of k(B)T, which means that the interface is extremely flexible and diffuse, and K(A) is close to zero, which means that the surface tension is zero. We argue that the CG model can be used if the headgroups are uncharged, dipolar, or effectively dipolar due to headgroup charge screening induced by counterion condensation.     

Polymerization of methacrylates. I. Polymerization reactivity of picryl methacrylate

The polymerization of picryl (PMA), 2,4-dinitrophenyl (2,4-DNMA),2,6-dinitrophenyl (2,6-DNMA), 2-methyl-4,6-dinitrophenyl (MDNMA), and 2,6-dimethylphenyl methacrylates (DMMA) was carried out in benzene at 60°C. PMA, 2,6-DNMA, and MDNMA did not undergo radical homopolymerization, while 2,4-DNMA and DMMA did. The results suggest that the growing radical readily attacks the oxygen atom of the nitro group at the 2 position of the terminal phenyl group due to the steric effect of the substituent at the 6 position, resulting in chain termination. PMA formed a charge-transfer complex with 2-naphthyl methacrylate (NMA). The stoichiometric composition of this complex was shown to be 1:1 molar complex. PMA was readily copolymerized with NMA. The amount of solvent affected the composition of the copolymer obtained at a given same mole fraction in feed. The results suggest that charge-transfer interaction between the ester groups affects the copolymerization mechanism.     

Interpenetrating polymer networks of poly(carbonate-urethane) and polymethyl methacrylate

We report the synthesis and characterization of interpentrating polymer networks (IPN's), pseudo INPN's, and blends of urethane-containing aliphatic polycarbonates (PCU) and polymethyl methacrylate (PMMA). The simultaneous full IPN's of PCU's and PMMA over the whole composition range have only one T_g, as determined by DSC and DMA which, together with transmission electron microscopy observations, suggest a single phase morphology even though the linear chain blends are completely immiscible. The full IPN's exhibit a maximum tensile strength at a certain composition, and superior solvent and heat resistance as compared to the pseudo IPN's, linear blends, and the pure crosslinked PCU's and PMMA.     

Polyepichlorohydrin polyurethane/poly(methyl methacrylate) interpenetrating polymer networks

Polyurethane (PU) based on polyepichlorohydrin/poly(methyl methacrylate) (PECH/PMMA) interpenetrating polymer networks (IPNs) was synthesized by a simultaneous method. The effects of composition, hydroxyl group number of PECH, NCO/OH ratio and crosslinking agent content in IPNs were investigated in detail. Some other glycols, such as poly(ethylene glycol), poly(propylene glycol) and hydroxyl-terminated polybutadiene, were also used to obtain PU/PMMA IPNs. The interpenetrating and fracture behaviors of the IPNs are explained briefly.     

Polymerization in a magnetic field. X. Solvent effect in poly(methyl methacrylate) synthesis

The effect of a continuous external magnetic field on the free radical polymerization of methyl methacrylate was studied with respect to some properties of the solvents used in the reaction. The studies were performed through dilatometric technique in and out a magnetic field of 0.25 T. Ten different solvents were used to underline the dependence between the magnetic field presence, the reaction medium, and the development of the polymerization process. The intervened magnetokinetic effects are attributed to the changes in the multiplicity of the radical pairs owing to the magnetic field influence. There is an interdependence among the viscosity and molar polarization of the solvents and the magnetic field effect. © 1996 John Wiley & Sons, Inc.     

Molecular dynamics simulations of a poly(p-phenylene) oligomer

We have studied the conformation and coefficient of thermal expansion in the poly(p-phenylene) oligomer p-sexiphenyl (C₃₆H₂₆) by molecular dynamics simulations. Studies of the backbone phenyl–phenyl torsion angle in a simulated p-sexiphenyl crystal at room temperature indicate the presence of torsional librations of approximately ±20°. Further analysis of the phenyl–phenyl backbone torsion angle in less closely packed regions of the simulated crystal (crystal ends) indicate the presence of 180° phenyl ring flips, in agreement with solid-state deuterium NMR data on poly(p-phenylene oligomers). The linear coefficient of thermal expansion was also calculated and found to be negative, in qualitative agreement with experimental data on rigid-rod compounds. © 1993 John Wiley & Sons, Inc.     

Thermodynamic characterization of interpenetrating polymer networks of poly (carbonate-urethane) and poly (methyl methacrylate)

Crosslinked poly (methyl methacrylate) and poly (carbonate-urethane) as well as full interpenetrating network on their base (IPN) were characterized by precise heat capacity measurements in the temperature interval 1.2–150 K. As judged by the positive sign of the excess Gibbs free energy in the whole temperature interval above 30 K, the apparently single-phase state of IPN is thermodynamically unstable; however, its kinetic stability is ensured by permanent chamical crosslinks prohibiting the incipient phase separation. © 1994 John Wiley & Sons, Inc.     

Fracture surface morphology and phase relationships of polystyrene/poly(methyl methacrylate) systems. I. Low-molecular-weight polystyrene in poly(methyl methacrylate)

Samples of low-molecular-weight polystyrene (PS) in poly(methyl methacrylate) (PMMA) were prepared by first dissolving PS in methyl methacrylate monomer and then polymerizing the monomer. Forty-three specimens of varying number-average molecular weight (2100–49,000) and composition (5–40 wt %) of PS were prepared, and the surface morphology and phase relationships studied by scanning electron microscopy. Four distinct types of phase relationships were observed: (i) a single phase consisting of PS dissolved in PMMA; (ii) PS dispersed in PMMA; (iii) PMMA dispersed in PS; and (iv) regions of PS dispersed in PMMA coexisting with regions of PMMA dispersed in PS. Values of the size and population density of the dispersed particles are reported. Finally, the size and distribution of the dispersed particles and the various types of phase relationships are discussed in terms of the ternary polystyrene/poly(methyl methacrylate)/methyl methacrylate phase diagram.     

Degradation of oligo(lactone) branches linked to poly(methacrylate) networks

Partially resorbable composite materials were developed for bone applications and were prepared by crosslinking copolymerization of oligo(lactone) macromonomers with selected comonomers in the presence of hydroxy apatite. The composites were incubated in aqueous solution of different pH and electrolyte content over a period up to 300 days. The percentage of the released hydroxycarboxylic acid from the oligo(lactone) branches was determined by the titration or by high performance liquid chromatography (HPLC). A catalytic influence of cholesterol esterase (CE) on the hydrolysis of the materials was observed. Determination of CE activity in the presence of the different composite materials showed that the half-life time of CE was the lowest in the presence of the material undergoing the fastest degradation. After degradation the composite materials have preserved their coherence and have shown higher glass transition temperatures in the dried state than the materials before degradation.     

Minimum free energy pathways and free energy profiles for conformational transitions based on atomistic molecular dynamics simulations

An efficient method for the calculation of minimum free energy pathways and free energy profiles for conformational transitions is presented. Short restricted perturbation-targeted molecular dynamics trajectories are used to generate an approximate free energy surface. Approximate reaction pathways for the conformational change are constructed from one-dimensional line segments on this surface using a Monte Carlo optimization. Accurate free energy profiles are then determined along the pathways by means of one-dimensional adaptive umbrella sampling simulations. The method is illustrated by its application to the alanine "dipeptide." Due to the low computational cost and memory demands, the method is expected to be useful for the treatment of large biomolecular systems.     

Amphiphilic networks. VII. Synthesis and characterization of pH-sensitive poly(sulfoethyl methacrylate)-1-polyisobutylene networks

A series of environmentally sensitive amphiphilic networks consisting of 2-sulfoethyl methacrylate (SEMA) chains linked by methacrylate-ditelechelic polyisobutylene (MA–PIB–MA) chains have been prepared and characterized. Network composition was determined after sequential solvent extraction by elemental analysis. These networks are two-phase microheterogeneous systems containing hydrophobic rubbery PIB domains (T_g ～ ?60°C) and hydrophilic poly(2-sulfoethyl methacrylate) domains (T_g ～ ?15°C). They exhibit large contact-angle hysteresis in water which is due to surface segmental mobility and microheterogeneity. By increasing the SEMA content of the networks the contact-angle hysteresis increases. This phenomenon is due to an increase in the advancing contact angle most likely caused by the migration of the nonpolar PIB domains toward the surface and concomitant decrease of the receding contact angle. These amphiphilic networks exhibit non-Fickian swelling in n-heptane, as well as in water, and show pH-sensitive swelling in aqueous media. They rapidly and reversibly swell and deswell in response to increasing or decreasing the pH of the media (cycling between pH = 2 and 12). © 1994 John Wiley & Sons, Inc.     

Determination of glass transition temperature of polyimides from atomistic molecular dynamics simulations and machine-learning algorithms

Glass transition temperature (T_g) plays an important role in controlling the mechanical and thermal properties of a polymer. Polyimides as an important category of engineering plastics have wide applications because of their superior heat resistance and mechanical strength. The capability of predicting T_g for a polyimide a priori is therefore highly desirable in order to expedite the design and discovery of new polyimide polymers with targeted properties and applications. Here we explore three different approaches to either compute T_g for a polyimide via all-atom molecular dynamics simulations or predict T_g via a mathematical model generated by using machine-learning algorithms to analyze existing data collected from the literature. Our simulations reveal that T_g can be determined from examining the diffusion coefficient of simple gas molecules in a polyimide as a function of temperature and the results are comparable to those derived from data on polymer density versus temperature and actually closer to the available experimental data. Furthermore, the predictive model of T_g derived with machine-learning algorithms can be used to estimate T_g successfully within an uncertainty of about 20 degrees, even for polyimides yet to be synthesized experimentally.     

Extracting parameters for base-pair level models of DNA from molecular dynamics simulations

 A method is described to extract a complete set of sequence-dependent energy parameters for a rigid base-pair model of DNA from molecular dynamics (MD) simulations. The method is properly consistent with equilibrium statistical mechanics and leads to effective inertia parameters for the base-pair units as well as stacking and stiffness parameters for the base-pair junctions. We give explicit formulas that yield a complete set of base-pair model parameters in terms of equilibrium averages that can be estimated from a time series generated in an MD simulation. The expressions to be averaged depend strongly both on the choice of coordinates used to describe rigid-body orientations and on the choice of strain measures at each junction. Received: 12 July 2000 / Accepted: 5 January 2001 / Published online: 3 May 2001     

Interpenetrating polymer networks of poly(dimethyl siloxane–urethane) and poly(methyl methacrylate)

Simultaneous IPNs of poly(dimethyl siloxane-urethane) (PDMSU)/poly(methyl methacrylate) (PMMA) and related isomers have been prepared by using new oligomers of bis(β-hydroxyethoxymethyl)poly(dimethyl siloxane)s (PDMS diols) and new crosslinkers biuret triisocyanate (BTI) and tris(β-hydroxylethoxymethyl dimethylsiloxy) phenylsilane (Si-triol). Their phase morphology have been characterized by DSC and SEM. The SEM phase domain size is decreased by increasing crosslink density of the PDMSU network. A single phase IPN of PDMSU/PMMA can be made at an M_c = 1000 and 80 wt % of PDMSU. All of the pseudo- or semi-IPNs and blends of PDMSU and PMMA were phase separated with phase domain sizes ranging from 0.2 to several micrometers. The full IPNs of PDMSU/PMMA have better thermal resistance compared to the blends of linear PDMSU and linear PMMA. © 1993 John Wiley & Sons, Inc.