TY - JOUR

T1 - Coarse-grained computations for a micellar system

AU - Kopelevich, Dmitry I.

AU - Panagiotopoulos, Athanassios Z.

AU - Kevrekidis, Yannis

PY - 2005

Y1 - 2005

N2 - We establish, through coarse-grained computation, a connection between traditional, continuum numerical algorithms (initial value problems as well as fixed point algorithms), and atomistic simulations of the Larson model of micelle formation. The procedure hinges on the (expected) evolution of a few slow, coarse-grained mesoscopic observables of the Monte Carlo simulation, and on (computational) time scale separation between these and the remaining "slaved," fast variables. Short bursts of appropriately initialized atomistic simulation are used to estimate the (coarse grained, deterministic) local dynamics of the evolution of the observables. These estimates are then in turn used to accelerate the evolution to computational stationarity through traditional continuum algorithms (forward Euler integration, Newton-Raphson fixed point computation). This "equation-free" framework, bypassing the derivation of explicit, closed equations for the observables (e.g., equations of state), may provide a computational bridge between direct atomistic/stochastic simulation and the analysis of its macroscopic, system-level consequences.

AB - We establish, through coarse-grained computation, a connection between traditional, continuum numerical algorithms (initial value problems as well as fixed point algorithms), and atomistic simulations of the Larson model of micelle formation. The procedure hinges on the (expected) evolution of a few slow, coarse-grained mesoscopic observables of the Monte Carlo simulation, and on (computational) time scale separation between these and the remaining "slaved," fast variables. Short bursts of appropriately initialized atomistic simulation are used to estimate the (coarse grained, deterministic) local dynamics of the evolution of the observables. These estimates are then in turn used to accelerate the evolution to computational stationarity through traditional continuum algorithms (forward Euler integration, Newton-Raphson fixed point computation). This "equation-free" framework, bypassing the derivation of explicit, closed equations for the observables (e.g., equations of state), may provide a computational bridge between direct atomistic/stochastic simulation and the analysis of its macroscopic, system-level consequences.

UR - http://www.scopus.com/inward/record.url?scp=22944466441&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=22944466441&partnerID=8YFLogxK

U2 - 10.1063/1.1839173

DO - 10.1063/1.1839173

M3 - Article

C2 - 15740298

AN - SCOPUS:22944466441

SN - 0021-9606

VL - 122

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

IS - 4

M1 - 044907

ER -