Detailed molecular dynamics simulation of the self-diffusion of n-alkane and cis-1,4 polyisoprene oligomer melts

The effect of chemical constitution on polyisoprene dynamics

Polyisoprene (PI) melts have been studied, with most reports focusing on systems with high 1,4-cis content. In contrast, 1,4-trans PI homopolymers or random copolymers have seldom been examined, despite a handful of investigations suggesting a distinct dynamic behavior. Herein, we employ all-atom simulations to investigate the effect of chemical architecture on the dynamics of cis and trans-PI homopolymers, as well as copolymers. We examine the thermodynamic, conformational, and structural properties of the polymers and validate the performance of the models. We probe chain dynamics, revealing that cis-PI presents accelerated translation and reorientation modes relative to trans as recorded by the mean square displacement of the chain center-of-mass as well as by the characteristic times of the lower modes in a Rouse analysis. Interestingly, progressing to higher modes, we observe a reversal with trans units exhibiting faster dynamics. This was further confirmed by calculations of local carbon-hydrogen vector reorientation dynamics, which offer a microscopic view of segmental mobility. To obtain insight into the simulation trajectories, we evaluate the intermediate incoherent scattering function that supports a temperature-dependent crossover in relative mobility that extends over separations beyond the Kuhn-length level. Finally, we analyzed the role of non-Gaussian displacements, which demonstrate that cis-PI exhibits increased heterogeneity in dynamics over short-timescales in contrast to trans-PI, where deviations persist over times extending to terminal dynamics. Our all-atom simulations provide a fundamental understanding of PI dynamics and the impact of microstructure while providing important data for the design and optimization of PI-based materials.

Accelerated Diffusion Following Deprotection in Chemically Amplified Resists

Polymeric chemically amplified resists (CARs) are critical materials for high-throughput lithographic processes. A photoactivated acid-anion catalyst changes the polymer's solubility via a deprotection reaction, which enables pattern development through selective dissolution. To capture observed reaction kinetics, reaction-diffusion models employ a catalyst diffusivity that is accelerated by reaction. However, the microscopic origin and factors contributing to this phenomena remain unclear. Herein, we employ detailed atomistic molecular dynamics simulations to examine the impact of protecting group removal and material relaxation on catalyst mobility. We report data on polymer density, catalyst dispersion, excess free volume, and segmental dynamics with increasing time/extent of deprotection. We then propose simple kinetic Monte Carlo algorithms that can describe both molecular dynamics simulations of deprotection reactions and experimental data.

Assessment of an anisotropic coarse-grained model for cis-1,4-polybutadiene: a bottom-up approach

The spherical representation usually utilized for the coarse-grained particles of soft matter systems is an assumption and pertinent studies have shown that both structural and dynamical properties can depend on anisotropic effects. On these grounds, we develop coarse-grained equations of motion which take into account explicitly the anisotropy of the beads. As a first step, this model incorporates only conservative terms. Inclusion of the dissipative and random terms is in principle possible but is beyond the scope of this study. The translational dynamics of the beads is tracked using the position and momentum of their center of mass, while their rotational dynamics is modeled by representing their orientation through the use of quaternions, similarly to the case of rigid bodies. The associated force and torque controlling the motion are derived from atomistic molecular dynamics (MD) simulations via a bottom-up approach and define a coarse-grained potential. The assumptions of the model are clearly stated and checked for a reference system of a cis-1,4-polybutadiene melt. In particular, the choice of the angular velocity as a slow variable is justified by comparing its dynamics to atomic vibrations. The accuracy of this approach to reproduce static structural features of the polymer melt is assessed.

Thermal Conductivities of Crosslinked Polyisoprene and Polybutadiene from Molecular Dynamics Simulations

For the first time, the thermal conductivities of vulcanized polybutadiene and polyisoprene have been investigated according to their degree of crosslinking. The C-C and C-S-S-C crosslink bridges, which can be obtained via vulcanization processes using peroxides and sulfur, respectively, are considered. The temperature dependence of the thermal conductivity of soft rubber derived from molecular dynamics (MD) simulations is in very good agreement with the experimental results. The contributions of bonded and non-bonded interactions in the MD simulations and their influence on the thermal conductivities of polyisoprene and polybutadiene are presented. The details are discussed in this paper.

In Situ Determination of Carbon Number Distributions of Mixtures of Linear Hydrocarbons Confined within Porous Media Using Pulsed Field Gradient NMR

Pulsed field gradient (PFG) NMR measurements, combined with a novel optimization method, are used to determine the composition of hydrocarbon mixtures of linear alkanes (C7-C16) in both the bulk liquid state and when imbibed within a porous medium of mean pore diameter 28.6 nm. The method predicts the average carbon number of a given mixture to an accuracy of ±1 carbon number and the mole fraction of a mixture component to within an average root-mean-square error of ±0.036 with just three calibration mixtures. Given that the method can be applied at any conditions of temperature and pressure at which the PFG NMR measurements are made, the method has the potential for application in characterizing hydrocarbon liquid mixtures inside porous media and at the operating conditions relevant to, for example, hydrocarbon recovery and heterogeneous catalysis.

Coarse-graining of polyisoprene melts using inverse Monte Carlo and local density potentials

Bottom-up coarse-graining of polymers is commonly performed by matching structural order parameters such as distribution of bond lengths, bending and dihedral angles, and pair distribution functions. In this study, we introduce the distribution of nearest-neighbors as an additional order parameter in the concept of local density potentials. We describe how the inverse-Monte Carlo method provides a framework for forcefield development that is capable of overcoming challenges associated with the parameterization of interaction terms in polymer systems. The technique is applied on polyisoprene melts as a prototype system. We demonstrate that while different forcefields can be developed that perform equally in terms of matching target distributions, the inclusion of nearest-neighbors provides a straightforward route to match both thermodynamic and conformational properties. We find that several temperature state points can also be addressed, provided that the forcefield is refined accordingly. Finally, we examine both the single-particle and the collective dynamics of the coarse-grain models, demonstrating that all forcefields present a similar acceleration relative to the atomistic systems.

Molecular dynamics study of natural rubber-fullerene composites: connecting microscopic properties to macroscopic behavior

Macroscopic and microscopic properties of fullerene (C₆₀)-cis-polyisoprene (cis-PI) composites at varying fullerene concentrations were investigated using atomistic molecular dynamics (MD) simulations over microsecond time scales. Results show that the introduction of fullerenes into a polymer matrix increases density, bulk modulus and heat capacity while thermal expansivity decreases. The presence of fullerenes slowed the diffusion of both C₆₀ and cis-PI. Moreover, increasing fullerene concentration results in ordering of the cis-PI chains at the cis-PI-fullerene interfaces and shrinking of bulk PI regions. Free energy calculations of fullerene dimerization suggest that fullerenes disperse at low and aggregate at high fullerene concentrations. Our multi-scaled analysis approach demonstrates the role of 'ordered' regions adjacent to the interface between cis-PI and fullerene in controlling the level of order and mobility of the cis-PI chains. The relationship between the microscopic behavior and the changes in mechanical and thermal properties are discussed. Our study is beneficial for further studies and development of advanced rubber technology for novel, cost-effective, material with very high stiffness and thermal endurance with optimizing conditions of filler contents.

Response of adhesive polymer interfaces to repeated mechanical loading and the spatial variation of diffusion coefficient and stresses in a deforming polymer film

Comprehensive molecular simulations are conducted to show that polymer crosslinks preserve the strength of solid-polymer (melt) interfaces when they are subjected to repeated mechanical loading. The spatial variation of the diffusion coefficient and local stresses is also investigated along the polymer thickness, during deformation. After each loading cycle, a reduction in entanglement strength is observed at the fracture site. The work of adhesion also decreases over consecutive loading cycles, when fracture is induced at the same site. Reduction in both, the work of adhesion and the entanglement strength, decreases as the crosslink density increases. Diffusion coefficient and stresses vary significantly and in a complex manner along the film thickness during the entire deformation process. These variations were due to peculiar configurations occurring at each instance of separation, which are analyzed and explained in this work. The variation of diffusion coefficient during deformation suggests that other dynamic properties, such as viscosity, also vary spatially during polymer deformation.

Slip Spring-Based Mesoscopic Simulations of Polymer Networks: Methodology and the Corresponding Computational Code

In previous work by the authors, a new methodology was developed for Brownian dynamics/kinetic Monte Carlo (BD/kMC) simulations of polymer melts. In this study, this methodology is extended for dynamical simulations of crosslinked polymer networks in a coarse-grained representation, wherein chains are modeled as sequences of beads, each bead encompassing a few Kuhn segments. In addition, the C++ code embodying these simulations, entitled Engine for Mesoscopic Simulations for Polymer Networks (EMSIPON) is described in detail. A crosslinked network of cis-1,4-polyisoprene is chosen as a test system. From the thermodynamic point of view, the system is fully described by a Helmholtz energy consisting of three explicit contributions: entropic springs, slip springs and non-bonded interactions. Entanglements between subchains in the network are represented by slip springs. The ends of the slip springs undergo thermally activated hops between adjacent beads along the chain backbones, which are tracked by kinetic Monte Carlo simulation. In addition, creation/destruction processes are included for the slip springs at dangling subchain ends. The Helmholtz energy of non-bonded interactions is derived from the Sanchez⁻Lacombe equation of state. The isothermal compressibility of the polymer network is predicted from equilibrium density fluctuations in very good agreement with the underlying equation of state and with experiment. Moreover, the methodology and the corresponding C++ code are applied to simulate elongational deformations of polymer rubbers. The shear stress relaxation modulus is predicted from equilibrium simulations of several microseconds of physical time in the undeformed state, as well as from stress-strain curves of the crosslinked polymer networks under deformation.

Coarse-grained simulations of cis- and trans-polybutadiene: A bottom-up approach

We apply the dissipative particle dynamics strategy proposed by Hijón et al. [Faraday Discuss. 144, 301-322 (2010)] and based on an exact derivation of the generalized Langevin equation to cis- and trans-1,4-polybutadiene. We prove that it is able to reproduce not only the structural but also the dynamical properties of these polymers without any fitting parameter. A systematic study of the effect of the level of coarse-graining is done on cis-1,4-polybutadiene. We show that as the level of coarse-graining increases, the dynamical properties are better and better reproduced while the structural properties deviate more and more from those calculated in molecular dynamics (MD) simulations. We suggest two reasons for this behavior: the Markovian approximation is better satisfied as the level of coarse-graining increases, while the pair-wise approximation neglects important contributions due to the relative orientation of the beads at large levels of coarse-graining. Finally, we highlight a possible limit of the Markovian approximation: the fact that in constrained simulations, in which the centers-of-mass of the beads are kept constant, the bead rotational dynamics become extremely slow.

Multiscale Molecular Simulations of Polymer-Matrix Nanocomposites: or What Molecular Simulations Have Taught us About the Fascinating Nanoworld

Following the substantial progress in molecular simulations of polymer-matrix nanocomposites, now is the time to reconsider this topic from a critical point of view. A comprehensive survey is reported herein providing an overview of classical molecular simulations, reviewing their major achievements in modeling polymer matrix nanocomposites, and identifying several open challenges. Molecular simulations at multiple length and time scales, working hand-in-hand with sensitive experiments, have enhanced our understanding of how nanofillers alter the structure, dynamics, thermodynamics, rheology and mechanical properties of the surrounding polymer matrices.

Predicting diffusion coefficients of chemicals in and through packaging materials

Most of the physicochemical properties in polymers such as activity and partition coefficients, diffusion coefficients, and their activation with temperature are accessible to direct calculations from first principles. Such predictions are particularly relevant for food packaging as they can be used (1) to demonstrate the compliance or safety of numerous polymer materials and of their constitutive substances (e.g. additives, residues…), when they are used: as containers, coatings, sealants, gaskets, printing inks, etc. (2) or to predict the indirect contamination of food by pollutants (e.g. from recycled polymers, storage ambiance…) (3) or to assess the plasticization of materials in contact by food constituents (e.g. fat matter, aroma…). This review article summarizes the classical and last mechanistic descriptions of diffusion in polymers and discusses the reliability of semi-empirical approaches used for compliance testing both in EU and US. It is concluded that simulation of diffusion in or through polymers is not limited to worst-case assumptions but could also be applied to real cases for risk assessment, designing packaging with low leaching risk or to synthesize plastic additives with low diffusion rates.

Temperature- and pressure-dependent densities, self-diffusion coefficients, and phase behavior of monoacid saturated triacylglycerides: toward molecular-level insights into processing

Using molecular-dynamics (MD) simulations the densities and self-diffusion coefficients of a range of liquid monoacid triacylglycerides (TAGs) have been studied as a function of temperature and, for the first time, pressure. While offset by their ambient properties, the response of the TAGs to temperature and pressure is qualitatively similar. Application of pressure was found to significantly increase densities and reduce diffusion of the TAG molecules, suggesting that it may have as much a role in processing and crystallizing TAGs as supercooling does. A solution of glycerol tripalmitate and glycerol trihexanoate was also studied, showing that application of pressure should lead to a significant decrease in the saturation point of the solution, which is an important consideration for processing TAGs. Different solid/liquid interfaces of glycerol tripalmitate have also been investigated. Although crystal growth could not be observed, dissolution of one interface was seen in the MD simulations. The results suggest that over moderate distances the melting of TAGs may be cooperative in nature, rather than involving dissolution of individual TAG molecules.