Phase behavior of polymer-nanorod composites: A comparative study using PRISM theory and molecular dynamics simulations

Well-dispersed composites of polymer and nanorods have many emerging applications and, therefore, are an important area of research. Polymer reference interaction site model (PRISM) theory and molecular dynamics simulations have become powerful tools in the study of the structure and phase behavior of polymer nanocomposites. In this work, we employ both PRISM theory and molecular dynamics simulations to determine the structure and spinodal phase diagram of 1% volume fraction of nanorods in a polymer melt. We make quantitative comparisons between the phase diagrams, which are reported as a function of nanorod aspect ratio and polymer-nanorod interactions. We find that both PRISM theory and molecular dynamics simulations predict the formation of contact aggregates at low polymer-nanorod attraction strength (gamma) and bridged aggregates at high polymer-nanorod attraction strength. They predict an entropic depletion-driven phase separation at low gamma and a bridging-driven spinodal phase separation at high gamma. The polymer and nanorods are found to form stable composites at intermediate values of the polymer-nanorod attraction strength. The fall of the bridging boundary and the gradual rise of the depletion boundary with the nanorod aspect ratio are predicted by both PRISM theory and molecular dynamics simulations. Hence, the miscible region narrows with increasing aspect ratio. The depletion boundaries predicted by theory and simulation are quite close. However, the respective bridging boundaries present a significant quantitative difference. Therefore, we find that theory and simulations qualitatively complement each other and display quantitative differences.