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.
Hydrocarbon polymers, foams and nanocomposites are increasingly being subjected to extreme environments. Molecular scale modeling of these materials offers insight into failure mechanisms and complex response. Prior classical molecular dynamics (MD) simulations of the principal shock Hugoniot for two hydrocarbon polymers, polyethylene (PE) and poly (4-methyl-1-pentene) (PMP) have shown good agreement with density functional theory (DFT) calculations and experiments conducted at Sandia National Laboratories. We extended these results to include low-density polymer foams using nonequilibrium MD techniques and found good quantitative agreement with experiment. Here, we have measured the local temperature during void collapse to investigate the formation of hot spots and their relationship to polymer dissociation in foams.