Heterogeneities in Polymer Structural and Dynamic Properties in Graphene and Graphene Oxide Nanocomposites: Molecular Dynamics Simulations

PolySMart: a general coarse-grained molecular dynamics polymerization scheme

The development of simulation methods to study the structure and dynamics of a macroscopically sized piece of polymer material is important as such methods can elucidate structure-property relationships. Several methods have been reported to construct initial structures for homo- and co-polymers; however, most of them are only useful for short linear polymers since one needs to pack and equilibrate the far-from-equilibrium initial structures, which is a tedious task for long or hyperbranched polymers and unfeasible for polymer networks. In this method article, we present PolySMart, i.e., an open-source python package, which can effectively produce fully equilibrated homo- and hetero-polymer melts and solutions with no limitation on the polymer topology and size, at a coarse-grained resolution and through a bottom-up approach. This python package is also capable of exploring the polymerization kinetics through its reactive scheme in realistic conditions so that it can model the multiple co-occurring polymerization reactions (with different reaction rates) as well as consecutive polymerizations under stoichiometric and non-stoichiometric conditions. Thus, the equilibrated polymer models are generated through correct polymerization kinetics. A benchmark and verification of the performance of the program for several realistic cases, i.e., for homo-polymers, co-polymers, and crosslinked networks, is given. We further discuss the capability of the program to contribute to the discovery and design of new polymer materials.

Effect of Electrostatic Interactions on the Interfacial Energy between Thermoplastic Polymers and Graphene Oxide: A Molecular Dynamics Study

In this study, the atomistic-scale mechanisms affecting the interfacial stability of a thermoplastic polymer/graphene oxide interface are investigated using molecular dynamics simulations. Different combinations of thermoplastic polymers (polyethersulfone (PES) and polyetherimide (PEI)) and graphene oxides modified with –O–, –OH, and –COOH are prepared. PES is found to be more strongly stabilized with modified/functionalized graphene oxide in the order of –COOH, –OH, –O–, which is opposite to the stability order of PEI. Our results suggest that these orders of stability are governed by a balance between the following two factors resulting from electrostatic interactions: (1) atoms with a strong charge bias attract each other, thereby stabilizing the interface; (2) the excluded-volume effect of the functional groups on graphene oxide destabilizes the interface by preventing π-π stacking of aromatic rings.

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.

Molecular dynamics simulations of hyperbranched poly(ethylene imine)-graphene oxide nanocomposites as dye adsorbents for water purification

Atomistically detailed molecular dynamics simulations were employed to study the adsorption capacity of graphene-oxide-based (GO) aqueous systems for the methylene blue (MB) dye in the presence of branched poly(ethylene imine) (BPEI) polymers. The polymeric component was either freely mixed or chemically attached to GO. The main focus was the elucidation of the effects originating from the presence of BPEI molecules in the association of MB with the formed GO complexes. The effect of temperature was also examined. It was found that the presence of the cationic BPEI molecules results in the formation of a distinct microenvironment characterized by a polymer-mediated interconnected morphology which promotes the development of larger-sized MB clusters. These clusters were found to form in the vicinity of the GO flakes, increasing thus the adsorption capacity of the dye molecules in the polymer-containing systems. Particularly in the system with the BPEI-functionalized GO flakes, a persistent percolated structure is formed, which results in a more restricted diffusion of the MB molecules, increasing thus significantly their residence time close to the GO surface. The clustering behavior of MB was found to be temperature-dependent in the BPEI-based models, providing useful information regarding the conditions for optimal adsorption performance of such membranes, in nanofiltration processes.

Interaction between Graphene-Based Materials and Small Ag, Cu, and CuO Clusters: A Molecular Dynamics Study

The excessive use of antibiotics has contributed to the rise in antibiotic-resistant bacteria, and thus, new antibacterial compounds must be developed. Composite materials based on graphene and its derivatives doped with metallic and metallic oxide nanoparticles, particularly Ag, Cu, and Cu oxides, hold great promise. These materials are often modified with polyethylene glycol (PEG) to improve their pharmacokinetic behavior and their solubility in biological media. In this work, we performed molecular dynamics (MD) simulations to study the interaction between small Ag, Cu, and CuO clusters and several graphene-based materials. These materials include pristine graphene (PG) and pristine graphene nanoplatelets (PGN) as well as PEGylated graphene oxide (GO_PEG) and PEGylated graphene oxide nanoplatelets (GO-PEG_N). We calculated the adsorption energies, mean equilibrium distances between the nanoparticles and graphene surfaces, and mean square displacement (MSD) of the nanoclusters. The results show that PEGylation favors the adsorption of the clusters on the graphene surfaces, causing an increase in adsorption energies and a decrease in both distances and MSD values. The strengthening of the interaction could be crucial to obtain effective antibacterial compounds.

Hydrophilic dangling chain interfacial segregation in polyurethane networks at aqueous interfaces and its underlying mechanisms: molecular dynamics simulations

Polymer networks with hydrophilic dangling chains are ideal candidates for many submerged applications, e.g., protein non-adhesive coatings with non-fouling behavior. The dangling chains segregate from the polymer network towards the water and form a brush-like structure at the interface. Several factors such as the polymer network structure, dangling chain length, and water/dangling chain interaction may all affect the interfacial performance of the polymer. Therefore, we employed a Martini based coarse-grained (CG) molecular dynamics (MD) simulation to elucidate the influences of the abovementioned parameters on dangling chain interfacial segregation. We built up several polyurethane (PU) networks based on poly(tetra methylene glycol) (PTMG), as a macrodiol, and methoxy poly(ethylene glycol) (mPEG), as a dangling chain, with varying molecular weights. We found out that the macrodiol/dangling chain length ratio considerably smaller than one impedes the migration of dangling chains towards the water interface, while the dangling chain hydrophilicity and length determine the polymer interfacial layer density/thickness. Then, we artificially changed the dangling chain affinity to water from an intermediate to a very attractive water/dangling chain interaction. We justified that a brush-like structure forms in two consecutive steps: first, a longitudinal, and then a lateral migration of dangling chains in water. The latter step results in a uniform interfacial layer over the polymer interface that mainly occurs in the case of the attractive water/dangling chain interaction.

Atomistic Molecular Dynamics Simulations of the Initial Crystallization Stage in an SWCNT-Polyetherimide Nanocomposite

Crystallization of all-aromatic heterocyclic polymers typically results in an improvement of their thermo-mechanical properties. Nucleation agents may be used to promote crystallization, and it is well known that the incorporation of nanoparticles, and in particular carbon-based nanofillers, may induce or accelerate crystallization through nucleation. The present study addresses the structural properties of polyetherimide-based nanocomposites and the initial stages of polyetherimide crystallization as a result of single-walled carbon nanotube (SWCNT) incorporation. We selected two amorphous thermoplastic polyetherimides ODPA-P3 and aBPDA-P3 based on 3,3′,4,4′-oxydiphthalic dianhydride (ODPA), 2,3′,3,4′-biphenyltetracarboxylic dianhydride (aBPDA) and diamine 1,4-bis[4-(4-aminophenoxy)phenoxy]benzene (P3) and simulated the onset of crystallization in the presence of SWCNTs using atomistic molecular dynamics. For ODPA-P3, we found that the planar phthalimide and phenylene moieties show pronounced ordering near the CNT (carbon nanotube) surface, which can be regarded as the initial stage of crystallization. We will discuss two possible mechanisms for ODPA-P3 crystallization in the presence of SWCNTs: the spatial confinement caused by the CNTs and π–π interactions at the CNT-polymer matrix interface. Based on our simulation results, we propose that ODPA-P3 crystallization is most likely initiated by favorable π–π interactions between the carbon nanofiller surface and the planar ODPA-P3 phthalimide and phenylene moieties.