Xiaozhen Yang

H-bond and conformations of donors and acceptors in model polyether based polyurethanes

Abstract The molecular mechanics (MM) method with COMPASS force field was used to study the H-bonds in the polyether based polyurethane model molecules. Availability of the calculation was firstly verified in comparison of some H-bonded model molecules, which were studied by using ab initio calculation, and those calculated by MM. Based on a urethane model molecule 1,3-dimethylcarbamate, which can be donor or acceptor and behaves in various conformations, it is reasonable to have a large number of H-bond interactions between various conformational donors and acceptors. For examining all the possible interaction patterns, we studied 57H-bond complexes. This systematic modeling covers well-known four types of interaction patterns, such as NH⋯OC (Type I), NH⋯O–CO (Type II), NH⋯NH (Type III), NH⋯COC (Type IV) in the system. Obtained H-bond energies were used to analyze the probabilities of the complexes. For the interaction within the hard segments, or Type I, Type II and Type III, a predominant H-bond complex has been found in the present study, which belongs to Type I. For the interaction between the hard segment and the soft segment, two conformations of Type IV were calculated to be existed.

A Replica Exchange Molecular Dynamics Simulation of a Single Polyethylene Chain: Temperature Dependence of Structural Properties and Chain Conformational Study at the Equilibrium Melting Temperature.

The conformational properties of a finite length polyethylene chain were explored over a wide range of temperatures using a replica exchange molecular dynamics simulation providing high quality simulation data representative for the equilibrium behavior of the chain molecule. The radial distribution function (RDF) and the structure factor S(q) of the chain as a function of temperature are analyzed in detail. The different characteristic peaks in the RDF and S(q) were assigned to specific distances in the chain and structural changes occurring with the temperature. In S(q), a peak characteristic for the order in the solid state was found and used to determine the equilibrium melting temperature. A detailed scaling analysis of the structure factor covering the full q range was performed according to the work of Hammouda. In the Θ region, a quantitative analysis of the full structure factor was done using the equivalent Kuhn chain, which enabled us to assign the Θ region of our chain and to demonstrate, in our particular case, the failure of the Gaussian chain approach. The chain conformational properties at the equilibrium melting temperature are discussed using conformational distribution functions, using the largest principal component of the radius of gyration and shape parameters as order parameters. We demonstrate that for the system studied here, the Landau free energy expression based on this conformational distribution information leads to erroneous conclusions concerning the thermodynamic transition behavior. Finally, we focus on the instantaneous conformational properties at the equilibrium melting temperature and give a detailed analysis of the conformational shapes using different shape parameters and a simulation snapshot. We show that the chain does not only take the lamellar rod-like and globular conformational shapes, typical of the solid and liquid states, but can also explore many other conformational states, including the toroidal conformational state. It is the first demonstration that a flexible molecule like PE can also take a toroidal conformational state, which is normally linked to stiffer chains.

Nascent Crystallization of a Growing Chain on a Catalyst Surface : A Nonequilibrium Molecular Dynamics Simulation Study

The growing chain molecular dynamics (GCMD) simulation method, a new nonequilibrium molecular dynamics code, is proposed to simulate the polymer chain aggregation behavior during polymerization on a catalyst surface. We found that the growing chain crystallizes on the surface in two stages: the nucleation stage and the crystal growth stage. In the first part of the nucleation period, the short polymerizing chain first absorbs on the surface and can be in either an ordered or disordered structure. Still in the nucleation period, when the chain reaches a degree of polymerization, about 100 bonds, the chain folds into a stable nucleus on the substrate with 3-5 stems. In the crystal growth stage where the polymerization also proceeds, we observed a stem elongation process in combination with a chain folding process. In the stem elongation step, the number of stems in the nucleus remains constant, and all the stems expand together to a length of ca. 5-25 ns. In the subsequent chain folding step, the stem length decreases about 20 bonds within a period of ca. 0.1-0.5 ns. During chain growth, the elongation process and the folding process occur in an alternating and repeated fashion. The crystallization mechanism of the polymerizing chain was discussed.

Analysis of multiple H-bond interactions in self-assembly between polyurethane with pendent carboxyl and poly(4-vinylpyridine).

As an extension study, FTIR and molecular simulation methods were combined in the present paper to analyze the H-bond interactions resulting from multiple donors and acceptors that have led to self-assembly based on segmented polyurethane with carboxyl (PUc) and poly(4-vinylpyridine) (P4VP) in our previous work. Of them, FTIR was used to analyze the H-bonding types and interactions as well as their changes before and after self-assembly; molecular mechanics (MM/COMPASS) was used to study the effect of possible conformations on the H-bonds involved and analyze the most probable H-bond patterns; quantum mechanics (QM/B3LYP) was used to help confirm the experimental FTIR band assignments and calculate the H-bond energy. It was found that two types of H-bonds exist, namely, COOH...P4VP (type I) and (OCO)NH...P4VP (type II), based on OH and NH as the strong donors in the interaction between PUc and P4VP. Strong evidence has been obtained for a type II H-bond, which is the specialty in PUc/P4VP assembly. The type I and type II H-bonding energies are -11.293 and -7.150 kcal/mol, respectively. The forming probability of the type I H-bond accounts for 95.87%, while that of the type II H-bond is 4.13%, showing the primary driving force for the assembly based on PUc and P4VP is still the H-bond between COOH and P4VP, yet the H-bonds based on NH and pyridyl in P4VP cannot be ignored.