6533b854fe1ef96bd12ae999

RESEARCH PRODUCT

Spinodal decomposition of polymer solutions: A parallelized molecular dynamics simulation

Kurt BinderLeonid YelashPeter VirnauWolfgang PaulMarcus Müller

subject

PhysicsMesoscopic physicsStructure formationThermodynamic equilibriumSpinodal decompositionMonte Carlo methodOrder (ring theory)02 engineering and technology021001 nanoscience & nanotechnology01 natural sciencesMolecular physicsMolecular dynamics0103 physical sciencesRelaxation (physics)Statistical physics010306 general physics0210 nano-technology

description

In simulations of phase separation kinetics, large length and time scales are involved due to the mesoscopic size of the polymer coils, and the structure formation on still larger scales of length and time. We apply a coarse-grained model of hexadecane dissolved in supercritical carbon dioxide, for which in previous work the equilibrium phase behavior has been established by Monte Carlo methods. Using parallelized simulations on a multiprocessor supercomputer, large scale molecular dynamics simulations of phase separation following pressure jumps are presented for systems containing $N=435\phantom{\rule{0.2em}{0ex}}136$ coarse-grained particles, which correspond to several millions of atoms in a box with linear dimension $447\phantom{\rule{0.3em}{0ex}}\mathrm{\AA{}}$. Even for large systems the phase separation can be observed up to the final, macroscopically segregated, equilibrium state. It is shown that in the segregation process the two order parameters of the system (density and concentration) are strongly coupled. The system does not follow the predicted growth law for the characteristic domain size $\mathcal{l}(t)\ensuremath{\propto}t$ in binary fluid mixtures for the range of times accessible in the simulation. Instead, it exhibits a distinctly slower growth, presumably due to the dynamic asymmetry of the constituents.

https://doi.org/10.1103/physreve.78.031801