Next Article in Journal
A Novel Surface Modification Strategy via Photopolymerized Poly-Sulfobetaine Methacrylate Coating to Prevent Bacterial Adhesion on Titanium Surfaces
Previous Article in Journal
Effect of Different Surface Treatments of Lithium Disilicate on the Adhesive Properties of Resin Cements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Percolated Network of Mixed Nanoparticles with Different Sizes in Polymer Nanocomposites: A Coarse-Grained Molecular Dynamics Simulation

1
State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, China
2
Key Laboratory of Beijing City on Preparation and Processing of Novel Polymer Materials, Beijing University of Chemical Technology, Beijing 100029, China
3
Beijing Engineering Research Center of Advanced Elastomers, Beijing University of Chemical Technology, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(12), 3301; https://doi.org/10.3390/ma14123301
Submission received: 13 May 2021 / Revised: 7 June 2021 / Accepted: 8 June 2021 / Published: 15 June 2021

Abstract

:
The size of real nanoparticles (NPs) is polydisperse which can influence the electrical property of polymer nanocomposites (PNCs). Here, we explored the percolated network of mixed NPs with different sizes (small NPs and big NPs) by adopting a molecular dynamics simulation. The simulated results reveal that the big NPs are adverse to building the percolated network compared to the small NPs. Thus, the percolation threshold becomes higher along with increasing the mixing ratio, which denotes the concentration ratio of big NPs to the total NPs. For a better understanding of it, the dispersion state and the number and the size of clusters are employed to analyze the percolated network, which can explain the percolation threshold well. Furthermore, by adopting the Sun’s theory (Macromolecules, 2009, 42, 459–463), small and big NPs exhibit a weak antagonistic effect in the simulation if their total concentration is fixed. On the one hand, the number of small NPs is larger than that of big NPs at the same concentration. In addition, one big NP can connect to more others than one small NP. These two contrast effects are responsible for it. Interestingly, the shear flow leads to more contact aggregation structure of NPs which is beneficial to build the new percolated networks. Especially, the big NPs play a more important role in forming the percolated network than small NPs. Consequently, the percolation threshold is reduced at a higher shear rate. In total, our research work provides a further understanding of how the mixed NPs with different sizes form the percolated network in polymer matrix.

1. Introduction

Conductive polymer nanocomposites (PNCs) are consistent with the polymer matrix and conductive particles (NPs) (graphene, carbon nanotube (CNT), carbon black (CB)), which have been applied in many fields (for example the electromagnetic interference shielding, sensor and conductors) [1,2]. According to the percolation theory, the percolated network is built if the volume fraction of NPs is beyond a critical percolation threshold [3]. As a result, the conductivity rises rapidly which makes the materials change from an insulator to a conductor. The concentration, the dispersion state of the conductive NPs as well as their shape will influence the distribution state of NPs, which indirectly affect the percolated network. It is proved that the high concentration of NPs is very necessary to build the percolated network, but it will damage the flexibility, elasticity and processing properties, which thus is necessary to lower the percolation threshold.
Currently, there are lots of works on improving the electrical conductivity of PNCs by tuning various experimental parameters. For instance, the NP surface is modified to improve their miscibility with chains which can reduce the percolation threshold [4]. Meanwhile, due to the branched morphology of high-structure CBs, they are helpful to form the percolated network than low-structure ones [5]. When the surface of CNTs is covered with the small and highly conductive particles, the percolation threshold can be further lowered [6]. Moreover, a percolation model is developed to predict the percolation threshold by inserting the aspect ratio of CNTs [7]. Hybrid NPs with different shapes are usually introduced to achieve a high conductivity in the polycarbonate matrix [8]. Furthermore, the addition of the large spherical silica inside the polymer matrix leads to a strong percolated network of CBs, which reduces the percolation threshold [2]. In addition, it is an effective method to modify the NP surface by grafted chains, which can tune their spatial distribution and optimize the electrical conductivity [9,10]. The external fields will change the percolated network which impacts the conductivity. It shows that the tensile field induces the nonuniform distribution of NPs, which forms more new percolated paths [11]. Moreover, the electrical property of PNCs in the shear flow can be lower or higher than that in the quiescent state, which is related to the initial network structure [12,13]. Interestingly, the electrical conductivity is improved by several magnitudes over time in the quiescent conditions, which is attributed to the NPs reaggregation [14]. By utilizing the change of the percolation networks under the strain/stress or pressure condition, the environment conditions can be monitored [15]. However, the percolated network is difficult to characterize accurately in experimental studies. Computer simulation is an important method which is good at characterizing the percolated network. By adopting the molecular dynamics simulation, the percolation threshold changes in an opposite direction with the polymer-NP interaction [16]. This is due to the appearance of the local bridging NPs via the polymer chains at a high interaction. By adopting a Monte Carlo simulation, the relationship between the percolation threshold and the aspect ratio of fibers can be described by an exponential function [17]. The index parameters varies from 2.0 to 2.6 for two-dimensional networks and from 1.7 to 2.0 for three-dimensional networks [18,19]. The percolation threshold exhibits a continuous decline with increasing the NP stiffness in the stable state while it reaches a minimum value in the shear flow [18,20]. The uniform distribution state improves the distance between NPs, which is not beneficial to form the percolated network than a partial contact structure [21]. At similar sizes, the rod or plate NPs are much easier to build the percolated network with than the spherical NPs [22]. As for grafted NPs, the distribution state and the phase behavior are two factors impacting the percolation threshold which is proved by model and theory [23,24]. Moreover, the shear field will change the initial percolated network and then another network will be reformed [21]. This process is affected by the shape and size of NPs and the shear strength [21,25]. Finally, the percolated network will be destroyed if the shear strength is beyond the certain rate [26].
Carbon black (CB) has excellent sensitivity to small strain and low cost, which is widely applied as a conductive filler [27,28]. Currently, the spherical NPs are assumed to be monodisperse in terms of their size. However, the real NPs exhibit various sizes, which will change the percolated network and the conductive property. Meanwhile, mixed NPs with different sizes may produce the synergistic or antagonistic effect which still remains unclear to our knowledge. To answer it, a molecular dynamics simulation is employed to explore the conductive behavior of the mixed NPs (small NPs and big NPs) filled PNCs. Especially, the respective roles of small NPs and big NPs in constructing the percolated network are discussed which is beneficial to uncover the conductive mechanism of PNCs.

2. Models

In this work, by adopting a coarse-grained model, first the initial configuration of PNCs is constructed where each polymer is consistent of 30 beads [29]. In addition, each system contains 1600 polymer chains. The m and σ denote the mass and diameter in each polymer bead, respectively. The size distribution of NPs is set to be the bimodal distribution, which is consistent with the theoretical or experimental works [30,31]. Thus, two types of spherical NPs with different sizes are adopted in the matrix. Considering the computational power, the diameters of small and big NPs are fixed to 1 σ and 4 σ , respectively [11,32].
The modified Lennard–Jones interaction potential is employed to model the interactions between different beads [29], given by
U i j ( r ) = { 4 ε i j [ ( σ r r E V ) 12 ( σ r r E V ) 6 ] r r E V < r c u t o f f 0 r r E V r c u t o f f
where r c u t o f f stands for the position where the interaction is truncated and shifted. The r E V denotes the excluded volume effect between the different beads. The interaction ε ij , the cutoff distance r c u t o f f and the r E V are presented in Table S1. These parameters can make the NPs relatively uniform distribution in the matrix, which can form the percolated network.
The chemical bonds within the polymer chains are simulated by the stiff finite extensible nonlinear elastic potential, given by
V F E N E   =   0.5 k R 0 2 l n [ 1 ( r R 0 ) 2 ]
where the k is set to be 30 ε σ 2 and R 0 is 1.5 σ , which are consistent with the previous works [29].
Since we did not focus on a specific polymer chain, the reduced units are adopted for all physical quantities, such as the energy units ε , the length units σ , the mass units m, and the time unit τ ( τ = ( m σ 2 / ε ) 1 / 2 ). At the beginning, the initial configuration of PNCs is constructed with a lower density in a simulation box. The simulation systems are then equilibrated with the isothermal-isobaric (NPT) ensemble for 2 × 105  τ with a timestep of δ t   =   0.001 τ . The P * and T * are set to be 0.0 and 1.0, respectively, which are controlled by the Nose–Hoover barostat and thermostat. The periodic boundary condition is applied during the whole simulation. It is noted that the polymer chains have been checked to experience fully relaxed for each system. The simulation process for calculating the percolated network has been described in the previous works [16,33], which is also shown in the supporting information. The diagram of a typical polymer nanocomposite is presented in Figure S1, which can help to understand the simulated systems. Then, the 10,000 configurations are obtained to calculate the conductive probability whose time interval is chosen to be 10 τ . The conductive probability is defined as the ratio of the number of the conductive configurations to that of all the configurations. Finally, the SLLOD methods are employed to simulate the shear flow, which is one widely used method for analyzing the shearing materials [34]. In addition, the special Lees–Edwards “sliding brick” boundary conditions [35] is adopted for the SLLOD equation. For more simulation details, refer to our previous works [16,33]. All MD runs are completed in the large scale atomic/molecular massively parallel simulator (LAMMPS) [36].

3. Results and Discussion

3.1. Ratio of Big Nanoparticles to the Total Nanoparticles

In general, the NPs are available in a wide variety of sizes in the real system, which will in turn influence the conductive behavior of PNCs. Thus, to understand the conductive behavior, we first explored the dependence of the conductive probability Λ on the α . Here, the α denotes the mixing ratio of big NPs, which is calculated by the concentration ratio of big NPs to the total NPs. Figure 1a presents the change of the conductive probability Λ with the concentration ( φ ) of NPs for different mixing ratio α . The concentration of NPs is defined to be volume ratio of all the NPs to the simulation box. As the increase of the φ , the percolated network slowly grows up and becomes completed. Then, the formation of the percolated network quickly improves the Λ from 0 to 1. Here, the percolation threshold φ c is defined as the volume fraction of NPs at Λ = 0.5, which is calculated by fitting the curves in Figure 1a with a hyperbolic tangent equation: Λ   =   0.5 + 0.5 tanh ( 2 ( φ φ c ) / d c ) where dc denotes the percolation width [37,38]. As shown in Figure 1b, the φ c exhibits a monotonous rise with the increase of the α which reflects that small NPs can build the percolated network at a relatively low concentration. This is because the number of small NPs is much larger than that of big NPs at the same concentration, which helps to build the percolated network.
To further understand the results, the distribution state of NPs is analyzed by the radial distribution function (RDF). The peak position at r = 1 σ for α     1.0 or r   =   4 σ for α   =   1.0 means the direct contact aggregation structure of NPs. The peak position at r = 2 σ for α   1.0 or r   =   5 σ for α   =   1.0 denotes the sandwiched structures of NPs. From Figure S2a, the RDF curves are similar for α < 1.0 due to more small NPs than big ones. Meanwhile, the low height of peaks in RDFs and the diagrams of NPs in Figure S2b can visually reflect a relatively uniform dispersion. The distribution behavior of particles can influence the conductive probability which however is more dependent on the percolated network. To analyze the percolated network, two important parameters are introduced, namely, the largest size C n and the number N c of clusters. The dependence of the C n and the N c on the φ for different α is presented in Figure 2.
It is found that the C n rises slowly at first, and then, the rising speed becomes larger when the φ is beyond a critical value. Corresponding to it, the N c exhibits a nonmonotonic change, which reaches the maximum at the moderate φ . This indicates that the initial small clusters emerge to become large clusters with increasing the φ . Meanwhile, N c is reduced as α goes up because of the reduced number of NPs. Thus, the C n is more sensitive and accurate than the N c to represent the percolated network. Furthermore, the diagrams of the largest clusters are shown in Figure 3 for various α and φ , which reflects the formation process of percolated networks. The largest cluster is isolated at low φ and then becomes larger with the φ . At last, there appears a percolated network which leads to the high Λ .
To describe more in depth their respective roles in constructing the percolated network, the ratio of big NPs in the largest cluster to the total big NPs is calculated in respect of the mixing ratio α in Figure 4.
Figure 4 also contains the concentration ratio α of big NPs to the whole NPs in the matrix (namely α ), which presents a better comparison. It is found that this ratio (red line) rises with the α which however is a weakly higher than the α (black line). This means that big NPs take part more in constructing the percolated network than small ones. This is mainly due to the large size of big NPs which can connect more other NPs. However, this effect is very weak which is due to the uniform dispersion of NPs. Furthermore, according to the Sun’ theory [20], the percolation threshold can be calculated according to the excluded volume theory, given by
φ S N φ 0 S N + φ B N φ 0 B N = 1
where φ S N and φ B N are the respective concentrations of small and big NPs in the ternary systems when the system reaches the percolation state. Namely, the total concentration of NPs is the percolation threshold which is shown in Figure 1b. φ 0 S N and φ 0 B N are the percolation thresholds of one kind of NP, respectively. If the simulated result is below Sun’s line, it reflects the synergy. Conversely, it represents the antagonism.
Figure 5 shows that the simulation values are slightly beyond that in Equation (3), which indicates a weak antagonistic effect. This is because on the one hand, the number of small NPs is larger than that of big NPs at the same concentration. On the other hand, one big NP can connect more other ones than one small NP. This can rationalize the obtained results. Furthermore, we analyzed the impact of the size polydispersity of NPs on the conductive probability Λ to describe more in depth their impact on building the percolated network. First, the corresponding averaged diameter (DNP) of NPs is calculated for different α in Figure S3. The DNP is calculated as the following equation:
D N P 3 ( N 1 + N 2 ) = D 1 3 N 1 + D 2 3 N 2
where N 1 and D 1 are the number and diameter of big NPs while N 2 and D 2 are the number and diameter of small NPs. The calculated Λ is presented in Figure 6a for the monodisperse NPs (namely, polydispersity index = 1.0). Similarly, the percolation threshold φ c is obtained which is shown in Figure 6b. It indicates that the φ c for monodisperse NPs is lower than that for polydisperse NPs. This indicates that the NPs with different sizes do not exhibit the synergistic effect which is consistent with the previous conclusion. However, it has reported a synergetic percolation for mixed fillers in experiments [39,40,41]. In these systems, the synergetic effect is realized by adding the additional second fillers into the first fillers which means that the total concentration of fillers increases. Figure 7 presents the change of the Λ with the φ B N at three fixed concentrations φ S N of small NPs as examples. The φ c is gradually reduced from the 27.3% to 24.5% with increasing the φ S N , which proves the experimental results. In total, the percolation threshold gradually rises by increasing the mixing ratio of big NPs. However, small and big NPs exhibit a weak antagonistic effect if their total concentration is fixed.

3.2. Shear Field

In this section, we aimed to investigate the effect of the shear flow on the percolated network. First, Figure 8a shows the conductive probability Λ with the α for different φ by fixing the shear rate γ ˙   = 0.1. The obtained percolation threshold φ c is put in Table S2 which slowly increases with increasing the α at γ ˙   = 0.1. Meanwhile, the φ c at γ   ˙ = 0.1 is lower than that at γ ˙   = 0.0. To understand it, Figure 8b,c presents the directional conductive probabilities Λ and Λ . The difference of φ c is small for Λ and Λ in Figure S4. This is because the spherical NPs is isotropic along the shear direction. Then, we analyzed how the distribution state of NPs is changed by the shear flow where the α   =   1.0 is chosen as an example in Figure S5. From the RDFs, the position at r   =   4 σ rises, reflecting that NPs appear more direct contact structures under the shear field. Meanwhile, one big NP can connect more other NPs than one small NP due to the size effect. Thus, both the direct contact aggregation structure and the size effect promote that big NPs participate more in forming the percolated network than small NPs. Namely, the concentration ratio of big NPs is higher at γ   ˙ = 0.1 than that at γ   ˙ = 0.0 which is presented in Figure 4. Furthermore, as shown in Figure S6, the largest cluster size C n gradually decreases with increasing the α which is in accord with φ c . Some diagrams of the largest clusters are presented in Figure 9 at γ ˙   = 0.0 and 0.1, which clearly presents the change of the percolated network. It is found that the C n is enhanced by the shear flow. In summary, more new percolated paths appear under the shear flow, which reduces the φ c . At last, the α = 0.75 is taken as an example to clarify the relationship between the Λ and the γ ˙ . The change of the Λ with the γ ˙ is calculated in Figure 10. The φ c is reduced from 13.0% to 11.3% with increasing the γ ˙ from 0.0 to 0.5. Then, both the Λ and Λ are calculated in Figure S7 which gradually rise with increasing the γ ˙ . This is attributed to more direct contact structures of NPs induced by the shear field. Furthermore, the C n is analyzed to characterize the percolated network in Figure 11, which is consistent with Λ . In total, the conductive property is improved with increasing the shear strength.

3.3. Discussion

When mapping this simulation systems to real polymers, the interaction energy ε is roughly 2.5–3.4 kJ·mol−1 for the common polymers [29]. The interfacial interaction is chosen to be 2.0 in this model, which means that the interaction energy is 5.0–6.8 kJ·mol−1. The average binding energy between polymer chains and silica particles is about 4.2–24 kJ·mol−1, which depends on the filler surface activity and chemical functionality [42,43]. Our simulated interfacial interaction is roughly within the experimental range. Meanwhile, the simulated persistence length is 0.68 σ while it is 0.35–0.76 nm for real polymers [44]. Thus, the σ  is roughly 1 nm. Furthermore, the diffusion coefficients (D) of NPs with different sizes are calculated which are 1.52 × 10−2 τ 1  and 1.02 × 10−3 τ 1 , respectively. The simulated shear strength  γ ˙ is 0.02–0.5 τ 1 which is larger than the experimental values. Thus, the range of the Peclet number ( P e = γ D ) is from 1.3 to 490, which can be comparable to the experimental value [45]. Moreover, the dispersion states and network of NPs have also explored by this model [46,47,48,49], which are consistent with the theoretical and experimental results [50,51]. Thus, our models are relatively reasonable.
The conductive property of PNCs has been investigated, which is determined by the size, shape, volume fraction of NPs, the interfacial interaction and so on. For example, the dispersion state of rod NPs first changes from the aggregation to dispersion and then forms the local bridging structure. However, the percolation threshold shows an anti N-type with increasing the interfacial interaction [16]. However, the percolation threshold is reduced by increasing the interfacial interaction for spherical NPs [52]. The uniform distribution state improves the distances between NPs, which is no more beneficial to form the percolated network than a partial contact structure [21]. In addition, the dependence of the fiber aspect ratio on the percolation threshold follows the exponential function [17]. Especially, the index parameters vary from 2.0 to 2.6 for two-dimensional networks and from 1.7 to 2.0 for three-dimensional networks [18,19]. Moreover, the shape of NPs influences the formation of the percolated network, which reveals that the rod or plate NPs are more effective to improve the conductive probability than the spherical NPs [22]. The high stiffness of NPs can reduce the percolation threshold. However, it is minimum under the shear flow [18,20]. In this work, we mainly investigated the effect of mixed NPs with different sizes on the conductive probability of PNCs. It is found that there is a weak antagonistic effect for two kinds of NPs with different sizes at the fixed concentration of NPs. Meanwhile, the shear field is beneficial to reduce the percolation threshold compared with the quiescent state, which is attributed to the change of the percolated network. Moreover, in experiments, the effects of the polymer-NP interaction [4], the structures and aspect ratio of NPs [5,7], the chemical grafting modification of NPs [9,10] and the hybrid NPs with different shapes [8] on the conductive property of PNCs have been investigated. Moreover, the experiment reported that the percolation threshold of PNCs declines with the decrease of the filler size [53,54], which is consistent with our simulation result. However, the effect of mixed NPs on the formation of the percolated network has not been investigated in experiments to our knowledge. This maybe because the monodisperse NPs are difficult to obtain in experiments where NPs are normally polydisperse. Meanwhile, the size of NPs is difficult to be accurately determined. It is noted that some details of our model are not different with those of experimental reports, such as the size of polymer chains and NPs and the distribution state of NPs. Meanwhile, the conductive probability is used to denote the conductive property rather than the electrical conductivity. However, we are confident that some aspects of our model capture the behavior of the experimental systems. In particular, we simulated the interactions between polymer chains and NPs, which are within the experimental range. In addition, the size distribution of NPs in the model is consistent with that in experiments. Meanwhile, the simulated shear field is the same as that of the experiments. Lastly, how the surface size of NPs affects the formation of the percolated network is also an interesting topic, which deserves to be investigated in our further work.

4. Conclusions

In this work, a molecular dynamics simulation is employed to explore the formation of percolated networks of mixed nanoparticles (NPs) with different sizes (small NPs and big NPs) in polymer nanocomposites (PNCs). It is found that the big NPs are adverse for building the percolated network. As the increase of the concentration ratio of big NPs to the whole NPs, the percolation threshold rises from 5.04% to 27.3%. Furthermore, the percolated network is characterized by analyzing the dispersion state, the largest size and the number of clusters, which rationalizes the percolation threshold. Meanwhile, a weak antagonistic effect is observed for small and big NPs at the fixed concentration of NPs. On one hand, the number of small NPs is larger than that of big NPs at the same concentration. Meanwhile, one big NP can connect more other NPs than one small NP. These two contrast effects can explain the results. Furthermore, the shear field is beneficial to improve the conductive probability compared with the quiescent state. This is because the NPs form more direct contact aggregation structures which can link other particles to build the new percolated network. Therefore, the big NPs participate more in forming the new percolated paths than small NPs. Finally, the percolation threshold is lowered with increasing the shear strength. In total, our work provides further understanding of the effect of the mixed NPs with different sizes on building the percolation network of PNCs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14123301/s1, Figure S1. (a) The radial distribution function for different mixing ratios ( α ) where the concentration φ of nanoparticles (NPs) is their percolation threshold. (b) Diagrams of NPs where the polymer chains are neglected for clarity, Figure S2. The averaged diameter (DNP) of nanoparticles for different mixing ratios α , Figure S3. The percolation threshold φ c for the parallel directional conductive probability Λ and perpendicular directional conductive probability Λ in respect of the mixing ratio α , Figure S4. The radial distribution function of nanoparticles for the mixing ratio α = 1.0, Figure S5. The largest size C n of cluster as a function of the concentration φ of nanoparticles for different mixing ratios ( α ), Figure S6. (a) The parallelly directional conductive probability Λ and (b) perpendicularly directional conductive probability Λ as a function of the concentration φ of nanoparticles for different shear rates ( γ ˙ ), Table S1 All the interaction parameter ε ij , the cutoff distance r c u t o f f and the r E V .

Author Contributions

Conceptualization, X.Z. and Y.N.; methodology, H.L.; software, H.W.; investigation, X.Z. and Y.N.; writing—original draft preparation, Y.G.; writing—review and editing, Y.G., L.Z. and X.Z.; supervision, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Beijing Natural Science Foundation, grant number 2212043, and National Natural Science Foundation of China, grant number 51973012.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

This work is supported by the Beijing Natural Science Foundation (2212043) and National Natural Science Foundation of China (51973012). The authors acknowledge the National Supercomputer Center in Guangzhou, Lvliang and Shenzhen.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alig, I.; Pötschke, P.; Lellinger, D.; Skipa, T.; Pegel, S.; Kasaliwal, G.R.; Villmow, T. Establishment, morphology and properties of carbon nanotube networks in polymer melts. Polymer 2012, 53, 4–28. [Google Scholar] [CrossRef]
  2. Bhagavatheswaran, E.S.; Parsekar, M.; Das, A.; Le, H.H.; Wiessner, S.; Stöckelhuber, K.W.; Schmaucks, G.; Heinrich, G. Construction of an Interconnected Nanostructured Carbon Black Network: Development of Highly Stretchable and Robust Elastomeric Conductors. J. Phys. Chem. C 2015, 119, 21723–21731. [Google Scholar] [CrossRef]
  3. Balberg, I. Recent developments in continuum percolation. Philos. Mag. B 1987, 56, 991–1003. [Google Scholar] [CrossRef]
  4. Bréchet, Y.; Cavaillé, J.Y.; Chabert, E.; Chazeau, L.; Dendievel, R.; Flandin, L.; Gauthier, C. Polymer Based Nanocomposites: Effect of Fill-Filler and Filler-Matrix Interactions. Adv. Eng. Mater. 2001, 3, 571–577. [Google Scholar] [CrossRef]
  5. Ren, D.; Zheng, S.; Huang, S.; Liu, Z.; Yang, M. Effect of the carbon black structure on the stability and efficiency of the conductive network in polyethylene composites. J. Appl. Polym. Sci. 2013, 129, 3382–3389. [Google Scholar] [CrossRef]
  6. Wang, Y.; Wang, X.; Cao, X.; Gong, S.; Xie, Z.; Li, T.; Wu, C.; Zhu, Z.; Li, Z. Effect of nano-scale Cu particles on the electrical property of CNT/polymer nanocomposites. Compos. Part A Appl. Sci. Manuf. 2021, 143, 106325. [Google Scholar] [CrossRef]
  7. White, S.I.; Mutiso, R.M.; Vora, P.; Jahnke, D.; Hsu, S.; Kikkawa, J.M.; Li, J.; Fischer, J.E.; Winey, K.I. Electrical Percolation Behavior in Silver Nanowire-Polystyrene Composites: Simulation and Experiment. Adv. Funct. Mater. 2010, 20, 2709–2716. [Google Scholar] [CrossRef]
  8. Hilarius, K.; Lellinger, D.; Alig, I.; Villmow, T.; Pegel, S.; Pötschke, P. Influence of shear deformation on the electrical and rheological properties of combined filler networks in polymer melts: Carbon nanotubes and carbon black in polycarbonate. Polymer 2013, 54, 5865–5874. [Google Scholar] [CrossRef]
  9. Srivastava, S.; Agarwal, P.; Archer, L.A. Tethered Nanoparticle–Polymer Composites: Phase Stability and Curvature. Langmuir 2012, 28, 6276–6281. [Google Scholar] [CrossRef] [PubMed]
  10. Chevigny, C.; Dalmas, F.; Di Cola, E.; Gigmes, D.; Bertin, D.; Boue, F.; Jestin, J. Polymer-Grafted-Nanoparticles Nanocomposites: Dispersion, Grafted Chain Conformation, and Rheological Behavior. Macromolecules 2011, 44, 122–133. [Google Scholar] [CrossRef]
  11. Ning, N.; Miao, C.; Zou, H.; Shao, Q.; Wang, S.; Zhang, L.; Tian, M. A new insight on the variation of the electric conductivity and conductive network of silver-coated glass particles/silicone elastomer composites under tensile strain. Compos. Sci. Technol. 2016, 136, 46–52. [Google Scholar] [CrossRef]
  12. Skipa, T.; Lellinger, D.; Böhm, W.; Saphiannikova, M.; Alig, I. Influence of shear deformation on carbon nanotube networks in polycarbonate melts: Interplay between build-up and destruction of agglomerates. Polymer 2010, 51, 201–210. [Google Scholar] [CrossRef]
  13. Alig, I.; Lellinger, D.; Engel, M.; Skipa, T.; Pötschke, P. Destruction and formation of a conductive carbon nanotube network in polymer melts: In-line experiments. Polymer 2008, 49, 1902–1909. [Google Scholar] [CrossRef]
  14. Alig, I.; Skipa, T.; Engel, M.; Lellinger, D.; Pegel, S.; Pötschke, P. Electrical conductivity recovery in carbon nanotube–polymer composites after transient shear. Phys. Status Solidi (b) 2007, 244, 4223–4226. [Google Scholar] [CrossRef]
  15. Chen, J.; Zhu, Y.; Huang, J.; Zhang, J.; Pan, D.; Zhou, J.; Ryu, J.E.; Umar, A.; Guo, Z. Advances in Responsively Conductive Polymer Composites and Sensing Applications. Polym. Rev. 2020, 61, 157–193. [Google Scholar] [CrossRef]
  16. Gao, Y.; Cao, D.; Liu, J.; Shen, J.; Wu, Y.; Zhang, L. Molecular dynamics simulation of the conductivity mechanism of nanorod filled polymer nanocomposites. Phys. Chem. Chem. Phys. 2015, 17, 22959–22968. [Google Scholar] [CrossRef]
  17. Ma, H.; Gao, X.-L. A three-dimensional Monte Carlo model for electrically conductive polymer matrix composites filled with curved fibers. Polymer 2008, 49, 4230–4238. [Google Scholar] [CrossRef]
  18. Li, C.; Thostenson, E.T.; Chou, T.-W. Effect of nanotube waviness on the electrical conductivity of carbon nanotube-based composites. Compos. Sci. Technol. 2008, 68, 1445–1452. [Google Scholar] [CrossRef]
  19. Dalmas, F.; Dendievel, R.; Chazeau, L.; Cavaillé, J.-Y.; Gauthier, C. Carbon nanotube-filled polymer composites. Numerical simulation of electrical conductivity in three-dimensional entangled fibrous networks. Acta Mater. 2006, 54, 2923–2931. [Google Scholar] [CrossRef]
  20. Sun, Y.; Bao, H.-D.; Guo, Z.-X.; Yu, J. Modeling of the Electrical Percolation of Mixed Carbon Fillers in Polymer-Based Composites. Macromolecules 2009, 42, 459–463. [Google Scholar] [CrossRef]
  21. Eken, A.; Tozzi, E.; Klingenberg, D.; Bauhofer, W. A simulation study on the effects of shear flow on the microstructure and electrical properties of carbon nanotube/polymer composites. Polymer 2011, 52, 5178–5185. [Google Scholar] [CrossRef]
  22. Kwon, S.; Cho, H.W.; Gwon, G.; Kim, H.; Sung, B.J. Effects of shape and flexibility of conductive fillers in nanocomposites on percolating network formation and electrical conductivity. Phys. Rev. E 2016, 93, 032501. [Google Scholar] [CrossRef] [PubMed]
  23. Pryamtisyn, V.; Ganesan, V.; Panagiotopoulos, A.; Liu, H.; Kumar, S.K. Modeling the anisotropic self-assembly of spherical polymer-grafted nanoparticles. J. Chem. Phys. 2009, 131, 221102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Jayaraman, A.; Schweizer, K.S. Effective Interactions and Self-Assembly of Hybrid Polymer Grafted Nanoparticles in a Homopolymer Matrix. Macromolecules 2009, 42, 8423–8434. [Google Scholar] [CrossRef]
  25. Eken, A.; Tozzi, E.; Klingenberg, D.; Bauhofer, W. Combined effects of nanotube aspect ratio and shear rate on the carbon nanotube/polymer composites. Polymer 2012, 53, 4493–4500. [Google Scholar] [CrossRef]
  26. Kwon, G.; Heo, Y.; Shin, K.; Sung, B.J. Electrical percolation networks of carbon nanotubes in a shear flow. Phys. Rev. E 2012, 85, 011143. [Google Scholar] [CrossRef]
  27. Huynh, M.T.T.; Cho, H.-B.; Suzuki, T.; Suematsu, H.; Nguyen, S.T.; Niihara, K.; Nakayama, T. Electrical property enhancement by controlled percolation structure of carbon black in polymer-based nanocomposites via nanosecond pulsed electric field. Compos. Sci. Technol. 2018, 154, 165–174. [Google Scholar] [CrossRef]
  28. Liu, X.; Krückel, J.; Zheng, G.; Schubert, D.W. Electrical conductivity behaviour of sheared poly(methyl methacrylate)/carbon black composites. Compos. Sci. Technol. 2014, 100, 99–104. [Google Scholar] [CrossRef]
  29. Kremer, K.; Grest, G.S. Dynamics of entangled linear polymer melts: A molecular-dynamics simulation. J. Chem. Phys. 1990, 92, 5057–5086. [Google Scholar] [CrossRef]
  30. Majidian, M.; Grimaldi, C.; Forró, L.; Magrez, A. Role of the particle size polydispersity in the electrical conductivity of carbon nanotube-epoxy composites. Sci. Rep. 2017, 7, 12553. [Google Scholar] [CrossRef] [Green Version]
  31. Meyer, H.; Van Der Schoot, P.; Schilling, T. Percolation in suspensions of polydisperse hard rods: Quasi universality and finite-size effects. J. Chem. Phys. 2015, 143, 044901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Li, H.; Wu, H.; Li, B.; Gao, Y.; Zhao, X.; Zhang, L. Molecular dynamics simulation of fracture mechanism in the double inter-penetrated cross-linked polymer. Polymer 2020, 199, 122571. [Google Scholar] [CrossRef]
  33. Wu, H.; Li, H.; Zhang, W.; Li, F.; Li, B.; Gao, Y.; Zhao, X.; Zhang, L. Percolation of polydisperse nanorods in polymer nanocomposites: Insights from molecular dynamics simulation. Compos. Sci. Technol. 2020, 196, 108208. [Google Scholar] [CrossRef]
  34. Tuckerman, M.E.; Mundy, C.; Balasubramanian, S.; Klein, M.L. Modified nonequilibrium molecular dynamics for fluid flows with energy conservation. J. Chem. Phys. 1997, 106, 5615–5621. [Google Scholar] [CrossRef]
  35. Lees, A.; Edwards, S.F. The computer study of transport processes under extreme conditions. J. Phys. C Solid State Phys. 1972, 5, 1921. [Google Scholar] [CrossRef]
  36. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  37. Chen, Y.; Wang, S.; Pan, F.; Zhang, J. A Numerical Study on Electrical Percolation of Polymer-Matrix Composites with Hybrid Fillers of Carbon Nanotubes and Carbon Black. J. Nanomater. 2014, 2014, 614797. [Google Scholar] [CrossRef] [Green Version]
  38. Nam, S.; Cho, H.W.; Lim, S.; Kim, D.; Kim, H.; Sung, B.J. Enhancement of Electrical and Thermomechanical Properties of Silver Nanowire Composites by the Introduction of Nonconductive Nanoparticles: Experiment and Simulation. ACS Nano 2013, 7, 851–856. [Google Scholar] [CrossRef]
  39. Ma, P.-C.; Liu, M.-Y.; Zhang, H.; Wang, S.-Q.; Wang, R.; Wang, K.; Wong, Y.K.; Tang, B.Z.; Hong, S.H.; Paik, K.W.; et al. Enhanced Electrical Conductivity of Nanocomposites Containing Hybrid Fillers of Carbon Nanotubes and Carbon Black. ACS Appl. Mater. Interfaces 2009, 1, 1090–1096. [Google Scholar] [CrossRef]
  40. Palza, H.; Garzon, C.; Arias, O. Modifying the electrical behaviour of polypropylene/carbon nanotube composites by adding a second nanoparticle and by annealing processes. Express Polymer Lett. 2012, 6, 639–646. [Google Scholar] [CrossRef] [Green Version]
  41. Socher, R.; Krause, B.; Hermasch, S.; Wursche, R.; Pötschke, P. Electrical and thermal properties of polyamide 12 composites with hybrid fillers systems of multiwalled carbon nanotubes and carbon black. Compos. Sci. Technol. 2011, 71, 1053–1059. [Google Scholar] [CrossRef] [Green Version]
  42. Hanson, D.E.; Hawley, M.; Houlton, R.; Chitanvis, K.; Rae, P.; Orler, E.B.; Wrobleski, D.A. Stress softening experiments in silica-filled poly-dimethylsiloxane provide insight into a mechanism for the Mullins effect. Polymer 2005, 46, 10989–10995. [Google Scholar] [CrossRef]
  43. Bueche, F. Molecular basis for the mullins effect. J. Appl. Polym. Sci. 1960, 4, 107–114. [Google Scholar] [CrossRef]
  44. Hiemenz, P.C.; Lodge, T.P. Polymer Chemistry; CRC Press: London, UK, 2007; p. 587. [Google Scholar]
  45. Larson, R.G. The Structure and Rheology of Complex Fluids; Oxford University Press: New York, NY, USA, 1999; p. 663. [Google Scholar]
  46. Gao, Y.; Hu, F.; Wu, Y.; Liu, J.; Zhang, L. Understanding the structural evolution under the oscillatory shear field to determine the viscoelastic behavior of nanorod filled polymer nanocomposites. Comput. Mater. Sci. 2018, 142, 192–199. [Google Scholar] [CrossRef]
  47. Shen, J.; Liu, J.; Li, H.; Gao, Y.; Li, X.; Wu, Y.; Zhang, L. Molecular dynamics simulations of the structural, mechanical and visco-elastic properties of polymer nanocomposites filled with grafted nanoparticles. Phys. Chem. Chem. Phys. 2015, 17, 7196–7207. [Google Scholar] [CrossRef]
  48. Shen, J.; Liu, J.; Gao, Y.; Li, X.; Zhang, L. Elucidating and tuning the strain-induced non-linear behavior of polymer nanocomposites: A detailed molecular dynamics simulation study. Soft Matter. 2014, 10, 5099–5113. [Google Scholar] [CrossRef]
  49. Liu, J.; Gao, Y.; Cao, D.; Zhang, L.; Guo, Z. Nanoparticle Dispersion and Aggregation in Polymer Nanocomposites: Insights from Molecular Dynamics Simulation. Langmuir 2011, 27, 7926–7933. [Google Scholar] [CrossRef] [PubMed]
  50. Meli, L.; Green, P.F. Aggregation and Coarsening of Ligand-Stabilized Gold Nanoparticles in Poly(methyl methacrylate) Thin Films. ACS Nano 2008, 2, 1305–1312. [Google Scholar] [CrossRef]
  51. Hooper, J.B.; Schweizer, K.S. Contact Aggregation, Bridging, and Steric Stabilization in Dense Polymer−Particle Mixtures. Macromolecules 2005, 38, 8858–8869. [Google Scholar] [CrossRef]
  52. Feng, Y.; Zou, H.; Tian, M.; Zhang, L.; Mi, J. Relationship between Dispersion and Conductivity of Polymer Nanocomposites: A Molecular Dynamics Study. J. Phys. Chem. B 2012, 116, 13081–13088. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, P.; Ding, T. Conductivity and piezoresistivity of conductive carbon black filled polymer composite. J. Appl. Polym. Sci. 2010, 116, 2035–2039. [Google Scholar] [CrossRef]
  54. Huang, Y.; Wang, W.; Zeng, X.; Guo, X.; Zhang, Y.; Liu, P.; Ma, Y.; Zhang, Y. Effects of the filler size on the electrical percolation threshold of carbon black-carbon nanotube-polymer composites. J. Appl. Polym. Sci. 2018, 135, 46517. [Google Scholar] [CrossRef]
Figure 1. (a) Conductive probability Λ with the concentration φ of nanoparticles and (b) the percolation threshold φ c in respect of the mixing ratio ( α ).
Figure 1. (a) Conductive probability Λ with the concentration φ of nanoparticles and (b) the percolation threshold φ c in respect of the mixing ratio ( α ).
Materials 14 03301 g001
Figure 2. (a) The largest size C n and (b) the number N c of clusters with respect to the concentration φ of nanoparticles for various mixing ratios ( α ).
Figure 2. (a) The largest size C n and (b) the number N c of clusters with respect to the concentration φ of nanoparticles for various mixing ratios ( α ).
Materials 14 03301 g002
Figure 3. Diagrams of the largest cluster at four concentrations φ of nanoparticles (NPs) for different mixing ratios ( α ).
Figure 3. Diagrams of the largest cluster at four concentrations φ of nanoparticles (NPs) for different mixing ratios ( α ).
Materials 14 03301 g003
Figure 4. Ratio of big nanoparticles (NPs) to the whole NPs in the largest cluster for two shear rates γ ˙ for different mixing ratios ( α ).
Figure 4. Ratio of big nanoparticles (NPs) to the whole NPs in the largest cluster for two shear rates γ ˙ for different mixing ratios ( α ).
Materials 14 03301 g004
Figure 5. The dependence of the concentration of small nanoparticles (NPs) ( φ S N ) on that of big NPs ( φ B N ) at the percolation threshold.
Figure 5. The dependence of the concentration of small nanoparticles (NPs) ( φ S N ) on that of big NPs ( φ B N ) at the percolation threshold.
Materials 14 03301 g005
Figure 6. (a) Conductive probability Λ with the concentration φ of monodisperse nanoparticles and (b) their percolation threshold φ c in respect of the mixing ratio (α).
Figure 6. (a) Conductive probability Λ with the concentration φ of monodisperse nanoparticles and (b) their percolation threshold φ c in respect of the mixing ratio (α).
Materials 14 03301 g006
Figure 7. Conductive probability Λ with the concentration φ of big nanoparticles (BN) for three fixed concentrations of small nanoparticles (SN).
Figure 7. Conductive probability Λ with the concentration φ of big nanoparticles (BN) for three fixed concentrations of small nanoparticles (SN).
Materials 14 03301 g007
Figure 8. (a) Conductive probability Λ , (b) the directional Λ parallel to the shear direction, and (c) the directional Λ perpendicular to the shear direction with the concentration φ of nanoparticles for different ratio ( α ).
Figure 8. (a) Conductive probability Λ , (b) the directional Λ parallel to the shear direction, and (c) the directional Λ perpendicular to the shear direction with the concentration φ of nanoparticles for different ratio ( α ).
Materials 14 03301 g008
Figure 9. Snapshots of the largest clusters at two shear rates γ   ˙ = 0.0 and 0.1 for different ratio ( α ).
Figure 9. Snapshots of the largest clusters at two shear rates γ   ˙ = 0.0 and 0.1 for different ratio ( α ).
Materials 14 03301 g009
Figure 10. Conductive probability Λ with the concentration φ of nanoparticles for different shear rates ( γ ˙ ).
Figure 10. Conductive probability Λ with the concentration φ of nanoparticles for different shear rates ( γ ˙ ).
Materials 14 03301 g010
Figure 11. The largest cluster size C n with the concentration φ of nanoparticles for different shear rates ( γ ˙ ).
Figure 11. The largest cluster size C n with the concentration φ of nanoparticles for different shear rates ( γ ˙ ).
Materials 14 03301 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, X.; Nie, Y.; Li, H.; Wu, H.; Gao, Y.; Zhang, L. Percolated Network of Mixed Nanoparticles with Different Sizes in Polymer Nanocomposites: A Coarse-Grained Molecular Dynamics Simulation. Materials 2021, 14, 3301. https://doi.org/10.3390/ma14123301

AMA Style

Zhao X, Nie Y, Li H, Wu H, Gao Y, Zhang L. Percolated Network of Mixed Nanoparticles with Different Sizes in Polymer Nanocomposites: A Coarse-Grained Molecular Dynamics Simulation. Materials. 2021; 14(12):3301. https://doi.org/10.3390/ma14123301

Chicago/Turabian Style

Zhao, Xiuying, Yun Nie, Haoxiang Li, Haoyu Wu, Yangyang Gao, and Liqun Zhang. 2021. "Percolated Network of Mixed Nanoparticles with Different Sizes in Polymer Nanocomposites: A Coarse-Grained Molecular Dynamics Simulation" Materials 14, no. 12: 3301. https://doi.org/10.3390/ma14123301

APA Style

Zhao, X., Nie, Y., Li, H., Wu, H., Gao, Y., & Zhang, L. (2021). Percolated Network of Mixed Nanoparticles with Different Sizes in Polymer Nanocomposites: A Coarse-Grained Molecular Dynamics Simulation. Materials, 14(12), 3301. https://doi.org/10.3390/ma14123301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop