Publications

Abstract not provided.

Reaching exceptionally long times up to 500 ns in equilibrium and nonequilibrium molecular dynamics simulations studies, we have attained a fundamental molecular understanding of the correlation of ionomer clusters structure and multiscale dynamics, providing new insight into one critical, long-standing challenge in ionic polymer physics. The cluster structure in melts of sulfonated polystyrene with Na+ and Mg2+ counterions are resolved and correlated with the dynamics on multiple length and time scales extracted from measurements of the dynamic structure factor and shear rheology. We find that as the morphology of the ionic clusters changes from ladderlike for Na+ to disordered structures for Mg2+, the dynamic structure factor is affected on the length scale corresponding to the ionic clusters. Rheology studies show that the viscosity for Mg2+ melts is higher than for Na+ ones for all shear rates, which is well correlated with the larger ionic clusters' size for the Mg2+ melts.

The conformation and dynamics of luminescent polymers collapsed into nanoparticles or polydots were studied using fully atomistic molecular dynamics (MD) simulations, providing a first insight into their internal dynamics. Controlling the conformation and dynamics of confined polymers is essential for realization of the full potential of polydots in nanomedicine and biotechnology. Specifically, the shape and internal dynamics of polydots that consist of highly rigid dialkyl p-phenylene ethynylene (PPE) are probed as a function of temperature. At room temperature, the polydots are spherical without any correlations between the aromatic rings on the PPE backbone. With increasing temperature, they expand and become slightly aspherical; however, the polymers remain confined. The coherent dynamic structure factor reveals that the internal motion of the polymer backbone is arrested, and the side chains dominate the internal dynamics of the polydots. These new soft nanoparticles retain their overall shape and dynamics over an extended temperature range, and their conformation is tunable via their degree of expansion.

Molecular dynamics simulations of a coarse-grained bead-spring model have been used to study the effects of molecular crowding on the accumulation of tension in the backbone of bottle-brush polymers tethered to a flat substrate. The number of bottle-brushes per unit surface area, Σ, as well as the lengths of the bottle-brush backbones Nbb (50 ≤ Nbb ≤ 200) and side chains Nsc (50 ≤ Nsc ≤ 200) were varied to determine how the dimensions and degree of crowding of bottle-brushes give rise to bond tension amplification along the backbone, especially near the substrate. From these simulations, we have identified three separate regimes of tension. For low Σ, the tension is due solely to intramolecular interactions and is dominated by the side chain repulsion that governs the lateral brush dimensions. With increasing Σ, the interactions between bottle-brush polymers induce compression of the side chains, transmitting increasing tension to the backbone. For large Σ, intermolecular side chain repulsion increases, forcing side chain extension and reorientation in the direction normal to the surface and transmitting considerable tension to the backbone.

Nanoparticles (NPs) grafted with organic layers form hybrids able to retain their unique properties through integration into the mesoscopic scale. The organic layer structure and response often determine the functionality of the hybrids on the mesoscopic length scale. Using molecular dynamics (MD) simulations, we probe the conformation of luminescent rigid polymers, dialkyl poly(p-phenylene ethynylene)s (PPE), end-grafted onto a silica nanoparticle in different solvents as the molecular weights and polymer coverages are varied. We find that, in contrast to NP-grafted flexible polymers, the chains are fully extended independent of the solvent. In toluene and decane, which are good solvents, the grafted PPEs chains assume a similar conformation to that observed in dilute solutions. In water, which is a poor solvent for the PPEs, the polymer chains form one large cluster but remain extended. The radial distribution of the chains around the core of the nanoparticle is homogeneous in good solvents, whereas in poor solvents clusters are formed independent of molecular weights and coverages. The clustering is distinctively different from the response of grafted flexible and semiflexible polymers.

The conformation of dinonyl poly para phenylene ethynylenes (PPEs) with carboxylate side chains, equilibrated in solvents of different quality have been studied using molecular dynamics simulations. PPEs are of interest because of their tunable electro-optical properties, chemical diversity, and functionality which are essential in wide range of applications. The polymer conformation determines the conjugation length and their assembly mode and affects electro-optical properties which are critical in current and potential uses. This study investigates the effect of carboxylate fraction on PPEs side chains on the conformation of chains in the dilute limit, in solvents of different quality. The dinonyl PPE chains are modeled atomistically, where the solvents are modeled both implicitly and explicitly. Dinonyl PPEs maintained a stretched out conformation up to a carboxylate fraction f of 0.7 in all solvents studied. The nonyl side chains are extended and oriented away from the PPE backbone in toluene and in implicit good solvent, whereas in water and implicit poor solvent, the nonyl side chains are collapsed toward the PPE backbone. Rotation around the aromatic ring is fast and no long range correlations are seen within the backbone.

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how the coarse-graining scale affects the measured dynamics. Iterative Boltzmann inversion is used to derive coarse-grained potentials with 2-6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large-scale dynamics. Using these models we simulate polyethylene melts for times over 500 μs to study the viscoelastic properties of well-entangled polymer melts.

Tailoring the nature of individual segments within ion containing block co-polymers is one critical design tool to achieve desired properties. The local structure including the size and distribution of the ionic blocks, as well as the long range correlations, are crucial for their transport ability. Here, we present molecular dynamics simulations on the effects of varying the concentrations of the ionizable groups on the conformations of pentablock ionomer that consist of a center block of ionic sulfonated styrene tethered to polyethylene and terminated by a bulky substituted styrene in dilute solutions. Sulfonation fractions f (0 ≤ f ≤ 0.55), spanning the range from ionomer to polyelectrolytes, were studied. Results for the equilibrium conformation of the chains in water and a 1:1 mixture of cyclohexane and heptane are compared to that in implicit poor solvents with dielectric constants ε = 1.0 and 77.73. In water, the pentablock collapses with the sulfonated groups on the outer surface. As f increases, the ionic, center block increasingly segregates from the hydrophobic regions. In the 1:1 mixture of cyclohexane and heptane, the flexible blocks swell, while the center ionic block collapses for f > 0. For f = 0, all blocks swell. In both implicit poor solvents, the pentablock collapses into a nearly spherical shape for all f. The sodium counterions disperse widely throughout the simulation cell for both water and ε = 77.73, whereas for ε = 1.0 and mixture of cyclohexane and heptane, the counterions largely condense onto the collapsed pentablock.

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions φ in the vicinity of the jamming transition at φc. Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal "Fickian" diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to φc. The time required to recover normal diffusion t∗ scales as (φ - φc)-0.5 and the long-time diffusivity D∞ ∼ (φ - φc)0.5, or D∞ ∼ 1/t∗. It is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both t∗ and D∞ in the limit φ → φc.

Ionic clusters control the structure, dynamics, and transport in soft matter. Incorporating a small fraction of ionizable groups in polymers substantially reduces the mobility of the macromolecules in melts. These ionic groups often associate into random clusters in melts, where the distribution and morphology of the clusters impact the transport in these materials. Here, using molecular dynamic simulations we demonstrate a clear correlation between cluster size and morphology with the polymer mobility in melts of sulfonated polystyrene. We show that in low dielectric media ladderlike clusters that are lower in energy compared with spherical assemblies are formed. Reducing the electrostatic interactions by enhancing the dielectric constant leads to morphological transformation from ladderlike clusters to globular assemblies. Decrease in electrostatic interaction significantly enhances the mobility of the polymer.

The high mechanical stiffness of single-nanoparticle-thick membranes is believed to result from the local structure of ligand coatings that mediate interactions between nanoparticles. These ligand structures are not directly observable experimentally. We use molecular dynamics simulations to observe variations in ligand structure and simultaneously measure variations in membrane mechanical properties. We have shown previously that ligand end group has a large impact on ligand structure and membrane mechanical properties. Here we introduce and apply quantitative molecular structure measures to these membranes and extend analysis to multiple nanoparticle core sizes and ligand lengths. Simulations of nanoparticle membranes with a nanoparticle core diameter of 4 or 6 nm, a ligand length of 11 or 17 methylenes, and either carboxyl (COOH) or methyl (CH3) ligand end groups are presented. In carboxyl-terminated ligand systems, structure and interactions are dominated by an end-to-end orientation of ligands. In methyl-terminated ligand systems large ordered ligand structures form, but nanoparticle interactions are dominated by disordered, partially interdigitated ligands. Core size and ligand length also affect both ligand arrangement within the membrane and the membrane's macroscopic mechanical response, but are secondary to the role of the ligand end group. Moreover, the particular end group (COOH or CH3) alters the nature of how ligand length, in turn, affects the membrane properties. The effect of core size does not depend on the ligand end group, with larger cores always leading to stiffer membranes. Asymmetry in the stress and ligand density is observed in membranes during preparation at a water-vapor interface, with the stress asymmetry persisting in all membranes after drying.

Abstract not provided.

Large-scale molecular dynamics simulations are used to study the internal relaxations of chains in nanoparticle (NP)/polymer composites. We examine the Rouse modes of the chains, a quantity that is closest in spirit to the self-intermediate scattering function, typically determined in an (incoherent) inelastic neutron scattering experiment. Our simulations show that for weakly interacting mixtures of NPs and polymers, the effective monomeric relaxation rates are faster than in a neat melt when the NPs are smaller than the entanglement mesh size. In this case, the NPs serve to reduce both the monomeric friction and the entanglements in the polymer melt, as in the case of a polymer-solvent system. However, for NPs larger than half the entanglement mesh size, the effective monomer relaxation is essentially unaffected for low NP concentrations. Even in this case, we observe a strong reduction in chain entanglements for larger NP loadings. Thus, the role of NPs is to always reduce the number of entanglements, with this effect only becoming pronounced for small NPs or for high concentrations of large NPs. Our studies of the relaxation of single chains resonate with recent neutron spin echo (NSE) experiments, which deduce a similar entanglement dilution effect.

In this study, the properties of mechanically stable single-nanoparticle (NP)-thick membranes have largely been studied at room temperature. How these membranes soften as nanoparticle ligands disorder with increasing temperature is unknown. Molecular dynamics simulations are used to probe the temperature dependence of the mechanical and nanostructural properties of nanoparticle membranes made of 6 nm diameter Au nanoparticles coated with dodecanethiol ligands and terminated with either methyl (CH₃) or carboxyl (COOH) terminal groups. For methyl-terminated ligands, interactions along the alkane chain provide mechanical stiffness, with a Young's modulus of 1.7 GPa at 300 K. For carboxyl-terminated chains, end-group interactions are significant, producing stiffer membranes at all temperatures, with a Young's modulus of 3.8 GPa at 300 K. For both end-group types, membrane stiffness is reduced to zero at about 400 K. Ligand structure and mechanical properties of membranes at 300 K that have been annealed at 400 K are comparable to samples that do not undergo thermal annealing.

The conformation of single chain pentablock ionomers (A-B-C-B-A) containing randomly sulfonated polystyrene in the center block, tethered to poly-ethylene-r-propylene end-capped by poly-t-butyl styrene is studied in dilute solutions by molecular dynamics simulations. Multi-block copolymers offer a means to tailor several properties into one molecule, taking advantage of their rich phase diagram together with unique properties of specific blocks. For this pentablock the ionic block facilitates transport while the A and B components are incorporated for mechanical stability. The present study investigates the confirmation of a single chain of pentablock ionomer of molecular weight M_w ~ 50,000 g/mol and sulfonated polystyrene of the same molecular weight as that of the center block for six sulfonation fractions f from f=0.0-0.55. For the sulfonated systems Na⁺ counterions are included. Results for the equilibrium conformation of the chains and the three blocks in water and 1:1 mixture of cyclohexane and n-heptane are compared to simulations in implicit poor solvents with dielectric constants ε =1.0 and 77.73. In water, the pentablock is collapsed with sulfonated groups on the outer surface. As the sulfonation fraction f increases, the ionic, center block is increasingly segregated from the hydrophobic regions. In the 1:1 mixture of cyclohexane and heptane both the flexible and end blocks are swollen while the center ionic block is collasped for f>0, while for f=0 all blocks are swollen. In both implicit poor solvents the pentablock is collapsed into a nearly spherical shape for all f. The sodium counterions are dispersed widely throughout the simulation cell for both water and ε =77.73 whereas for ε =1.0 the counterions are largely condensed on the collapsed pentablock.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Structured polymers offer a means to tailor transport pathways within mechanically stable manifolds. The building block of such a membrane is examined, namely a single large pentablock co-polymer that consists of a center block of a randomly sulfonated polystyrene, designed for transport, tethered to poly-Ethylene-R-Propylene and end-Capped by poly-T-Butyl styrene, for mechanical stability, using molecular dynamics simulations. The polymer structure in a cyclohexane-Heptane mixture, a technologically viable solvent, and in water, a poor solvent for all segments and a ubiquitous substance is extracted. In all solvents the pentablock collapsed into nearly spherical aggregates where the ionic block is segregated. In hydrophobic solvents, the ionic block resides in the center, surrounded by swollen intermix of flexible and end blocks. In water all blocks are collapsed with the sulfonated block residing on the surface. Our results demonstrate that solvents drive different local nano-Segregation, providing a gateway to assemble membranes with controlled topology.

Abstract not provided.

Large-scale molecular simulations are performed to investigate tensile failure of polymer interfaces as a function of welding time t. Changes in the tensile stress, mode of failure and interfacial fracture energy GI are correlated to changes in the interfacial entanglements as determined from primitive path analysis. Bulk polymers fail through craze formation, followed by craze breakdown through chain scission. At small t welded interfaces are not strong enough to support craze formation and fail at small strains through chain pullout at the interface. Once chains have formed an average of about one entanglement across the interface, a stable craze is formed throughout the sample. The failure stress of the craze rises with welding time and the mode of craze breakdown changes from chain pullout to chain scission as the interface approaches bulk strength. The interfacial fracture energy GI is calculated by coupling the simulation results to a continuum fracture mechanics model. As in experiment, GI increases as t1/2 before saturating at the average bulk fracture energy Gb. As in previous simulations of shear strength, saturation coincides with the recovery of the bulk entanglement density. Before saturation, GI is proportional to the areal density of interfacial entanglements. Immiscibiltiy limits interdiffusion and thus suppresses entanglements at the interface. Even small degrees of immiscibility reduce interfacial entanglements enough that failure occurs by chain pullout and GI 蠐 Gb.

We report classical atomistic molecular dynamics simulations of alkanethiol-coated gold nanoparticles solvated in water and decane, as well as at water/vapor interfaces. The structure of the coatings is analyzed as a function of various functional end groups, including amine and carboxyl groups in various ionization states. We study both neutral and charged end groups for two different chain lengths (9 and 17 carbons). For the charged end groups, we simulated both mono- and divalent counterions. For the longer alkanes, we find significant local bundling of chains on the nanoparticle surface, which results in highly asymmetric coatings. In general, the charged end groups attenuate this effect by enhancing the water solubility of the nanoparticles. On the basis of the coating structures and density profiles, we can qualitatively infer the overall solubility of the nanoparticles. This asymmetry in the alkanethiol coatings is likely to have a significant effect on aggregation behavior. Our simulations elucidate the mechanism by which modulating the end group charge state can be used to control coating structure and therefore nanoparticle solubility and aggregation behavior.

We use molecular dynamics simulations of the Kremer–Grest (KG) bead–spring model of polymer chains of length between 10 and 500, and a closely related analogue that allows for chain crossing, to clearly delineate the effects of entanglements on the length-scale-dependent chain relaxation in polymer melts. We analyze the resulting trajectories using the Rouse modes of the chains and find that entanglements strongly affect these modes. The relaxation rates of the chains show two limiting effective monomeric frictions, with the local modes experiencing much lower effective friction than the longer modes. The monomeric relaxation rates of longer modes vary approximately inversely with chain length due to kinetic confinement effects. The time-dependent relaxation of Rouse modes has a stretched exponential character with a minimum of stretching exponent in the vicinity of the entanglement chain length. None of these trends are found in models that allow for chain crossing. As a result, these facts, in combination, argue for the confined motion of chains for time scales between the entanglement time and their ultimate free diffusion.

We present a review and critique of several methods for the simulation of the dynamics of colloidal suspensions at the mesoscale. We focus particularly on simulation techniques for hydrodynamic interactions, including implicit solvents (Fast Lubrication Dynamics, an approximation to Stokesian Dynamics) and explicit/particle-based solvents (Multi-Particle Collision Dynamics and Dissipative Particle Dynamics). Several variants of each method are compared quantitatively for the canonical system of monodisperse hard spheres, with a particular focus on diffusion characteristics, as well as shear rheology and microstructure. In all cases, we attempt to match the relevant properties of a well-characterized solvent, which turns out to be challenging for the explicit solvent models. Reasonable quantitative agreement is observed among all methods, but overall the Fast Lubrication Dynamics technique shows the best accuracy and performance. We also devote significant discussion to the extension of these methods to more complex situations of interest in industrial applications, including models for non-Newtonian solvent rheology, non-spherical particles, drying and curing of solvent and flows in complex geometries. This work identifies research challenges and motivates future efforts to develop techniques for quantitative, predictive simulations of industrially relevant colloidal suspension processes.

Self-healing of polymer films often takes place as the molecules diffuse across a damaged region, above their melting temperature. Using molecular dynamics simulations we probe the healing of polymer films and compare the results with those obtained for thermal welding of homopolymer slabs. These two processes differ from each other in their interfacial structure since damage leads to increased polydispersity and more short chains. A polymer sample was cut into two separate films that were then held together in the melt state. The recovery of the damaged film was followed as time elapsed and polymer molecules diffused across the interface. The mass uptake and formation of entanglements, as obtained from primitive path analysis, are extracted and correlated with the interfacial strength obtained from shear simulations. We find that the diffusion across the interface is significantly faster in the damaged film compared to welding because of the presence of short chains. Though interfacial entanglements increase more rapidly for the damaged films, a large fraction of these entanglements are near chain ends. As a result, the interfacial strength of the healing film increases more slowly than for welding. For both healing and welding, the interfacial strength saturates as the bulk entanglement density is recovered across the interface. However, the saturation strength of the damaged film is below the bulk strength for the polymer sample. At saturation, cut chains remain near the healing interface. They are less entangled and as a result they mechanically weaken the interface. Chain stiffness increases the density of entanglements, which increases the strength of the interface. Our results show that a few entanglements across the interface are sufficient to resist interfacial chain pullout and enhance the mechanical strength. © 2014 American Physical Society.