How Thick Is the Polymer Interphase in Nanocomposites? Probing It by Local Stress Anisotropy and Gas Solubility

Molecular Dynamic and Dissipative Particle Dynamic Simulation on the Miscibility of NR/CR Blends

Natural rubber (NR) exhibits good elasticity, flexural resistance, wear resistance, and excellent mechanical properties, and it has been widely used in aerospace, transportation, medical, and health fields. For NR, however, the resistance to thermal-oxidation and ozone aging is fairly poor. Although aging properties of NR can be significantly improved with the incorporation of chloroprene rubber (CR) according to some references, the miscibility between NR and CR, the morphologies of the binary blends, and so on are revealed ambiguously. In this work, molecular dynamics simulation (MD) and dissipative particle dynamics (DPD) simulation were carried out to predict the compatibility between natural rubber and chloroprene rubber in view of Flory-Huggins parameters. The morphologies of the blends were obtained with the use of the DPD method. The simulation results were furtherly examined by means of Fourier transform infrared spectroscopy (FT-IR) and dynamic mechanical analysis (DMA). It was found that the miscibility between NR and CR is poor. Nevertheless, the miscibility could be improved when the content of CR is 50% or 90%. In addition, spinodal decomposition with a critical temperature of 390 K would take place according to the phase diagram. Microphase structure such as spherical, lamellar, and bicontinuous phases can be found with different contents of CR in the blends with the results of morphologies analysis.

Atomistic-scale analysis of the deformation and failure of polypropylene composites reinforced by functionalized silica nanoparticles

Interfacial adhesion between polymer matrix and reinforcing silica nanoparticles plays an important role in strengthening polypropylene (PP) composite. To improve the adhesion strength, the surface of silica nanoparticles can be modified by grafted functional molecules. Using atomistic simulations, we examined the effect of functionalization of silica nanoparticles by hexamethyldisilazane (HMDS) and octyltriethoxysilane (OTES) molecules on the deformation and failure of silica-reinforced PP composite. We found that the ultimate tensile strength (UTS) of PP composite functionalized by OTES (28 MPa) is higher than that of HMDS (25 MPa), which is in turn higher than that passivated only by hydrogen (22 MPa). To understand the underlying mechanistic origin, we calculated the adhesive energy and interfacial strength of the interphase region, and found that both the adhesive energy and interfacial strength are the highest for the silica nanoparticles functionalized by OTES molecules, while both are the lowest by hydrogen. The ultimate failure of the polymer composite is initiated by the cavitation in the interphase region with the lowest mass density, and this cavitation failure mode is common for all the examined PP composites, but the cavitation position is dependent on the tail length of the functional molecules. The present work provides interesting insights into the deformation and cavitation failure mechanisms of the silica-reinforced PP composites, and the findings can be used as useful guidelines in selecting chemical agents for surface treatment of silica nanoparticles.

Thermoplastic bio-nanocomposites: From measurement of fundamental properties to practical application

Although the discovery of plastic has revolutionized materials used in many industries and by consumers, their non-biodegradable nature has led to one of the greatest problems of our times: plastic waste in the environment. Bioplastics which are biobased and biodegradable, have been suggested as alternatives for their fossil based counterparts, but their properties often do not meet the requirements that standard plastics would, and are in clear need of improvement. One way to do so is by the addition of nanoparticles which, when homogeneously dispersed, have been reported to result in great improvements. However, in practice, homogenous distribution of nanoparticles is not that trivial due to their tendency to aggregate, also after addition to the polymer matrix. Although theoretical frameworks to prevent this process are available, we feel that the options explored in practice are often rather trial and error in nature. For that reason, we review the theories available, aiming to facilitate the design of the nanocomposites for a sustainable future. We first discuss thermodynamic frameworks which revolve around nanoparticle aggregation. To minimize nanoparticle aggregation, the nanoparticle and polymer can be selected in such a way that they have similar polar and dispersive surface energies. The second part is dedicated to nanocomposite processing, where kinetic effects act on the nanocomposite material therewith influencing its final morphology, although it is good to point out that other factors such as reaggregation also affect the final nanocomposite morphology. The third section is dedicated to how nanoparticles affect the polymer matrix to which they are added. We describe how interactions at an atomic scale, result in the formation of an interphasial region which ultimately leads to changed bulk material properties. From these three sections, we conclude that three parameters are often overlooked when designing nanocomposites, namely the surface energies of the nanoparticles and polymers, the aggregation bond energy or strength, and the interphase region. Therefore, in the fourth section we provide an overview of techniques to identify these three parameters. We finish with a summery and outlook for the design of bio nanocomposites, where we bring all insights from the previous four sections together.

Efficient Hybrid Particle-Field Coarse-Grained Model of Polymer Filler Interactions: Multiscale Hierarchical Structure of Carbon Black Particles in Contact with Polyethylene

In the present study, we propose, validate, and give first applications for large-scale systems of coarse-grained models suitable for filler/polymer interfaces based on carbon black (CB) and polyethylene (PE). The computational efficiency of the proposed approach, based on hybrid particle-field models (hPF), allows large-scale simulations of CB primary particles of realistic size (∼20 nm) embedded in PE melts. The molecular detailed models, here introduced, allow a microscopic description of the bound layer, through the analysis of the conformational behavior of PE chains adsorbed on different surface sites of CB primary particles, where the conformational behavior of adsorbed chains is different from models based on flat infinite surfaces. On the basis of the features of the systems, an optimized version of OCCAM code for large-scale (up to more than 8 million of beads) parallel runs is proposed and benchmarked. The computational efficiency of the proposed approach opens the possibility of a computational screening of the bound layer, involving the optimal combination of surface chemistry, size, and shape of CB aggregates and the molecular weight distribution of the polymers achieving an important tool to address the polymer/fillers interface and interphase engineering in the polymer industry.

Interface and Interphase in Polymer Nanocomposites with Bare and Core-Shell Gold Nanoparticles

Metal nanoparticles are used to modify/enhance the properties of a polymer matrix for a broad range of applications in bio-nanotechnology. Here, we study the properties of polymer/gold nanoparticle (NP) nanocomposites through atomistic molecular dynamics, MD, simulations. We probe the structural, conformational and dynamical properties of polymer chains at the vicinity of a gold (Au) NP and a functionalized (core/shell) Au NP, and compare them against the behavior of bulk polyethylene (PE). The bare Au NPs were constructed via a systematic methodology starting from ab-initio calculations and an atomistic Wulff construction algorithm resulting in the crystal shape with the minimum surface energy. For the functionalized NPs the interactions between gold atoms and chemically adsorbed functional groups change their shape. As a model polymer matrix we consider polyethylene of different molecular lengths, from the oligomer to unentangled Rouse like systems. The PE/Au interaction is parametrized via DFT calculations. By computing the different properties the concept of the interface, and the interphase as well, in polymer nanocomposites with metal NPs are critically examined. Results concerning polymer density profiles, bond order parameter, segmental and terminal dynamics show clearly that the size of the interface/interphase, depends on the actual property under study. In addition, the anchored polymeric chains change the behavior/properties, and especially the chain density profile and the dynamics, of the polymer chain at the vicinity of the Au NP.

The interfacial structure and dynamics in a polymer nanocomposite containing small attractive nanoparticles: a full atomistic molecular dynamics simulation study

We study the interfacial structure and dynamics of a polymer nanocomposite (PNC) composed of octaaminophenyl polyhedral oligomeric silsesquioxane (OAPS) and poly(2-vinylpyridine) (P2VP) by performing full atomistic molecular dynamics simulations. There are eight aminophenyl groups grafted on the surface of the OAPS particle and the particle has a size comparable to the Kuhn segment of P2VP. These aminophenyl groups can form hydrogen bonds (HBs) with pyridine rings from surrounding P2VP chains. We found that OAPS can form ∼2 HBs on average with surrounding polymer chains. The effect of the HBs is investigated in detail by either switching on or off these HBs in our simulation. By analyzing the interfacial static packing structure and dynamic properties, we demonstrate that the system has an ∼1 nm interface width, similar to the OAPS particle size. We also found that HBs can prevent the further penetration of polymers into the inner zone (grafting layer) of the OAPS, and therefore keep the P2VP chains in the outer layer (>1 nm), remaining bulk-like, which is well consistent with experimental results. In addition, we found that NP diffusion is coupled to the absorbed polymer chains, which also dramatically slows down the diffusion of polymer segments in return. The core-shell model in which the NP and absorbed polymers diffuse as a single object is validated here at the full atomistic level. These results provide atomistic insights into the unique structure and dynamics in the small attractive NP-polymer interfacial region. We hope these results will be helpful for the understanding of peculiar phenomena in attractive polymer nanocomposites containing small NPs.

Decoupling subsurface inhomogeneities: a 3D finite element approach for contact nanomechanical measurements

Novel material properties can be attained when embedding three-dimensional (3D) nanoparticles (NPs) in a variety of polymeric matrices. These inhomogeneities influence the bulk mechanical response due to the local high modulus mismatch between the particles and the matrix. The degree of the mechanical mismatch that is seen near a composite surface depends on the geometry/shape and spatial location and orientation of the particle with respect to the external contact loading. Isolating each particle's contribution to the surrounding elastic field can be numerically discerned but is experimentally complex, as there are limited direct characterization approaches available at the nanoscale. Atomic force microscopy (AFM) instrumentation is one such method that can quantify subsurface particle stiffness effects on nanocomposites with a resolution of a few nanometers. This work studies the spatial and geometrical effects of subsurface silver NPs on the local composite stiffness of a polystyrene matrix using 3D finite element (FE) models to interpret contact resonance (CR) AFM measurements. The present FE-AFM findings suggest both particle shape and particle orientation have a significant role in the degree of uniformity of the stiffness distribution in the embedding matrix. The applied CR-AFM technique shows that the NP geometry can be clearly distinguished when such inhomogeneities are relatively close, 17 nm, to a free surface whereas material-interface measurements at deeper subsurfaces are obscured by experimental noise. This work demonstrates that (i) numerical solutions can assist in qualitatively elucidating nanoinstrumentation stiffness profiles in terms of particle shape and orientation and (ii) CR-AFM measurements can quantify the influence of particle geometry and orientation on the surface nanomechanics of nanocomposite materials.

Solvation Effects on the Permeation and Aging Performance of PIM-1-Based MMMs for Gas Separation

Membranes are particularly attractive for lowering the energy intensity of separations as they eliminate phase changes. While many tantalizing polymers are known, limitations in selectivity and stability slightly preclude further development. Mixed-matrix membranes may address these shortcomings. Key to their realization is the intimate mixing between the polymer and the additive to eliminate nonselective transport, improve selectivity, and resist physical aging. Polymers of intrinsic microporosity (PIMs) have inherently promising gas transport properties. Here, we show that porous additives can improve transport and resist aging in PIM-1. We develop a simple, low-cost, and scalable hyper-cross-linked polymer (poly-dichloroxylene, pDCX), which was hydroxylated to form an intimate mixture with the polar PIM-1. Solvent variation allowed control of physical aging rates and improved selectivity for smaller gases. This detailed study has allowed many interactions within mixed matrix membranes to be directly elucidated and presents a practical means to stabilize porous polymers for separation applications.

Insights into Nano-Scale Physical and Mechanical Properties of Epoxy/Boehmite Nanocomposite Using Different AFM Modes

Understanding the interaction between nanoparticles and the matrix and the properties of interphase is crucial to predict the macroscopic properties of a nanocomposite system. Here, we investigate the interaction between boehmite nanoparticles (BNPs) and epoxy using different atomic force microscopy (AFM) approaches. We demonstrate benefits of using multifrequency intermodulation AFM (ImAFM) to obtain information about conservative, dissipative and van der Waals tip-surface forces and probing local properties of nanoparticles, matrix and the interphase. We utilize scanning kelvin probe microscopy (SKPM) to probe surface potential as a tool to visualize material contrast with a physical parameter, which is independent from the mechanics of the surface. Combining the information from ImAFM stiffness and SKPM surface potential results in a precise characterization of interfacial region, demonstrating that the interphase is softer than epoxy and boehmite nanoparticles. Further, we investigated the effect of boehmite nanoparticles on the bulk properties of epoxy matrix. ImAFM stiffness maps revealed the significant stiffening effect of boehmite nanoparticles on anhydride-cured epoxy matrix. The energy dissipation of epoxy matrix locally measured by ImAFM shows a considerable increase compared to that of neat epoxy. These measurements suggest a substantial alteration of epoxy structure induced by the presence of boehmite.

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.

Dramatic Increase in Polymer Glass Transition Temperature under Extreme Nanoconfinement in Weakly Interacting Nanoparticle Films

Properties of polymers in polymer nanocomposites and nanopores have been shown to deviate from their respective bulk properties due to physical confinement as well as polymer-particle interfacial interactions. However, separating the confinement effects from the interfacial effects under extreme nanoconfinement is experimentally challenging. Capillary rise infiltration enables polymer infiltration into nanoparticle (NP) packings, thereby confining polymers within extremely small pores and dramatically increasing the interfacial area, providing a good system to systematically distinguish the role of each effect on polymer properties. In this study, we investigate the effect of spatial confinement on the glass transition temperature ( T_g) of polystyrene (PS) infiltrated into SiO₂ NP films. The degree of confinement is tuned by varying the molecular weight of polymers, the size of NPs (diameters between 11 and 100 nm, producing 3-30 nm average pore sizes), and the fill-fraction of PS in the NP films. We show that in these dense NP packings the T_g of confined PS, which interacts weakly with SiO₂ NPs, significantly increases with decreasing pore size such that for the two molecular weights of PS studied the T_g increases by up to 50 K in 11 nm NP packings, while T_g is close to the bulk T_g in 100 nm NP packings. Interestingly, as the fill-fraction of PS is decreased, resulting in the accumulation of the polymer in the contacts between nanoparticles, hence an increased specific interfacial area, the T_g further increases relative to the fully filled films by another 5-8 K, indicating the strong role of geometrical confinement as opposed to the interfacial effects on the measured T_g values.

Segregation of liquid crystal mixtures in topological defects

The structure and physical properties of liquid crystal (LC) mixtures are a function of composition, and small changes can have pronounced effects on observables, such as phase-transition temperatures. Traditionally, LC mixtures have been assumed to be compositionally homogenous. The results of chemically detailed simulations presented here show that this is not the case; pronounced deviations of the local order from that observed in the bulk at defects and interfaces lead to significant compositional segregation effects. More specifically, two disclination lines are stabilized in this work by introducing into a nematic liquid crystal mixture a cylindrical body that exhibits perpendicular anchoring. It is found that the local composition deviates considerably from that of the bulk at the interface with the cylinder and in the defects, thereby suggesting new assembly and synthetic strategies that may capitalize on the unusual molecular environment provided by liquid crystal mixtures.

Liquid crystal mixtures are used in commercial applications and their composition affects their properties. Here Rahimi et al. use atomistic simulations to show that defects influence the molecular arrangement of the mixture components leading to a deviation of the local order from that of the bulk.

Investigation of interphase effects in silica-polystyrene nanocomposites based on a hybrid molecular-dynamics-finite-element simulation framework

A recently developed hybrid method is employed to study the mechanical behavior of silica-polystyrene nanocomposites (NCs) under uniaxial elongation. The hybrid method couples a particle domain to a continuum domain. The region of physical interest, i.e., the interphase around a nanoparticle (NP), is treated at molecular resolution, while the surrounding elastic continuum is handled with a finite-element approach. In the present paper we analyze the polymer behavior in the neighborhood of one or two nanoparticle(s) at molecular resolution. The coarse-grained hybrid method allows us to simulate a large polymer matrix region surrounding the nanoparticles. We consider NCs with dilute concentration of NPs embedded in an atactic polystyrene matrix formed by 300 chains with 200 monomer beads. The overall orientation of polymer segments relative to the deformation direction is determined in the neighborhood of the nanoparticle to investigate the polymer response to this perturbation. Calculations of strainlike quantities give insight into the deformation behavior of a system with two NPs and show that the applied strain and the nanoparticle distance have significant influence on the deformation behavior. Finally, we investigate to what extent a continuum-based description may account for the specific effects occurring in the interphase between the polymer matrix and the NPs.

Effect of Interfacial Bonding on Interphase Properties in SiO2/Epoxy Nanocomposite: A Molecular Dynamics Simulation Study

Atomistic molecular dynamics simulations have been performed to explore the effect of interfacial bonding on the interphase properties of a nanocomposite system that consists of a silica nanoparticle and the highly cross-linked epoxy matrix. For the structural properties, results show that interfacial covalent bonding can broaden the interphase region by increasing the radial effect range of fluctuated mass density and oriented chains, as well as strengthen the interphase region by improving the thermal stability of interfacial van der Waals excluded volume and reducing the proportion of cis conformers of epoxy segments. The improved thermal stability of the interphase region in the covalently bonded model results in an increase of ∼21 K in the glass transition temperature (Tg) compared to that of the pure epoxy. It is also found that interfacial covalent bonding mainly restricts the volume thermal expansion of the model at temperatures near or larger than Tg. Furthermore, investigations from mean-square displacement and fraction of immobile atoms point out that interfacial covalent and noncovalent bonding induces lower and higher mobility of interphase atoms than that of the pure epoxy, respectively. The obtained critical interfacial bonding ratio when the interphase and matrix atoms have the same mobility is 5.8%. These results demonstrate that the glass transitions of the interphase and matrix will be asynchronous when the interfacial bonding ratio is not 5.8%. Specifically, the interphase region will trigger the glass transition of the matrix when the ratio is larger than 5.8%, whereas it restrains the glass transition of the matrix when the ratio is smaller than 5.8%.

Understanding n-Octane Behavior near Graphene with Scaled Solvent-Solute Attractions

We employ molecular dynamics simulations of n-octane near a layered graphene surface to study the related phenomena of solvation, density fluctuations, wettability, and structure and dynamics of n-octane molecules in the inhomogeneous interfacial environment. That solvation in bulk n-octane displays a lengthscale-dependent crossover similar to that of hydrophobic solvation in water is known. Here we show that, near an extended graphene interface having attractive interactions with n-octane, lengthscale-dependent solvation is similar to that in the bulk and displays a small to large crossover. However, as the n-octane-graphene interactions are reduced to make the surface increasingly solvophobic, the crossover behavior is modulated and essentially absent near the most solvophobic surfaces, similar to that in water near hydrophobic interfaces. We show that the macroscopic measure of wettability, namely, the contact angle, characterizes n-octane-graphene coupling over a limited range of attractions. In contrast, molecular measures such as the free energy of cavity formation or the local compressibility in the interfacial region provide an effective measure of this coupling over a broader range of attractions. Finally, as n-octane-graphene attractions are increased, the n-octane liquid displays a wetting transition and corresponding change from sigmoidal to layered density profile. Analysis of the local structure shows that n-octane molecules prefer approximately linear conformations and surface-parallel orientations near the graphene surface, and their translational dynamics slow down with increasing n-octane-graphene attractions. Our study highlights molecular scale behavior of n-octane molecules that is relevant to understanding nanoparticle-solvent coupling in composite materials with enhanced mechanical or thermal properties.

Interfacial thermal transport and structural preferences in carbon nanotube–polyamide-6,6 nanocomposites: how important are chemical functionalization effects?

We investigate the influence of chemically functionalized CNTs on the structural properties of the surrounding polyamide-6,6 matrix as well as the interfacial thermal conductivity of polymer–CNT nanocomposites.