Viscoelastic Properties of Polymer-Grafted Nanoparticle Composites from Molecular Dynamics Simulations

Promoting Cross-Link Reaction of Polymers by the Matrix-Filler Interface Effect: Role of Coupling Agents and Intermediate Linkers

The matrix-filler interface effect plays an important role in determining the structural stability and mechanical properties of polymer composite materials, which remain ambiguous and need to be studied. The network-forming dynamics of poly(3,3-bis (azidomethyl) oxetane-tetrahydrofuran) (PBT) at the ammonium perchlorate (AP) surface was studied by using atomistic molecular dynamics simulation, considering the additives of curing agent toluene diisocyanate (TDI), cross-linker trimethylolpropane (TMP), and coupling agent triethanolamine (TEA). The presence of the AP surface promotes chain cross-link reaction, which is attributed to the increased production of intermediate linkers formed by TDI, TMP, and TEA. The intermediate linker has three reactive sites that can react with PBT main chains to form a cross-linked structure. Owing to the strong interaction with the AP surface, the coupling agent TEA plays a dominant role in forming the intermediate linker. At the early stage of network forming (reaction ratio r < 30%), the AP surface adsorbs TEA, which leads to a maximum contact density to PBT. As r increases to 60%, the density of intermediate linkers near the AP surface reaches a maximum value. Consequently, the chain cross-link reactions between the intermediate linker and PBT main chains are enhanced as r > 60%. This work explains the micromechanism of the promotion of chain cross-link reaction by the interface effect and provides important insights on designing polymer materials with high mechanical properties.

Characterizing the shear response of polymer-grafted nanoparticles

Grafting polymer chains to the surface of nanoparticles overcomes the challenge of nanoparticle dispersion within nanocomposites and establishes high-volume fractions that are found to enable enhanced material mechanical properties. This study utilizes coarse-grained molecular dynamics simulations to quantify how the shear modulus of polymer-grafted nanoparticle (PGN) systems in their glassy state depends on parameters such as strain rate, nanoparticle size, grafting density, and chain length. The results are interpreted through further analysis of the dynamics of chain conformations and volume fraction arguments. The volume fraction of nanoparticles is found to be the most influential variable in deciding the shear modulus of PGN systems. A simple rule of mixture is utilized to express the monotonic dependence of shear modulus on the volume fraction of nanoparticles. Due to the reinforcing effect of nanoparticles, shortening the grafted chains results in a higher shear modulus in PGNs, which is not seen in linear systems. These results offer timely insight into calibrating molecular design parameters for achieving the desired mechanical properties in PGNs.

Modeling the viscoelastic relaxation dynamics of soft particles via molecular dynamics simulation-informed multi-dimensional transition-state theory

Viscoelastic soft colloidal particles have been widely explored in mechanical, chemical, pharmaceutical and other engineering applications due to their unique combination of viscosity and elasticity. The characteristic viscoelastic relaxation time shows an Arrhenius-type (or super-Arrhenius due to temperature-dependent transition attempts) thermally-activated behavior, but a holistic explanation from the relevant transition-state theory remains elusive. In this paper, the viscoelastic relaxation times of Lennard-Jones soft colloidal particle systems, including a single particle type system and a binary particle mixture based on the Kob-Andersen model, are determined using molecular dynamics (MD) simulations as the benchmark. First, the particle systems show a non-Maxwellian behavior after comparing the MD-predicted viscoelastic relaxation time and dynamic moduli (storage and loss modulus) to the classic Maxwell viscoelastic model and the recent particle local connectivity theory. Surprisingly, neither the Maxwell relaxation time τ_Maxwell (obtained from the static shear viscosity η and the high-frequency shear modulus G_∞) nor the particle local connectivity lifetime τ_LC can capture the super-Arrhenius temperature-dependent behavior in the MD-predicted relaxation time τ_MD. Then, the particle dissociation and association transition kinetics, fractal dimensions of the particle systems, and neighbor particle structure (obtained from the radial distribution functions) are shown to collectively determine the viscoelastic relaxation time. These factors are embedded into a new multi-dimensional transition kinetics model to directly estimate the viscoelastic relaxation time τ_Model, which is found to agree with the MD-predicted τ_MD remarkably well. This work highlights the microscopic origin of viscoelastic relaxation dynamics of soft colloidal particles, and theoretically connects rheological dynamics and transition kinetics in soft matters.

Tensile and Viscoelastic Behavior in Nacre-Inspired Nanocomposites: A Coarse-Grained Molecular Dynamics Study

Organisms hold an extraordinarily evolutionary advantage in forming complex, hierarchical structures across different length scales that exhibit superior mechanical properties. Mimicking these structures for synthesizing high-performance materials has long held a fascination and has seen rapid growth in the recent past thanks to high-resolution microscopy, design, synthesis, and testing methodologies. Among the class of natural materials, nacre, found in mollusk shells, exhibits remarkably high mechanical strength and toughness. The highly organized "brick and mortar" structure at different length scales is a basis for excellent mechanical properties and the capability to dissipate energy and propagation in nacre. Here, we employ large-scale atomistic coarse-grained molecular dynamics simulations to study the mechanical and viscoelastic behavior of nacre-like microstructures. Uniaxial tension and oscillatory shear simulations were performed to gain insight into the role of complex structure-property relationships. Specifically, the role played by the effect of microstructure (arrangement of the crystalline domain) and polymer-crystal interactions on the mechanical and viscoelastic behavior is elucidated. The tensile property of the nanocomposite was seen to be sensitive to the microstructure, with a staggered arrangement of the crystalline tablets giving rise to a 20-30% higher modulus and lower tensile strength compared to a columnar arrangement. Importantly, the staggered microstructure is shown to have a highly tunable mechanical behavior with respect to the polymer-crystal interactions. The underlying reasons for the mechanical behavior are explained by showing the effect of polymer chain mobility and orientation and the load-carrying capacity for the constituents. Viscoelastic responses in terms of the storage and loss moduli and loss tangent are studied over three decades in frequency and again highlight the differences brought about by the microstructure. We show that our coarse-grained models offer promising insights into the design of novel biomimetic structures for structural applications.

Screening confinement of entanglements: Role of a self-propelling end inducing ballistic chain reptation

Synthetic and natural nanomaterials with self-propelling mechanisms continue to be explored to boost chain mobility beyond normal reptation in the crowded environments of entangled chains. Here we employ scaling theory and numerical simulations to demonstrate that activating one chain end of a singular or isolated chain boosts entanglement-constrained chain reptation from the one-dimensional diffusive mobility as described by the de Gennes-Edwards-Doi model to ballistic motion along the entanglement tube contour. The active chain is effectively screened from the constraint of entanglements on length scales exceeding the tube size.

Tunable Orientation and Assembly of Polymer-Grafted Nanocubes at Fluid-Fluid Interfaces

Self-assembly of faceted nanoparticles is a promising route for fabricating nanomaterials; however, achieving low-dimensional assemblies of particles with tunable orientations is challenging. Here, we demonstrate that trapping surface-functionalized faceted nanoparticles at fluid-fluid interfaces is a viable approach for controlling particle orientation and facilitating their assembly into unique one- and two-dimensional superstructures. Using molecular dynamics simulations of polymer-grafted nanocubes in a polymer bilayer along with a particle-orientation classification method we developed, we show that the nanocubes can be induced into face-up, edge-up, or vertex-up orientations by tuning the graft density and differences in their miscibility with the two polymer layers. The orientational preference of the nanocubes is found to be governed by an interplay between the interfacial area occluded by the particle, the difference in interactions of the grafts with the two layers, and the stretching and intercalation of grafts at the interface. The resulting orientationally constrained nanocubes are then shown to assemble into a variety of unusual architectures, such as rectilinear strings, close-packed sheets, bilayer ribbons, and perforated sheets, which are difficult to obtain using other assembly methods. Our work thus demonstrates a versatile strategy for assembling freestanding arrays of faceted nanoparticles with possible applications in plasmonics, optics, catalysis, and membranes, where precise control over particle orientation and position is required.

Simulation of the Coronal Dynamics of Polymer-Grafted Nanoparticles

Polymer-grafted nanoparticles (PGNPs) are an important component of many advanced materials. The interplay between the nanoparticle surface curvature and spatial confinement by neighboring chains produces a complex set of structural and dynamical behaviors in the polymer corona surrounding the nanoparticle. For example, experiments have shown that the inner portion of the corona is more stretched and relaxes more slowly than the outer region. Here, we perform systematic core-modified dissipative particle dynamics (CM-DPD) simulations and analyze the relaxation dynamics using proper orthogonal decomposition (POD) of the monomer coordinates. We find that grafted chains relax more slowly than free chains and that the relaxation time of the grafted chains scales inversely with the confinement strength. For PGNPs in a polymer melt, the relaxation processes are always Rouse-like. However, we observe either Zimm-like or Rouse-like dynamics for PGNPs in solution depending on the confinement strength.

Fracture in Silica/Butadiene Rubber: A Molecular Dynamics View of Design-Property Relationships

Despite intense investigation, the mechanisms governing the mechanical reinforcement of polymers by dispersed nanoparticles have only been partially clarified. This is especially true for the ultimate properties of the nanocomposites, which depend on their resistance to fracture at large deformations. In this work, we adopt molecular dynamics simulations to investigate the mechanical properties of silica/polybutadiene rubber, using a quasi-atomistic model that allows a meaningful description of bond breaking and fracture over relatively large length scales. The behavior of large nanocomposite models is explored systematically by tuning the cross-linking, grafting densities, and nanoparticle concentration. The simulated stress-strain curves are interpreted by monitoring the breaking of chemical bonds and the formation of voids, up to complete rupture of the systems. We find that some chemical bonds, and particularly the S-S linkages at the rubber-nanoparticle interface, start breaking well before the appearance of macroscopic features of fracture and yield.

Computational Design for the Damping Characteristics of Poly(ether ether ketone)

To investigate the damping characteristics of poly(ether ether ketone) (PEEK), various potential modifications of the molecular structure, including sulfonate groups, hydroxyl groups, amino groups, carboxyl groups, methyl groups, fluorines, and benzene rings, were considered. It was found that these functional groups can mediate both the storage and loss modulus of PEEK derivatives, and the loss factors of PEEK derivatives are sensitive to the content and type of functional groups, indicating an ideal designability of energy dissipation performance. The reciprocating process of H-bonds and C_π-H bonds breaking and reforming during material deformation and the available free volume in the material are critical to the energy dissipation capacities in polymers.

Examining the self-assembly of patchy alkane-grafted silica nanoparticles using molecular simulation

In this work, molecular dynamics simulations are used to examine the self-assembly of anisotropically coated "patchy" nanoparticles. Specifically, we use a coarse-grained model to examine silica nanoparticles coated with alkane chains, where the poles of the grafted nanoparticle are bare, resulting in strongly attractive patches. Through a systematic screening process, the patchy nanoparticles are found to form dispersed, string-like, and aggregated phases, dependent on the combination of alkane chain length, coating chain density, and the fractional coated surface area. Correlation analysis is used to identify the ability of various particle descriptors to predict bulk phase behavior from more computationally efficient single grafted nanoparticle simulations and demonstrates that the solvent-accessible surface area of the nanoparticle core is a key predictor of bulk phase behavior. The results of this work enhance our knowledge of the phase space of patchy nanoparticles and provide a powerful approach for future screening of these materials.

Polymer Nanocomposite-based Coatings for Corrosion Protection

Corrosion of metals induces enormous loss of material performance and increase of cost, which has been a common and intractable issue that needs to be addressed urgently. Coating technology has been acknowledged to be the most economic and efficient approach to retard the metal corrosion. For several decades, polymers have been recognized as an effective anticorrosion coating material in both industries and scientific communities, as they demonstrate good barrier properties, ease of altering properties and massive production. Nanomaterials show distinctively different physical and chemical properties compared with their bulk counterparts, which have been considered as highly promising functional materials in various applications, impacting virtually all the fields of science and technologies. Recently, the introduction of nanomaterials with various properties into polymer matrix to form a polymer nanocomposite has been devoted to improve anticorrosive ability of polymer coatings. In this review article, we highlight the recent advances and synopsis of these high-performance polymer nanocomposites as anticorrosive coating materials. We expect that this work could be helpful for the researchers who are interested in the development of functional nanomaterials and advanced corrosion protection technology.

Viscoelasticity of Nanosheet-Filled Polymer Composites: Three Regimes in the Enhancement of Moduli

We employed nonequilibrium molecular dynamics simulations to elucidate the viscoelastic properties of nanosheet (NS)-filled polymer composites. The effects of NS loadings and NS-polymer interaction on viscoelasticity were examined. The simulation results show that the NS-filled polymer composites exhibit an enhanced storage modulus and loss modulus as the NSs are loaded. There are three regimes of the enhanced process based on the NS loadings. At lower NS loadings, the motion of polymers slows down owing to the interaction between NSs and polymers, and the polymer chains generally follow the Rouse dynamics. As the NS loadings increase, the polymer chains are confined between NSs, leading to a substantial increment in dynamic moduli. At higher NS loadings, a transient network is formed, which strengthens the dynamic moduli further. In addition, the attractive NS-polymer interaction can improve the dispersion of NSs and increase the storage and loss moduli. The present work could provide essential information for designing high-performance hybrid polymeric materials.

A Novel Approach to Atomistic Molecular Dynamics Simulation of Phenolic Resins Using Symthons

Materials science is beginning to adopt computational simulation to eliminate laboratory trial and error campaigns—much like the pharmaceutical industry of 40 years ago. To further computational materials discovery, new methodology must be developed that enables rapid and accurate testing on accessible computational hardware. To this end, the authors utilise a novel methodology concept of intermediate molecules as a starting point, for which they propose the term ‘symthon’ (The term ‘Symthon’ is being used as a simulation equivalent of the synthon, popularised by Dr Stuart Warren in ‘Organic Synthesis: The Disconnection Approach’, OUP: Oxford, 1983.) rather than conventional monomers. The use of symthons eliminates the initial monomer bonding phase, reducing the number of iterations required in the simulation, thereby reducing the runtime. A novel approach to molecular dynamics, with an NVT (Canonical) ensemble and variable unit cell geometry, was used to generate structures with differing physical and thermal properties. Additional script methods were designed and tested, which enabled a high degree of cure in all sampled structures. This simulation has been trialled on large-scale atomistic models of phenolic resins, based on a range of stoichiometric ratios of formaldehyde and phenol. Density and glass transition temperature values were produced, and found to be in good agreement with empirical data and other simulated values in the literature. The runtime of the simulation was a key consideration in script design; cured models can be produced in under 24 h on modest hardware. The use of symthons has been shown as a viable methodology to reduce simulation runtime whilst generating accurate models.

Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

The nature of alterations to dynamics and vitrification in the nanoscale vicinity of interfaces-commonly referred to as "nanoconfinement" effects on the glass transition-has been an open question for a quarter century. We first analyze experimental and simulation results over the last decade to construct an overall phenomenological picture. Key features include the following: after a metrology- and chemistry-dependent onset, near-interface relaxation times obey a fractional power law decoupling relation with bulk relaxation; relaxation times vary in a double-exponential manner with distance from the interface, with an intrinsic dynamical length scale appearing to saturate at low temperatures; the activation barrier and vitrification temperature T_g approach bulk behavior in a spatially exponential manner; and all these behaviors depend quantitatively on the nature of the interface. We demonstrate that the thickness dependence of film-averaged T_g for individual systems provides a poor basis for discrimination between different theories, and thus we assess their merits based on the above dynamical gradient properties. Entropy-based theories appear to exhibit significant inconsistencies with the phenomenology. Diverse free-volume-motivated theories vary in their agreement with observations, with approaches invoking cooperative motion exhibiting the most promise. The elastically cooperative nonlinear Langevin equation theory appears to capture the largest portion of the phenomenology, although important aspects remain to be addressed. A full theoretical understanding requires improved confrontation with simulations and experiments that probe spatially heterogeneous dynamics within the accessible 1-ps to 1-year time window, minimal use of adjustable parameters, and recognition of the rich quantitative dependence on chemistry and interface.

Orientational phase behavior of polymer-grafted nanocubes

Surface functionalization of nanoparticles with polymer grafts was recently shown to be a viable strategy for controlling the relative orientation of shaped nanoparticles in their higher-order assemblies. In this study, we investigated in silico the orientational phase behavior of coplanar polymer-grafted nanocubes confined in a thin film. We first used Monte Carlo simulations to compute the two-particle interaction free-energy landscape of the nanocubes and identify their globally stable configurations. The nanocubes were found to exhibit four stable phases: those with edge-edge and face-face orientations, and those exhibiting partially overlapped slanted and parallel faces previously assumed to be metastable. Moreover, the edge-edge configuration originally thought to involve kissing edges instead displayed partly overlapping edges, where the extent of the overlap depends on the attachment positions of the grafts. We next formulated analytical scaling expressions for the free energies of the identified configurations, which were used for constructing a comprehensive phase diagram of nanocube orientation in a multidimensional parameter space comprising of the size and interaction strength of the nanocubes and the Kuhn length and surface density of the grafts. The morphology of the phase diagram was shown to arise from an interplay between polymer- and surface-mediated interactions, especially differences in their scalings with respect to nanocube size and grafting density across the four phases. The phase diagram provided insights into tuning these interactions through the various parameters of the system for achieving target configurations. Overall, this work provides a framework for predicting and engineering interparticle configurations, with possible applications in plasmonic nanocomposites where control over particle orientation is critical.

Structure-Property Relationships via Recovery Rheology in Viscoelastic Materials

The recoverable strain is shown to correlate to the temporal evolution of microstructure via time-resolved small-angle neutron scattering and dynamic shear rheology. Investigating two distinct polymeric materials of wormlike micelles and fibrin network, we demonstrate that, in addition to the nonlinear structure-property relationships, the shear and normal stress evolution is dictated by the recoverable strain. A distinct sequence of physical processes under large amplitude oscillatory shear (LAOS) is identified that clearly contains information regarding both the steady-state flow curve and the linear-regime frequency sweep, contrary to most interpretations that LAOS responses are either distinct from or somehow intermediate between the two cases. This work provides a physically motivated and straightforward path to further explore the structure-property relationships of viscoelastic materials under dynamic flow conditions.

Interfacial Assembly of Tunable Anisotropic Nanoparticle Architectures

We propose a strategy for assembling spherical nanoparticles (NPs) into anisotropic architectures in a polymer matrix. The approach takes advantage of the interfacial tension between two mutually immiscible polymers forming a bilayer and differences in the compatibility of the two polymer layers with polymer grafts on particles to trap NPs within two-dimensional planes parallel to the interface. The ability to precisely tune the location of the entrapment planes via the NP grafting density, and to trap multiple interacting particles within distinct planes, can then be used to assemble NPs into unconventional arrangements near the interface. We carry out molecular dynamics simulations of polymer-grafted NPs in a polymer bilayer to demonstrate the viability of the proposed approach in both trapping NPs at tunable distances from the interface and assembling them into a variety of unusual nanostructures. We illustrate the assembly of NP clusters, such as dimers with tunable tilt relative to the interface and trimers with tunable bending angle, as well as anisotropic macroscopic phases, including serpentine and branched structures, ridged hexagonal monolayers, and square-ordered bilayers. We also develop a theoretical model to predict the preferred positions and free energies of NPs trapped at or near the interface that could help guide the design of polymer-grafted NPs for achieving target NP architectures. Overall, this work suggests that interfacial assembly of NPs could be a promising approach for fabricating next-generation polymer nanocomposites with potential applications in plasmonics, electronics, optics, and catalysis where precise arrangement of polymer-embedded NPs is required for function.

A Transferable, Multi-Resolution Coarse-Grained Model for Amorphous Silica Nanoparticles

Despite the ubiquity of nanoparticles in modern materials research, computational scientists are often forced to choose between simulations featuring detailed models of only a few nanoparticles or simplified models with many nanoparticles. Herein, we present a coarse-grained model for amorphous silica nanoparticles with parameters derived via potential matching to atomistic nanoparticle data, thus enabling large-scale simulations of realistic models of silica nanoparticles. Interaction parameters are optimized to match a range of nanoparticle diameters in order to increase transferability with nanoparticle size. Analytical functions are determined such that interaction parameters can be obtained for nanoparticles with arbitrary coarse-grained fidelity. The procedure is shown to be extensible to the derivation of cross-interaction parameters between coarse-grained nanoparticles and other moieties and validated for systems of grafted nanoparticles. The optimization procedure used is available as an open-source Python package and should be readily extensible to models of non-silica nanoparticles.

Enhanced Mechanical Properties of Polymer Nanocomposites Using Dopamine-Modified Polymers at Nanoparticle Surfaces in Very Low Molecular Weight Polymers

While incorporation of nanoparticles in a polymer matrix generally enhances the physical properties, effective control of the nanoparticle/polymer interface is often challenging. Here, we report a dramatic enhancement of the mechanical properties of polymer nanocomposites (PNCs) using a simple physical grafting method. The PNC consists of low molecular weight poly(ethylene glycol) (PEG) and silica nanoparticles whose surfaces are modified with dopamine-modified PEG (DOPA-mPEG) brush polymers. With DOPA-mPEG grafting, the nanoparticle surface can be readily altered, and the shear modulus of the PNC is increased by a factor of 10⁵ at an appropriate surface grafting density. The detailed microstructure and mechanical properties are examined with small-angle X-ray scattering (SAXS) and oscillatory rheometry experiments. The attractive interactions between particles induced by DOPA-mPEG grafting dramatically improve the mechanical properties of PNCs even in an unentangled polymer matrix, which shows a much higher shear modulus than that of a highly entangled polymer matrix.

Interfacial Properties and Hopping Diffusion of Small Nanoparticle in Polymer/Nanoparticle Composite with Attractive Interaction on Side Group

The diffusion dynamics of fullerene (C 60 ) in unentangled linear atactic polystyrene (PS) and polypropylene (PP) melts and the structure and dynamic properties of polymers in interface area are investigated by performing all-atom molecular dynamics simulations. The comparison of the results in two systems emphasises the influence of local interactions exerted by polymer side group on the diffusion dynamics of the nanoparticle. In the normal diffusive regime at long time scales, the displacement distribution function (DDF) follows a Gaussian distribution in PP system, indicating a normal diffusion of C 60 . However, we observe multiple peaks in the DDF curve for C 60 diffusing in PS melt, which indicates a diffusion mechanism of hopping of C 60 . The attractive interaction between C 60 and phenyl ring side groups are found to be responsible for the observed hopping diffusion. In addition, we find that the C 60 is dynamically coupled with a subsection of a tetramer on PS chain, which has a similar size with C 60 . The phenyl ring on PS chain backbone tends to have a parallel configuration in the vicinity of C 60 surface, therefore neighbouring phenyl rings can form chelation effect on the C 60 surface. Consequently, the rotational dynamics of phenyl ring and the translational diffusion of styrene monomers are found to be slowed down in this interface area. We hope our results can be helpful for understanding of the influence of the local interactions on the nanoparticle diffusion dynamics and interfacial properties in polymer/nanoparticle composites.

Effects of chemically heterogeneous nanoparticles on polymer dynamics: insights from molecular dynamics simulations

The dispersion of solid nanoparticles within polymeric materials is widely used to enhance their performance. Many scientific and technological aspects of the resulting polymer nanocomposites have been studied, but the role of the structural and chemical heterogeneity of the nanoparticles has just started to be appreciated. For example, simulations of polymer films on planar heterogeneous surfaces revealed unexpected, non-monotonic activation energy to diffusion on varying the surface composition. Motivated by these intriguing results, here we simulate via molecular dynamics a different, fully three-dimensional system, in which the heterogeneous nanoparticles are incorporated in a polymer melt. The nanoparticles are roughly spherical assemblies of strongly and weakly attractive sites, in fractions of f and 1 - f, respectively. We show that the polymer diffusion is still characterized by a non-monotonic dependence of the activation energy on f. The comparison with the case of homogeneous nanoparticles clarifies that the effect of the heterogeneity increases on approaching the polymer glass transition.

Molecular dynamics simulation of the viscoelasticity of polymer nanocomposites under oscillatory shear: effect of interfacial chemical coupling

In this work by adopting coarse-grained molecular dynamics simulation, we focus our attention on investigating the effect of the chemical coupling between polymer and nanoparticles (NPs) on the viscoelastic properties of polymer nanocomposites (PNCs). Firstly we examine the effect of the interfacial chemical coupling on the non-linear behavior, such as the change of the storage moduli, the loss moduli and the loss factor as a function of the strain amplitude. Besides the reinforcing effect contributed by the interfacial chemical interaction, a much smaller loss factor is also observed attributed to less molecular friction and dissipation. Meanwhile, the effects of temperature, frequency, and the interfacial physical interaction between NPs and polymers on the viscoelastic properties are also probed. To uncover the structural and dynamic effect of the interfacial chemical coupling, we calculate the radial distribution function of polymer chains around NPs, the content of the polymer beads in the first layer of the interfacial region under quiescent and dynamic conditions, the incoherent intermediate dynamic structure factor of the polymer beads, which are chemically or physically tethered to the NPs, and all the polymer beads of the system, the quantitative comparison of the mean relaxation time for different interfacial chemical coupling, and the mean-square displacement of the polymer chains. Lastly we analyze the change of the interfacial energy such as the physical and chemical energies during oscillatory shear. Through these analyses, we conclude that with the increase of the interfacial chemical coupling, the change extent of the interfacial physical interaction versus the periodic strain decreases, attributed to a much smaller adsorption–desorption reversible process. This can rationalize the much weaker non-linear behavior or the “Payne effect”. Based on these results, we anticipate that a better molecular-level understanding is provided on the effect of the interfacial coupling on the viscoelastic properties of PNCs.

In this work by adopting coarse-grained molecular dynamics simulation, we focus attention on investigating the effect of the chemical coupling between polymer and nanoparticles (NPs) on the viscoelastic properties of polymer nanocomposites (PNCs).

Modeling the Effect of Polymer Chain Stiffness on the Behavior of Polymer Nanocomposites

Due to their central role in industrial formulations spanning from food packaging to smart coatings, polymer nanocomposites have been the object of remarkable attention over the last two decades. Incorporating nanoparticles (NPs) into a polymer matrix modifies the conformation and mobility of the polymer chains at the NP-polymer interface and can potentially provide materials with enhanced properties as compared to pristine polymers. To this end, it is crucial to predict and control the ability of NPs to diffuse and achieve a good dispersion in the polymer matrix. Understanding how to control the NPs' dispersion is a challenging task controlled by the delicate balance between enthalpic and entropic contributions, such as NP-polymer interaction, NP size and shape, and polymer chain conformation. By performing molecular dynamics (MD) simulations, we investigate the effect of polymer chains' stiffness on the mobility of spherical NPs that establish weak or strong interactions with the polymer. Our results show a sound dependence of the NPs' diffusivity on the long-range order of the polymer melt, which undergoes an isotropic-to-nematic phase transition upon increasing chain stiffness. This phase transition induces a dynamical anisotropy in the nematic phase, with the NPs preferentially diffusing along the nematic director rather than in the directions perpendicular to it. Not only does this tendency determine the NPs' mobility and degree of dispersion in the polymer matrix, but it also influences the resistance to flow of the polymer nanocomposite when a shear is applied. In particular, to assess the role of the chains' conformation on the macroscopic response of our model PNC, we employ reverse nonequilibrium MD to calculate the zero-shear viscosity in both the isotropic and nematic phases, and unveil a plasticizing effect at increasing chain stiffness when the shear is applied along the nematic axis.

Multiscale Molecular Simulations of Polymer-Matrix Nanocomposites: or What Molecular Simulations Have Taught us About the Fascinating Nanoworld

Following the substantial progress in molecular simulations of polymer-matrix nanocomposites, now is the time to reconsider this topic from a critical point of view. A comprehensive survey is reported herein providing an overview of classical molecular simulations, reviewing their major achievements in modeling polymer matrix nanocomposites, and identifying several open challenges. Molecular simulations at multiple length and time scales, working hand-in-hand with sensitive experiments, have enhanced our understanding of how nanofillers alter the structure, dynamics, thermodynamics, rheology and mechanical properties of the surrounding polymer matrices.

Fish DNA-modified clays: Towards highly flame retardant polymer nanocomposite with improved interfacial and mechanical performance

Deoxyribonucleic Acid (DNA) has been recently found to be an efficient renewable and environmentally-friendly flame retardant. In this work, for the first time, we have used waste DNA from fishing industry to modify clay structure in order to increase the clay interactions with epoxy resin and take benefit of its additional thermal property effect on thermo-physical properties of epoxy-clay nanocomposites. Intercalation of DNA within the clay layers was accomplished in a one-step approach confirmed by FT-IR, XPS, TGA, and XRD analyses, indicating that d-space of clay layers was expanded from ~1.2 nm for pristine clay to ~1.9 nm for clay modified with DNA (d-clay). Compared to epoxy nanocomposite containing 2.5%wt of Nanomer I.28E organoclay (m-clay), it was found that at 2.5%wt d-clay loading, significant enhancements of ~14%, ~6% and ~26% in tensile strength, tensile modulus, and fracture toughness of epoxy nanocomposite can be achieved, respectively. Effect of DNA as clay modifier on thermal performance of epoxy nanocomposite containing 2.5%wt d-clay was evaluated using TGA and cone calorimetry analysis, revealing significant decreases of ~4000 kJ/m² and ~78 kW/m² in total heat release and peak of heat release rate, respectively, in comparison to that containing 2.5%wt of m-clay.

Influence of Morphology on the Mechanical Properties of Polymer Nanocomposites Filled with Uniform or Patchy Nanoparticles

In this work we perform molecular-dynamics simulations, both on the coarse-grained and the chemistry-specific levels, to study the influence of morphology on the mechanical properties of polymer nanocomposites (PNCs) filled with uniform spherical nanoparticles (which means without chemical modification) and patchy spherical nanoparticles (with discrete, attractive interaction sites at prescribed locations on the particle surface). Through the coarse-grained model, the nonlinear decrease of the elastic modulus (G') and the maximum of the viscous modulus (G″) around the shear strain of 10% is clearly reproduced. By turning to the polybutadiene model, we examine the effect of the shear amplitude and the interaction strength among uniform NPs on the aggregation kinetics. Interestingly, the change of the G' as a function of the aggregation time exhibited a maximum value at intermediate time attributed to the formation of a polymer-bridged filler network in the case of strong interaction between NPs. By imposing a dynamic periodic shear, we probe the change of the G' as a function of the strain amplitude while varying the interaction strength between uniform NPs and its weight fraction. A continuous filler network is developed at a moderate shear amplitude, which is critically related to the interaction strength between NPs and the weight fraction of the fillers. In addition, we study the self-assembly of the patchy NPs, which form the typical chain-like and sheet-like structures. For the first time, the effect of these self-assembled structures on the viscoelastic and stress-strain behavior of PNCs is compared. In general, in the coarse-grained model we focus on the size effect of the rough NPs on the Payne effect, while some other parameters such as the dynamic shear flow, the interaction strength between NPs, the weight fraction, and the chemically heterogeneous surface of the NPs are explored for the chemistry-specific model.

Simulated Permeation and Characterization of PEGylated Gold Nanoparticles in a Lipid Bilayer System

PEGylated gold nanoparticles are considered suitable nanocarriers for use in biomedical applications and targeted drug delivery systems. In our previous investigation with the alkanethiol-functionalized gold nanoparticle, we found that permeation across a protein-free phospholipid membrane resulted in damaging effects of lipid displacement and water and ion leakage. In the present study, we carry out a series of coarse-grained molecular simulations to explore permeation of lipid bilayer systems by a PEGylated gold nanoparticle, especially at the bulk-liquid-lipid interface as well as the interface between the two lipid leaflets. Initially, we examine molecular-level details of a PEGylated gold nanoparticle (constructed from cycled annealing) in water and find a distribution of ligand configurations (from mushroom to brush states) present in nanoparticles with medium to high surface coverage. We also find that the characteristic properties of the PEGylated gold nanoparticle do not change when it is placed in a salt solution. In our permeation studies, we investigate events of water and ion penetration as well as lipid translocation while varying the ligand length, nanoparticle surface coverage, and ion concentration gradient of our system. Results from our studies show the following: (1) The number of water molecules in the interior of the membrane during ligand-coated nanoparticle permeation increases with PEGn-SH surface coverage, ligand length, and permeation velocity but is not sensitive to the ion concentration gradient. (2) Lipid molecules do not leave the membrane; instead they complete trans-bilayer lipid flip-flop with longer ligands and higher surface coverages. (3) The lack of formation of stable water pores prevents ion translocation. (4) The PEGylated nanoparticle causes less damage to the membrane overall due to favorable interactions with the lipid headgroups which may explain why experimentalists observe endocytosis of PEGylated nanocarriers in vivo.

Entanglements in polymer nanocomposites containing spherical nanoparticles

We investigate the polymer packing around nanoparticles and polymer/nanoparticle topological constraints (entanglements) in nanocomposites containing spherical nanoparticles in comparison to pure polymer melts using molecular dynamics (MD) simulations. The polymer-nanoparticle attraction leads to good dispersion of nanoparticles. We observe an increase in the number of topological constraints (decrease of total entanglement length Ne with nanoparticle loading in the polymer matrix) in nanocomposites due to nanoparticles, as evidenced by larger contour lengths of the primitive paths. An increase of the nanoparticle radius reduces the polymer-particle entanglements. These studies demonstrate that the interaction between polymers and nanoparticles does not affect the total entanglement length because in nanocomposites with small nanoparticles, the polymer-nanoparticles topological constraints dominate.

Computationally Guided Assembly of Oriented Nanocubes by Modulating Grafted Polymer-Surface Interactions

The bottom-up fabrication of ordered and oriented colloidal nanoparticle assemblies is critical for engineering functional nanomaterials beyond conventional polymer-particle composites. Here, we probe the influence of polymer surface ligands on the self-orientation of shaped metal nanoparticles for the formation of nanojunctions. We examine how polymer graft-surface interactions dictate Ag nanocube orientation into either edge-edge or face-face nanojunctions. Specifically, we investigate the effect of end-functionalized polymer grafts on nanocube assembly outcomes, such as interparticle angle and interparticle distance. Our assembly results can be directly mapped onto our theoretical phase diagrams for nanocube orientation, enabling correlation of experimental variables (such as graft length and metal binding strength) with computational parameters. These results represent an important step toward unifying modeling and experimental approaches to understanding nanoparticle-polymer self-assembly.