Coarse-grained kinetic Monte Carlo models: Complex lattices, multicomponent systems, and homogenization at the stochastic level

On-lattice kinetic Monte Carlo (KMC) simulations have extensively been applied to numerous systems. However, their applicability is severely limited to relatively short time and length scales. Recently, the coarse-grained MC (CGMC) method was introduced to greatly expand the reach of the lattice KMC technique. Herein, we extend the previous spatial CGMC methods to multicomponent species and/or site types. The underlying theory is derived and numerical examples are presented to demonstrate the method. Furthermore, we introduce the concept of homogenization at the stochastic level over all site types of a spatially coarse-grained cell. Homogenization provides a novel coarsening of the number of processes, an important aspect for complex problems plagued by the existence of numerous microscopic processes (combinatorial complexity). As expected, the homogenized CGMC method outperforms the traditional KMC method on computational cost while retaining good accuracy.

The role of molecular interactions and interfaces in diffusion: Permeation through single-crystal and polycrystalline microporous membranes

In this second paper of a two part series, we investigate the implications of the interfacial phenomenon, caused by adsorbate-adsorbate interactions coupled with the difference in adsorbate density between the zeolite and the gas phase, upon benzene permeation through single-crystal and polycrystalline microporous NaX membranes. The high flux predicted for thin single-crystal membranes reveals that substantially enhanced flux should be expected in submicron films. Simulations also indicate that the standard local equilibrium assumption made for larger scale membranes is inapplicable at the submicron scale associated with nanometer size grains of thin and/or polycrystalline membranes. Apparent activation energies predicted for benzene permeation through NaX membranes via kinetic Monte Carlo (KMC) simulations are in good agreement with laboratory experiments. The simulations also uncover temperature-dependent flux pathways leading to non-Arrhenius behavior observed experimentally. The failure of the Darken approximation, especially in the presence of the interfacial phenomenon, leads to a substantial overprediction of the flux. Simulations of polycrystalline membranes suggest that this same interfacial phenomenon leads to resistance that can reduce flux by an order of a magnitude with only moderate polycrystallinity.

Molecular sieve valves driven by adsorbate-adsorbate interactions: hysteresis in permeation of microporous membranes

A recently derived mesoscopic framework describing activated micropore diffusion is employed to explore system criticality in microporous membranes under nonequilibrium conditions. Rapid exploration of parameter space, possible with this continuum framework, elucidates a novel temperature-induced ignition and extinction of the molecular flux under a macroscopic gradient in pressure (chemical potential). Deviation from equilibrium like phase behavior (i.e., shifting and narrowing of phase envelopes and double hysteresis) derives from asymmetry of the coupled boundaries of the nonequilibrium membrane. We confirm this new phase behavior, akin to "opening" and "closing" of a molecular valve, via gradient kinetic Monte Carlo simulations of thin one-dimensional and three-dimensional systems. The heat of adsorption, strength of adsorbate-adsorbate intermolecular forces, and chemical potential gradient are all shown to control 'valve' actuation, suggesting potential implications in chemical sensing and novel diffusion control.