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 -