Nanoparticle Geometrical Effect on Structure, Dynamics and Anisotropic Viscosity of Polyethylene Nanocomposites

The mechanical property and microscopic deformation mechanism of nanoparticle-contained graphene foam materials under uniaxial compression

Nanoparticle-contained graphene foams have found more and more practical applications in recent years, which desperately requires a deep understanding on basic mechanics of this hybrid material. In this paper, the microscopic deformation mechanism and mechanical properties of such a hybrid material under uniaxial compression, that are inevitably encountered in applications and further affect its functions, are systematically studied by the coarse-grained molecular dynamics simulation method. Two major factors of the size and volume fraction of nanoparticles are considered. It is found that the constitutive relation of nanoparticle filled graphene foam materials consists of three parts: the elastic deformation stage, deformation with inner re-organization and the final compaction stage, which is much similar to the experimental measurement of pristine graphene foam materials. Interestingly, both the initial and intermediate modulus of such a hybrid material is significantly affected by the size and volume fraction of nanoparticles, due to their influences on the microstructural evolution. The experimentally observed 'spacer effect' of such a hybrid material is well re-produced and further found to be particle-size sensitive. With the increase of nanoparticle size, the micro deformation mechanism will change from nanoparticles trapped in the graphene sheet, slipping on the graphene sheet, to aggregation outside the graphene sheet. Beyond a critical relative particle size 0.26, the graphene-sheet-dominated deformation mode changes to be a nanoparticle-dominated one. The final microstructure after compression of the hybrid system converges to two stable configurations of the 'sandwiched' and 'randomly-stacked' one. The results should be helpful not only to understand the micro mechanism of such a hybrid material in different applications, but also to the design of advanced composites and devices based on porous materials mixed with particles.

Local Structural Fingerprints of Nanoparticle-Bound Polymer Layers

The molecular structure of bound layers at attractive polymer-nanoparticle interfaces strongly influences the properties of nanocomposites. Thus, a unifying theoretical framework that can provide insights into the correlations between the molecular structure of the bound layers, their thermodynamics, and macroscopic properties is highly desirable. In this work, molecular dynamics simulations were used in combination with local fingerprint analysis of configurational entropy and interaction energy at the segmental scale, with the goal to establish such physical grounds. The thickness of bound polymer layers is found to be independent of the polymer chain length, as deduced from density oscillations at the surface of a nanotube. The local configurational entropy of layers is estimated from pair correlations in equilibrium structures. By plotting mean layer entropy vs internal energy on a phase diagram, a one-to-one equivalence is established between the local structures of layers and their thermodynamic properties. Moreover, a gradient in local dynamics of segments in bound layers under equilibrium is observed normal to the nanoparticle surface. The relaxation times of individual layers show correspondence to their phase diagram fingerprints, thus suggesting that a unified perspective can be envisioned for such materials built on the grounds of locally heterogeneous interfaces.

Designing high thermal conductivity of cross-linked epoxy resin via molecular dynamics simulations

Coarse-grained (CG) non-equilibrium molecular dynamics simulation was used to study the thermal conductivity of a cross-linked network composed of epoxy resin (E51) and polyether amine (PEA). By probing the mechanism of heat transfer in the cross-linked epoxy resin, we systematically explored the effects of the crosslinking degree, chain length and multi-functional groups of the curing agent on the thermal conductivity behavior. Our results indicate that the thermal conductivity is mainly dependent on the chain length and the functional groups of the curing agent. A shorter chain length and a curing agent with more functional groups contribute to higher thermal conductivity, while the crosslinking degree has a negligible effect. Moreover, it is revealed that the thermal conductivity is manipulated by the non-bonding interaction energy (E_pair) and the vibrational density. In general, our work could provide some guidelines for the design and fabrication of a cross-linked epoxy network with high thermal conductivity.

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.

Graphene nanosheets as reinforcement and cell-instructive material in soft tissue scaffolds

Mechanical strength of polymeric scaffolds deteriorates quickly in the physiological mileu. This can be minimized by reinforcing the polymeric matrix with graphene, a planar two-dimensional material with unique physicochemical and biological properties. Association between the sheet and polymer chains offers a range of porosity commensurate with tissue requirements. Besides, studies suggest that corrugated structure of graphene offers desirable bio-mechanical cues for tissue regeneration. This review covers three important aspects of graphene-polymer composites, (a) the opportunity on reinforcing the polymer matrix with graphene, (b) challenges associated with limited aqueous processability of graphene, and (c) physiological signaling in the presence of graphene. Among numerous graphene materials, our discussion is limited to graphene oxide (GO) and reduced graphene oxide (rGO) nanosheets. Challenges associated with limited dispersity of hydrophobic sheets within the polymeric matrix have been discussed at molecular level.

Slip-spring simulations of different constraint release environments for linear polymer chains

The constraint release (CR) mechanism has important effects on polymer relaxation and the chains will show different relaxation behaviour in conditions of monodisperse, bidisperse and other topological environments. By comparing relaxation data of linear polyisoprene (PI) chains dissolved in very long matrix and monodisperse melts, Matsumiya et al. showed that CR mechanism accelerates both dielectric and viscoelastic relaxation (Matsumiya et al. 2013 Macromolecules 46, 6067. (doi:10.1021/ma400606n)). In this work, the experimental data reported by Matsumiya et al. are reproduced using the single slip-spring (SSp) model and the CR accelerating effects on both dielectric and viscoelastic relaxation are validated by simulations. This effect on viscoelastic relaxation is more pronounced. The coincidence for end-to-end relaxation and the viscoelastic relaxation has also been checked using probe version SSp model. A variant of SSp with each entanglement assigning a characteristic lifetime is also proposed to simulate various CR environment flexibly. Using this lifetime version SSp model, the correct relaxation function can be obtained with equal numbers of entanglement destructions by CR and reptation/contour length fluctuation (CLF) for monodisperse melts. Good agreement with published experiment data is also obtained for bidisperse melts, which validates the ability to correctly describe the CR environment of the lifetime version model.

Modeling of Entangled Polymer Diffusion in Melts and Nanocomposites: A Review

This review concerns modeling studies of the fundamental problem of entangled (reptational) homopolymer diffusion in melts and nanocomposite materials in comparison to experiments. In polymer melts, the developed united atom and multibead spring models predict an exponent of the molecular weight dependence to the polymer diffusion very similar to experiments and the tube reptation model. There are rather unexplored parameters that can influence polymer diffusion such as polymer semiflexibility or polydispersity, leading to a different exponent. Models with soft potentials or slip-springs can estimate accurately the tube model predictions in polymer melts enabling us to reach larger length scales and simulate well entangled polymers. However, in polymer nanocomposites, reptational polymer diffusion is more complicated due to nanoparticle fillers size, loading, geometry and polymer-nanoparticle interactions.

Designing Superlattice Structure via Self-Assembly of One-Component Polymer-Grafted Nanoparticles

The control of the self-assembly of the nanocrystals into ordered structures has been extensively investigated, but fewer efforts have been devoted to studying one-component polymer-grafted nanoparticles (OPNPs). Herein, through coarse-grained molecular dynamics simulation, we design a novel nanoparticle (NP) grafted with polymer chains, focusing on its self-assembled structures. First, we examine the effects of length and density of grafted polymer chains by calculating the radial distribution function between NPs, as well as through direct visualization. We observe a monotonic change of the arranged morphology of grafted-NPs as a function of the density of grafted polymer chains, which indicates that the increase of the grafting density contributes to the order of the morphology. Meanwhile, we find that much longer grafted polymer chains worsen the regularity of the morphology. Then, we probe the influence of the stiffness of grafted polymer chains (denoted by K ranging from 0 to 500) on the order of grafted-NPs, finding that the order of the structure exhibits a nonmonotonic behavior as a function of K at moderate grafting density. For high grafting density, the order of the morphology is initially enhanced and becomes saturated as a function of K. For the effect of K on the stress-strain behavior, the system with the lowest order demonstrates the most remarkable reinforced mechanical behavior for both low and high grafting density. Last, we establish the phase diagram by varying the stiffness and density of the grafted polymer chains, which contains the amorphous, ordered, and superlattice structures, respectively. In general, our simulated results provide guidelines to tailor the self-assembly of the OPNPs by taking advantage of the length, density, and stiffness of grafted polymer chains.

Molecular dynamics simulation of thermo-mechanical behaviour of elastomer cross-linked via multifunctional zwitterions

Coarse-grained (CG) molecular dynamics simulations have been employed to study the thermo-mechanical response of a physically cross-linked network composed of zwitterionic moieties and fully flexible elastomeric polymer chains.

Morphology, Thermal, Mechanical Properties and Rheological Behavior of Biodegradable Poly(butylene succinate)/poly(lactic acid) In-Situ Submicrofibrillar Composites

In this study, biodegradable poly(butylene succinate)/poly(lactic acid) (PBS/PLA) in-situ submicrofibrillar composites with various PLA content were successfully produced by a triple-screw extruder followed by a hot stretching−cold drawing−compression molding process. This study aimed to investigate the effects of dispersed PLA submicro-fibrils on the thermal, mechanical and rheological properties of PBS/PLA composites. Morphological observations demonstrated that the PLA phases are fibrillated to submicro-fibrils in the PBS/PLA composites, and all the PLA submicro-fibrils produced seem to have a uniform diameter of about 200nm. As rheological measurements revealed, at low frequencies, the storage modulus (G’) of PBS/PLA composites has been increased by more than four orders of magnitude with the inclusion of high concentrations (15 wt % and 20 wt %) of PLA submicro-fibrils, which indicates a significant improvement in the elastic responses of PBS melt. Dynamic Mechanical Analysis (DMA) results showed that the glass transition temperature (T_g) of PBS phase slightly shifted to the higher temperature after the inclusion of PLA. DSC experiments proved that fiber morphology of PLA has obvious heterogeneous nucleation effect on the crystallization of PBS. The tensile properties of the PBS/PLA in-situ submicrofibrillar composites are also improved compared to neat PBS.

Designing the Slide-Ring Polymer Network with both Good Mechanical and Damping Properties via Molecular Dynamics Simulation

Through coarse-grained molecular dynamics simulation, we have successfully designed the chemically cross-linked (fixed junction) and the slide-ring (SR) systems. Firstly, we examine the dynamic properties such as the mean-square displacement, the bond, and the end-to-end autocorrelation functions as a function of the cross-linking density, consistently pointing out that the SR system exhibits much lower mobility compared with the fixed junction one at the same cross-linking density. This is further validated by a relatively higher glass transition temperature for the SR system compared with that of the fixed junction one. Then, we calculated the effect of the cross-linking density on the stretch-recovery behavior for the SR and fixed junction systems. Although the chain orientation of the SR system is higher than that of the fixed-junction system, the tensile stress is smaller than the latter. We infer that much greater chain sliding can occur during the stretch, because the movable ring structure homogeneously sustains the external force of the SR system, which, therefore, leads to much larger permanent set and higher hysteresis during the recovery process compared with the fixed-junction one. Based on the stretch-recovery behavior for various cross-linking densities, we obtain the change of the hysteresis loss, which is larger for the SR system than that of the fixed junction system. Lastly, we note that the relatively bigger compressive stress for the SR system results from the aggregation of the rigid rings compared with the fixed junction system. In general, compared with the traditionally cross-linked system, a deep molecular-level insight into the slide-ring polymer network is offered and thus is believed to provide some guidance to the design and preparation of the slide-ring polymer network with both good mechanical and damping properties.

Effect of copolymer sequence on structure and relaxation times near a nanoparticle surface

We simulate a simple nanocomposite consisting of a single spherical nanoparticle surrounded by coarse-grained polymer chains. The polymers are composed of two different monomer types that differ only in their interaction strengths with the nanoparticle. We examine the effect of adjusting copolymer sequence on the structure as well as the end-to-end vector autocorrelation, bond vector autocorrelation, and self-intermediate scattering function relaxation times as a function of distance from the nanoparticle surface. We show how the range and magnitude of the interphase of slowed dynamics surrounding the nanoparticle depend strongly on sequence blockiness. We find that, depending on block length, blocky copolymers can have faster or slower dynamics than a random copolymer. Certain blocky copolymer sequences lead to relaxation times near the nanoparticle surface that are slower than those of either homopolymer system. Thus, tuning copolymer sequence could allow for significant control over the nanocomposite behavior.

Transport Diffusion of Light Gases in Polyethylene Using Atomistic Simulations

We explore the temperature dependence of the self-, corrected-, and transport-diffusivities of CO₂, CH₄, and N₂ in a polyethylene (PE) polymer membrane through equilibrium molecular dynamics simulations. We also investigate the morphology of the polymer membrane based on the intermolecular radial distribution function, free volume, and pore size distribution analysis. The results indicate the existence of 1.5-3 Å diameter pores in the PE membrane, and with the increase in the temperature, the polymer swells linearly with changing slope at 450 K in the absence of gas and exponentially in the presence of gas. The gas adsorption isotherms extracted via a two-step methodology, considering the dynamics and structural transitions in the polymer matrix upon gas adsorption, were fitted using a "two-mode sorption" model. Our results suggest that CO₂ adsorbs strongly, whereas N₂ shows weak adsorption in PE. The results demonstrate that CO₂ is more soluble, whereas N₂ is least soluble. Further, it is found that an increase in the temperature negatively impacts the solubility of CO₂ and CH₄ but positively for N₂; this reverse solubility behavior is due to increased availability of pores accessible to N₂, which are kinetically closed at the lowest temperatures. The reported self-diffusivities of the gases from our simulations are on the order of 10^-6 cm²/s, consistent with the experimental evidence, whereas transport-diffusivities are 2 orders of magnitude higher than self-diffusivities. Furthermore, the temperature dependence of the self-diffusivity follows Arrhenius behavior, whereas the transport-diffusivity follows non-Arrhenius behavior having different activation energies in low and high temperature regions. Also, it is seen that loading has little effect on the self- and corrected-diffusion coefficients of all gases in the PE membrane.

Molecular insights into the effect of graphene packing on mechanical behaviors of graphene reinforced cis-1,4-polybutadiene polymer nanocomposites

We adopt molecular dynamics simulation to study the graphene packing patterns on chain structure, dynamics, uniaxial tension and visco-elastic behaviors.

Epitaxial crystallization of precisely bromine-substituted polyethylene induced by carbon nanotubes and graphene

Epitaxial crystallization of precisely bromine-substituted polyethylene induced by carbon nanotubes and graphene.

Tailoring the dispersion of nanoparticles and the mechanical behavior of polymer nanocomposites by designing the chain architecture

The dispersion of nanoparticles with different polymer–nanoparticle interaction strengths and chain architectures.

Multiscale simulations of PS-SiO2 nanocomposites: from melt to glassy state

The interaction energetics, molecular packing, entanglement network properties, segmental dynamics, and elastic constants of atactic polystyrene-amorphous silica nanocomposites in the molten and the glassy state are studied via molecular simulations using two interconnected levels of representation: (a) a coarse-grained one, wherein each polystyrene repeat unit is mapped onto a single "superatom" and the silica nanoparticle is viewed as a solid sphere. Equilibration at all length scales at this level is achieved via connectivity-altering Monte Carlo simulations. (b) A united-atom (UA) level, wherein the polymer chains are represented in terms of a united-atom forcefield and the silica nanoparticle is represented in terms of a simplified, fully atomistic model. Initial configurations for UA molecular dynamics (MD) simulations are obtained by reverse mapping well-equilibrated coarse-grained configurations. By analysing microcanonical UA MD trajectories, the polymer density profile is studied and the polymer is found to exhibit layering in the vicinity of the nanoparticle surface. An estimate of the enthalpy of mixing between polymer and nanoparticles, derived from the UA simulations, compares favourably against available experimental values. The dynamical behaviour of polystyrene (in neat and filled melt systems) is characterized in terms of bond orientation and dihedral angle time autocorrelation functions. At low concentration in the molten polymer matrix, silica nanoparticles are found to cause a slight deceleration of the segmental dynamics close to their surface compared to the bulk polymer. Well-equilibrated coarse-grained long-chain configurations are reduced to entanglement networks via topological analysis with the CReTA algorithm, yielding a slightly lower density of entanglements in the filled than in the neat systems. UA melt configurations are glassified by MD cooling. The elastic moduli of the resulting glassy nanocomposites are computed through an analysis of strain fluctuations in the undeformed state and through explicit mechanical deformation by MD, showing a stiffening of the polymer in the presence of nanoparticles. UA simulation results for the elastic constants are compared to continuum micromechanical calculations invoked in homogenization models of the overall mechanical behaviour of heterogeneous materials. They can be interpreted in terms of the presence of an "interphase" of approximate thickness 2 nm around the nanoparticles, with elastic constants intermediate between those of the filler and the matrix.

A Comparison of the Elastic Properties of Graphene- and Fullerene-Reinforced Polymer Composites: The Role of Filler Morphology and Size

Nanoscale carbon-based fillers are known to significantly alter the mechanical and electrical properties of polymers even at relatively low loadings. We report results from extensive molecular-dynamics simulations of mechanical testing of model polymer (high-density polyethylene) nanocomposites reinforced by nanocarbon fillers consisting of graphene flakes and fullerenes. By systematically varying filler concentration, morphology, and size, we identify clear trends in composite stiffness with reinforcement. To within statistical error, spherical fullerenes provide a nearly size-independent level of reinforcement. In contrast, two-dimensional graphene flakes induce a strongly size-dependent response: we find that flakes with radii in the 2–4 nm range provide appreciable enhancement in stiffness, which scales linearly with flake radius. Thus, with flakes approaching typical experimental sizes (~0.1–1 μm), we expect graphene fillers to provide substantial reinforcement, which also is much greater than what could be achieved with fullerene fillers. We identify the atomic-scale features responsible for this size- and morphology-dependent response, notably, ordering and densification of polymer chains at the filler–matrix interface, thereby providing insights into avenues for further control and enhancement of the mechanical properties of polymer nanocomposites.

Reversible wrinkles of monolayer graphene on a polymer substrate: toward stretchable and flexible electronics

The surface instability of monolayer graphene supported by a soft (polymer) substrate under equal-biaxial compression has been explored through large scale coarse-grained molecular simulations. Regardless of the interfacial adhesion strength between the graphene and the substrate, herringbone wrinkles have always been observed due to their lowest energy status, compared with the checkerboard, hexagonal, triangular and one dimensional sinusoidal modes. Moreover, the graphene-polymer substrate interaction energy has a negligible effect on the critical strain for the onset of these wrinkles. Yet, if the graphene is bonded to a rigid (non-deformable) substrate, the critical strain increases with increasing graphene-substrate interfacial strength. The surface wrinkles of graphene are delayed and suppressed by the strong bonding of graphene to the rigid substrate. Besides, only localized folds and crumples have been observed on the surface of graphene, when graphene-substrate interaction energy is strong enough. All these observations signal that the deformability (stiffness) of the substrate plays an essential role in determining the morphology of supported graphene under compression. In addition, when a flat graphene is attached on a highly pre-strained (50%) polymer substrate, wrinkles will be formed on its surface during the relaxation of pre-strain within the polymer substrate. The wrinkled graphene could be stretched up to 50% without fracture, accompanied by the diminishing of surface wrinkles. Therefore, it opens a new avenue to enhance the stretchability of graphene materials, and enables the future applications of graphene and other 2D materials in stretchable and flexible electronics.

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.

Effect of nanoparticles on fibril formation and mechanical performance of olefinic block copolymer (OBC)/polypropylene (PP) microfibrillar composites

Adding nanoparticles provides an alternative strategy to tune the morphological evolution and mechanical performances of OBC/PP microfibrillar composites.

Poly(lactic acid)-Based in Situ Microfibrillar Composites with Enhanced Crystallization Kinetics, Mechanical Properties, Rheological Behavior, and Foaming Ability

Melt blending is one of the most promising techniques for eliminating poly(lactic acid)'s (PLA) numerous drawbacks. However, success in a typical melt blending process is usually achieved through the inclusion of high concentrations of a second polymeric phase which can compromise PLA's green nature. In a pioneering study, we introduce the production of in situ microfibrillar PLA/polyamide-6 (PA6) blends as a cost-effective and efficient technique for improving PLA's properties while minimizing the required PA6 content. Predominantly biobased products, with only 3 wt % of in situ generated PA6 microfibrils (diameter ≈200 nm), were shown to have dramatically improved crystallization kinetics, mechanical properties, melt elasticity and strength, and foaming-ability compared with PLA. Crucially, the microfibrillar blends were produced using an environmentally friendly and cost-effective process. Both of these qualities are essential in guarantying the viability of the proposed technique for overcoming the obstacles associated with the vast commercialization of PLA.

Pulling-force-induced elongation and alignment effects on entanglement and knotting characteristics of linear polymers in a melt

We employ a primitive path (PP) algorithm and the Gauss linking integral to study the degree of entanglement and knotting characteristics of linear polymer model chains in a melt under the action of a constant pulling force applied to selected chain ends. Our results for the amount of entanglement, the linking number, the average crossing number, the writhe of the chains and their PPs and the writhe of the entanglement strands all suggest a different response at the length scale of entanglement strands than that of the chains themselves and of the corresponding PPs. Our findings indicate that the chains first stretch at the level of entanglement strands and next the PP (tube) gets oriented with the "flow." These two phases of the extension and alignment of the chains coincide with two phases related to the disentanglement of the chains. Soon after the onset of external force the PPs attain a more entangled conformation, and the number of nontrivially linked end-to-end closed chains increases. Next, the chains disentangle continuously to attain an almost unentangled conformation. Using the linking matrix of the chains in the melt, we furthermore show that these phases are accompanied by a different scaling of the homogeneity of the global entanglement in the system. The homogeneity of the end-to-end closed chains first increases to a maximum and then decreases slowly to a value characterizing a completely unlinked system.