Interphase Structure in Silica–Polystyrene Nanocomposites: A Coarse-Grained Molecular Dynamics Study

Morphology-Driven Nanofiller Size Measurement Integrated with Micromechanical Finite Element Analysis for Quantifying Interphase in Polymer Nanocomposites

This study focused on an innovative practical method using computer vision for particle size measurement, which serves as a key precursor for predicting the elastic modulus of polymer nanocomposites. This approach involved the morphological segmentation of the nanodispersed phase. It aimed, for the first time, to address the impractical conditions resulting from the assumption of idealized single-particle sizes in a monodispersed system during modeling. Subsequently, a micromechanical finite element framework was employed to determine the interphase thickness and modulus in ultrahigh molecular weight polyethylene/nanozeolite composites, following the quantification of nanoparticle sizes. The size measurement approach relied on morphological images extracted from scanning electron microscopy micrographs of impact-fractured surfaces. To compute the interphase thickness, experimental data was fitted to an interphase-inclusive upper-bound Hashin-Shtrikman model, with the measured average particle size per composition serving as a crucial input. Subsequently, the interphase elastic modulus was computed based on its thickness, employing a hybrid modified-Hashin-Hansen and Maxwell model. These estimated interfacial variables were then utilized as inputs for the finite element model to determine the tensile modulus. A comparison between the model results and measured data revealed a maximum discrepancy of 3.29%, indicating the effectiveness of the methodology employed in quantifying interfacial properties.

Numerical investigations of translocation characteristics of nano-silica lunar dust across pulmonary surfactant monolayer

The interactions between nano-silica lunar dust (NSLD) on the moon surface and pulmonary surfactant (PS) monolayer will pose risks to astronaut health in future manned lunar exploration missions, but the specifics of these interactions are unknown. This study investigates them using the coarse-grained molecular dynamics method considering different sizes (5, 10, and 15 nm) and shapes (sphere, ellipsoid, and cube), with special focus on the unique morphology of NSLDs with bugles. The key findings are as follows: (1) The 10 nm and 15 nm NSLDs embed in the PS monolayer through the major sphere of spherical-type, major ellipsoid of ellipsoidal-type, or one edge of cubic-type NSLDs upon contact the PS monolayer. (2) Adsorbed NSLDs cause a higher Sz value (ASz > 0.84), while embedded NSLDs cause a lower Sz value (0.47 < ASz < 0.83) that decreases with an increase in the number of bulges. (3) The embedding process absorbs 50-342 dipalmitoylphosphatidylcholine (DPPC) molecules, reducing the PS monolayer area by 0.21%-6.05%. NSLDs with bulges absorb approximately 9-126 additional DPPC molecules and cause a 0.05%-3.22% reduction in the PS monolayer area compared to NSLDs without bulges. (4) NSLDs move obliquely or vertically within the PS monolayer, displaying two distinct stages with varying velocities. Their movement direction and speed are influenced by the increasing complexity of NSLD with more bulges on them. In general, larger NSLDs with sharper shapes and increasing complex morphology of more bulges cause more significant damages to the PS monolayer. These findings have implications for safeguarding astronaut health in future manned lunar exploration missions.

Molecular dynamics simulation of the impact of the surface topology of carbon black on the mechanical properties of elastomer nanocomposites

Carbon black has always played a pivotal role in reinforcing elastomers because it remarkably improves the mechanical properties. The reinforcing effect of carbon black is influenced by its grades, which mainly depend on the difference in the structure of the carbon black particles. Despite many traditional experiments on the performance of carbon black composites, there has been less emphasis on reinforcement mechanisms due to the challenges associated with unraveling the intermolecular interactions. In this paper, a coarse grained molecular dynamics simulation was employed to examine the relationship between the morphology of the carbon black particles and the mechanical properties of the elastomer nanocomposites. Specifically, three different morphological carbon black nanoparticle models, including the smooth particle model, rough particle model, and the rough ellipsoid model, were constructed first. We then focused on investigating the changes of the mechanical properties by systematically varying the filling fraction of the carbon black particles, and the strength of the interfacial interaction between the filler and the rubber. The results indicated that the surface roughness and the filler's shape had a significant impact on the mechanical properties of the filled rubber models. The mechanical enhancement effect of the rough ellipsoidal carbon black is around 50-400% higher than that of the smooth carbon black, and the stronger the interfacial interactions, the more pronounced the enhancement. In addition, the rough ellipsoid filled system has low hysteresis, low permanent deformation, and high fatigue resistance. In general, this work explores the strengthening mechanism of carbon black on the elastomer at the molecular level and generates new insight into the design and fabrication of novel reinforcing fillers.

Approaches and Perspective of Coarse-Grained Modeling and Simulation for Polymer-Nanoparticle Hybrid Systems

Molecular modeling and simulations have emerged as effective and indispensable tools to characterize polymeric systems. They provide fundamental and essential insights to design a product of the required properties and to improve the understanding of a phenomenon at the molecular level for a particular system. The polymer-nanoparticle hybrids are materials with outstanding properties and correspondingly large applications whose study has benefited from this new paradigm. However, despite the significant expansion of modern day computational powers, investigation of the long time and large length scale phenomenon in polymeric and polymer-nanoparticle systems is still a challenging task to complete through all-atom molecular dynamics (AA-MD) simulations. To circumvent this problem, a variety of coarse-grained (CG) models have been proposed, ranging from the generic CG models for qualitative properties predictions to more realistic chemically specific CG models for quantitative properties predictions. These CG models have already delivered some success stories in the study of several spatial and temporal evolutions of many processes. Some of these studies were beyond the feasibility of traditional atomistic resolution models due to either the size or the time constraints. This review captures the different types of popular CG approaches that are utilized in the investigation of the microscopic behavior of polymer-nanoparticle hybrid systems. The rationale of this article is to furnish an overview of the popular CG approaches and their applications, to review several important and most recent developments, and to delineate the perspectives on future directions in the field.

Tuning the polymer thermal conductivity through structural modification induced by MoS2 bilayers

The present work refers to a physical and structural study of nanoconfined polymers in polymer-MoS₂ nanocomposites as a function of MoS₂-MoS₂ interlayer distance. We applied reverse nonequilibrium molecular dynamics (RNEMD) simulations to investigate the thermal conductivity (λ) of polyamide oligomers confined by MoS₂ bilayers. The results of this study indicate that thermal conductivity of polymer can be considerably enhanced when polymer chains are confined by MoS₂ sheets, this behavior is more pronounced by charged surfaces. The presence of MoS₂ surfaces leads to elongation as well as preferential alignment of polymer chains parallel to the MoS₂ surfaces, which in turn results in higher order and denser packing of polymer content and hence larger thermal conductivity in comparison to the bulk polymer. Additionally, the analysis of the number of hydrogen bonds (HBs) in confined polymer chains suggests that a combined effect of the mentioned structural modification and enlarged values of HBs may cooperatively contribute to high polymer thermal conductivity, facilitating phonon transport. The results reported here suggest a significant way to design confined polymer-MoS₂ composites for significantly improving thermal conductivity for a wide variety of applications.

Versatilely Manipulating the Mechanical Properties of Polymer Nanocomposites by Incorporating Porous Fillers: A Molecular Dynamics Simulation Study

Polymer nanocomposites (PNCs) have been attracting myriad scientific and technological attention due to their promising mechanical and functional properties. However, there remains a need for an efficient method that can further strengthen the mechanical performance of PNCs. Here, we propose a strategy to design and fabricate novel PNCs by incorporating porous fillers (PFs) such as metal-organic frameworks with ultrahigh specific surface areas and tunable nanospaces to polymer matrices via coarse-grained molecular dynamics simulations. Three important parameters─the polymer chain stiffness (k), the interaction strength between the PF center and the end functional groups of polymer chains (ε_center end), and the PF weight fraction (w)─are systematically examined. First, attributed to the penetration of polymer chains into PFs at a strong ε_center end, the dimension of polymer chains such as the radius of gyration and the end-to-end distance increases greatly as a function of k compared to the case of the neat polymer system. The penetration of polymer chains is validated by characterizing the radial distribution function between end functional groups and filler centers, as well as the visualization of the snapshots. Also, the dispersion state of PFs tends to be good because of the chain penetration. Then, the glass transition temperature ratio of PNCs to that of the neat systems exhibits a maximum in the case of k = 5ε, indicating that the strongest interlocking between polymer chains and PFs occurs at intermediate chain stiffness. The polymer chain dynamics of PNCs decreases to a plateau at k = 5ε and then becomes stable, and the relative mobility to that of the neat system as well presents the same variation trend. Furthermore, the mechanical property under uniaxial deformation is thoroughly studied, and intermediates k, ε_center end, and w can bring about the best mechanical property. This is because of the robust penetration and interaction, which is confirmed by calculating the stress of every component of PNCs with and without end functional groups and PF centers as well as the nonbonded interaction energy change between different components. Finally, the optimal condition (k = 5.36ε, ε_center end = 5.29ε, and w = 6.54%) to design the PNC with superior mechanical behavior is predicted by Gaussian process regression, an active machine learning (ML) method. Overall, incorporating PFs greatly enhances the entanglements and interactions between polymer chains and nanofillers and brings effective mechanical reinforcements with lower filler weight fractions. We anticipate that this will provide new routes to the design of mechanically reinforced PNCs.

Prediction of the structure and mechanical properties of polycaprolactone-silica nanocomposites and the interphase region by molecular dynamics simulations: the effect of PEGylation

Polymer/silica (PS) nanocomposites are, among numerous combinations of inorganic/organic nanocomposites, one of the most important materials reported in the literature and have been employed in a wide variety of applications. Due to this great interest in the scientific and industry community, knowledge about their physiochemistry allows for a better understanding of their development and improvement. One area of interest found in biopolymers is silica, where silica nanoparticles can be used to increase their mechanical properties and give them higher opportunities to replace synthetic plastics. With this aim in mind, molecular dynamics (MD) simulations were used to predict the structure and mechanical properties of the interphase region and nanocomposite systems of polycaprolactone (PCL), a common poly(hydroxy acid) type biopolymer, reinforced with silica nanoparticles. Two types of nanoparticles were studied to assess the effect of PEGylation: hydroxyl (ungrafted) and polyethylene glycol (PEG) (grafted or PEGylated) functionalized silica. The interaction energy between the nanoparticle and the polymeric matrix was determined, showing an increase of the affinity between each component due to the PEGylation of the nanoparticle. Through the analysis of polymer density profiles, the structure and thickness of the interphase region were determined, and it was observed that PEGylation increased the interphase thickness from 10.80 Å to 13.04 Å while it decreased the peak and average polymer density of the interphase region. Using compressed and expanded molecular models of the neat PCL polymer, the mechanical properties of the interphase region were related to its density through an interpolation model, and mechanical property profiles were obtained, from which the average values of the Young's modulus, Poisson's ratio and shear modulus of the interphase region were calculated. Finally, the mechanical properties of the nanocomposites were determined by molecular mechanics simulations, showing that the silica nanoparticles increased the stiffness of the composite system to about 7-8% with respect to that of the neat polymer, having a 2.09% weight of bare silica or 2.82% weight of PEGylated silica. PEGylation did not show an additional effect on the overall mechanical properties. A mean field micromechanics model (Mori-Tanaka) corroborated the properties calculated for the interphase region using MD simulations. It was concluded that the PEGylation of the nanoparticle improved the affinity, and thus the dispersion, of the silica nanoparticles towards the PCL matrix, but with no further increase in the mechanical properties of the composite.

Addressing Surface Effects at the Particle-Continuum Interface in a Molecular Dynamics and Finite Elements Coupled Multiscale Simulation Technique

Atomistic-to-continuum coupling methods are used to unravel molecular mechanisms of polymers and polymer composites. These multiscale techniques advantageously combine the computational efficiency of continuum approaches while keeping the accuracy of particle-based methods. The Capriccio method [Pfaller et al. Comput. Methods Appl. Mech. Eng. 2013, 260, 109-129.] is a well-proven multiscale technique, which connects finite elements (FE) with molecular dynamics (MD) in a partitioned-domain approach. A vital aspect of these multiscale methods is to provide physically sound boundary conditions to the particle domain suppressing any interface effects at the domain boundary occurring due to the coupling. These interfacial coupling artifacts still pose a significant problem, especially for amorphous polymers due to their highly irregular microstructure. We solve this problem by extending the particle-continuum interface by a layer of passive atoms which move with the outer continuum, thereby providing the missing interactions with a surrounding polymer bulk to the inner particle region. This solution allows us to successfully reproduce structural and mechanical properties obtained under conventional periodic boundary conditions, like density, stress, Young's modulus, and Poisson's ratio. Furthermore, we demonstrate the application of a nonaffine deformation by means of a simple bending test. In general, our revised method provides a framework to apply complex deformations for molecular scientists, while it allows the engineering community to examine challenging phenomena such as fracture behavior at a molecular level.

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.

Heterogeneous Dynamics of Polymer Melts Exerted by Chain Loops Anchored on the Substrate: Insights from Molecular Dynamics Simulation

Understanding polymer-substrate interfacial dynamics at the molecular level is crucial for tailoring the properties of polymer ultrathin films (PUFs). Herein, through coarse-grained molecular dynamics simulation, the effect of length (N_loop) and rigidity (K_loop) of loop chains on the dynamics of linear chains is systematically explored, in which the loop chains are adsorbed on a solid substrate and the linear chains are covered on the loop chains. It is found that there is an optimal K_loop, which strongly confines the motion of the linear chains. Meanwhile, compared to increasing the rigidity of the loop chains, increasing the length of the loop chains can more effectively confine the motion of the linear chains. More interestingly, we observe that the mismatch of the length (ΔN) and rigidity (ΔK) between the loop and linear chains leads to dynamic asymmetry (ΔD_c). The relationship between the ΔN, ΔK, and ΔD_c are found to follow the mathematical expression of ΔD_c ∼ (ΔN)^α(ΔK)^β, in which the values of α and β are around 4.58 and 0.83, separately. Remarkably, using the Gaussian process regression model, we construct a master curve of diffusion coefficient on the segmental and chain length scales of the linear chains as a function of N_loop and K_loop, which is further validated by our simulated prediction. In general, this work provides a fundamental understanding of polymer interfacial dynamics at the molecular level, enlightening some rational principles for manipulating the physical properties of PUFs.

Uncertainty Quantification Guided Parameter Selection in a Fully Coupled Molecular Dynamics-Finite Element Model of the Mechanical Behavior of Polymers

The objective of investigating macroscopic polymer properties with a low computing cost and a high resolution has led to the development of efficient hybrid simulation tools. Systems generated from such simulation tools can fail in service if the effect of uncertainty of model inputs on its outputs is not accounted for. This work focuses on quantifying the effect of parametric uncertainty in our coarse-grained molecular dynamics-finite element coupling approach using uncertainty quantification. We consider uniaxial deformation simulations of a polystyrene sample at T = 100 K in our study. Parametric uncertainty is assumed to originate from parameters in the molecular dynamics model with a nonperiodic boundary (the force constant between polymer beads and anchor points, the number of anchor points, and the size of the surrounding dissipative particle dynamics domain) and a parameter to blend the energies of particles and continuum (weighting factor). Key issues that arise in uncertainty quantification are discussed on the basis of the quantities of interest including mass density, end-to-end distance, and radial distribution function. This work reveals the influence of key input parameters on the properties of polymer structure and facilitates the determination of those parameters in the application of this hybrid molecular dynamics-finite element approach.

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.

Optimizing the heterogeneous network structure to achieve polymer nanocomposites with excellent mechanical properties

Designing and optimizing the polymer network structure at the molecular level to manipulate its mechanical properties are of great scientific significance. Although heterogeneous multi-network structures have been extensively investigated, little effort has been devoted to investigating heterogeneous single-networks with a well-defined interface. Herein, through coarse-grained molecular dynamics simulation, we successfully fabricated a heterogeneous single-network, which was divided into several regions with different crosslink densities. Firstly, we found that there is an optimal crosslink density ratio between high and low crosslink density regions to obtain the best stress-strain behavior. Secondly, the effect of the regularity of the network topology (by changing the distribution of two-phase regions) on mechanical properties was also studied. It was clearly observed that the polymer network showed better elastic response and mechanical properties as the distribution of two-phase regions became uniform. Finally, we investigated the effect of the selective distribution of nanoparticles (NPs) on mechanical properties by introducing NPs into a pre-designed multiphase network. Results showed that the selective distribution of NPs in the high crosslink density region had a more significant effect on the mechanical reinforcement. Generally, our simulated results may provide some guidelines to design polymer network structures to achieve high-performance polymer nanocomposites with excellent mechanical properties.

Computer simulation of zwitterionic polymer brush grafted silica nanoparticles to modify polyvinylidene fluoride membrane

Dissipative particle dynamics (DPD) simulations was adopted to investigate the modification of polyvinylidene fluoride (PVDF) membrane by adding zwitterionic polymer brush poly(sulfobetaine methacrylate)- tetraethyl orthosilicate (PSBMA-TEOS) grafted silicon nanoparticles (SNPs) to the casting solution. The effects of polymer concentration and grafting architecture (PSBMA length and SNPs grafting ratio) on membrane morphology are discussed. When the polymer concentration reaches 40%, part of the SNPs is embedded in the membrane; the optimal polymer concentration is around 25-30%. In the SNPs system with the grafting ratio of 1, some SNPs are eluted into solution during phase separation. Compared with different grafting architectures, M_8-5, M_10-5 and M_12-5 system (M_x-y, where x represents the length of the zwitterionic polymer brush and y represents the grafting ratio of the silica nanoparticles) exhibited stable membrane morphologies. This work can provide guidance for the design and modification of organic-inorganic composite membrane and help understand the distribution of modified materials on the membrane surface.

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.

Role of Interface Structure and Chain Dynamics on the Diverging Glass Transition Behavior of SSBR-SiO2-PIL Elastomers

Intermolecular interactions between the constituents of a polymer nanocomposite at the polymer-particle interface strongly affect the segmental mobility of polymer chains, correlated with their glass transition behavior, and are responsible for the improved dynamical viscoelastic properties. In this work, we emphasized on the evolution of characteristic interfaces and their dynamics in silica (SiO2 NP)-reinforced, solution-polymerized, styrene butadiene rubber (SSBR) composites, whose relative prevalence varied with the phosphonium ionic liquid (PIL) volume fraction, used as an interfacial modifier. The molecular origins of such interfaces were examined through systematic dielectric spectroscopy, molecular dynamics (MD) simulations, and dynamic-mechanical analyses. The PIL facilitated H-bonding, cation-π, surface-phenyl, and van der Waals interfacial interactions between SSBR and SiO2 NP, thereby regulating the polymer chain dynamics, orientation, and mean-square displacement. Specifically, the mass density profiles from MD simulations revealed the dynamic gradient of polymer chains in the interfacial region as a function of radial distance from the center of mass of the SiO2 NP surface. The results showed a structuring effect to result in well-resolved density peaks at specific radial distances with the tangential orientation of styrene monomers in the vicinity of the SiO2 NP surface. These domino effects highlighted strong interfacial interactions to have an indispensable effect on the viscoelastic performance and thermal motion of SSBR molecular chains, leading to a higher glass transition temperature (T g) by ∼15 K, validating the experimental data. More importantly, our results gave new insights into the fundamental understanding of the fact that the strength of intermolecular interactions induced by PIL at the polymer-particle interface is the key to control the α-relaxation dynamics and T g optimization, desired for specific applications.

Material identification for improving the strength of silica/SBR interface using MD simulations

Comprehensive molecular dynamics simulations are conducted to identify material modifications which can improve strength and reduce hysteresis losses at the nanointerfaces formed between silica, silane coupling agent (SCA) and styrene-butadiene rubber (SBR), all of which are important ingredients of green tyres. Improving strength and reducing hysteresis losses at such interfaces are expected to reduce rolling resistance (RR), consequently lowering greenhouse emissions. Various modifications considered in this work include a variety of SBR blends, several SCA and surface occupancies of SCA on the silica surface. To tackle a large number of combinations possible and identify modifications which may improve the nature of the interfaces, a hierarchical computational framework is developed. The reduced sample space of such material modifications may be more amenable to comprehensive and computationally or experimentally expensive studies. It was found that some amino-based SCA in combination with certain blends of SBR can improve the interfaces strength and lower hysteresis losses, when compared to the commonly used bis[3-(triethoxysilyl)propyl]tetrasulfide (TESPT), which is a sulphur-based SCA.

Mechanical and Self-Healing Behavior of Matrix-Free Polymer Nanocomposites Constructed via Grafted Graphene Nanosheets

Through molecular dynamics (MD) simulation, the structure and mechanical properties of matrix-free polymer nanocomposites (PNCs) constructed via polymer-grafted graphene nanosheets are studied. The dispersion of graphene sheets is characterized by the radial distribution function (RDF) between graphene sheets. We observe that a longer polymer chain length L_g leads to a relatively better dispersion state attributed to the formation of a better brick-mud structure, effectively screening the van der Waals interactions between sheets. By tuning the interaction strength ε_end-end between end functional groups of grafted chains, we construct physical networks with various robustness characterized by the formation of the fractal clusters at high ε_end-end values. The effects of ε_end-end and L_g on the mechanical properties are examined, and the enhancement of the stress-strain behavior is observed with the increase of ε_end-end and L_g. Structural evolution during deformation is quantified by calculating the orientation of the graphene sheets and their distribution, the stress decomposition, and the size of the clusters formed between end groups and their distribution. Then, we briefly study the effects of time and temperature on the self-healing behavior of these unique PNCs in the rubbery state. Lastly, the self-healing kinetics is quantitatively analyzed. In general, this work can provide some rational guidelines to design and fabricate matrix-free PNCs with both excellent mechanical and self-healing properties.

Multi-scale modeling of the polymer-filler interaction

We report mesoscopic simulations of the interaction between a silica nanoparticle and cis-1,4-polybutadiene chains with realistic coarse-(CG) grained models. The CG models are obtained with a bottom-up Bayesian method based on trajectory matching of atomistic configurations of the system. We then investigate the structural properties of the interfacial region as a function of the grafting density and polymer chain length. We take advantage of the realistic CG models to explore the dynamics of the nanoparticle over a period of 10 microseconds. We show that the dynamics of the nanoparticle is affected by the grafting density and the polymer chain length of the grafted chains.

Extensive CGMD Simulations of Atactic PS Providing Pseudo Experimental Data to Calibrate Nonlinear Inelastic Continuum Mechanical Constitutive Laws

In this contribution, we present a characterization methodology to obtain pseudo experimental deformation data from CG MD simulations of polymers as an inevitable prerequisite to choose and calibrate continuum mechanical constitutive laws. Without restriction of generality, we employ a well established CG model of atactic polystyrene as exemplary model system and simulate its mechanical behavior under various uniaxial tension and compression load cases. To demonstrate the applicability of the obtained data, we exemplarily calibrate a viscoelastic continuum mechanical constitutive law. We conclude our contribution by a thorough discussion of the findings obtained in the numerical pseudo experiments and give an outline of subsequent research activities. Thus, this work contributes to the field of multiscale simulation methods and adds a specific application to the body of knowledge of CG MD simulations.

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.

EFFECT OF THE NANOFILLER SHAPE ON THE CONDUCTIVE NETWORK FORMATION OF POLYMER NANOCOMPOSITES VIA A COARSE-GRAINED SIMULATION

It is very important to improve the electrical conductivity of polymer nanocomposites, which can widen their application. The effect of the nanofiller shape on the relationship between the nanofiller microstructure and the conductive probability of the nanofiller filled polymer nanocomposites (PNCs) has been investigated in detail by employing a coarse-grained molecular dynamics simulation. Four kinds of nanofiller shapes are considered: rod filler, Y filler, X filler, and sphere filler. First, the mean square radius of gyration gradually decreases from rod filler, Y filler, X filler, to sphere filler, which reflects the highest aspect ratio for rod filler. Meanwhile, the dispersion state of the nanofiller is relatively uniform in the matrix. The conductive probability (denoted by the formation probability of the conductive network) is adopted to stand for the conductive property. The results show that the conductive probability gradually decreases from rod filler, Y filler, X filler, to sphere filler, which is attributed to their gradually decreased size. In summary, the nanofiller shape affects the electric conductive property of PNCs.

Miscibility and Nanoparticle Diffusion in Ionic Nanocomposites

We investigate the effect of various spherical nanoparticles in a polymer matrix on dispersion, chain dimensions and entanglements for ionic nanocomposites at dilute and high nanoparticle loading by means of molecular dynamics simulations. The nanoparticle dispersion can be achieved in oligomer matrices due to the presence of electrostatic interactions. We show that the overall configuration of ionic oligomer chains, as characterized by their radii of gyration, can be perturbed at dilute nanoparticle loading by the presence of charged nanoparticles. In addition, the nanoparticle's diffusivity is reduced due to the electrostatic interactions, in comparison to conventional nanocomposites where the electrostatic interaction is absent. The charged nanoparticles are found to move by a hopping mechanism.

Interface Characterization of Epoxy Resin Nanocomposites: A Molecular Dynamics Approach

In polymer nanocomposites, the interface region between the matrix and the fillers has been identified as a key interaction region that strongly determines the properties of the final material. Determining its structure is crucial from several points of view, from modeling (i.e., properties prediction) to materials science (i.e., understanding properties/structure relationships). In the presented paper, a method for characterizing the interface region of polymer nanocomposites is described using molecular dynamics (MD) simulations. In particular, the structure of the polymer within the interface region together with its dimension in terms of thickness were analyzed through density profiles. Epoxy resin nanocomposites based on diglycidyl ether of bisphenol A (DGEBA) were studied using this approach, and the interface region with triple walled carbon nanotubes (TWCNT) and carbon fibers (CF) was characterized. The effect of carbon nanotube diameter, type of hardener, and effect of epoxy resin cross-linking degree on interface thickness were analyzed using MD models. From this analysis no general rule on the effect of these parameters on the interface thickness could be established, since in some cases overlapping effects between the analyzed parameters were observed, and each specific case needs to be analyzed independently in detail. Results show that the diameter has an impact on interface thickness, but this effect is affected by the cross-linking degree of the epoxy resin. The type of hardener also has a certain influence on the interface thickness.

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.

Perspective: Outstanding theoretical questions in polymer-nanoparticle hybrids

This topical review discusses the theoretical progress made in the field of polymer nanocomposites, i.e., hybrid materials created by mixing (typically inorganic) nanoparticles (NPs) with organic polymers. It primarily focuses on the outstanding issues in this field and is structured around five separate topics: (i) the synthesis of functionalized nanoparticles; (ii) their phase behavior when mixed with a homopolymer matrix and their assembly into well-defined superstructures; (iii) the role of processing on the structures realized by these hybrid materials and the role of the mobilities of the different constituents; (iv) the role of external fields (electric, magnetic) in the active assembly of the NPs; and (v) the engineering properties that result and the factors that control them. While the most is known about topic (ii), we believe that significant progress needs to be made in the other four topics before the practical promise offered by these materials can be realized. This review delineates the most pressing issues on these topics and poses specific questions that we believe need to be addressed in the immediate future.

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.

Tailoring the mechanical properties by molecular integration of flexible and stiff polymer networks

Designing a multiple-network structure at the molecular level to tailor the mechanical properties of polymeric materials is of great scientific and technological importance. Through the coarse-grained molecular dynamics simulation, we successfully construct an interpenetrating polymer network (IPN) composed of a flexible polymer network and a stiff polymer network. First, we find that there is an optimal chain stiffness for a single network (SN) to achieve the best stress-strain behavior. Then we turn to study the mechanical behaviors of IPNs. The result shows that the stress-strain behaviors of the IPNs appreciably exceed the sum of that of the corresponding single flexible and stiff network, which highlights the advantage of the IPN structure. By systematically varying the stiffness of the stiff polymer network of the IPNs, optimal stiffness also exists to achieve the best performance. We attribute this to a much larger contribution of the non-bonded interaction energy. Last, the effect of the component concentration ratio is probed. With the increase of the concentration of the flexible network, the stress-strain behavior of the IPNs is gradually enhanced, while an optimized concentration (around 60% molar ration) of the stiff network occurs, which could result from the dominant role of the enthalpy rather than the entropy. In general, our work is expected to provide some guidelines to better tailor the mechanical properties of the IPNs made of a flexible network and a stiff network, by manipulating the stiffness of the stiff polymer network and the component concentration ratio.

Molecular insight into the Mullins effect: irreversible disentanglement of polymer chains revealed by molecular dynamics simulations

The debate regarding the possible molecular origins of the Mullins effect has been ongoing since its discovery. Molecular dynamics (MD) simulations were carried out to elucidate the underlying mechanism of the Mullins effect. For the first time, the key characteristics associated with the Mullins effect, including (a) the majority of stress softening occurring in the first stretch, (b) continuous softening with stress increase, (c) a permanent set, and (d) recovery with heat treatment, are captured by molecular modeling. It is discovered that the irreversible disentanglement of polymer chains is physically sufficient to interpret these key characteristics, providing molecular evidence for this long-controversial issue. Our results also reveal that filled polymers exhibit three distinct regimes, i.e., the polymer matrix, the interface, and the filler. When subjected to external strain, the polymer matrix suffers from excess deformation, indicating strong heterogeneity within the filled polymer, which offers molecular insight for the formulation of physics-based constitutive relations for filled polymers.

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

The selectivity of nanoparticles for polydispersed ligand chains during the grafting-to process: a computer simulation study

By constructing a grafting-to reaction model of polydispersed polymer chains to bind onto nanoparticles (NPs), we elucidate the changes of grafting density, polydispersity index and chain length distribution of grafted ligand chains as a dependence of the feeding polymer chains.

Simulational insights into the mechanical response of prestretched double network filled elastomers

This paper deals with molecular-dynamics simulations of the mechanical properties of prestretched double network filled elastomers. To this end, we firstly validated the accuracy of this method, and affirmed that the produced stress-strain characteristics were qualitatively consistent with Lesser's experimental results on the prestretched tri-block copolymers with a competitive double network. Secondly, we investigated the effect of the crosslinking network ratio on the mechanical properties of the prestretched double network homopolymers under uniaxial tension. We found that the prestretched double network contributes greatly to the enhanced tensile stress and ultimate strength at fracture, as well as to the lower permanent set (the residual strain) and dynamic hysteresis loss, both parallel and perpendicular to the prestretching direction. Notably, though, an anisotropic behavior was observed: in the parallel direction, both the first and the second crosslinked networks bore the external force; whereas in the perpendicular direction, only the second crosslinked network was relevantly effective. Finally, the polymer nanocomposites with a prestretched double network exhibited tensile mechanical properties similar to those of the studied homopolymers with prestretched double networks. Summing up the results, it can be concluded that the incorporation of prestretched double networks with a specified crosslinking network ratio seems to be a promising method for manipulating the mechanical properties of elastomers and their nanocomposites, as well as for introducing anisotropy in their mechanical response.