Molecular Dynamics Simulation of a Silica Nanoparticle in Oligomeric Poly(methyl methacrylate): A Model System for Studying the Interphase Thickness in a Polymer–Nanocomposite via Different Properties

Nanoscopically Confined 1,4-Polybutadiene Melts: Exploring Confinement by Alumina Nanorod and Nanopore Systems

We present molecular dynamics simulations of a chemically realistic model of 1,4-polybutadiene (PBD) in contact with curved alumina surfaces. We contrast the behavior of PBD infiltrated into alumina pores with a curvature radius of about three times the radius of gyration of the chains to its behavior next to a melt dispersed alumina rod of equal absolute curvature. These confinement types represent situations occurring in polymer melts loaded with nanoparticles due to nanoparticle aggregation. While there are observable differences in structure and dynamics due to the different types of geometric confinement, the main effects stem from the strong attraction of PBD to the alumina surfaces. This strong attraction leads to a deformation of the chains in contact to the surfaces. We focus on temperatures well above the bulk glass transition temperature, but even at these high temperatures, the layers next to the alumina surfaces show glass-like relaxation behavior. We analyze the signature of this glassy behavior for neutron scattering or nuclear magnetic resonances experiments.

Computational investigation of the effects of polymer grafting on the effective interaction between silica nanoparticles in water

Understanding and control of the effective interaction between nanoscale building blocks (colloids or nanoparticles) dispersed in a solvent is an important prerequisite for the development of bottom-up design strategies for soft functional materials. Here, we have employed all-atom molecular dynamics simulations to investigate the impact of polymer grafting on the solvent-mediated effective interaction between the silica nanoparticles (Si-NPs) in water, and in turn, on its bulk structural and thermodynamic properties. We found that the nature of the short grafting polymers [characterized by their interaction with water (hydrophobicity or hydrophilicity) and molecular weight] has a profound effect on the range and strength of the effective interaction between the Si-NPs. The hydrophobic polymer [such as polyethylene (PE)]-grafting of Si-NP gives rise to a more attractive interaction between the Si-NPs compared to the hydrophilic polymer [such as polyethylene glycol (PEG)] and non-grafted cases. This study further provides fundamental insights into the molecular origin of the observed behavior of the effective pair interactions between the grafted Si-NPs. For PE-grafted Si-NPs, the confined water (water inside the cavity formed by a pair of Si-NPs) undergoes a partial dewetting transition on approaching below a critical inter-particle separation leading to a stronger attractive interaction. Furthermore, we report that the effective attraction between the PE-grafted Si-NPs can be reliably controlled by changing the grafting PE density. We have also investigated the bulk structural and thermodynamic behavior of the coarse-grained Si-NP system where the particles interact via effective interaction in the absence of water. We believe that the insights gained from this work are important prerequisites for formulating rational bottom-up design strategies for functional materials where nano- (or, colloidal) particles are the building blocks.

Model Approach to Thermal Conductivity in Hybrid Graphene-Polymer Nanocomposites

The thermal conductivity of epoxy nanocomposites filled with self-assembled hybrid nanoparticles composed of multilayered graphene nanoplatelets and anatase nanoparticles was described using an analytical model based on the effective medium approximation with a reasonable amount of input data. The proposed effective thickness approach allowed for the simplification of the thermal conductivity simulations in hybrid graphene@anatase TiO2 nanosheets by including the phenomenological thermal boundary resistance. The sensitivity of the modeled thermal conductivity to the geometrical and material parameters of filling particles and the host polymer matrix, filler's mass concentration, self-assembling degree, and Kapitza thermal boundary resistances at emerging interfaces was numerically evaluated. A fair agreement of the calculated and measured room-temperature thermal conductivity was obtained.

Molecular Origin of the Reinforcement Effect and Its Strain-Rate Dependence in Polymer Nanocomposite Glass

We investigate the molecular origin of mechanical reinforcement in a polymer nanocomposite (PNC) under a glass state via molecular dynamics simulations. The strength of the PNC system is found to be reinforced mainly via reduced plastic deformations of the nanoparticle neighborhood (NN). Such a reinforcement effect is found to decay with an increase in the strain rate. The Arrhenius-Eyring relation is used to analyze its origin. The amplitude of the reinforcement is found to be determined by the difference between the energy barrier (ΔE) for the activation of NN and the work (W) done by the applied stress to conquer that barrier. A larger strain rate is found to result in a larger W and, hence, a weaker reinforcement effect. Such a strain-rate dependence is verified in the experimental tensile tests of a poly(vinyl alcohol)/SiO2 composite system. These results not only provide a new understanding of the molecular origin of the reinforcement effect in the PNC system, but also pave the way for a better design of the PNC material properties.

The Role of Polymer Chain Stiffness and Guest Nanoparticle Loading in Improving the Glass Transition Temperature of Polymer Nanocomposites

The impact of polymer chain stiffness characterized by the bending modulus (k_θ) on the glass transition temperature (T_g) of pure polymer systems, as well as polymer nanocomposites (PNCs), is investigated using molecular dynamics simulations. At small k_θ values, the pure polymer system and respective PNCs are in an amorphous state, whereas at large k_θ values, both systems are in a semicrystalline state with a glass transition at low temperature. For the pure polymer system, T_g initially increases with k_θ and does not change obviously at large k_θ. However, the T_g of PNCs shows interesting behaviors with the increasing volume fraction of nanoparticles (f_NP) at different k_θ values. T_g tends to increase with f_NP at small k_θ, whereas it becomes suppressed at large k_θ.

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.

Supramolecular Self-Assembly of Dipalmitoylphosphatidylcholine and Carbon Nanotubes: A Dissipative Particle Dynamics Simulation Study

Dissipative particle dynamics simulations were performed to investigate the self-assembly of dipalmitoylphosphatidylcholine (DPPC) as a model lipid membrane on the surface of carbon nanotubes (CNTs). The influence of surface curvature of CNTs on self-assembly was investigated by performing simulations on solutions of DPPC in water in contact with CNTs of different diameters: CNT (10, 10), CNT (14, 14), CNT (20, 20), and CNT (34, 34). DPPC solutions with a wide range of concentrations were chosen to allow for formation of lipid structures of various surface densities, ranging from a submonolayer to a well-organized monolayer and a CNT covered with a lipid monolayer immersed in a planar lipid bilayer. Our results are indicative of a sequence of phase-ordering processes for DPPC on the surface of CNTs. At low surface coverages, the majority of hydrocarbon tail groups of DPPC are in contact with the CNT surface. Increasing the surface coverage leads to the formation of hemimicellar aggregates, and at high surface coverages close to the saturation limit, an organized lipid monolayer self-assembles. An examination of the mechanism of self-assembly reveals a two-step mechanism. The first step involves densification of DPPC on the CNT surface. Here, the lipid molecules do not adopt the order of the target phase (lipid monolayer on the CNT surface). In the second step, when the lipid density on the CNT surface is above a threshold value (close to saturation), the lipid molecules reorient themselves to form an organized monolayer around the tube. Here, the DPPC molecules adopt stretched conformations normal to the surface, the end hydrocarbon groups adsorb on the surface, and the head groups occupy the outermost part of the monolayer. The saturation density and the degree of lipid ordering on the CNT surface depend on the surface curvature. The saturation density increases with increased surface curvature, and better-ordered structures are formed on less curved surfaces.

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.

Infrared Linear Dichroism for the Analysis of Molecular Orientation in Polymers and in Polymer Composites

The mechanical properties of polymeric materials are strongly affected by molecular orientation occurring under processing conditions. Infrared dichroism is particularly well suited for characterizing polymer chain orientation at a molecular level. The usefulness of this technique has been demonstrated through various applications in homopolymers, semi-crystalline polymers, copolymers, polymer blends, as well as in polymer composites. Determination of molecular orientation can be carried out in the mid- and near-infrared ranges and very small dichroic effects can be detected with the use of a photoelastic modulator. Chain orientation in polymer composites is seen to increase with the filler content in the case of a strong interface between the two phases, making possible a quantification of the degree of bonding between the host polymeric matrix and the incorporated inclusions.

Polymerization and Structure of Opposing Polymer Brushes Studied by Computer Simulations

A model of the polymerization process during the formation of a pair of polymer brushes was designed and investigated. The obtained system consisted of two impenetrable parallel surfaces with the same number of chains grafted on both surfaces. Coarse-grained chains embedded in nodes of a face-centered cubic lattice with excluded volume interactions were obtained by a ‘grafted from’ procedure. The structure of synthesized macromolecular systems was also studied. Monte Carlo simulations using the dynamic lattice liquid model were employed using dedicated parallel machine ARUZ in a large size and time scale. The parameters of the polymerization process were found to be crucial for the proper structure of the brush. It was found that for high grafting densities, chains were increasingly compressed, and there is surprisingly little interpenetration of chains from opposite surfaces. It was predicted and confirmed that in a polydisperse sample, the longer chains have unique configurations consisting of a stretched stem and a coiled crown.

A High-Similarity Modeling Method for Low-Porosity Porous Material and Its Application in Bearing Cage Self-Lubrication Simulation

The porous oil-containing cage achieves the storage, spillage, and suction of lubricating oil by its micro-pore structure, thus ensuring the self-lubricating performance of the bearing. Carrying out fast and accurate modeling of the cage microscopic pore structure is the key to the analysis of the self-lubricating mechanism of bearings. In response to the issues where current modeling methods of porous materials have a low similarity of pore distribution, morphology, structure, and size characteristics, and the transition of pore surfaces is sharp, this paper proposed a modeling method of a highly similar micro-pore structure based on the idea of median filtering, the quartet structure generation set (QSGS), and the slice method. By extracting and analyzing the pore characteristics of the porous model and comparing them with the experimental results of CT scanning, the advantages of the modeling method in terms of morphology and pore connectivity were verified. Finally, by carrying out simulation analysis of the centrifugal force of oil splashing and capillary oil absorption on the constructed model by combining the parameters of porous structures such as porosity and tortuosity, the advantages of the modeling method in the construction of the porous model and multi-physical field analysis were further verified.

Dynamics of Opposing Polymer Brushes: A Computer Simulation Study

Opposing polymer brush systems were synthesized and investigated by molecular modeling. Chains were restricted to a face-centered cubic lattice with the excluded volume interactions only. The system was confined between two parallel impenetrable walls, with the same number of chains grafted to each surface. The dynamic properties of such systems were studied by Monte Carlo simulations based on the dynamic lattice liquid model and using a highly efficient parallel machine ARUZ, which enabled the study of large systems and long timescales. The influence of the surface density and mean polymer length on the system dynamic was discussed. The self-diffusion coefficient of the solvent depended strongly on the degree of polymerization and on the polymer concentration. It was also shown that it is possible to capture changes in solvent mobility that can be attributed to the regions of different polymer densities.

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.

Insights from modeling into structure, entanglements, and dynamics in attractive polymer nanocomposites

Conformations, entanglements and dynamics in attractive polymer nanocomposites are investigated in this work by means of coarse-grained molecular dynamics simulation, for both weak and strong confinements, in the presence of nanoparticles (NPs) at NP volume fractions φ up to 60%. We show that the behavior of the apparent tube diameter dapp in such nanocomposites can be greatly different from nanocomposites with nonattractive interactions. We find that this effect originates, based on a mean field argument, from the geometric confinement length dgeo at strong confinement (large φ) and not from the bound polymer layer on NPs (interparticle distance ID <2Rg) as proposed recently based on experimental measurements. Close to the NP surface, the entangled polymer mobility is reduced in attractive nanocomposites but still faster than the NP mobility for volume fractions beyond 20%. Furthermore, entangled polymer dynamics is hindered dramatically by the strong confinement created by NPs. For the first time using simulations, we show that the entangled polymer conformation, characterized by the polymer radius of gyration Rg and form factor, remains basically unperturbed by the presence of NPs up to the highest volume fractions studied, in agreement with various experiments on attractive nanocomposites. As a side-result we demonstrate that the loose concept of ID can be made a microscopically well defined quantity using the mean pore size of the NP arrangement.

Multiscale modelling of heterogeneous fillers in polymer composites: the case of polyisoprene and carbon black

The dispersion of inorganic particles within polymeric materials is an extensively used method to enhance their mechanical properties. One of the major challenges in the simulation of polymer composites is to model the uneven surface of the fillers which strongly affects the dynamics of the adsorbed polymers and consequently the macroscopic mechanical properties of the final composite. Here we propose a new multiscale approach that, using experimental adsorption data, constructs the filler surface to statistically reproduce the surface defects. We use this approach to analyse the structure and dynamics of highly entangled polyisoprene melt in contact with different realistic carbon black samples. We show that the presence of the heterogeneous surface has a negligible influence on the structure of the polymer chains but a major effect on their dynamics and the surface wettability.

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.

Dispersibility of Poly(vinyl acetate) Modified Silica Nanoparticles in Carbon Dioxide with Several Cosolvents

The dispersibility and stabilization of silica nanoparticles with surface-capped poly(vinyl acetate) (PVAc) chains are examined in carbon dioxide with four different cosolvents. Three surface coverages of silica-PVAc were formed by using different weight ratios of the silica and PVAc. The dispersibilities of three silica-PVAc nanoparticles in CO₂ with the four cosolvents were tested in a rotatable high-pressure variable-volume view cell. The effects of surface coverage, cosolvent type, pressure, and particle concentration on dispersion were investigated. Results show that, in the experimental pressure range (5.5 to 20 MPa), the pressure has no significant effect on the dispersion of nanoparticles, and the cosolvent is the key factor in dispersing silica-PVAc particles in CO₂. 1-Butanol is an adequate cosolvent to disperse silica-PVAc in CO₂ with any coverage of PVAc on the surface of the particles when the concentration of particles is smaller than 0.31 wt %. Ethanol can only improve the dispersibility of particles with a high surface coverage of PVAc when the concentration of particles is smaller than 0.14 wt %. 1-Hexanol and ethyl acetate cannot disperse the particles in CO₂ with any coverage of PVAc. Molecular dynamics simulations were carried out to study the nanoparticle-CO₂-cosolvent dispersions. Results suggest that 1-butanol has a good solubility in the CO₂ condensed phase and can effectively absorb onto the nanoparticle surface, which help to prevent the formation of nanoparticle aggregation. The precipitation of nanoparticles in the nanoparticle/1-hexanol/CO₂ and nanoparticle/ethyl acetate/CO₂ systems is attributed to the relatively low solubility of CO₂ in 1-hexanol and ethyl acetate. The precipitation of nanoparticles in the nanoparticle/ethanol/CO₂ system is the result of less hindrance of ethanol molecules to the aggregation of nanoparticles.

The Structure of Dipolar Polymer Brushes and Their Interaction in the Melt. Impact of Chain Stiffness

By using the numerical lattice Scheutjens-Fleer self-consistent field (SF-SCF) method we have studied the effect of the restricted flexibility of grafted chains on the structure and mutual interaction of two opposing planar conventional and A-type dipolar brushes. Brushes are immersed in the solvent consisting of chains similar to the grafted ones. The increase of the chain rigidity enhances the segregation of grafted chains in a A-type dipolar brush into two populations: backfolded chains with terminal monomers near the grafting surface and chains with the ends at the brush periphery. The fraction of backfolded chains grows by an increase of the Kuhn segment length. It is shown that two opposite A-type dipolar brushes from semi-rigid chains are attracted to each other at short distances. The attraction becomes more pronounced and begins at larger distances for more rigid chains with the same brush characteristics: polymerization degree, grafting density, and dipole moments of monomer units. This attraction is connected with the dipole-dipole interactions between chains of oncoming brushes with oppositely directed dipoles penetrating deeply into each other upon contact. This effect of the chain rigidity is opposite to that for conventional brushes without dipoles in the chains. For such brushes, an increase in the chain rigidity leads to the enhanced repulsion between them.

Assembly of Polymer-Grafted Nanoparticles in Polymer Matrices

Leibler pioneered the idea that long enough matrix polymers of length P will spontaneously dewet a chemically identical polymer layer, comprising chains of length N, densely end-grafted to a flat surface ("brush"). This entropically driven idea is routinely used to explain experiments in which 10-20 nm diameter nanoparticles (NPs) densely grafted with polymer chains are found to phase separate from chemically identical melts for P/N ≳4. At lower grafting densities, these effects are also thought to underpin the self-assembly of grafted NPs into a variety of structures. To explore the validity of this picture, we conducted large-scale molecular dynamics simulations of grafted NPs in a chemically identical polymer melt. For the NPs we consider, in the ≈5 nm diameter range, we find no phase separation even for P/N = 10 in the absence of attractions. Instead, we find behavior that more closely parallels experiments when all of the chain monomers are equally attractive to each other but repel the NPs. Our results thus imply that experimental situations investigated to date are dominated by the surfactancy of the NPs, which is driven by the chemical mismatch between the inorganic core and the organic ligands (the graft and free chains are chemically identical). Entropic effects, that is, the translational entropy of the NPs and the matrix, the entropy of mixing of the grafts and the matrix, and the conformational entropy of the chains appear to thus have a second-order effect even in the context of these model systems.

Molecular Dynamics Simulation of the Ionic Liquid 1-n-Butyl-3-Methylimidazolium Methylsulfate [Bmim][MeSO4 ]: Interfacial Properties at the Silica and Vacuum Interfaces

Molecular dynamics simulations are done to investigate the structure and dynamics of a thin [Bmim][MeO₄ ] film in contact with a hydroxylated silica surface on one side and with vacuum on the other. An examination of the microscopic structure of ionic liquid (IL) film shows that strong layered anionic/cationic structures are formed at both interfaces. At the silica interface, the imidazolium rings are closer to the silica surface (compared to anions) and are coplanar with it. At the vacuum interface, the charged imidazolium ring more concentrates in the interior of the film, but the butyl side chain stretches out toward the vacuum interface. While there exists an excess concentration of the cations at the silica interface, at the vacuum interface an excess concentration of anions (dissolved in the butyl chain) is found. The influence of the interface on the dynamical properties is shown to depend on their time scales. A short-time dynamical property, such as hydrogen bond formation is not noticeably perturbed at the interface. In contrary, long-time properties such as ion-pair formation/rupture and translation of ions across the film are largely decelerated at the silica interface but are accelerate at the vacuum interface. Our findings indicate that the structural relaxation time of ion-pairs, is comparable to diffusion time scale in the IL film. Therefore, ion-pairs are not stable species; the IL is composed of short-lived ion-pairs and freely diffusing ions. However, the structural relaxation times of ion-pairs is still long enough (comparable to the time scale of diffusion) to conclude that correlated motions of counterions influence the macroscopic properties of IL, such as diffusion and ionic conductivity. In this respect, we have shown that correcting the Nernst-Einstein equation for the joint translation of ion-pairs considerably improves the accuracy of calculated ionic conductivities.

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.

Optimizing PEGylation of TiO2 Nanocrystals through a Combined Experimental and Computational Study

PEGylation of metal oxide nanoparticles is the common approach to improve their biocompatibility and in vivo circulation time. In this work, we present a combined experimental and theoretical study to determine the operating condition that guarantee very high grafting densities, which are desirable in any biomedical application. Moreover, we present an insightful conformational analysis spanning different coverage regimes and increasing polymer chain lengths. Based on 13C NMR measurements and molecular dynamics simulations, we show that classical and popular models of polymer conformation on surfaces fail in determining the mushroom-to-brush transition point and prove that it actually takes place only at rather high grafting density values.

Temperature dependence of the interfacial bonding characteristics of silica/styrene butadiene rubber composites: a molecular dynamics simulation study

Based on our previous studies on the modification of in-chain styrene butadiene rubber (SBR) using 3-mercaptopropionic acid as well as its composites filled with silica, we further constructed two types of models (amorphous and layered) to investigate the temperature dependence of the interfacial bonding characteristics of silica/SBR composites via molecular dynamics (MD) simulation. The competing effects of rubber–rubber interactions and filler–rubber interactions were identified, and the relationship between the competing effects and the temperature was determined. Besides this, the effect of temperature on the mobility and distribution of SBR chains on the surface of silica was investigated. It was found that the stronger the interfacial interactions, the less sensitive the motion of SBR chains to temperature. Finally, the number and length of hydrogen bonds as a function of temperature were analyzed. These simulated results deepened the understanding of interface temperature dependence of the silica/SBR composites and gave a molecular level explanation for the existence of an optimum modifier content (14.2 wt%) that is temperature independent.

Temperature dependence of the interface between silica and styrene butadiene rubber modified by 3-mercaptopropionic acid was investigated by molecular dynamics simulation.

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.

Integral Equation Prediction of the Structure of Alternating Copolymer Nanocomposites near a Substrate

The packing structure and phase behavior of polymer-nanoparticle mixtures under confinement play an important role in developing strategies for rational design of nanomaterials. However, understanding the microscopic dispersion and aggregation mechanism of polymer nanocomposites is a great challenge through experimental techniques. In this work, the microscopic structure of alternating copolymer nanocomposites (ACNs) near a substrate is investigated systematically through extension of the inhomogeneous polymer reference interaction site model (PRISM) theory. In order to characterize the flexibility and internal chain stiffness of copolymers, a semiflexible chain model is introduced to describe the intramolecular correlations between different monomers. Based on the bridge functionals derived from the fluids density functional theory, the modified hypernetted chain closure is integrated with the PRISM equation to form a full theoretical framework to capture the density distributions of ACNs. The influence of the particle volume fraction, nanoparticle diameter, and adsorption strengths between different interaction sites on the packing structure of ACNs under confinement is analyzed and discussed in detail. With the increase of the particle volume fraction, the size asymmetry between nanoparticles and copolymer monomers can greatly influence the density profiles of ACNs near a substrate. Increasing the nanoparticle diameter, the density distribution of nanoparticles experiences a process from absorbing onto the solid surface to segregating from the wall to larger distances. With increasing the adsorption strength between copolymers and nanoparticles, the density distribution of nanoparticles decreases, which is similar to the case of nanoparticles containing attractive interactions. All these characteristics of ACNs show that the current inhomogeneous PRISM theory can give a detailed description of the packing behavior of different segments. Predictive approaches could be desired and developed for design control of alternating copolymer nanocomposites under confinement.

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.

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.

On the calculation of the potential of mean force between atomistic nanoparticles

We study the potential of mean force (PMF) between atomistic silica and gold nanoparticles in the vacuum by using molecular dynamics simulations. Such an investigation is devised in order to fully characterize the effective interactions between atomistic nanoparticles, a crucial step to describe the PMF in high-density coarse-grained polymer nanocomposites. In our study, we first investigate the behavior of silica nanoparticles, considering cases corresponding to different particle sizes and assessing results against an analytic theory developed by Hamaker for a system of Lennard-Jones interacting particles (H.C. Hamaker, Physica A 4, 1058 (1937)). Once validated the procedure, we calculate effective interactions between gold nanoparticles, which are considered both bare and coated with polyethylene chains, in order to investigate the effects of the grafting density [Formula: see text] on the PMF. Upon performing atomistic molecular dynamics simulations, it turns out that silica nanoparticles experience similar interactions regardless of the particle size, the most remarkable difference being a peak in the PMF due to surface interactions, clearly apparent for the larger size. As for bare gold nanoparticles, they are slightly interacting, the strength of the effective force increasing for the coated cases. The profile of the resulting PMF resembles a Lennard-Jones potential for intermediate [Formula: see text], becoming progressively more repulsive for high [Formula: see text] and low interparticle separations.

Vertically aligned double wall carbon nanotube arrays adsorbent for pure and mixture adsorption of H2S, ethylbenzene and carbon monoxide, grand canonical Monte Carlo simulation

In this study, pure and ternary adsorption of hydrogen sulfide (H₂S), ethylbenzene (EB), and carbon monoxide (CO) on different arrays of zigzag double wall carbon nanotube was investigated using grand canonical Monte Carlo simulations. The internal diameters of nanotube were fixed at 2r = 50.17 Å while nanotube wall distances were different values from d = 0 Å to d = 150 Å. Pure simulation results indicated that adsorption quantity of H₂S and EB in low pressure ranges of P = 1.9 bar to P = 3.1 bar was at least 100% more than CO adsorption quantities. At high pressure ranges of P = 23.1 bar to P = 38.2 bar H₂S adsorption was greater than EB and CO by about 200 molecules per unit cell (UC) at low nanotube distances. This was related to smaller kinetic diameter and greater dipole moment of H₂S compared to EB and CO. At higher nanotube distance the effect of size however disappears and all three gases approach to adsorption quantity of about 800 molecules/UC. Graphical representation of adsorption areas showed that H₂S and CO form multilayer adsorption around nanotube inner and outer walls while EB fill the whole space uniformly without any congestion around the walls. Ternary adsorption results EB/CO and H₂S/CO selectivity are greater than EB/H₂S selectivity. In addition, at smaller nanotube distances H₂S/CO selectivity is generally higher than EB/CO selectivity, which at higher nanotube distance the order becomes revers suggesting that size dependent effects on adsorption vanishes. Isosteric heat of adsorption shows that the order of EB > H₂S > CO suggesting that ethylbenzene interaction with nanotube arrays was strongest. Although H₂S has a greater dipole moment and smaller molecular dimension, EB adsorption at higher nanotube distance is greater than H₂S by at least 50% probably because EB is less volatile.

Electrical, Mechanical and Electromechanical Properties of Graphene-Thermoset Polymer Composites Produced Using Acetone-DMF Solvents

Recently, graphene-polymer composites gained a central role in advanced stress and strain sensing. A fundamental step in the production of epoxy-composites filled with graphene nanoplatelets (GNPs) consists in the exfoliation and dispersion of expanded graphite in a proper solvent, in the mixing of the resulting GNP suspension with the polymer matrix, and in the final removal of the solvent from the composite before curing through evaporation. The effects of traces of residual solvent on polymer curing process are usually overlooked, even if it has been found that even a small amount of residual solvent can affect the mechanical properties of the final composite. In this paper, we show that residual traces of N,N'-Dimethylformamide (DMF) in vinylester epoxy composites can induce relevant variations of the electrical, mechanical and electromechanical properties of the cured GNP-composite. To this purpose, a complete analysis of the morphological and structural characteristics of the composite samples produced using different solvent mixtures (combining acetone and DMF) is performed. Moreover, electrical, mechanical and electromechanical properties of the produced composites are assessed. In particular, the effect on the piezoresistive response of the use of DMF in the solvent mixture is analyzed using an experimental strain dependent percolation law to fit the measured electromechanical data. It is shown that the composites realized using a higher amount of DMF are characterized by a higher electrical conductivity and by a strong reduction of Young's Modulus.