Bridging simulation length scales with cellular automata
John J. Karnes, Esteban D. Gadea, Shakkira Erimban, Ignacio J. Bombau, Valeria Molinero
Abstract
Multiscale simulation requires coupling physics models that operate at different characteristic length and time scales, because no single method spans the range needed for most problems of interest. Moving to a coarser-grained representation unlocks longer length and time scales, but it discards the microscopic interactions that build morphology. A fine-grained model can be initialized from an arbitrary packing and left to self-assemble into a physically meaningful structure; a lower-resolution model cannot, and must inherit its starting configuration from a higher-fidelity simulation. The length scales accessible to the coarse-grained model are therefore set not by the coarse-grained method itself, but by the largest fine-grained configuration that can be affordably equilibrated. A representative example is the scale-up from particle-based molecular dynamics (MD) to a lattice-based representation such as kinetic Monte Carlo (kMC). In this work, we present a cellular automata (CA) approach for generating arbitrarily large lattice starting configurations. CA is a natural fit for this task: short-ranged local rules drive the evolution of a lattice, and their repeated application gives rise to emergent long-range order, thematically mirroring how short-ranged interactions in MD produce self-assembled morphology. We use a configuration from a higher-fidelity simulation as a training set and learn the CA rules from it via logistic regression. As a proof of principle, we develop these rules for a hydrated anion exchange membrane (AEM), generate new starting configurations, and benchmark their performance in mesoscale kMC simulations against an MD-derived "ground truth." We then demonstrate the ability to generate substantially larger lattices and show that their behavior in kMC is consistent with that of the smaller CA benchmark configurations.
Create a lesson
Related papers
A Gaussian process coarse-grained potential for Na-montmorillonite
Yalda Pedram, Yaoting Zhang, Laurent Brochard et al.
First-principles theory of phonon renormalization from nonlinear electron-phonon interactions
Florian Kluibenschedl, Matthew Houtput, Jacques Tempere et al.
Spin-Lattice Dynamics and Interactions in Magnonic Spinels
Hari Paudyal, Yuri Suzuki, Michael E. Flatté et al.
Magnon-Phonon Dynamics in Multidimensional Antiferromagnetic Oxides
Yogendra Limbu, Michael E. Flatté, Durga Paudyal
Strain-Induced Metal-to-Insulator Transition in Antiferromagnetic SrCrO3 Thin Films
S. Jöhr, A. Carta, J. Moreno et al.
Tuning the Coercive Field in Ferroelectric Hf0.5Zr0.5O2-Al2O3 Heterostructures via Interfacial Charge Dynamics
Marshall B. Frye, Chanyoung Kim, Jeong-Woo Sun et al.