Publications

Srivastava, Ishan S.; Lechman, Jeremy B.; Bolintineanu, Dan S.; Silbert, Leonardo E.; Grest, Gary S.

Abstract not provided.

Srivastava, Ishan S.; Lechman, Jeremy B.; Silbert, Leionardo E.; Grest, Gary S.

Abstract not provided.

Grest, Gary S.; Peters, Brandon L.; Salerno, K.M.; Ge, Ting G.; Perahia, Dvora P.

Abstract not provided.

Lane, James M.; Srivastava, Ishan S.; Peters, Brandon L.; Fan, Hongyou F.; Grest, Gary S.; Salerno, K.M.

Abstract not provided.

Physical Review E

Koski, Jason K.; Moore, Stan G.; Grest, Gary S.; Stevens, Mark J.

The role of an external field on capillary waves at the liquid-vapor interface of a dipolar fluid is investigated using molecular dynamics simulations. For fields parallel to the interface, the interfacial width squared increases linearly with respect to the logarithm of the size of the interface across all field strengths tested. The value of the slope decreases with increasing field strength, indicating that the field dampens the capillary waves. With the inclusion of the parallel field, the surface stiffness increases with increasing field strength faster than the surface tension. For fields perpendicular to the interface, the interfacial width squared is linear with respect to the logarithm of the size of the interface for small field strengths, and the surface stiffness is less than the surface tension. Above a critical field strength that decreases as the size of the interface increases, the interface becomes unstable due to the increased amplitude of the capillary waves.

npj Computational Materials

Meng, Dong; Kumar, Sanat K.; Grest, Gary S.; Mahynski, Nathan A.; Panagiotopoulos, Athanassios Z.

A large body of experimental work has established that athermal colloid/polymer mixtures undergo a sequence of transitions from a disordered fluid state to a colloidal crystal to a second disordered phase with increasing polymer concentration. These transitions are driven by polymer-mediated interparticle attraction, which is a function of both the polymer density and size. It has been posited that the disordered state at high polymer density is a consequence of strong interparticle attractions that kinetically inhibit the formation of the colloidal crystal, i.e., the formation of a non-equilibrium gel phase interferes with crystallization. Here we use molecular dynamics simulations and density functional theory on polymers and nanoparticles (NPs) of comparable size and show that the crystal-disordered phase coexistence at high polymer density for sufficiently long chains corresponds to an equilibrium thermodynamic phase transition. While the crystal is, indeed, stabilized at intermediate polymer density by polymer-induced intercolloid attractions, it is destabilized at higher densities because long chains lose significant configurational entropy when they are forced to occupy all of the crystal voids. Our results are in quantitative agreement with existing experimental data and show that, at least in the nanoparticle limit of sufficiently small colloidal particles, the crystal phase only has a modest range of thermodynamic stability.

ACS Macro Letters

Peters, Brandon L.; Pike, Darin Q.; Rubinstein, Michael; Grest, Gary S.

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.

Grest, Gary S.

Abstract not provided.

Langmuir

Aryal, Dipak; Agrawal, Anupriya; Perahia, Dvora; Grest, Gary S.

Controlling the structure and dynamics of thin films of ionizable polymers at water interfaces is critical to their many applications. As the chemical diversity within one polymer is increased, controlling the structure and dynamics of the polymer, which is a key to their use, becomes a challenge. Here molecular dynamics simulations (MD) are used to obtain molecular insight into the structure and dynamics of thin films of one such macromolecule at the interface with water. The polymer consists of an ABCBA topology with randomly sulfonated polystyrene (C), tethered symmetrically to flexible poly(ethylene-r-propylene) blocks (B), and end-capped by a poly(t-butylstyrene) block (A). The compositions of the interfacial and bulk regions of thin films of the ABCBA polymers are followed as a function of exposure time to water. We find that interfacial rearrangements take place where buried ionic segments migrate toward the water interface. The hydrophobic blocks collapse and rearrange to minimize their exposure to water. The water that initially drives interfacial reengagements breaks the ionic clusters within the film, forming a dynamic hydrophilic internal network within the hydrophobic segments.

Macromolecules

Aryal, Dipak; Agrawal, Anupriya; Perahia, Dvora; Grest, Gary S.

Structure and dynamics of melts of copolymers with an ABCBA topology, where C is an ionizable block, have been studied by fully atomistic molecular dynamics (MD) simulations. Introducing an ionizable block for functionality adds a significant element to the coupled set of interactions that determine the structure and dynamics of the macromolecule. The polymer consists of a randomly sulfonated polystyrene C block tethered to a flexible poly(ethylene-r-propylene) bridge B and end-capped with poly(tert-butylstyrene) A. The chemical structure and topology of these polymers constitute a model for incorporation of ionic blocks within a framework that provides tactility and mechanical stability. Here we resolve the structure and dynamics of a structured polymer on the nanoscale constrained by ionic clusters. We find that the melts form intertwined networks of the A and C blocks independent of the degree of sulfonation of the C block with no long-range order. The cluster cohesiveness and morphology affect both macroscopic translational motion and segmental dynamics of all the blocks.

Srivastava, Ishan S.; Bolintineanu, Dan S.; Silbert, Leonardo E.; Lechman, Jeremy B.; Grest, Gary S.

Abstract not provided.

Grest, Gary S.

Coupled length and time scales determine the dynamic behavior of polymers and polymer nanocomposites and underlie their unique properties. To resolve the properties over large time and length scales it is imperative to develop coarse grained models which retain the atomistic specificity. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details a nd 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 polyethylene as a model system, we probe how the coarse - graining scale affects the measured dynamics with different number methylene group s per coarse - grained beads. Using these models we simulate polyethylene melts for times over 500 ms to study the viscoelastic properties of well - entangled polymer melts and large nanoparticle assembly as the nanoparticles are driven close enough to form nanostructures.

Peters, Brandon L.; Grest, Gary S.

Abstract not provided.

Journal of Chemical Physics

Wijesinghe, Sidath; Maskey, Sabina; Perahia, Dvora; Grest, Gary S.

Long-lived soft nanoparticles, formed by conjugated polymers, constitute a new class of far-from-equilibrium responsive structures for nano-medicine. Tethering ionizable groups to the polymers enables functionality. However concurrently, the ionic groups perturb the delicate balance of interactions that governs these particles. Using fully atomistic molecular dynamics simulations, this study probed the effects of charged groups tethered to poly para phenylene ethynylene substituted by alkyl groups on the polymer conformation and dynamics in confined geometry. We find that the ionizable groups affect the entire shape of the polydots and impact the conformation and dynamics of the polymer.

Journal of Chemical Theory and Computation

Peters, Brandon L.; Salerno, K.M.; Agrawal, Anupriya; Perahia, Dvora; Grest, Gary S.

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.

Grest, Gary S.

Abstract not provided.

Grest, Gary S.

Abstract not provided.

Nature Communications

Li, Binsong; Bian, Kaifu B.; Lane, J.M.; Salerno, K.M.; Grest, Gary S.; Ao, Tommy A.; Hickman, Randy J.; Wise, Jack L.; Wang, Zhongwu; Fan, Hongyou F.

Gold nanostructured materials exhibit important size- and shape-dependent properties that enable a wide variety of applications in photocatalysis, nanoelectronics and phototherapy. Here we show the use of superfast dynamic compression to synthesize extended gold nanostructures, such as nanorods, nanowires and nanosheets, with nanosecond coalescence times. Using a pulsed power generator, we ramp compress spherical gold nanoparticle arrays to pressures of tens of GPa, demonstrating pressure-driven assembly beyond the quasi-static regime of the diamond anvil cell. Our dynamic magnetic ramp compression approach produces smooth, shockless (that is, isentropic) one-dimensional loading with low-temperature states suitable for nanostructure synthesis. Transmission electron microscopy clearly establishes that various gold architectures are formed through compressive mesoscale coalescences of spherical gold nanoparticles, which is further confirmed by in-situ synchrotron X-ray studies and large-scale simulation. This nanofabrication approach applies magnetically driven uniaxial ramp compression to mimic established embossing and imprinting processes, but at ultra-short (nanosecond) timescales.

Peters, Brandon L.; Grest, Gary S.

Abstract not provided.

Lane, James M.; Salerno, K.M.; Grest, Gary S.; Fan, Hongyou F.

Abstract not provided.

Grest, Gary S.; Koski, Jason K.

Abstract not provided.

Macromolecules

Ge, Ting; Kalathi, Jagannathan T.; Halverson, Jonathan D.; Grest, Gary S.; Rubinstein, Michael

The motion of nanoparticles (NPs) in entangled melts of linear polymers and nonconcatenated ring polymers are compared by large-scale molecular dynamics simulations. The comparison provides a paradigm for the effects of polymer architecture on the dynamical coupling between NPs and polymers in nanocomposites. Strongly suppressed motion of NPs with diameter d larger than the entanglement spacing a is observed in a melt of linear polymers before the onset of Fickian NP diffusion. This strong suppression of NP motion occurs progressively as d exceeds a and is related to the hopping diffusion of NPs in the entanglement network. In contrast to the NP motion in linear polymers, the motion of NPs with d > a in ring polymers is not as strongly suppressed prior to Fickian diffusion. The diffusion coefficient D decreases with increasing d much slower in entangled rings than in entangled linear chains. NP motion in entangled nonconcatenated ring polymers is understood through a scaling analysis of the coupling between NP motion and the self-similar entangled dynamics of ring polymers.

Nanoscale

Aryal, Dipak; Grest, Gary S.; Perahia, Dvora

Directing the formation of nanostructures that serve as building blocks of membranes presents an immense step towards engineering controlled polymeric ion transport systems. Using the exquisite atomic detail captured by molecular dynamics simulations, we follow the assembly of a co-polymer that consists of polystyrene sulfonate tethered symmetrically to hydrophobic blocks, realizing a new type of long lived solvent-responsive soft nanoparticle.

Poster Sessions 2017 - Core Programming Area at the 2017 AIChE Annual Meeting

Ge, Ting; Rubinstein, Michael; Grest, Gary S.

Abstract not provided.

Journal of Chemical Physics

Etampawala, Thusitha N.; Aryal, Dipak; Osti, Naresh C.; He, Lilin; Heller, William T.; Willis, Carl L.; Grest, Gary S.; Perahia, Dvora

The self-assembly of multiblock copolymers in solutions is controlled by a delicate balance between inherent phase segregation due to incompatibility of the blocks and the interaction of the individual blocks with the solvent. The current study elucidates the association of pentablock copolymers in a mixture of selective solvents which are good for the hydrophobic segments and poor for the hydrophilic blocks using small angle neutron scattering (SANS). The pentablock consists of a center block of randomly sulfonated polystyrene, designed for transport, tethered to poly-ethylene-r-propylene and end-capped by poly-t-butyl styrene, for mechanical stability. We find that the pentablock forms ellipsoidal core-shell micelles with the sulfonated polystyrene in the core and Gaussian decaying chains of swollen poly-ethylene-r-propylene and poly-t-butyl styrene tertiary in the corona. With increasing solution concentration, the size of the micelle, the thickness of the corona, and the aggregation number increase, while the solvent fraction in the core decreases. In dilute solution the micelle increases in size as the temperature is increased, however, temperature effects dissipate with increasing solution concentration.