Coarse‐grained A‐graft‐B model of poly(lactic acid) for molecular dynamics simulations

A mesoscopic model of poly(lactic acid) is developed where the polymer is represented as an A‐graft‐B chain with monomer units consisting of two covalently connected beads. A coarse‐graining algorithm is proposed to convert an atomistic model of PLA into a coarse‐grained one. The developed model is based on atomistic simulations of oligolactides to take into account terminal groups correctly. It was used for coarse‐grained simulations of polylactide. Gyration radii and end to end distances of polymer chains as well as the density of the polymer melt are shown to be in a good agreement with those obtained from atomistic simulations. The thermal expansion coefficients of the OLA melts calculated using the coarse‐grained model are in reasonable agreement with those obtained from all‐atom molecular dynamics. The model provides a 17‐fold speedup compared with atomistic calculations. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 604–612     

Electrohydrodynamics of spherical polyampholyte‐grafted nanoparticles: Multiscale simulations by coupling of molecular dynamics and lattice‐boltzmann method

We perform multiscale simulations based on the coupling of molecular dynamics and lattice‐Boltzmann (LB) method to study the electrohydrodynamics of a polyampholyte‐grafted spherical nanoparticle. The long‐range hydrodynamic interactions are modeled by coupling the movement of particles to a LB fluid. Our results indicate that the net‐neutral soft particle moves with a nonzero mobility under applied electric fields. We systematically explore the effects of different parameters, including the chain length, grafting density, electric field, and charge sequence, on the structures of the polymer layer and the electrophoretic mobility of the soft particle. It shows that the mobility of nanoparticles has remarkable dependence on these parameters. We find that the deformation of the polyampholyte chains and the ion distribution play dominant roles in modulating the electrokinetic behavior of the polyampholyte‐grafted particle. The enhancement or reduction in the accumulation of counterions around monomers can be attributed to the polymer layer structure and the conformational transition of the chains in the electric field. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55, 1435–1447     

Hybrid particle‐field molecular dynamics simulations: Parallelization and benchmarks

The parallel implementation of a recently developed hybrid scheme for molecular dynamics (MD) simulations (Milano and Kawakatsu, J Chem Phys 2009, 130, 214106) where self‐consistent field theory (SCF) and particle models are combined is described. Because of the peculiar formulation of the hybrid method, considering single particles interacting with density fields, the most computationally expensive part of the hybrid particle‐field MD simulation can be efficiently parallelized using a straightforward particle decomposition algorithm. Benchmarks of simulations, including comparisons of serial MD and MD‐SCF program profiles, serial MD‐SCF and parallel MD‐SCF program profiles, and parallel benchmarks compared with efficient MD program GROMACS 4.5.4 are tested and reported. The results of benchmarks indicate that the proposed parallelization scheme is very efficient and opens the way to molecular simulations of large scale systems with reasonable computational costs. © 2012 Wiley Periodicals, Inc.     

Molecular dynamics simulations of cation‐controlled aggregation in fluorene–triarylamine copolymers

Systematic molecular dynamics (MD) simulations are applied, combining the sulfonated fluorene–triarylamine copolymer anion poly((N‐(4‐butylphenyl)‐diphenylamine)‐alt‐fluorene 9,9'di‐n‐propane sulfonate) with a wide range of monovalent cations of different cation field strengths. The resulting MD trajectories are analyzed for each of the conjugated polyelectrolyte focusing on ion aggregation, and its influence on the static and dynamic polymer morphology as well as on the charge carrier mobilities that are relevant for the use of such hole transporting materials in organic solar cells and organic light emitting devices. A pronounced variation of the degree of cation clustering with the cation field strength is found to control the polymer morphology, cation mobility and thereby the time evolution of the Coulomb energy landscape for hole transport. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 965–974     

Molecular dynamics simulation of poly(3‐hexylthiophene) helical structure In Vacuo and in amorphous polymer surrounding

The stability of poly(3‐hexylthiophene) (P3HT) helical structure has been investigated in vacuo and in amorphous polymer surrounding via molecular dynamics‐based simulations at temperatures below and above the P3HT melting point. The results show that the helical chain remains stable at room temperature both in vacuo and in amorphous surrounding, and promptly loses its structure at elevated temperatures. However, the amorphous surrounding inhibits the destruction of the helix at higher temperatures. In addition, it is shown that the electrostatic interactions do not significantly affect the stability of the helical structure. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 2448–2456     

Electric‐field‐induced molecular alignment of side‐chain liquid‐crystalline polyacetylenes containing biphenyl mesogens

Electric‐field‐induced molecular alignments of side‐chain liquid‐crystalline polyacetylenes [? {HC?C[(CH₂)_mOCO‐biph‐OC₇H₁₅]}? , where biph is 4,4′‐biphenylyl and m is 3 (PA3EO7) or 9 (PA9EO7)] were studied with X‐ray diffraction and polarized optical microscopy. An orientation as high as 0.84 was obtained for PA9EO7. Furthermore, the molecular orientation of PA9EO7 was achieved within a temperature range between the isotropic‐to‐smectic A transition temperature and 115 °C, and this suggested that the orientational packing was affected by the thermal fluctuation of the isotropic liquid and the mobility of the mesogenic moieties. The maximum achievable orientation for PA9EO7 was much greater than that for PA3EO7. This was the first time that the electric‐field‐induced molecular orientation of a side‐chain liquid‐crystalline polymer with a stiff backbone was studied. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1333–1341, 2004     

Synthesis,solid‐state,and charge‐transport properties of conjugated polythiophene‐S,S‐dioxides

An alkylated semiconducting polymer comprising alternating bithiophene‐[all]‐S,S‐dioxide and aromatic monothiophene units in the polymer backbone was synthesized with the intent of modifying the energy gap and lowest unoccupied molecular orbital for use as a stable n‐type semiconductor. Films spun from this semiconducting polymer were characterized utilizing X‐ray scattering, near edge X‐ray absorption fine structure spectroscopy, ultraviolet photoelectron spectroscopy, and thin‐film field effect transistors to determine how oxidation of the thiophene ring systems impacts the structural and electronic properties of the polymer. The thiophene‐S,S‐dioxide polymers have lower optical and electrical band gaps than corresponding thiophene polymers. X‐ray scattering results indicate that the polymers are well ordered with the π–π stacking distances increased by 0.4 Å relative to analogous thiophene polymers. The electrical stability of these polymers is poor in transistors with a drop in the field effect mobility by approximately one order of magnitude upon addition of just 5% of the thiophene‐S,S‐dioxide unit in a copolymer with thiophene. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2013     

Short‐time elasticity of polymer melts: Tobolsky conjecture and heterogeneous local stiffness

An extended molecular‐dynamics study of the short‐time “glassy” elasticity exhibited by a polymer melt of linear fully‐flexible chains above the glass transition is presented. The focus is on the infinite‐frequency shear modulus G_∞ manifested in the picosecond time scale and the relaxed plateau G_p reached at later times and terminated by the structural relaxation. The local stiffness of the interactions with the first neighbors of each monomer exhibits marked distribution with average value given by G_∞. In particular, the neighborhood of the end monomers of each chain are softer than the one of the inner monomers, so that G_∞ increases with the chain length. G_p is not affected by the chain length and is largely set by the nonbonding interactions, thus confirming for polymer melts the conjecture formulated by Tobolsky for glassy polymers. © 2015 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 1401–1407     

Molecular dynamics simulations of PNIPAM‐co‐PEGMA copolymer hydrophilic to hydrophobic transition in NaCl solution

Classical molecular dynamics simulations were carried out to investigate the hydrophilic to hydrophobic transition of PNIPAM‐co‐PEGMA close to its lower critical solution temperature (LCST) in 1 M NaCl solution. PNIPAM‐co‐PEGMA is a copolymer of poly(N‐isopropylacrylamide) (PNIPAM) and poly(ethylene glycol) methacrylate (PEGMA). The copolymer consists of 38 monomer units of NIPAM with two PEGMA chains attached to the PNIAPM backbone. The PNIPAM‐co‐PEGMA was observed to go through the hydrophilic?hydrophobic conformational change for simulations at temperature slightly above its LCST. Na⁺ ions were found to bind strongly and directly with amide O, even more strongly with the O atoms on PEGAMS chains, whereas Cl^? ions only exhibit weak interaction with the polymer. Significantly a novel caged stable metal‐organic complex involving a Na⁺ ion coordinated by six O atoms from the copolymer was observed after the PNIPAM‐co‐PEGMA copolymer went through conformational transition to form a hydrophobic folded structure. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2011     

Synthesis of comb‐like polystyrene with poly(N‐phenyl maleimide‐alt‐p‐chloromethyl styrene) as macroinitiator

The copolymerization of N‐phenyl maleimide and p‐chloromethyl styrene via reversible addition–fragmentation chain transfer (RAFT) process with AIBN as initiator and 2‐(ethoxycarbonyl)prop‐2‐yl dithiobenzoate as RAFT agent produced copolymers with alternating structure, controlled molecular weights, and narrow molecular weight distributions. Using poly(N‐phenyl maleimide‐alt‐p‐chloromethyl styrene) as the macroinitiator for atom transfer radical polymerization of styrene in the presence of CuCl/2,2′‐bipyridine, well‐defined comb‐like polymers with one graft chain for every two monomer units of backbone polymer were obtained. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 2069–2075, 2006     

Facile one‐pot/one‐step technique for preparation of side‐chain functionalized polymers: Combination of SET‐RAFT polymerization of azide vinyl monomer and click chemistry

An azido‐containing functional monomer, 11‐azido‐undecanoyl methacrylate, was successfully polymerized via ambient temperature single electron transfer initiation and propagation through the reversible addition–fragmentation chain transfer (SET‐RAFT) method. The polymerization behavior possessed the characteristics of “living”/controlled radical polymerization. The kinetic plot was first order, and the molecular weight of the polymer increased linearly with the monomer conversion while keeping the relatively narrow molecular weight distribution (M_w/M_n ≤ 1.22). The complete retention of azido group of the resulting polymer was confirmed by ¹H NMR and FTIR analysis. Retention of chain functionality was confirmed by chain extension with methyl methacrylate to yield a diblock copolymer. Furthermore, the side‐chain functionalized polymer could be prepared by one‐pot/one‐step technique, which is combination of SET‐RAFT and “click chemistry” methods. © 2011 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2012     

Equilibrium statistics of weakly slip‐linked Gaussian polymer chains

We calculate the free energy and the pressure of a weakly slip‐linked Gaussian polymer chains. We show that the equilibrium statistics of a slip‐linked system is different from one of the corresponding ideal chain system without any constraints by slip‐links. It is shown that the pressure of a slip‐linked system decreases compared with the ideal system, which implies that slip‐linked chains spontaneously form aggregated cluster‐like compact structures. These are qualitatively consistent with previous theoretical analyses or multichain simulations. We also show that repulsive potentials between chains, which have been phenomenologically utilized in simulations, can cancel the artificial pressure decrease. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2011     

Anisotropic self‐assembly and gelation in aqueous methylcellulose—theory and modeling

Recent experimental studies demonstrated that the aqueous methylcellulose (MC) polymer chains in water can form nanoscale fibrils (diameter ～14 nm, persistence length ～60 nm), and those fibrils can organize into networks at higher temperatures and/or concentrations, forming the commonly observed gel. Here we propose that the fibrils are one‐dimensional self‐assemblies of stacked, fused polymer rings that are formed at elevated temperatures due to the changing nature of the MC‐water hydrogen bonding. This mechanism is analogous to the coil‐helix transition in polypeptides, although it is not clear whether the MC fibrils possess chirality. We perform coarse‐grained molecular simulations of MC chain structure at temperatures both above and below the hypothesized coil‐to‐ring transition, with CG forcefield tuned by atomistic molecular dynamics simulations, and observe the expected conformational change. We then develop a statistical mechanical theory to predict the fibril self‐assembly, gelation and rheology as function of temperature and concentration. The findings are in reasonable agreement with experimental data and could be generalized to other carbohydrate polymers. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 1624–1636     

Anionic polymerization of N,N‐dialkyl‐2‐methylenebut‐3‐enamide: Effect of alkyl substituents on polymerization behavior

Anionic polymerizations of three 1,3‐butadiene derivatives containing different N,N‐dialkyl amide functions, N,N‐diisopropylamide (DiPA), piperidineamide (PiA), and cis‐2,6‐dimethylpiperidineamide (DMPA) were performed under various conditions, and their polymerization behavior was compared with that of N,N‐diethylamide analogue (DEA), which was previously reported. When polymerization of DiPA was performed at ?78 °C with potassium counter ion, only trace amounts of oligomers were formed, whereas polymers with a narrow molecular weight distribution were obtained in moderate yield when DiPA was polymerized at 0 °C in the presence of LiCl. Decrease in molecular weight and broadening of molecular weight distribution were observed when polymerization was performed at a higher temperature of 20 °C, presumably because of the effect of ceiling temperature. In the case of DMPA, no polymer was formed at 0 °C and polymers with relatively broad molecular weight distributions (M_w/M_n = 1.2) were obtained at 20 °C. The polymerization rate of PiA was much faster than that of the other monomers, and poly(PiA) was obtained in high yield even at ?78 °C in 24 h. The microstructure of the resulting polymers were exclusively 1,4‐ for poly(DMPA), whereas 20–30% of the 1,2‐structure was contained in poly(DiPA) and poly(PiA). © 2010 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 48: 3714–3721, 2010     

Influence of molecular weight on the short‐channel effect in polymer‐based field‐effect transistors

In this study, we demonstrate how the intrinsic properties of a polymer can influence the electrical characteristics of organic field‐effect transistors (OFETs). OFETs fabricated with three batches of poly[2‐methoxy,5‐(3′,7′‐dimethyl‐octyloxy)]‐p‐phenylene vinylene (MDMO‐PPV) were investigated. The properties of the polymers were initially investigated using Fourier transform infrared spectroscopy (FTIR), impedance spectroscopy (IS), gel permeation chromotography (GPC), and cyclic voltammetry (CV), respectively. The structure and purity of the polymer batches were found to be very comparable, but the molecular weight (M_n and M_w) and polydispersity (PDI = M_w/M_n), varied between the samples and the HOMO and LUMO levels of the polymers were found to depend on the molecular weight properties. OFETs were then fabricated with the polymers and electrically characterized. It was observed that the channel current and the field‐effect mobility increase with increasing polymer molecular weight. The output characteristics of the transistors, on the other hand, were found to depend on the PDI of the polymer. Saturation of the channel current occurs at higher source–drain voltages and short‐channel behavior was observed to start at longer channel lengths for polymers with a higher PDI. This behavior is observed to be thickness dependent, and the short‐channel behavior was more pronounced for thicker MDMO‐PPV films. These results are explained in terms of influences of chain packing and ordering and high bulk currents on the FET output and transistor parameters. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 50: 117–124, 2012     

Enhancing the electroluminescent efficiency of poly{5‐methoxy‐2‐[(2′‐ethyl‐hexyl)‐oxy]‐p‐phenylenevinylene}/poly(2,3‐diphenyl‐5‐octyl‐p‐phenylenevinylene) polyblends by the introduction of chemical bonding through a thermal treatment

A significant improvement in the electroluminescence (EL) properties was observed for a poly{5‐methoxy‐2‐[(2′‐ethyl‐hexyl)‐oxy]‐p‐phenylenevinylene} (MEH–PPV)/poly(2,3‐diphenyl‐5‐octyl‐p‐phenylenevinylene) (DPO–PPV) blend after a thermal treatment at 200 °C for 2 h in vacuo to furnish the chemical bonding between polymer chains. ¹H NMR spectroscopy and two‐photon excitation microscopy revealed that the chemical bonding turned the immiscible polyblend into a system more like a block copolymer with a vertically segregated morphology. Because both the lowest unoccupied molecular orbital and highest occupied molecular orbital levels of MEH–PPV in the wetting layer were higher than those of DPO–PPV in the upper layer, the heterojunction between the two layers of the polymers fit the category of so‐called type II heterojunctions. As a result, the turn‐on voltage of the polymer light‐emitting diode prepared with the thermally treated polyblend decreased to ～0.6 V, and the EL emission intensities and quantum efficiencies increased to about 4 times those of the untreated polyblend. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 62–69, 2006     

Dependence of the luminescent properties and chain length of poly[2‐methoxy‐5‐(2′‐ethyl‐hexyloxy)‐1,4‐phenylene vinylene] on the formation of cis‐vinylene bonds during Gilch polymerization

The presence of cis‐vinylene bonds in Gilch‐polymerized poly[2‐methoxy‐5‐(2′‐ethyl‐hexyloxy)‐1,4‐phenylene vinylene] is reported. Through fractionation, species with a weight‐average molecular weight of less than 37,000 exhibited an abnormal blueshift of photoluminescence spectra in toluene solutions, and this was attributed to the presence of cis‐vinylene bonds, as verified by NMR spectroscopy. Surprisingly, the fractionated species (～1 wt %) with a weight‐average molecular weight of 5000 were mostly linked by the cis‐vinylene bonds. The concentration decreased with the molecular weight until a molecular weight of 37,000 was reached; at that point, the polymer chains contained mainly trans‐vinylene bonds. Obviously, the formation of cis‐vinylene bonds strongly inhibited the growth of polymer chains during Gilch polymerization. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 43: 2520–2526, 2005     

Synthesis of methyl‐substituted azobenzene–carbazole conjugated copolymers with photoinduced structural changes

Azobenzene switches its structure instantaneously by reversible trans‐to‐cis and cis‐to‐trans photoisomerization with light irradiations. Dynamic change in polymer structure is expected via introducing an azobenzene unit into the main chain. In this study, a set of methyl‐substituted azobenzene–carbazole conjugated copolymers is synthesized by the Suzuki–Miyaura coupling method. Introduction of methyl substituents to the azobenzene unit of the monomer, and polymerization in a high‐boiling solvent improve the molecular weight of the polymer. Decrease of effective conjugation length due to the twisted structure of the main chain allows progress of photoisomerization. The microstructure of the polymer was determined with grazing incidence X‐ray diffraction (GIXD) measurements using synchrotron radiation. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2019, 57, 1756–1764     

End‐group fidelity in nitroxide‐mediated living free‐radical polymerizations

In this study, new nitroxides based on the 2,2,5‐trimethyl‐4‐phenyl‐3‐azahexane‐3‐oxy skeleton were used to examine chain‐end control during the preparation of polystyrene and poly(t‐butyl acrylate) under living free‐radical conditions. Alkoxyamine‐based initiators with a chromophore attached to either the initiating fragment or the mediating nitroxide fragment were prepared, and the extent of the incorporation of the chromophores at either the initiating end or the propagating chain end was determined. In contrast to 2,2,6,6‐tetramethyl piperidinoxy (TEMPO), the incorporation of the initiating and terminating fragment into the polymer chain was extremely high. For both poly(t‐butyl acrylate) and polystyrene with molecular weights less than or equal to 70,000, incorporations at the initiating end of greater than 97% were observed. At the terminating chain end, incorporations of greater than 95% were obtained for molecular weights less than or equal to 50,000. The level of incorporation tended to decrease slightly at higher molecular weights because of the loss of the alkoxyamine propagating unit, which had important consequences for block copolymer formation. These results clearly show that these new α‐H nitroxides could control the polymerization of vinyl monomers such as styrene and t‐butyl acrylate to an extremely high degree, comparable to anionic and atom transfer radical polymerization procedures. © 2000 John Wiley & Sons, Inc. J Polym Sci A: Polym Chem 38: 4749–4763, 2000