Predicting χ of Polymer Blends Using Atomistic Morphing Simulations

The Flory Huggins interaction parameter chi measures the compatibility of different species in mixtures and governs their phase behavior. We have previously used molecular dynamics (MD) simulations and thermodynamic integration along the path of transformation of one species to another (morphing), to determine chi in coarse-grained bead-spring models of polymer blends. In this work, we use united-atom (UA) MD simulations and morphing to calculate chi for real polymer blends: (1) poly(ethylene) and poly(ethylene oxide), (2) poly(styrene) and poly(2-vinyl pyridine), (3) poly(isoprene) and saturated poly(isoprene), and (4) poly(styrene) and poly(alpha-methyl styrene). These examples require different schemes for transforming chains: changing Lennard Jones parameters and partial charges (case 1 and 2), transforming double bonds to single bonds (case 3), and making atoms disappear (case 4). For the first three blends, chi predictions agree reasonably with experiments but are sensitive to the choice of force field parameters. For poly(styrene)/poly(alpha-methyl styrene), we reach the limits of the morphing method.

Automated summary

The study used molecular dynamics simulations and morphing method to calculate the interaction parameter chi for real polymer blends, with different transformation schemes needed for different polymers. Predictions for the first three blends agree reasonably with experiments but are sensitive to force field parameters.