We demonstrate how an iterative method for potential inversion from distribution functions developed for simple liquid systems can be generalized to polymer systems. It uses the differences in the potentials of mean force between the distribution functions generated from a guessed potential and the true distribution functions to improve the effective potential successively. The optimization algorithm is very powerful: convergence is reached for every trial function in few iterations. As an extensive test case we coarse-grained an atomistic all-atom model of polyisoprene (PI) using a 13:1 reduction of the degrees of freedom. This procedure was performed for PI solutions as well as for a PI melt. Comparisons of the obtained force fields are drawn. They prove that it is not possible to use a single force field for different concentration regimes.
We present a detailed study of a new, optimized coarse‐grained (CG) model of polystyrene (PS) and compare it with a recently published one (Harmandaris et al., Macromolecules 2006, 39, 6708). By implementing a different mapping scheme, the new model, augmented with softer nonbonded interactions, better reproduces the local chain conformations and melt packing observed in atomistic simulations of atactic PS. Both models properly predict the bonded distributions and are capable of simulating different tacticities without needing sidegroups. Both CG models fit dynamic data from long atomistic simulations after determining the scale factor for the simulation time. Together with a rigorous back‐mapping procedure from the mesoscopic to atomistic description, this opens up a very feasible way for generating very long atomistic trajectories.magnified image
We develop coarse-grained force fields for poly (vinyl alcohol) and poly (acrylic acid) oligomers. In both cases, one monomer is mapped onto a coarsegrained bead. The new force fields are designed to match structural properties such as radial distribution functions of various kinds derived by atomistic simulations of these polymers. The mapping is therefore constructed in a way to take into account as much atomistic information as possible. On the technical side, our approach consists of a simplex algorithm which is used to optimize automatically non-bonded parameters as well as bonded parameters. Besides their similar conformation (only the functional side group differs), poly (acrylic acid) was chosen to be in aqueous solution in contrast to a poly (vinyl alcohol) melt. For poly (vinyl alcohol) a non-optimized bond angle potential turns out to be sufficient in connection with a special, optimized non-bonded potential. No torsional potential has to be applied here. For poly (acrylic acid), we show that each peak of the radial distribution function is usually dominated by some specific model parameter(s). Optimization of the bond angle parameters is essential. The coarse-grained forcefield reproduces the radius of gyration R G of the atomistic model. As a first application, we use the force field to simulate longer chains and compare the hydrodynamic radius R H with experimental data.
The aim of this study is to understand deeper the thermal diffusion transport process (Ludwig-Soret effect) at the microscopic level. For that purpose, the recently developed reverse nonequilibrium molecular dynamics method was used to calculate Soret coefficients of various systems in a systematic fashion. We studied binary Lennard-Jones (LJ) fluids near the triple point (of one of the components) in which we separately changed the ratio of one of the LJ parameters mass, atomic diameter and interaction strength while keeping all other parameters fixed and identical. We observed that the magnitude of the Soret coefficient depends on all three ratios. Concerning its sign we found that heavier species, smaller species and species with higher interaction strengths tend to accumulate in the cold region whereas the other ones (lighter, bigger or weaker bound) migrate to the hot region of our simulation cell. Additionally, the superposition of the influence of the various parameters was investigated as well as more realistic mixtures. We found that in the experimentally relevant parameter range the contributions are nearly additive and that the mass ratio often is the dominating factor.
We implemented a coarse-graining procedure to construct mesoscopic models of complex molecules. The final aim is to obtain better results on properties depending on slow modes of the molecules. Therefore the number of particles considered in molecular dynamics simulations is reduced while conserving as many properties of the original substance as possible. We address the problem of finding nonbonded interaction parameters which reproduce structural properties from experiment or atomistic simulations. The approach consists of optimizing automatically nonbonded parameters using the simplex algorithm to fit structural properties like the radial distribution function as target functions. Moreover, any mix of structural and thermodynamic properties can be included in the target function. Different spherically symmetric inter-particle potentials are discussed. Besides demonstrating the method for Lennard-Jones liquids, it is applied to several more complex molecular liquids such as diphenyl carbonate, tetrahydrofurane, and monomers of poly(isoprene).
This work addresses the question of the scaling behavior of polyelectrolytes in solution for a realistic prototype. We show results of a combined experimental (light scattering) and theoretical (computer simulations) investigation of structural properties of the sodium salt of poly (acrylic acid) (Na-PAA). Experimentally, we determined the molecular weight (MW) and the hydrodynamic radius (RH) by static and dynamic light scattering for six different Na-PAA samples in aqueous NaCl-containing solution (0.1–1 mol/l) of polydispersity DP between 1.5 and 1.8. On the computational side, three different variants of a newly developed mesoscopic force field for Na-PAA were employed to determine RH for monodisperse systems of the same MW as in the experiments. The force field effectively incorporates atomistic information and one coarse-grained bead corresponds to one PAA monomer. We find that RH matches with the experimental data for all investigated samples. The effective scaling exponent for RH is found to be around 0.55, which is well below its asymptotic value for good solvents. Additionally, data for the radius of gyration (RG) are presented.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.