Effect of Filler Particle Size on the Properties of Model Nanocomposites

Enhancing Multifunctionality: Optimal Properties of Iron-Oxide-Reinforced Polyvinylidene Difluoride Unveiled Through Full Atom Molecular Dynamics Simulations

Nanocomposites made of magnetite (Fe3O4) nanoparticles (NP)s with different surface chemistry and polyvinyl difluoride (PVDF) polymer were investigated using full atom molecular dynamics (MD) simulation. NPs with hydroxyl (OH), hexanoic, and oleic acid terminations were considered in this study. The effect of each surface chemistry was investigated in terms of the mechanical properties, the distribution of the internal energy around the NP, and the chain polarization gradient from the interface to the bulk. From this investigation, we find that oleic acid termination, although the most popular, is less favorable for interfacial interaction and local polarization. The OH-terminated NP results in the best configuration for the properties investigated. The hexanoic acid-grafted NP presents a good compromise. Hydrogen bonding governs the induced response of the nanocomposites. Although the hexanoic acid grafted NP presents less hydrogen bonding than the OH-terminated case, the conformation of the hexanoic acid acts as a mobility flow inhibitor, leading to a performance comparable to that of the OH-terminated NP composite. This work led to investigating routes to make nanocomposite materials with optimized properties. These results shed light on the multiple combinations offered by nanocomposites that go beyond the conventional effects of size.

Revealing the Role of Chain Conformations on the Origin of the Mechanical Reinforcement in Glassy Polymer Nanocomposites

Understanding the mechanism of mechanical reinforcement in glassy polymer nanocomposites is of paramount importance for their tailored design. Here, we present a detailed investigation, via atomistic simulation, of the coupling between density, structure, and conformations of polymer chains with respect to their role in mechanical reinforcement. Probing the properties at the molecular level reveals that the effective mass density as well as the rigidity of the matrix region changes with filler volume fraction, while that of the interphase remains constant. The origin of the mechanical reinforcement is attributed to the heterogeneous chain conformations in the vicinity of the nanoparticles, involving a 2-fold mechanism. In the low-loading regime, the reinforcement comes mainly from a thin, single-molecule, 2D-like layer of adsorbed polymer segments on the nanoparticle, whereas in the high-loading regime, the reinforcement is dominated by the coupling between train and bridge conformations; the latter involves segments connecting neighboring nanoparticles.

Dynamics of Polymer Chains in Poly(ethylene oxide)/Silica Nanocomposites via a Combined Computational and Experimental Approach

The dynamics of polymer chains in poly(ethylene oxide)/silica (PEO/SiO₂) nanoparticle nanohybrids have been investigated via a combined computational and experimental approach involving atomistic molecular dynamics simulations and dielectric relaxation spectroscopy (DRS) measurements. The complementarity of the approaches allows us to study systems with different polymer molecular weights, nanoparticle radii, and compositions across a broad range of temperatures. We study the effects of spatial confinement, which is induced by the nanoparticles, and chain adsorption on the polymer's structure and dynamics. The investigation of the static properties of the nanocomposites via detailed atomistic simulations revealed a heterogeneous polymer density layer at the vicinity of the PEO/SiO₂ interface that exhibited an intense maximum close to the inorganic surface, whereas the bulk density was reached for distances ∼1-1.2 nm away from the nanoparticle. For small volume fractions of nanoparticles, the polymer dynamics, probed by the atomistic simulations of low-molecular-weight chains at high temperatures, are consistent with the presence of a thin adsorbed layer that exhibits slow dynamics, with the dynamics far away from the nanoparticle being similar to those in the bulk. However, for high volume fractions of nanoparticles (strong confinement), the dynamics of all polymer chains were predicted slower than that in the bulk. On the other hand, similar dynamics were found experimentally for both the local β-process and the segmental dynamics for high-molecular-weight systems measured at temperatures below the melting temperature of the polymer, which were probed by DRS. These differences can be attributed to various parameters, including systems of different molecular weights and nanoparticle states of dispersion, the different temperature range studied by the different methods, the potential presence of a reduced-mobility PEO/SiO₂ interfacial layer that does not contribute to the dielectric spectrum, and the presence of amorphous-crystalline interfaces in the experimental samples that may lead to a different dynamical behaviors of the PEO chains.

Nanoparticle Shape Influence over Poly(lactic acid) Barrier Properties by Molecular Dynamics Simulations

Climate change is leading us to search for new materials that allow a more sustainable environmental situation in the long term. Poly(lactic acid) (PLA) has been proposed as a substitute for traditional plastics due to its high biodegradability. Various components have been added to improve their mechanical, thermal, and barrier properties. The modification of the PLA barrier properties by introducing nanoparticles with different shapes is an important aspect to control the molecular diffusion of oxygen and other gas compounds. In this work, we have described changes in oxygen diffusion by introducing nanoparticles of different shapes through molecular dynamics simulations. Our model illustrates that the existence of curved surfaces and the deposition of PLA around them by short chains generate small holes where oxygen accumulates, forming clusters and reducing their mobility. From the several considered shapes, the sphere is the most suitable structure to improve the barrier properties of the PLA.

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.

Coupling between Polymer Conformations and Dynamics Near Amorphous Silica Surfaces: A Direct Insight from Atomistic Simulations

The dynamics of polymer chains in the polymer/solid interphase region have been a point of debate in recent years. Its understanding is the first step towards the description and the prediction of the properties of a wide family of commercially used polymeric-based nanostructured materials. Here, we present a detailed investigation of the conformational and dynamical features of unentangled and mildly entangled cis-1,4-polybutadiene melts in the vicinity of amorphous silica surface via atomistic simulations. Accounting for the roughness of the surface, we analyze the properties of the polymer chains as a function of their distance from the silica slab, their conformations and the chain molecular weight. Unlike the case of perfectly flat and homogeneous surfaces, the monomeric translational motion parallel to the surface was affected by the presence of the silica slab up to distances comparable with the extension of the density fluctuations. In addition, the intramolecular dynamical heterogeneities in adsorbed chains were revealed by linking the conformations and the structure of the adsorbed chains with their dynamical properties. Strong dynamical heterogeneities within the adsorbed layer are found, with the chains possessing longer sequences of adsorbed segments ("trains") exhibiting slower dynamics than the adsorbed chains with short ones. Our results suggest that, apart from the density-dynamics correlation, the configurational entropy plays an important role in the dynamical response of the polymers confined between the silica slabs.

Structure and Dynamics of Highly Attractive Polymer Nanocomposites in the Semi-Dilute Regime: The Role of Interfacial Domains and Bridging Chains

Detailed molecular dynamics (MD) simulations are employed to study how the presence of adsorbed domains and nanoparticle bridging chains affect the structural, conformational, thermodynamic, and dynamic properties of attractive polymer nanocomposite melts in the semi-dilute regime. As a model system we have chosen an unentangled poly(ethylene glycol) (PEG) matrix containing amorphous spherical silica nanoparticles with different diameters and at different concentrations. Emphasis is placed on properties such as the polymer mass density profile around nanoparticles, the compressibility of the system, the mean squared end-to-end distance of PEG chains, their orientational and diffusive dynamics, the single chain form factor, and the scattering functions. Our analysis reveals a significant impact of the adsorbed, interfacial polymer on the microscopic dynamic and conformational properties of the nanocomposite, especially under conditions favoring higher surface-to-volume ratios (e.g., for small nanoparticle sizes at fixed nanoparticle loading, or for higher silica concentrations). Simultaneously, adsorbed polymer chains adopt graft-like conformations, a feature that allows them to considerably extend away from the nanoparticle surface to form bridges with other nanoparticles. These bridges drive the formation of a nanoparticle network whose strength (number of tie chains per nanoparticle) increases substantially with increasing concentration of the polymer matrix in nanoparticles, or with decreasing nanoparticle size at fixed nanoparticle concentration. The presence of hydroxyl groups at the ends of PEG chains plays a key role in the formation of the network. If hydroxyl groups are substituted by methoxy ones, the simulations reveal that the number of bridging chains per nanoparticle decreases dramatically, thus the network formed is less dense and less strong mechanically, and has a smaller impact on the properties of the nanocomposite. Our simulations predict further that the isothermal compressibility and thermal expansion coefficient of PEG-silica nanocomposites are significantly lower than those of pure PEG, with their values decreasing practically linear with increasing concentration of the nanocomposite in nanoparticles.

Anomalous Terminal Shear Viscosity Behavior of Polycarbonate Nanocomposites Containing Grafted Nanosilica Particles

Viscosity controls an important issue in polymer processing. This paper reports on the terminal viscosity behavior of a polymer melt containing grafted nanosilica particles. The melt viscosity behavior of the nanocomposites was found to depend on the interaction between the polymer matrix and the nanoparticle surface. In the case of polycarbonate (PC) nanocomposites, the viscosity decreases by approximately 25% at concentrations below 0.7 vol% of nanosilica, followed by an increase at higher concentrations. Chemical analysis shows that the decrease in viscosity can be attributed to in situ grafting of PC on the nanosilica surface, leading to a lower entanglement density around the nanoparticle. The thickness of the graft layer was found to be of the order of the tube diameter, with the disentangled zone being approximately equal to the radius of gyration (Rg) polymer chain. Furthermore, it is shown that the grafting has an effect on the motion of the PC chains at all timescales. Finally, the viscosity behavior in the PC nanocomposites was found to be independent of the molar mass of PC. The PC data are compared with polystyrene nanocomposites, for which the interaction between the polymer and nanoparticles is absent. The results outlined in this paper can be utilized for applications with low shear processing conditions, e.g., rotomolding, 3D printing, and multilayer co-extrusion.

Multiscale Analysis of Elastic Properties of Nano-Reinforced Materials Exhibiting Surface Effects. Application for Determination of Effective Shear Modulus

This work concerns a multiscale analysis of nano-reinforced heterogeneous materials. Such materials exhibit surface effects that must be taken into account in the homogenization procedure. In this study, a coherent imperfect interface model was employed to characterize the jumps of mechanical properties through the interface region between the matrix and the nanofillers. As the hypothesis of scale separation was adopted, a generalized self-consistent micromechanical scheme was employed for the determination of the homogenized elastic moduli. An explicit calculation for the determination of effective shear modulus is presented, together with a numerical application illustrating the surface effect. It is shown that the coherent imperfect interface model is capable of exploring the surface effect in nano-reinforced materials, as demonstrated experimentally in the literature.

Molecular structure and multi-body potential of mean force in silica-polystyrene nanocomposites

We perform a systematic application of the hybrid particle-field molecular dynamics technique [Milano, et al., J. Chem. Phys., 2009, 130, 214106] to study interfacial properties and potential of mean force (PMF) for separating nanoparticles (NPs) in a melt. Specifically, we consider Silica NPs bare or grafted with Polystyrene chains, aiming to shed light on the interactions among free and grafted chains affecting the dispersion of NPs in the nanocomposite. The proposed hybrid models show good performances in catching the local structure of the chains, and in particular their density profiles, documenting the existence of the "wet-brush-to-dry-brush" transition. By using these models, the PMF between pairs of ungrafted and grafted NPs in Polystyrene matrix are calculated. Moreover, we estimate the three-particle contribution to the total PMF and its role in regulating the phase separation on the nanometer scale. In particular, the multi-particle contribution to the PMF is able to give an explanation of the complex experimental morphologies observed at low grafting densities. More in general, we propose this approach and the models utilized here for a molecular understanding of specific systems and the impact of the chemical nature of the systems on the composite final properties.

Novel three-dimensional interphase characterisation of polymer nanocomposites using nanoscaled topography

Mechanical properties of polymer nanocomposites depend primarily on nanointerphases as transitional zones between nanoparticles and surrounding matrices. Due to the difficulty in the quantitative characterisation of nanointerphases, previous literature generally deemed such interphases as one-dimensional uniform zones around nanoparticles by assumption for analytical or theoretical modelling. We hereby have demonstrated for the first time direct three-dimensional topography and physical measurement of nanophase mechanical properties between nanodiameter bamboo charcoals (NBCs) and poly (vinyl alcohol) (PVA) in polymer nanocomposites. Topographical features, nanomechanical properties and dimensions of nanointerphases were systematically determined via peak force quantitative nanomechanical tapping mode. Significantly different mechanical properties of nanointerphases were revealed as opposed to those of individual NBCs and PVA matrices. Non-uniform irregular three-dimensional structures and shapes of nanointerphases are manifested around individual NBCs, which can be greatly influenced by nanoparticle size and roughness, and nanoparticle dispersion and distribution. Elastic moduli of nanointerphases were experimentally determined in range from 25.32 ± 3.4 to 66.3 ± 3.2 GPa. Additionally, it is clearly shown that the interphase modulus strongly depends on interphase surface area and interphase volume. Different NBC distribution patterns from fully to partially embedded nanoparticles are proven to yield a remarkable reduction in elastic moduli of nanointerphases.

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.

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.

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%.

New insights into reinforcement mechanism of nanoclay-filled isoprene rubber during uniaxial deformation by in situ synchrotron X-ray diffraction

The existence of a denser network domain formed by incorporation of filler and its vital role in determining the strain-induced crystallization behavior of nanocomposites is proved by in situ synchrotron X-ray diffraction characterization.

Bound Polymer Layer in Nanocomposites

There has been considerable interest in characterizing the polymer layer that is effectively irreversibly bound to nanoparticles (NPs) because it is thought to underpin the unusual thermomechanical properties of polymer nanocomposites (PNC). We study PNCs formed by mixing silica nanoparticles (NPs) with poly-2-vinylpyridine (P2VP) and compare the bound layer thickness δ determined by three different methods. We show that the thickness obtained by thermogravimetric analysis (TGA) and assuming that the bound layer has a density corresponding to a dense melt clearly underestimates the real bound layer thickness. A more realistic extent of the bound layer is obtained by in situ measurements of the interaction pair potential between NPs in PNCs via analysis of TEM micrographs; we verify these estimates using Dynamic Light Scattering (DLS) in θ solvent. Our results confirm the existence of long-ranged interactions between NPs corresponding roughly in size to the radius of gyration of the bound chains.