Erik Nies

Dynamics of the crystal to plastic crystal transition in the hydrogen bonded N-isopropylpropionamide

N-Isopropylpropionamide (NiPPA), which can self-associate via hydrogen bonds, was found to undergo a solid-solid transition as identified by DSC and X-ray diffraction. Below the melting temperature of 51 °C NIPPA adopts a plastic crystalline state with a tetragonal unit cell until it transforms into an ordered crystal with a monoclinic structure at temperatures ≤10 °C. Dielectric spectroscopy was used to characterize the dynamics of the system, determining the activation parameters for the plastic to crystalline phase transition. The activation enthalpy is relatively high, as expected for a system that involves hydrogen bonds. However, most of the activation energy as the plastic phase assumes a more crystalline state is due to the activation entropy, suggesting that the increased cooperativity observed in the relaxation processes is due to a steric locking of the molecules.

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.

A parallel algorithm for step- and chain-growth polymerization in molecular dynamics.

Classical Molecular Dynamics (MD) simulations provide insight into the properties of many soft-matter systems. In some situations, it is interesting to model the creation of chemical bonds, a process that is not part of the MD framework. In this context, we propose a parallel algorithm for step- and chain-growth polymerization that is based on a generic reaction scheme, works at a given intrinsic rate and produces continuous trajectories. We present an implementation in the ESPResSo++ simulation software and compare it with the corresponding feature in LAMMPS. For chain growth, our results are compared to the existing simulation literature. For step growth, a rate equation is proposed for the evolution of the crosslinker population that compares well to the simulations for low crosslinker functionality or for short times.

Polymer networks by molecular dynamics simulation: Formation, thermal, structural and mechanical properties

A molecular dynamics simulation method is presented and used in the study of the formation of polymer networks. We study the formation of networks representing the methylene repeating units as united atoms. The network formation is accomplished by cross-linking polymer chains with dedicated functional end groups. The simulations reveal that during the cross-linking process, initially branched molecules are formed before the gel point; approaching the gel point, larger branched entities are formed through integration of smaller branched molecules, and at the gel point a network spanning the simulation box is obtained; beyond the gel point the network continues to grow through the addition of the remaining molecules of the sol phase onto the gel (the network); the final completion of the reaction occurs by intra-network connection of dangling ends onto unsaturated cross-linkers. The conformational properties of the strands in the undeformed network are found to be very similar with the conformational properties of the chains before cross-linking. The uniaxial deformation of the formed networks is investigated and the modulus determined from the stress-strain curves shows reciprocal scaling with the precursor chain length for networks formed from sufficiently large precursor chains (N ≥ 20).

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.

On the Thermal Expansion of Water and the Phase Behavior of Macromolecules in Aqueous Solution

Water is crucial for the existence of life as we know it, and many have wondered what makes water so special. Here we point out the analogies between the pressure-temperature dependence of the isobaric thermal expansion of water (αp) and the pressure-temperature phase behavior of macromolecules in aqueous solutions. We suggest that αp could be the key to understand why water is to be the so-called ‘matrix of life’.