Multi-scale modeling of polyimides

Reverse mapping method for complex polymer systems

We present a comprehensive approach for reverse mapping of complex polymer systems in which the connectivity is created by the simulation of chemical reactions at the coarse-grained scale. Within the work, we use a recently developed generic adaptive reverse mapping procedure that we adapt to handle the varying connectivity structure resulting from the chemical reactions. The method is independent of the coarse-grained and fine-grained force-fields by design and relies only on a single control parameter. We employ the approach to reverse map four different systems: a three-component epoxy network, a trimethylol melamine network, a hyperbranched polymer, and a polyethylene terephthalate. In the case of the epoxy network, we use two fine-scale representations: fully atomistic and united-atom. Whereas the trimethylol melamine network the hyperbranched polymer and the polyethylene terephthalate are reverse mapped to the fully atomistic description. After the reverse mapping, we examine the fine-grained structure by comparing the radial distribution functions with respect to the control parameter. Moreover, in the case of the epoxy we perform tensile-test experiments and examine the resulting Young's modulus. In all cases, we show how the properties of the reverse mapped systems depend on the control parameters. In general, we see that the results are relatively insensitive to the control parameter and the resulting atomistic systems are stable. Only for the trimethylol melamine network, we notice chemically incorrect conformations when the reverse mapping is performed too fast. We provide a remedy for this issue. © 2017 Wiley Periodicals, Inc.

Molecular Structuring and Percolation Transition in Hydrated Sulfonated Poly(ether ether ketone) Membranes

The extent of phase separation and water percolation in sulfonated membranes are the key to their performance in fuel cells. Toward this, the effect of hydration on the morphology and transport characteristics of sulfonated poly(ether ether ketone), sPEEK, membrane is investigated using atomistic molecular dynamics simulation at various hydration levels(λ: number of water molecules per sulfonate group). The evolution of local morphology is investigated using structural correlations and minimum pair distances. Transport properties are probed using mean squared displacements and diffusion coefficients. The water-sulfonate interaction in sPEEK is found to be stronger than that in Nafion, as observed in experiments. Analyses indicate the presence of narrow connected path of water and hydronium at λ = 4 and large domains, spanning half the simulation box, at λ = 15. The behavior of membrane water remains far from bulk as indicated by its diffusion coefficient. The persistence of small isolated water clusters demonstrates the extent of phase separation in sPEEK to be lesser than that in Nafion. Analyses at molecular and collective levels suggest the occurrence of a percolation transition between λ = 8 and 10, which leads to a connected network of water channels in the membrane, thereby boosting the hydronium mobility.

Size and chain length effects on structural behaviors of biphenylcyclohexane-based liquid crystal nanoclusters by a coarse-grained model

Size and chain length effects on structural behaviors of liquid crystal nanoclusters were examined by a coarse-grained model and the configurational-bias Monte Carlo (CBMC) simulation. The nanoclusters investigated in this study are composed of the biphenylcyclohexane-based BCH5H liquid crystal molecule and its derivatives. Results of the study show that the average energy decreases (i.e., more negative) as the cluster size (i.e., the number of molecules) increases. With the increasing cluster size, the equilibrium conformation of the nanocluster changes gradually from a pipe-like structure (for the smaller systems) to a ball-like cluster (for the larger systems). The order parameter of the system reduces with the transition of the equilibrium conformation. Regarding the chain length effect, the pipe-like equilibrium conformation (for the smaller systems) was observed more close to a pipe as the length of the tail alkyl chain of the derivatives extended. However, due to the flexibility of the tail alkyl chain, the pipe conformation of the system deflects slightly about its cyclohexyl group as the tail extends further.

Single-chain and aggregation properties of semiconducting polymer solutions investigated by coarse-grained langevin dynamics simulation

A coarse-grained (CG) model and Langevin dynamics scheme are proposed to investigate the material properties in dilute solution of a model semiconducting conjugated polymer, poly(2-methoxy-5-(2'-ethylhexyloxy)-1,4-phenylenevinylene) (MEH-PPV). While the intra- and intermolecular potentials for the CG particle (currently, a monomer unit) were determined from the molecular dynamics (MD) simulation of a united atomistic model, fluctuation-dissipation forces arising from the treatment of a solvent field were self-consistently constructed from the measured particle diffusivity in a given solvent (i.e., chloroform or toluene) through the atomistic MD simulation. It is shown that the resultant Langevin dynamics simulation, which is substantially more efficient than the counterpart MD simulation of the same CG model, is able to capture the dynamic (such as center-of-mass diffusivity) as well as the structural (such as radius of gyration) features of the investigated polymer solutions. Essential material properties that can now be directly studied include the following: Scaling exponents for estimating the exact solvent qualities were, for the first time, determined for the two solvent systems investigated; the persistence length obtained was also noted to be in excellent agreement with early experimental estimations. Preliminary observations on the supramolecular aggregation properties were in good agreement with the general observations from a wide range of recent experiments, and shed light on the essential impact of solvent quality on the supramolecular aggregation structures.

Developing a coarse-grained force field for the diblock copolymer poly(styrene-b-butadiene) from atomistic simulation

We have developed a coarse-grained force field for the poly(styrene-b-butadiene) diblock copolymer. We describe the computational methods and discuss how they were applied to develop a coarse-grained force field for this diblock copolymer from the atomistic simulation. The new force field contains three different bonds, four angles, five dihedral angles, and three nonbonded terms. We successfully tested this coarse-grained model against the chain properties, including static and dynamic properties, derived from the atomistic simulation; the results suggest that the coarse-grained force field is an effective model.