Local Dynamics of Syndiotactic Poly(methyl methacrylate) Using Molecular Dynamics Simulation

Confinement Effects in Dynamics of Ionic Liquids with Polymer-Grafted Nanoparticles

Ionic liquid mixed with poly(methyl methacrylate)-grafted nanoparticle aggregates at low particle concentrations was shown to exhibit different dynamics and ionic conductivity than that of pure ionic liquid in our previous studies. In this work, we report on the quasi-elastic neutron scattering results on ionic liquid containing polymer-grafted nanoparticles at the higher particle concentration. The diffusivity of imidazolium (HMIM + ) cations of 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (HMIM-TFSI) in the presence of poly(methyl methacrylate)-grafted iron oxide nanoparticles and the ionic conductivity of solutions were discussed through the confinement. Analysis of the elastic incoherent structure factor suggested the confinement radius decreased with the addition of grafted particles in HMIM-TFSI/solvent mixture, indicating the confinement that is induced by the high concentration of grafted particles, shrinks the HMIM-TFSI restricted volume. We further conjecture that this enhanced diffusivity occurs as a result of the local ordering of cations within aggregates of poly(methyl methacrylate)-grafted particles.

Disentangling Self-Atomic Motions in Polyisobutylene by Molecular Dynamics Simulations

We present fully atomistic molecular dynamics simulations on polyisobutylene (PIB) in a wide temperature range above the glass transition. The cell is validated by direct comparison of magnitudes computed from the simulation and measured by neutron scattering on protonated samples reported in previous works. Once the reliability of the simulation is assured, we exploit the information in the atomic trajectories to characterize the dynamics of the different kinds of atoms in PIB. All of them, including main-chain carbons, show a crossover from Gaussian to non-Gaussian behavior in the intermediate scattering function that can be described in terms of the anomalous jump diffusion model. The full characterization of the methyl-group hydrogen motions requires accounting for rotational motions. We show that the usually assumed statistically independence of rotational and segmental motions fails in this case. We apply the rotational rate distribution model to correlation functions calculated for the relative positions of methyl-group hydrogens with respect to the carbon atom at which they are linked. The contributions to the vibrational density of states are also discussed. We conclude that methyl-group rotations are coupled with the main-chain dynamics. Finally, we revise in the light of the simulations the hypothesis and conclusions made in previously reported neutron scattering investigations on protonated samples trying to address the origin of the dielectric β-process.

Molecular Dynamics Study of the Swelling of Poly(methyl methacrylate) in Supercritical Carbon Dioxide

The swelling of a poly (methyl methacrylate) in supercritical carbon dioxide was studied by means of full atomistic classical molecular dynamics simulation. In order to characterize the polymer swelling, we calculated various properties related to the density, structure, and dynamics of polymer chains as a function of the simulation time, temperature, and pressure. In addition, we compared the properties of the macromolecular chains in supercritical CO₂ with the properties of the corresponding bulk system at the same temperature and atmospheric pressure. It was shown that diffusion of CO₂ molecules into the polymer led to a significant increase in the chain mobility and distances between them. Analysis of diffusion coefficients of CO₂ molecules inside and outside the poly(methyl methacrylate) sample has shown that carbon dioxide actively interacts with the functional groups of poly (methyl methacrylate). Joint analysis of the radial distribution functions obtained from classical molecular dynamics and of the averaging interatomic distances from Car-Parrinello molecular dynamics allows us to make a conclusion about the possibility of formation of weak hydrogen bonds between the carbon dioxide oxygen atom and the hydrogen atoms of the polymer methyl groups.

Insights into the structural dynamics of poly lactic-co-glycolic acid at terahertz frequencies

The mechanical properties of an amorphous copolymer are directly related to the dynamic processes occurring at the molecular level.

Structure and dynamics of stereo-regular poly(methyl-methacrylate) melts through atomistic molecular dynamics simulations

Poly(methyl-methacrylate), PMMA, is a disubstituted vinyl polymer whose properties depend significantly on its tacticity. Here we present a detailed study of the structure, conformation, and dynamics of syndiotactic, atactic, and isotactic PMMA melts at various temperatures (580, 550, 520, and 490 K) via all-atom molecular dynamics simulations. The calculated volumetric properties are close to experimental data. The T_g and chain dimensions of PMMA model systems are found to depend strongly on tacticity in agreement with experimental findings. The backbone bonds in trans (t), diads in tt, and inter-diads in t|t torsional states are the most populated for all PMMA stereo-chemistries and their fractions increase with the number of syndiotactic sequences. Also, the effective torsional barrier heights for the backbone, ester side group, and α-methyl group are larger for syndiotactic PMMA compared to the isotactic one. The structure of the PMMA chains is studied by computing the intra- and inter-chain static structure factors, S(q), and the radial pair distribution functions. In the first peak of S(q), both intra- and inter-chain components contribute, whereas the second and third peaks mainly come from inter- and intra-chain parts, respectively. For all PMMA stereo-isomers a clear tendency of ester-methyl groups to aggregate is observed. The local dynamics are studied by analyzing torsional autocorrelation functions for various dihedral angles. A wide spectrum of correlation times and different activation energies are observed for the motions of different parts of PMMA chains. The stereo-chemistry affects the backbone, ester side group, and α-methyl motions, whereas the ester-methyl rotation remains unaffected. The dynamic heterogeneity of the PMMA chains is also studied in detail for the different stereo-chemistries via the temperature dependence of the stretching exponent. Furthermore, the reorientational dynamics at the chain level and translational dynamics for monomer and chain centers-of-mass are analyzed.

Local Structure Contributions to Surface Tension of a Stereoregular Polymer

We have used all-atom molecular dynamics (MD) simulations to calculate the surface tension of melt poly(methyl methacrylate) (PMMA) as a function of tacticity. Computation of surface tension using the Kirkwood-Buff approach required hundreds of nanoseconds of equilibration. The computed slopes of surface tension versus temperature are in very good agreement with reported experimental values. Using a rigorous treatment of the true interface, which takes into account the molecular roughness, we find that isotactic PMMA, in comparison to syndiotactic and atactic PMMA, shows a larger surface concentration of polar ester-methyl and carbonyl groups on the surface versus nonpolar α-methyl groups. A mechanistic hypothesis based on the helical nature of the isotactic PMMA chains, their relative flexibility, and their reported conformational energies is proposed to explain the trends in composition near the surface. We highlight here how surface composition and surface tension are controlled by both polarity and steric constraints imposed by tacticity.

Interaction potential models for bulk ZnS, ZnS nanoparticle, and ZnS nanoparticle-PMMA from first-principles

An ab initio derived transferable polarizable force-field has been developed for Zinc sulphide (ZnS) nanoparticle (NP) and ZnS NP-PMMA nanocomposite. The structure and elastic constants of bulk ZnS using the new force-field are within a few percent of experimental observables. The new force-field show remarkable ability to reproduce structures and nucleation energies of nanoclusters (Zn1S1-Zn12S12) as validated with that of the density functional theory calculations. A qualitative agreement of the radial distribution functions of Zn-O, in a ZnS nanocluster-PMMA system, obtained using molecular mechanics molecular dynamics (MD) and ab initio MD (AIMD) simulations indicates that the ZnS-PMMA interaction through Zn-O bonding is explained satisfactorily by our force-field.

Effects of nanoscopic-confinement on polymer dynamics

The static and dynamic behavior of polymers in confinement close to interfaces can be very different from that in the bulk. Among the various geometries, intercalated nanocomposites, in which polymer films of ∼1 nm thickness reside between the parallel inorganic surfaces of layered silicates in a well-ordered multilayer, offer a unique avenue for the investigation of the effects of nanoconfinement on polymer structure and dynamics by utilizing conventional analytical techniques and macroscopic specimens. In this article, we provide a review of research activities mainly in our laboratory on polymer dynamics under severe confinement utilizing different polymer systems: polar and non-polar polymers were mixed with hydrophilic or organophilic silicates, respectively, whereas hyperbranched polymers were studied in an attempt to probe the effect of polymer-surface interactions by altering the number and the kinds of functional groups in the periphery of the branched polymers. The polymer dynamics was probed by quasielastic neutron scattering and dielectric relaxation spectroscopy and was compared with that of the polymers in the bulk. In all cases, very local sub-Tg processes related to the motion of side and/or end groups as well as the segmental α-relaxation were identified with distinct differences recorded between the bulk and the confined systems. Confinement was found not to affect the very local motion in the case of the linear chains whereas it made it easier for hyperbranched polymers due to modifications of the hydrogen bond network. The segmental relaxation in confinement becomes faster than that in the bulk, exhibits Arrhenius temperature dependence and is observed even below the bulk Tg due to reduced cooperativity in the confined systems.

Generation of Well-Relaxed All-Atom Models of Large Molecular Weight Polymer Melts: A Hybrid Particle-Continuum Approach Based on Particle-Field Molecular Dynamics Simulations

A procedure based on Molecular Dynamics (MD) simulations employing soft potentials derived from self-consistent field (SCF) theory (named MD-SCF) able to generate well-relaxed all-atom structures of polymer melts is proposed. All-atom structures having structural correlations indistinguishable from ones obtained by long MD relaxations have been obtained for poly(methyl methacrylate) (PMMA) and poly(ethylene oxide) (PEO) melts. The proposed procedure leads to computational costs mainly related on system size rather than to the chain length. Several advantages of the proposed procedure over current coarse-graining/reverse mapping strategies are apparent. No parametrization is needed to generate relaxed structures of different polymers at different scales or resolutions. There is no need for special algorithms or back-mapping schemes to change the resolution of the models. This characteristic makes the procedure general and its extension to other polymer architectures straightforward. A similar procedure can be easily extended to the generation of all-atom structures of block copolymer melts and polymer nanocomposites.

Molecular structure of poly(methyl methacrylate) surface II: Effect of stereoregularity examined through all-atom molecular dynamics

Utilizing all-atom molecular dynamics (MD), we have analyzed the effect of tacticity and temperature on the surface structure of poly(methyl methacrylate) (PMMA) at the polymer-vacuum interface. We quantify these effects primarily through orientation, measured as the tilt with respect to the surface normal, and the surface number densities of the α-methyl, ester-methyl, carbonyl, and backbone methylene groups. Molecular structure on the surface is a complex interplay between orientation and number densities and is challenging to capture through sum frequency generation (SFG) spectroscopy alone. Independent quantification of the number density and orientation of chemical groups through all-atom MD presents a comprehensive model of stereoregular PMMA on the surface. SFG analysis presented in part I of this joint publication measures the orientation of molecules that are in agreement with MD results. We observe the ester-methyl groups as preferentially oriented, irrespective of tacticity, followed by the α-methyl and carbonyl groups. SFG spectroscopy also points to ester-methyl being dominant on the surface. The backbone methylene groups show a very broad angular distribution, centered along the surface plane. The surface number density ratios of ester-methyl to α-methyl groups show syndiotactic PMMA having the lowest value. Isotactic PMMA has the highest ratios of ester- to α-methyl. These subtle trends in the relative angular orientation and number densities that influence the variation of surface structure with tacticity are highlighted in this article. A more planar conformation of the syndiotactic PMMA along the surface (x-y plane) can be visualized through the trajectories from all-atom MD. Results from conformation tensor calculations for chains with any of their segments contributing to the surface validate the visual observation.