Optimizing the accuracy of lattice Monte Carlo algorithms for simulating diffusion.

Phys Rev E Stat Nonlin Soft Matter Phys

Department of Physics, University of Ottawa, 150 Louis-Pasteur, Ottawa, Ontario, Canada K1N 6N5.

Published: January 2012

The behavior of a lattice Monte Carlo (LMC) algorithm (if it is designed correctly) must approach that of the continuum system that it is designed to simulate as the time step and the mesh step tend to zero. However, we show for an algorithm for unbiased particle diffusion that if one of these two parameters remains fixed, the accuracy of the algorithm is optimal for a certain finite value of the other parameter. In one dimension, the optimal algorithm with moves to the two nearest neighbor sites reproduces the correct second and fourth moments (and minimizes the error for the higher moments at large times) of the particle distribution and preserves the first two moments of the first-passage time distributions. In two and three dimensions, the same level of accuracy requires simultaneous moves along two axes ("diagonal" moves). Such moves attempting to cross an impenetrable boundary should be projected along the boundary, rather than simply rejected. We also treat the case of absorbing boundaries. We discuss the relation between optimally accurate LMC algorithms and a particular case of lattice Boltzmann (LB) algorithms for simulating diffusion and compare the computational efficiency of optimal LMC and optimal LB algorithms.

Download full-text PDF

Source
http://dx.doi.org/10.1103/PhysRevE.85.016709DOI Listing

Publication Analysis

Top Keywords

lattice monte
8
monte carlo
8
algorithms simulating
8
simulating diffusion
8
optimizing accuracy
4
accuracy lattice
4
algorithms
4
carlo algorithms
4
diffusion behavior
4
behavior lattice
4

Similar Publications

A coarse-grained model of a two-dimensional colloidal suspension was designed. The model was athermal and, in addition, a lattice approximation was introduced. It consisted of solvent (monomer) molecules, dimer molecules, and immobile impenetrable obstacles that introduced additional heterogeneity into the system.

View Article and Find Full Text PDF

Magnetoplumbites are one of the most broadly studied families of hexagonal ferrites, typically with high magnetic ordering temperatures, making them excellent candidates for permanent magnets. However, magnetic frustration is rarely observed in magnetoplumbites. Herein, the discovery, synthesis, and characterization of the first Mn-based magnetoplumbite, as well as the first magnetoplumbite involving pnictogens (Sb), ASbMnO (A = K or Rb) are reported.

View Article and Find Full Text PDF

Background: A novel focal lattice-tip catheter allowing the delivery of either pulsed field (PF) or radiofrequency (RF) energy has recently received regulatory approval. The technology features a proprietary 3-dimensional electroanatomic mapping system.

Objective: We describe the first real-world and multicenter experience.

View Article and Find Full Text PDF

Statistics of random walks in geologically relevant conductance fields.

Phys Rev E

November 2024

Sonny Astani Department of Civil and Environmental Engineering, University of Southern California, Los Angeles, California 90089-2531, USA.

Statistics of diffusion, modeled by random walks, such as the mean number of distinct sites visited S(t) at time t, the mean probability P_{0}(t) of being at the origin of the walk, and the mean-squared displacements 〈R^{2}(t)〉 of the random walkers have been studied extensively in the past in both regular lattices and such disordered media as percolation clusters and other fractal structures, and universal power laws for such quantities have been derived. S(t) provides insight into reaction properties of geological formations, while P_{0}(t) is directly linked with the problem of back diffusion in remediation of groundwater aquifers. In all such studies, it was assumed that the conductances of the bonds that connect nearest-neighbor sites of the lattices are equal.

View Article and Find Full Text PDF

Polymorphs commonly exist for various materials, enabling phase engineering for diverse material properties. While the crystal structures of different polymorphs can, in principle, be experimentally characterized, interpreting and understanding complex crystal structures can be very challenging. Using Ga_{2}O_{3} as a prototype, here we show that the crystal structure of γ-Ga_{2}O_{3} has long been misinterpreted from either theory or experiment.

View Article and Find Full Text PDF

Want AI Summaries of new PubMed Abstracts delivered to your In-box?

Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!