Low temperature shape relaxation of 2-d islands by edge diffusion

We present a precise microscopic description of the limiting step for low temperature shape relaxation of two dimensional islands in which activated diffusion of particles along the boundary is the only mechanism of transport allowed. In particular, we are able to explain why the system is driven irreversibly towards equilibrium. Based on this description, we present a scheme for calculating the duration of the limiting step at each stage of the relaxation process. Finally, we calculate numerically the total relaxation time as predicted by our results and compare it with simulations of the relaxation process.Comment: 11 pages, 5 figures, to appear in Phys. Rev.

Continuum description of profile scaling in nanostructure decay

The relaxation of axisymmetric crystal surfaces with a single facet below the roughening transition is studied via a continuum approach that accounts for step energy g_1 and step-step interaction energy g_3>0. For diffusion-limited kinetics, free-boundary and boundary-layer theories are used for self-similar shapes close to the growing facet. For long times and g_3/g_1 < 1, (a) a universal equation is derived for the shape profile, (b) the layer thickness varies as (g_3/g_1)^{1/3}, (c) distinct solutions are found for different g_3/_1, and (d) for conical shapes, the profile peak scales as (g_3/g_1)^{-1/6}. These results compare favorably with kinetic simulations.Comment: 4 pages including 3 figure

Slow fluctuations in enhanced Raman scattering and surface roughness relaxation

We propose an explanation for the recently measured slow fluctuations and ``blinking'' in the surface enhanced Raman scattering (SERS) spectrum of single molecules adsorbed on a silver colloidal particle. We suggest that these fluctuations may be related to the dynamic relaxation of the surface roughness on the nanometer scale and show that there are two classes of roughness with qualitatively different dynamics. The predictions agree with measurements of surface roughness relaxation. Using a theoretical model for the kinetics of surface roughness relaxation in the presence of charges and optical electrical fields, we predict that the high-frequency electromagnetic field increases both the effective surface tension and the surface diffusion constant and thus accelerates the surface smoothing kinetics and time scale of the Raman fluctuations in manner that is linear with the laser power intensity, while the addition of salt retards the surface relaxation kinetics and increases the time scale of the fluctuations. These predictions are in qualitative agreement with the Raman experiments

Relaxation kinetics in two-dimensional structures

We have studied the approach to equilibrium of islands and pores in two dimensions. The two-regime scenario observed when islands evolve according to a set of particular rules, namely relaxation by steps at low temperature and smooth at high temperature, is generalized to a wide class of kinetic models and the two kinds of structures. Scaling laws for equilibration times are analytically derived and confirmed by kinetic Monte Carlo simulations.Comment: 6 pages, 7 figures, 1 tabl

Area-preserving dynamics of a long slender finger by curvature: a test case for the globally conserved phase ordering

A long and slender finger can serve as a simple ``test bed'' for different phase ordering models. In this work, the globally-conserved, interface-controlled dynamics of a long finger is investigated, analytically and numerically, in two dimensions. An important limit is considered when the finger dynamics are reducible to the area-preserving motion by curvature. A free boundary problem for the finger shape is formulated. An asymptotic perturbation theory is developed that uses the finger aspect ratio as a small parameter. The leading-order approximation is a modification of ``the Mullins finger" (a well-known analytic solution) which width is allowed to slowly vary with time. This time dependence is described, in the leading order, by an exponential law with the characteristic time proportional to the (constant) finger area. The subleading terms of the asymptotic theory are also calculated. Finally, the finger dynamics is investigated numerically, employing the Ginzburg-Landau equation with a global conservation law. The theory is in a very good agreement with the numerical solution.Comment: 8 pages, 4 figures, Latex; corrected typo

A Hybrid Monte Carlo Method for Surface Growth Simulations

We introduce an algorithm for treating growth on surfaces which combines important features of continuum methods (such as the level-set method) and Kinetic Monte Carlo (KMC) simulations. We treat the motion of adatoms in continuum theory, but attach them to islands one atom at a time. The technique is borrowed from the Dielectric Breakdown Model. Our method allows us to give a realistic account of fluctuations in island shape, which is lacking in deterministic continuum treatments and which is an important physical effect. Our method should be most important for problems close to equilibrium where KMC becomes impractically slow.Comment: 4 pages, 5 figure

Dynamics of an Unbounded Interface Between Ordered Phases

We investigate the evolution of a single unbounded interface between ordered phases in two-dimensional Ising ferromagnets that are endowed with single-spin-flip zero-temperature Glauber dynamics. We examine specifically the cases where the interface initially has either one or two corners. In both examples, the interface evolves to a limiting self-similar form. We apply the continuum time-dependent Ginzburg-Landau equation and a microscopic approach to calculate the interface shape. For the single corner system, we also discuss a correspondence between the interface and the Young tableau that represents the partition of the integers.Comment: 9 pages, 11 figures, 2-column revtex4 format. V2: references added and discussion section expanded slightly. Final version for PRE. V3: A few small additional editorial change

Surface Instability of Icicles

Quantitatively-unexplained stationary waves or ridges often encircle icicles. Such waves form when roughly 0.1 mm-thick layers of water flow down the icicle. These waves typically have a wavelength of 1cm approximately independent of external temperature, icicle thickness, and the volumetric rate of water flow. In this paper we show that these waves can not be obtained by naive Mullins-Sekerka instability, but are caused by a quite new surface instability related to the thermal diffusion and hydrodynamic effect of thin water flow.Comment: 11 pages, 5 figures, Late

Crossover effects in the Wolf-Villain model of epitaxial growth in 1+1 and 2+1 dimensions

A simple model of epitaxial growth proposed by Wolf and Villain is investigated using extensive computer simulations. We find an unexpectedly complex crossover behavior of the original model in both 1+1 and 2+1 dimensions. A crossover from the effective growth exponent $\beta_{\rm eff}\!\approx\!0.37$ to $\beta_{\rm eff}\!\approx\!0.33$ is observed in 1+1 dimensions, whereas additional crossovers, which we believe are to the scaling behavior of an Edwards--Wilkinson type, are observed in both 1+1 and 2+1 dimensions. Anomalous scaling due to power--law growth of the average step height is found in 1+1 D, and also at short time and length scales in 2+1~D. The roughness exponents $\zeta_{\rm eff}^{\rm c}$ obtained from the height--height correlation functions in 1+1~D ($\approx\!3/4$) and 2+1~D ($\approx\!2/3$) cannot be simultaneously explained by any of the continuum equations proposed so far to describe epitaxial growth.Comment: 11 pages, REVTeX 3.0, IC-DDV-93-00

Coarsening of Surface Structures in Unstable Epitaxial Growth

We study unstable epitaxy on singular surfaces using continuum equations with a prescribed slope-dependent surface current. We derive scaling relations for the late stage of growth, where power law coarsening of the mound morphology is observed. For the lateral size of mounds we obtain $\xi \sim t^{1/z}$ with $z \geq 4$. An analytic treatment within a self-consistent mean-field approximation predicts multiscaling of the height-height correlation function, while the direct numerical solution of the continuum equation shows conventional scaling with z=4, independent of the shape of the surface current.Comment: 15 pages, Latex. Submitted to PR