S. C. Ying

Phase diagram and commensurate-incommensurate transitions in the phase field crystal model with an external pinning potential.

We study the phase diagram and the commensurate-incommensurate transitions in a phase field model of a two-dimensional crystal lattice in the presence of an external pinning potential. The model allows for both elastic and plastic deformations and provides a continuum description of lattice systems, such as for adsorbed atomic layers or two-dimensional vortex lattices. Analytically, a mode expansion analysis is used to determine the ground states and the commensurate-incommensurate transitions in the model as a function of the strength of the pinning potential and the lattice mismatch parameter. Numerical minimization of the corresponding free energy shows reasonable agreement with the analytical predictions and provides details on the topological defects in the transition region. We find that for small mismatch the transition is of first order, and it remains so for the largest values of mismatch studied here. Our results are consistent with results of simulations for atomistic models of adsorbed overlayers.

Quantum Diffusion of H/Ni(111) through a Monte Carlo Wave Function Formalism

We consider a quantum system coupled to a dissipative background with many degrees of freedom using the Monte Carlo wave function method. Instead of dealing with a density matrix which can be very highly dimensional, the method consists of integrating a stochastic Schrödinger equation with a non-Hermitian damping term in the evolution operator, and with random quantum jumps. The method is applied to the diffusion of hydrogen on the Ni(111) surface below 100 K. We show that the recent experimental diffusion data for this system can be understood through an interband activation process, followed by quantum tunneling.

A dynamical mean field theory for the study of surface diffusion constants

Abstract We present a combined analytical and numerical approach based on the Mori projection operator formalism and Monte Carlo simulations to study surface diffusion within the lattice-gas model. In the present theory, the average jump rate and the susceptibility factor appearing are evaluated through Monte Carlo simulations, while the memory functions are approximated by the known results for a Langmuir gas model. This leads to a dynamical mean field theory (DMF) for collective diffusion, while approximate correlation effects beyond DMF are included for tracer diffusion. We apply our formalism to three very different strongly interacting systems and compare the results of the new approach with those of usual Monte Carlo simulations. We find that the combined approach works very well for collective diffusion, whereas for tracer diffusion the influence of interactions on the memory effects is more prominent.

Transverse thermal depinning and nonlinear sliding friction of an adsorbed monolayer.

We study the response of an adsorbed monolayer under a driving force as a model of sliding friction phenomena between two crystalline surfaces with a boundary lubrication layer. Using Langevin-dynamics simulation, we determine the nonlinear response in the direction transverse to a high symmetry direction along which the layer is already sliding. We find that below a finite transition temperature there exist a critical depinning force and hysteresis effects in the transverse response in the dynamical state when the adlayer is sliding smoothly along the longitudinal direction.

Dynamical transitions and sliding friction in the two-dimensional Frenkel-Kontorova model

The nonlinear response of an adsorbed layer on a periodic substrate to an external force is studied via a two dimensional uniaxial Frenkel-Kontorova model. The nonequlibrium properties of the model are simulated by Brownian molecular dynamics. Dynamical phase transitions between pinned solid, sliding commensurate and incommensurate solids and hysteresis effects are found that are qualitatively similar to the results for a Lennard-Jones model, thus demonstrating the universal nature of these features.

Comment on “Surface diffusion near the points corresponding to continuous phase transitions” [J. Chem. Phys. 109, 3197 (1998)]

It is well known that unlike static equilibrium properties, kinetic quantities in Monte Carlo simulations are very sensitive to the details of the algorithm used for the microscopic transition rates. This is particularly true near the critical region where fluctuations are pronounced. We demonstrate that when diffusion of oxygen adatoms near the order–disorder transition of a lattice-gas model of the O/W(110) model system is studied, the transition rates must be chosen carefully. In particular, we show that the choice by Uebing and Zhdanov [J. Chem. Phys. 109, 3197 (1998)] is inappropriate for the study of critical effects in diffusion.

Non-equilibrium surface diffusion in the OW(110) system

In this Letter, we present results of an extensive Monte Carlo study of the O/W(110) system under non-equilibrium conditions. We study the mean square displacements and long wavelength density fluctuations of adatoms. From these quantities, we define effective and time-dependent values for the collective and tracer diffusion mobilities. These mobilities reduce to the usual diffusion constants when equilibrium is reached. We discuss our results in view of existing experimental measurements of effective diffusion barriers, and the difficulties associated with interpreting non-equilibrium data.

Nonlinear driven response of a phase-field crystal in a periodic pinning potential.

We study numerically the phase diagram and the response under a driving force of the phase field crystal model for pinned lattice systems introduced recently for both one- and two-dimensional systems. The model describes the lattice system as a continuous density field in the presence of a periodic pinning potential, allowing for both elastic and plastic deformations of the lattice. We first present results for phase diagrams of the model in the absence of a driving force. The nonlinear response to a driving force on an initially pinned commensurate phase is then studied via overdamped dynamic equations of motion for different values of mismatch and pinning strengths. For large pinning strength the driven depinning transitions are continuous, and the sliding velocity varies with the force from the threshold with power-law exponents in agreement with analytical predictions. Transverse depinning transitions in the moving state are also found in two dimensions. Surprisingly, for sufficiently weak pinning potential we find a discontinuous depinning transition with hysteresis even in one dimension under overdamped dynamics. We also characterize structural changes of the system in some detail close to the depinning transition.

Velocity correlations and memory functions in surface diffusion

Abstract We introduce a generalized velocity correlation function and the corresponding memory function for the study of tracer diffusion of interacting particles on surfaces. These functions can be obtained by using discrete particle displacement variables either from experimental or simulation data. We study the behavior of these functions for diverse systems and find that in most cases, for both functions, there is an intermediate power-law decay ∝ t − x which spans about two orders of magnitude in time. The exponent x depends on the range and the strength of interactions. For the Langmuir gas with on-site exclusions only, x ≈2. For other interacting systems we find that, when strong attractive interactions are present, x tends to decrease from two to unity, while strong repulsion may cause x to be larger than two. These results demonstrate that the velocity correlation function can be used to gage the importance of interaction effects for surface diffusion.

Two-dimensional lattice hamiltonian for surface reconstruction

Abstract We demonstrate that the strong coupling surface structural phase transition as observed on the (001) face of bcc group transition metals can be studied with an effective two-dimensional lattice Hamiltonian. This Hamiltonian provides an accurate description of the zero temperature properties and is convenient for a detailed study of the finite temperature phase diagram and critical fluctuation properties of the structural phase transition.