Effect of Chain Length Dispersity on the Mobility of Entangled Polymers.

While nearly all theoretical and computational studies of entangled polymer melts have focused on uniform samples, polymer synthesis routes always result in some dispersity, albeit narrow, of distribution of molecular weights (D_{M}=M_{w}/M_{n}∼1.02-1.04). Here, the effects of dispersity on chain mobility are studied for entangled, disperse melts using a coarse-grained model for polyethylene. Polymer melts with chain lengths set to follow a Schulz-Zimm distribution for the same average M_{w}=36 kg/mol with D_{M}=1.0 to 1.16, were studied for times of 600-800 µs using molecular dynamics simulations. This time frame is longer than the time required to reach the diffusive regime. We find that dispersity in this range does not affect the entanglement time or tube diameter. However, while there is negligible difference in the average mobility of chains for the uniform distribution D_{M}=1.0 and D_{M}=1.02, the shortest chains move significantly faster than the longest ones offering a constraint release pathway for the melts for larger D_{M}.

A multi-chain polymer slip-spring model with fluctuating number of entanglements: Density fluctuations, confinement, and phase separation.

Coarse grained simulation approaches provide powerful tools for the prediction of the equilibrium properties of polymeric systems. Recent efforts have sought to develop coarse-graining strategies capable of predicting the non-equilibrium behavior of entangled polymeric materials. Slip-link and slip-spring models, in particular, have been shown to be capable of reproducing several key aspects of the linear response and rheology of polymer melts. In this work, we extend a previously proposed multi-chain slip-spring model in a way that correctly incorporates the effects of the fluctuating environment in which polymer segments are immersed. The model is used to obtain the equation of state associated with the slip-springs, and the results are compared to those of related numerical approaches and an approximate analytical expression. The model is also used to examine a polymer melt confined into a thin film, where an inhomogeneous distribution of polymer segments is observed, and the corresponding inhomogeneities associated with density fluctuations are reflected on the spatial slip-spring distribution.

A multichain polymer slip-spring model with fluctuating number of entanglements for linear and nonlinear rheology.

A theoretically informed entangled polymer simulation approach is presented for description of the linear and non-linear rheology of entangled polymer melts. The approach relies on a many-chain representation and introduces the topological effects that arise from the non-crossability of molecules through effective fluctuating interactions, mediated by slip-springs, between neighboring pairs of macromolecules. The total number of slip-springs is not preserved but, instead, it is controlled through a chemical potential that determines the average molecular weight between entanglements. The behavior of the model is discussed in the context of a recent theory for description of homogeneous materials, and its relevance is established by comparing its predictions to experimental linear and non-linear rheology data for a series of well-characterized linear polyisoprene melts. The results are shown to be in quantitative agreement with experiment and suggest that the proposed formalism may also be used to describe the dynamics of inhomogeneous systems, such as composites and copolymers. Importantly, the fundamental connection made here between our many-chain model and the well-established, thermodynamically consistent single-chain mean-field models provides a path to systematic coarse-graining for prediction of polymer rheology in structurally homogeneous and heterogeneous materials.

Fully atomistic simulations of the response of silica nanoparticle coatings to alkane solvents.

Molecular dynamics simulations are used to study the effect of passivating ligands of varying lengths grafted to a nanoparticle and placed in various alkane solvents. Average height and density profiles for methyl-terminated alkoxylsilane ligands (-O-Si(OH)(2)(CH(2))(n)CH(3), with n = 9, 17, and 35) attached to a 5-nm-diameter amorphous silica nanoparticle with coverages of between 1.0 and 3.0 chains/nm(2) are presented for explicitly modeled, short-chain hydrocarbon solvents and for implicit good and poor solvents. Three linear solvents, C(10)H(22) (decane), C(24)H(50), and C(48)H(96), and a branched solvent, squalene, were studied. An implicit poor solvent captured the effect of the longest chain length solvent at lower temperatures, while its temperature dependence was similar to that of the branched solvent squalene. In contrast, an implicit good solvent produced coating structures that were far more extended than those found in any of the explicit solvents tested and showed little dependence on temperature. Coatings equilibrated in explicit solvents were more compact in longer-chain solvents because of autophobic dewetting. Changes in the coating density profiles were more pronounced as the solvent chain length was increased from decane to C(24)H(50) than from C(24)H(50) to C(48)H(98) for all coatings. The response of coatings in squalene was not significantly different from that of the linear chain of equal mass. Significant interpenetration of the solvent chains with the brush coating was observed only for the shortest grafted chains in decane. In all cases, the methyl terminal group was not confined to the coating edge but was found throughout the entire coating volume, from the core to the outermost shell. Increasing the temperature from 300 to 500 K led to greater average brush heights, but the dependence was weak.

Toward Bottom-Up Understanding of Transport in Concentrated Battery Electrolytes.

Bottom-up understanding of transport describes how molecular changes alter species concentrations and electrolyte voltage drops in operating batteries. Such an understanding is essential to predictively design electrolytes for desired transport behavior. We herein advocate building a structure-property-performance relationship as a systematic approach to accurate bottom-up understanding. To ensure generalization across salt concentrations as well as different electrolyte types and cell configurations, the property-performance relation must be described using Newman's concentrated solution theory. It uses Stefan-Maxwell diffusivity, ij , to describe the role of molecular motions at the continuum scale. The key challenge is to connect ij to the structure. We discuss existing methods for making such a connection, their peculiarities, and future directions to advance our understanding of electrolyte transport.

Pushing and Pulling: A Dual pH Trigger Controlled by Varying the Alkyl Tail Length in Heme Coordinating Peptide Amphiphiles.

Some organisms in nature that undergo anaerobic respiration utilize 1D nanoscale arrays of densely packed cytochromes containing the molecule heme. The assemblies can be mimicked with 1D nanoscale fibrils composed of peptide amphiphiles designed to coordinate heme in dense arrays. To create such materials and assemblies, it is critical to understand the assembly process and what controls the various aspects of hierarchical assembly. MD simulations suggest that shorter alkyl chains on the peptide lead to more dynamic structures than the peptides with longer chains that yield kinetically trapped states. The hydration parameters manifest themselves experimentally through the observation of a dual pH trigger, which controls the peptide assembly rate, the heme binding affinity, and heme organization kinetics. Great strides in understanding the relative complexity of the self-assembly process in relation to incorporating a functional moiety like heme opens up many possibilities in developing abiotic assemblies for bioelectronic devices and assemblies.

Free energy calculations of the functional selectivity of 5-HT2B G protein-coupled receptor.

G Protein-Coupled Receptors (GPCRs) mediate intracellular signaling in response to extracellular ligand binding and are the target of one-third of approved drugs. Ligand binding modulates the GPCR molecular free energy landscape by preferentially stabilizing active or inactive conformations that dictate intracellular protein recruitment and downstream signaling. We perform enhanced sampling molecular dynamics simulations to recover the free energy surfaces of a thermostable mutant of the GPCR serotonin receptor 5-HT2B in the unliganded form and bound to a lysergic acid diethylamide (LSD) agonist and lisuride antagonist. LSD binding imparts a ∼110 kJ/mol driving force for conformational rearrangement into an active state. The lisuride-bound form is structurally similar to the apo form and only ∼24 kJ/mol more stable. This work quantifies ligand-induced conformational specificity and functional selectivity of 5-HT2B and presents a platform for high-throughput virtual screening of ligands and rational engineering of the ligand-bound molecular free energy landscape.

Polymers at Liquid/Vapor Interface.

Polymers confined to the liquid/vapor interface are studied using molecular dynamics simulations. We show that for polymers which are weakly immiscible with the solvent, the density profile perpendicular to the liquid/vapor interface is strongly asymmetric. On the vapor side of the interface, the density distribution falls off as a Gaussian with a decay length on the order of the bead diameter, whereas on the liquid side, the density profile decays as a simple exponential. This result differs from that of a polymer absorbed from a good solvent with the density profile decaying as a power law. As the surface coverage increases, the average end-to-end distance and chain mobility systematically decreases toward that of the homopolymer melt.

Coarse-Grained Modeling of Polyethylene Melts: Effect on Dynamics.

The distinctive viscoelastic behavior of polymers results from a coupled interplay of motion on multiple length and time scales. Capturing the broad time and length scales of polymer motion remains a challenge. Using polyethylene (PE) as a model macromolecule, we construct coarse-grained (CG) models of PE with three to six methyl groups per CG bead and probe two critical aspects of the technique: pressure corrections required after iterative Boltzmann inversion (IBI) to generate CG potentials that match the pressure of reference fully atomistic melt simulations and the transferability of CG potentials across temperatures. While IBI produces nonbonded pair potentials that give excellent agreement between the atomistic and CG pair correlation functions, the resulting pressure for the CG models is large compared with the pressure of the atomistic system. We find that correcting the potential to match the reference pressure leads to nonbonded interactions with much deeper minima and slightly smaller effective bead diameter. However, simulations with potentials generated by IBI and pressure-corrected IBI result in similar mean-square displacements (MSDs) and stress autocorrelation functions G(t) for PE melts. While the time rescaling factor required to match CG and atomistic models is the same for pressure- and non-pressure-corrected CG models, it strongly depends on temperature. Transferability was investigated by comparing the MSDs and stress autocorrelation functions for potentials developed at different temperatures.