Deformation and Flow in Dense Granular Matter from Stress-based Discrete Element Simulations
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Physical Review E
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
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
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.
Abstract not provided.
Langmuir
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
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.
Abstract not provided.
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.
Abstract not provided.
Journal of Chemical Physics
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
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.
Abstract not provided.
Abstract not provided.
Nature Communications
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.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Macromolecules
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
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
Abstract not provided.
Journal of Chemical Physics
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.
Journal of Chemical Physics
The coupled dynamics of entangled polymers, which span broad time and length scales, govern their unique viscoelastic properties. To follow chain mobility by numerical simulations from the intermediate Rouse and reptation regimes to the late time diffusive regime, highly coarse grained models with purely repulsive interactions between monomers are widely used since they are computationally the most efficient. Here using large scale molecular dynamics simulations, the effect of including the attractive interaction between monomers on the dynamics of entangled polymer melts is explored for the first time over a wide temperature range. Attractive interactions have little effect on the local packing for all temperatures T and on the chain mobility for T higher than about twice the glass transition Tg. These results, across a broad range of molecular weight, show that to study the dynamics of entangled polymer melts, the interactions can be treated as pure repulsive, confirming a posteriori the validity of previous studies and opening the way to new large scale numerical simulations.
European Physical Journal. Special Topics
Polymer dynamics creates distinctive viscoelastic behavior as a result of a coupled interplay of motion at the atomic length scale and motion of the entire macromolecule. Capturing the broad time and length scales of polymeric motion however, remains a challenge. Using linear polyethylene as a model system, we probe the effects of the degree of coarse graining on polymer dynamics. Coarse-grained (CG) potentials are derived using iterative Boltzmann inversion with λ methylene groups per CG bead (denoted CGλ) with λ = 2,3,4 and 6 from a fully-atomistic polyethylene melt simulation. By rescaling time in the CG models by a factor α, the chain mobility for the atomistic and CG models match. We show that independent of the degree of coarse graining, all measured static and dynamic properties are essentially the same once the dynamic scaling factor α and a non-crossing constraint for the CG6 model are included. The speedup of the CG4 model is about 3 times that of the CG3 model and is comparable to that of the CG6 model. Furthermore, using these CG models we were able to reach times of over 500 μs, allowing us to measure a number of quantities, including the stress relaxation function, plateau modulus and shear viscosity, and compare directly to experiment.
Journal of Chemical Physics
The crazing behavior of polymer nanocomposites formed by blending polymer grafted nanoparticles with an entangled polymer melt is studied by molecular dynamics simulations. We focus on the three key differences in the crazing behavior of a composite relative to the pure homopolymer matrix, namely, a lower yield stress, a smaller extension ratio, and a grafted chain length dependent failure stress. The yield behavior is found to be mostly controlled by the local nanoparticle-grafted polymer interfacial energy, with the grafted polymer-polymer matrix interfacial structure being of little to no relevance. Increasing the attraction between nanoparticle core and the grafted polymer inhibits void nucleation and leads to a higher yield stress. In the craze growth regime, the presence of "grafted chain" sections of ≈100 monomers alters the mechanical response of composite samples, giving rise to smaller extension ratios and higher drawing stresses than for the homopolymer matrix. The dominant failure mechanism of composite samples depends strongly on the length of the grafted chains, with disentanglement being the dominant mechanism for short chains, while bond breaking is the failure mode for chain lengths >10Ne, where Ne is the entanglement length.
ACS Macro Letters
Large-scale molecular dynamics simulations are used to study the dispersion of nanoparticles (NPs) in a polymer film during solvent evaporation. As the solvent evaporates, a dense polymer-rich skin layer forms at the liquid/vapor interface, which is either NP rich or poor depending on the strength of the NP/polymer interaction. When the NPs are strongly wet by the polymer, the NPs accumulate at the interface and form layers. However, when the NPs are only partially wet by the polymer, most NPs are uniformly distributed in the bulk of the polymer film, with the dense skin layer serving as a barrier to prevent the NPs from moving to the interface. Our results point to a possible route to employ less favorable NP/polymer interactions and fast solvent evaporation to uniformly disperse NPs in a polymer film, contrary to the common belief that strong NP/polymer attractions are needed to make NPs well dispersed in polymer nanocomposites.
Abstract not provided.
Abstract not provided.
Physical Review Letters
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.
Abstract not provided.
Macromolecules
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.
Macromolecules
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.
Langmuir
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.
Journal of Polymer Science, Part B: Polymer Physics
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.
Physical Review Letters
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.
Physical Review E - Statistical, Nonlinear, and Soft Matter Physics
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.
Physical Review E - Statistical, Nonlinear, and Soft Matter Physics
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.
Soft Matter
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.
Sandia journal manuscript; Not yet accepted for publication
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 Mw ~ 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.
Faraday Discussions
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
Physical Review Letters
Abstract not provided.
Macromolecular Theory and Simulations
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.
Nature Communications
Abstract not provided.
Macromolecules
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.