Site determination and thermally assisted tunneling in homogenous nucleation

A combined low-temperature scanning tunneling microscopy and density functional theory study on the binding and diffusion of copper monomers, dimers, and trimers adsorbed on Cu(111) is presented. Whereas atoms in trimers are found in fcc sites only, monomers as well as atoms in dimers can occupy the stable fcc as well as the metastable hcp site. In fact the dimer fcc-hcp configuration was found to be only 1.3 meV less favorable with respect to the fcc-fcc configuration. This enables a confined intra-cell dimer motion, which at temperatures below 5 K is dominated by thermally assisted tunneling.Comment: 4 pages, 4 figure

Diffusion rates of Cu adatoms on Cu(111) in the presence of an adisland nucleated at FCC or HCP sites

The surface diffusion of Cu adatoms in the presence of an adisland at FCC or HCP sites on Cu(111) is studied using the EAM potential derived by Mishin {\it et al.} [Phys. Rev. B {\bf 63} 224106 (2001)]. The diffusion rates along straight (with close-packed edges) steps with (100) and (111)-type microfacets (resp. step A and step B) are first investigated using the transition state theory in the harmonic approximation. It is found that the classical limit beyond which the diffusion rates follow an Arrhenius law is reached above the Debye temperature. The Vineyard attempt frequencies and the (static) energy barriers are reported. Then a comparison is made with the results of more realistic classical molecular dynamic simulations which also exhibit an Arrhenius-like behavior. It is concluded that the corresponding energy barriers are completely consistent with the static ones within the statistical errors and that the diffusion barrier along step B is significantly larger than along step A. In contrast the prefactors are very different from the Vineyard frequencies. They increase with the static energy barrier in agreement with the Meyer-Neldel compensation rule and this increase is well approximated by the law proposed by Boisvert {\it et al.} [Phys. Rev. Lett. {\bf 75} 469 (1995)]. As a consequence, the remaining part of this work is devoted to the determination of static energy barriers for a large number of diffusion events that can occur in the presence of an adisland. In particular, it is found that the corner crossing diffusion process for triangular adislands is markedly different for the two types of borders (A or B). From this set of results the diffusion rates of the most important atomic displacements can be predicted and used as input in Kinetic Monte-Carlo simulations

Dynamics of the rotational degrees of freedom in a supercooled liquid of diatomic molecules

Using molecular dynamics computer simulations, we investigate the dynamics of the rotational degrees of freedom in a supercooled system composed of rigid, diatomic molecules. The interaction between the molecules is given by the sum of interaction-site potentials of the Lennard-Jones type. In agreement with mode-coupling theory (MCT), we find that the relaxation times of the orientational time correlation functions C_1^(s), C_2^(s) and C_1 show at low temperatures a power-law with the same critical temperature T_c, and which is also identical to the critical temperature for the translational degrees of freedom. In contrast to MCT we find, however, that for these correlators the time-temperature superposition principle does not hold well and that also the critical exponent gamma depends on the correlator. We also study the temperature dependence of the rotational diffusion constant D_r and demonstrate that at high temperatures D_r is proportional to the translational diffusion constant D and that when the system starts to become supercooled the former shows an Arrhenius behavior whereas the latter exhibits a power-law dependence. We discuss the origin for the difference in the temperature dependence of D (or the relaxation times of C_l^(s) and D_r. Finally we present results which show that at low temperatures 180 degree flips of the molecule are an important component of the relaxation dynamics for the orientational degrees of freedom.Comment: 17 pages of RevTex, 12 figure

Potential Energy Landscape Equation of State

Depth, number, and shape of the basins of the potential energy landscape are the key ingredients of the inherent structure thermodynamic formalism introduced by Stillinger and Weber [F. H. Stillinger and T. A. Weber, Phys. Rev. A 25, 978 (1982)]. Within this formalism, an equation of state based only on the volume dependence of these landscape properties is derived. Vibrational and configurational contributions to pressure are sorted out in a transparent way. Predictions are successfully compared with data from extensive molecular dynamics simulations of a simple model for the fragile liquid orthoterphenyl.Comment: RevTeX4, 4 pages, 5 figure

Pressure induced structural and dynamical changes in liquid Si. An ab-initio study

The static and dynamic properties of liquid Si at high-pressure have been studied using the orbital free ab-initio molecular dynamics method. Four thermodynamic states at pressures 4, 8, 14 and 23 GPa are considered. The calculated static structure shows qualitative agreement with the available experimental data. We analize the remarkable structural changes occurring between 8 and 14 GPa along with its effect on several dynamic properties.Comment: 10 pages, 11 figures. Accepted for publication in Journal of Physics: Condensed Matte

Structural Relaxation and Mode Coupling in a Simple Liquid: Depolarized Light Scattering in Benzene

We have measured depolarized light scattering in liquid benzene over the whole accessible temperature range and over four decades in frequency. Between 40 and 180 GHz we find a susceptibility peak due to structural relaxation. This peak shows stretching and time-temperature scaling as known from $\alpha$ relaxation in glass-forming materials. A simple mode-coupling model provides consistent fits of the entire data set. We conclude that structural relaxation in simple liquids and $\alpha$ relaxation in glass-forming materials are physically the same. A deeper understanding of simple liquids is reached by applying concepts that were originally developed in the context of glass-transition research.Comment: submitted to New J. Phy

Testing Mode-Coupling Theory for a Supercooled Binary Lennard-Jones Mixture II: Intermediate Scattering Function and Dynamic Susceptibility

We have performed a molecular dynamics computer simulation of a supercooled binary Lennard-Jones system in order to compare the dynamical behavior of this system with the predictions of the idealized version of mode-coupling theory (MCT). By scaling the time $t$ by the temperature dependent $\alpha$-relaxation time $\tau(T)$, we find that in the $\alpha$-relaxation regime $F(q,t)$ and $F_s(q,t)$, the coherent and incoherent intermediate scattering functions, for different temperatures each follows a $q$-dependent master curve as a function of scaled time. We show that during the early part of the $\alpha$-relaxation, which is equivalent to the late part of the $\beta$-relaxation, these master curves are well approximated by the master curve predicted by MCT for the $\beta$-relaxation. This part is also fitted well by a power-law, the so-called von Schweidler law. We show that the effective exponent $b'$ of this power-law depends on the wave vector $q$ if $q$ is varied over a large range. The early part of the $\beta$-relaxation regime does not show the critical decay predicted by MCT. The $q$-dependence of the nonergodicity parameter for $F_{s}(q,t)$ and $F(q,t)$ are in qualitative agreement with MCT. On the time scale of the late $\alpha$-relaxation the correlation functions show a Kohlrausch-Williams-Watt behavior (KWW). The KWW exponent $\beta$ is significantly different from the effective von Schweidler exponent $b'$. At low temperatures the $\alpha$-relaxation time $\tau(T)$ shows a power-law behavior with a critical temperature that is the same as the one found previously for the diffusion constant [Phys. Rev. Lett. {\bf 73}, 1376 (1994)]. The critical exponent of this power-law and the von Schweidler exponent $b'$ fulfill the connection proposed by MCT between these two quantities. We also show that theComment: 28 Pages of REVTEX, Figures available from W. Ko

Molecular dynamics simulation of the fragile glass former ortho-terphenyl: a flexible molecule model

We present a realistic model of the fragile glass former orthoterphenyl and the results of extensive molecular dynamics simulations in which we investigated its basic static and dynamic properties. In this model the internal molecular interactions between the three rigid phenyl rings are described by a set of force constants, including harmonic and anharmonic terms; the interactions among different molecules are described by Lennard-Jones site-site potentials. Self-diffusion properties are discussed in detail together with the temperature and momentum dependencies of the self-intermediate scattering function. The simulation data are compared with existing experimental results and with the main predictions of the Mode Coupling Theory.Comment: 20 pages and 28 postscript figure