Development and Application of a Coarse-Grained Model for PNIPAM by Iterative Boltzmann Inversion and Its Combination with Lattice Boltzmann Hydrodynamics.

The polymer poly(N-isopropylacrylamide) (PNIPAM) is studied using a novel combination of multiscale modeling methodologies. We develop an iterative Boltzmann inversion potential of concentrated PNIPAM solutions and combine it with lattice Boltzmann as a Navier-Stokes equation solver for the solvent. We study in detail the influence of the methodology on statics and dynamics of the system. The combination is successful and significantly simpler and faster than other mapping techniques for polymer solution while keeping the correct hydrodynamics. The model can semiquantitatively describe the correct phase behavior and polymer dynamics.

Nonergodicity parameters of confined hard-sphere glasses.

Within a recently developed mode-coupling theory for fluids confined to a slit we elaborate numerical results for the long-time limits of suitably generalized intermediate scattering functions. The theory requires as input the density profile perpendicular to the plates, which we obtain from density functional theory within the fundamental-measure framework, as well as symmetry-adapted static structure factors, which can be calculated relying on the inhomogeneous Percus-Yevick closure. Our calculations for the nonergodicity parameters for both the collective as well as for the self motion are in qualitative agreement with our extensive event-driven molecular dynamics simulations for the intermediate scattering functions for slightly polydisperse hard-sphere systems at high packing fraction. We show that the variation of the nonergodicity parameters as a function of the wavenumber correlates with the in-plane static structure factors, while subtle effects become apparent in the structure factors and relaxation times of higher mode indices. A criterion to predict the multiple reentrant from the variation of the in-plane static structure is presented.

Round Robin Study: Molecular Simulation of Thermodynamic Properties from Models with Internal Degrees of Freedom.

Thermodynamic properties are often modeled by classical force fields which describe the interactions on the atomistic scale. Molecular simulations are used for retrieving thermodynamic data from such models, and many simulation techniques and computer codes are available for that purpose. In the present round robin study, the following fundamental question is addressed: Will different user groups working with different simulation codes obtain coinciding results within the statistical uncertainty of their data? A set of 24 simple simulation tasks is defined and solved by five user groups working with eight molecular simulation codes: DL_POLY, GROMACS, IMC, LAMMPS, ms2, NAMD, Tinker, and TOWHEE. Each task consists of the definition of (1) a pure fluid that is described by a force field and (2) the conditions under which that property is to be determined. The fluids are four simple alkanes: ethane, propane, n-butane, and iso-butane. All force fields consider internal degrees of freedom: OPLS, TraPPE, and a modified OPLS version with bond stretching vibrations. Density and potential energy are determined as a function of temperature and pressure on a grid which is specified such that all states are liquid. The user groups worked independently and reported their results to a central instance. The full set of results was disclosed to all user groups only at the end of the study. During the study, the central instance gave only qualitative feedback. The results reveal the challenges of carrying out molecular simulations. Several iterations were needed to eliminate gross errors. For most simulation tasks, the remaining deviations between the results of the different groups are acceptable from a practical standpoint, but they are often outside of the statistical errors of the individual simulation data. However, there are also cases where the deviations are unacceptable. This study highlights similarities between computer experiments and laboratory experiments, which are both subject not only to statistical error but also to systematic error.

Direct Phase Equilibrium Simulations of NIPAM Oligomers in Water.

NIPAM (N-isopropylacrylamide)-based polymers in water show many interesting properties in experiments, including a lower critical solution temperature (LCST) at 305 K and a conformational transition of single chains at the same temperature. The results of many simulation studies suggest that standard force fields are able to describe the conformational transition and the phase equilibrium well. We show by performing long molecular dynamics simulations of the direct liquid-liquid phase equilibrium of NIPAM trimers in water that there is no LCST in the expected temperature range for any of the force fields under study. The results show further that the relaxation times of single-chain simulations are considerably longer than anticipated. Conformational transitions of single polymers can therefore not necessarily be used as surrogates for a real phase transition.

Glass transition in confined geometry.

Extending mode-coupling theory, we elaborate a microscopic theory for the glass transition of liquids confined between two parallel flat hard walls. The theory contains the standard mode-coupling theory equations in bulk and in two dimensions as limiting cases and requires as input solely the equilibrium density profile and the structure factors of the fluid in confinement. We evaluate the phase diagram as a function of the distance of the plates for the case of a hard sphere fluid and obtain an oscillatory behavior of the glass transition line as a result of the structural changes related to layering.