Multiscale Kinetic Monte-Carlo for Simulating Epitaxial Growth

We present a fast Monte-Carlo algorithm for simulating epitaxial surface growth, based on the continuous-time Monte-Carlo algorithm of Bortz, Kalos and Lebowitz. When simulating realistic growth regimes, much computational time is consumed by the relatively fast dynamics of the adatoms. Continuum and continuum-discrete hybrid methods have been developed to approach this issue; however in many situations, the density of adatoms is too low to efficiently and accurately simulate as a continuum. To solve the problem of fast adatom dynamics, we allow adatoms to take larger steps, effectively reducing the number of transitions required. We achieve nearly a factor of ten speed up, for growth at moderate temperatures and large D/F.Comment: 7 pages, 6 figures; revised text, accepted by PR

An Alternative Method to Deduce Bubble Dynamics in Single Bubble Sonoluminescence Experiments

In this paper we present an experimental approach that allows to deduce the important dynamical parameters of single sonoluminescing bubbles (pressure amplitude, ambient radius, radius-time curve) The technique is based on a few previously confirmed theoretical assumptions and requires the knowledge of quantities such as the amplitude of the electric excitation and the phase of the flashes in the acoustic period. These quantities are easily measurable by a digital oscilloscope, avoiding the cost of expensive lasers, or ultrafast cameras of previous methods. We show the technique on a particular example and compare the results with conventional Mie scattering. We find that within the experimental uncertainties these two techniques provide similar results.Comment: 8 pages, 5 figures, submitted to Phys. Rev.

ADI splitting schemes for a fourth-order nonlinear partial differential equation from image processing

We present directional operator splitting schemes for the numerical solution of a fourth-order, nonlinear partial differential evolution equation which arises in image processing. This equation constitutes the H−1-gradient flow of the total variation and represents a prototype of higher-order equations of similar type which are popular in imaging for denoising, deblurring and inpainting problems. The efficient numerical solution of this equation is very challenging due to the stiffness of most numerical schemes. We show that the combination of directional splitting schemes with implicit time-stepping provides a stable and computationally cheap numerical realisation of the equation

by kinetic Monte Carlo

Computation of strained epitaxial growth in three dimension