Logarithmic islanding in submonolayer epitaxial growth

We investigate submonolayer epitaxial growth with a fixed monomer flux and irreversible aggregation of adatom islands due to an effective island diffusion, with a diffusivity for an mass k island proportional to . For , there is a steady state, while for , continuously evolving logarithmic islanding occurs in which the island density grows extremely slowly, as . In the latter regime, the island size distribution exhibits complex, but universal, multiple-scale mass dependence which we account for theoretically. Received: 3 June 1998 / Accepted: 13 July 1998     

Diffusion of adsorbates on random alloy surfaces

In this work the diffusion of non-interacting adsorbates on a random AB alloy surface is considered. For this purpose a simple cubic (sc), body-centered cubic (bcc) or face-centered cubic (fcc) auxiliary metal lattice is introduced. The auxiliary lattice is truncated parallel to its (100) plane in such a way that the fourfold hollow positions of the metal surface form a regular net of adsorption sites with square symmetry. The adsorption energy of each adsorption site is determined by its own environment, i.e. by the numbers of direct A or B neighbors. The Monte-Carlo method has been utilized to simulate surface diffusion of adsorbates on such energetically heterogeneous alloy surfaces and to calculate the tracer, jump and chemical diffusion coefficients. The chemical diffusion coefficient was calculated via two different approaches: the fluctuation and the Kubo-Green method. The influence of energetical heterogeneities on the surface diffusion is largely pronounced at low temperatures and low surface coverages, where most of the adatoms are trapped by deep adsorption sites. It was found that at low temperatures the sequential occupation of the different types of adsorption sites can be observed. Received: 24 October 1997 / Accepted: 17 December 1997     

The nonlinear Fokker-Planck equation with state-dependent diffusion - a nonextensive maximum entropy approach

Nonlinear Fokker-Planck equations (e.g., the diffusion equation for porous medium) are important candidates for describing anomalous diffusion in a variety of systems. In this paper we introduce such nonlinear Fokker-Planck equations with general state-dependent diffusion, thus significantly generalizing the case of constant diffusion which has been discussed previously. An approximate maximum entropy (MaxEnt) approach based on the Tsallis nonextensive entropy is developed for the study of these equations. The MaxEnt solutions are shown to preserve the functional relation between the time derivative of the entropy and the time dependent solution. In some particular important cases of diffusion with power-law multiplicative noise, our MaxEnt scheme provides exact time dependent solutions. We also prove that the stationary solutions of the nonlinear Fokker-Planck equation with diffusion of the (generalized) Stratonovich type exhibit the Tsallis MaxEnt form. Received 26 February 1999     

Study of surface-bulk mass transport and phase transformation in nano-TiO2 using hyperfine interaction technique

Phase transition from anatase to rutile for the 70 nm TiO₂ crystallite has been investigated by the time differential perturbed angular correlation (TDPAC) technique. The study involved the annealing of the TiO₂ nanocrystals, adsorbed with the nuclear probe (¹⁸¹Hf/¹⁸¹Ta) at trace level, at different temperatures for different durations. The TDPAC measurement was also supported by XRD measurement where the width of the peaks increases with the increase in annealing temperature indicating a crystal growth. The samples annealed up to 823 K for 4 h showed no phase transition, except for the growth of the crystallites. However, it showed phase transition at the same temperature (823 K), when annealed for longer duration, indicating the slower kinetics of the phase transition process. Further the sample, when annealed at 1123 K for 4 h, showed phase transition. It has also been observed that the ¹⁸¹Hf tracer, adsorbed on 70 nm anatase TiO₂, diffuses from surface to bulk during the phase transition process and the extent of diffusion in anatase differs from that in rutile phase. However, surface to bulk mass-transfer is found to play a significant role in the phase transition process.     

Growth mechanisms for wire-like epitaxial gold silicide islands on Si(1 1 0) surfaces

Epitaxial islands grown on various substrates are usually strained because of differences in lattice constants of the materials of the island and the substrate. Shape transition in the growth of strained islands has been proposed as a mechanism for strain relief and a way to form self-organized quantum wires. Shape transition usually leads to an elongated island growth. However, an elongated island growth may also be due to an anisotropic diffusion of material, the anisotropy being imposed by the symmetry of the substrate surface. In the present example, growth of gold silicide wire-like nanostructures on a Si(1 1 0) surface has been investigated by photoemission electron microscopy (PEEM). Growth of elongated unidirectional gold silicide islands, with an aspect ratio as large as 12:1, has been observed by PEEM following gold deposition on the Si substrate and subsequent annealing at the Au-Si eutectic temperature. Distribution of the width and the length of the gold silicide islands as a function of island area shows a feature similar to that for the shape transition. However, detailed investigations reveal that the elongated growth of gold silicide islands is rather mainly due to anisotropic diffusion of gold due to the twofold symmetry of the (1 1 0) surface of the Si substrate.     

Monte-Carlo analysis of tracer diffusion mechanisms in YBa2Cu3O6+c

A method for estimating, via the Monte-Carlo simulation, the most often realized diffusion mechanisms in 2D ordered structures is presented. Taking as an example the diffusion of oxygen ions in high temperature superconductor we propose several diffusion mechanisms and show to what extent they depend on the temperature and concentration of the diffusing particles. Our results are compared with the ones proposed earlier on the basis of energy arguments. We find also additional trajectories, different from those earlier reported in that system. Received 9 November 1998     

Influence of isotopic substitution on the diffusion and thermal diffusion coefficient of binary liquids

The mutual mass diffusion coefficient (D) and the thermal diffusion coefficient ( D _T) of the liquids acetone, benzene, benzene-d ₁, benzene-d ₃, benzene-d ₅, benzene-d ₆, benzene- ¹³C₆, n-hexane, toluene, 1, 2, 3, 4-tetrahydronaphtalene, isobutylbenzene, and 1, 6-dibromohexane in protonated and perdeuterated cyclohexane have been measured with a transient holographic grating technique at a temperature of 25 ^°C. The mass diffusion coefficient shows a pronounced concentration dependence. Perdeuteration of cyclohexane only leads to marginal changes of the mass diffusion coefficient. The Stokes-Einstein equation describes the limiting tracer diffusion coefficients well if the solute molecule is smaller than the solvent. It is not capable to describe the small isotope effect of a few percent. On the other hand, the isotope effect, which is independent of concentration, is in agreement with the Enskog theory, that does not provide the absolute value of the mass diffusion coefficient of the liquid mixtures. The thermal diffusion coefficient of all the binary mixtures shows a moderate and almost linear concentration dependence. Its isotope effect, which is the change of D _T upon deuteration of cyclohexane, varies with mole fraction. The thermophoretic force acting on any tracer molecule in cyclohexane changes by the same amount when cyclohexane is perdeuterated, irrespective of the magnitude of the thermophoretic force before deuteration. This change of the thermophoretic force is equal but of opposite sign to the difference between the thermophoretic forces acting on cyclohexane and perdeuterated cyclohexane as tracers in any of the above liquids.     

Kinetics of shape equilibration for two dimensional islands

We study the relaxation to equilibrium of two dimensional islands containing up to 20 000 atoms by Kinetic Monte Carlo simulations. We find that the commonly assumed relaxation mechanism - curvature-driven relaxation via atom diffusion - cannot explain the results obtained at low temperatures, where the island edges consist in large facets. Specifically, our simulations show that the exponent characterizing the dependence of the equilibration time on the island size is different at high and low temperatures, in contradiction with the above cited assumptions. Instead, we propose that - at low temperatures - the relaxation is limited by the nucleation of new atomic rows on the large facets: this allows us to explain both the activation energy and the island size dependence of the equilibration time. Received 7 December 1998 and Received in final form 18 March 1999     

DMRG studies on interchain coupling model for quasi-one-dimensional organic magnet

We introduce a solid-on-solid growth process which evolves by random deposition of dimers, surface diffusion, and evaporation of monomers from the edges of plateaus. It is shown that the model exhibits a robust transition from a smooth to a rough phase. The roughening transition is driven by an absorbing phase transition at the bottom layer of the interface, which displays the same type of critical behavior as the pair contact process with diffusion 2A↦3A, 2A↦. Received 14 October 2002 Published online 14 February 2003 RID="a" ID="a"e-mail: Haye.Hinrichsen@physik.uni-wuppertal.de     

Electron magnetic resonance investigation of gadolinium diffusion in zircon powders

The electron magnetic resonance (EMR) technique was used to investigate the diffusion of gadolinium in zircon (ZrSiO₄) powders. The EMR absorption intensity was measured for several annealing times and three different temperatures of isothermal annealing: 1273, 1323 and 1373 K. The activation energy for diffusion, calculated from the experimental data using a theoretical model based on the Fick equation, was found to be E_A=506±5 kJ mol⁻¹. This value is close to the ones for the diffusion of Gd in UO₂ and CeO₂, but much larger than for the diffusion of gadolinium in a compound with the same crystal structure as zircon, YVO₄. This is attributed to a difference in the relative sizes of the ions involved in the diffusion process.     

Diffusion and rheology in a model of glassy materials

We study self-diffusion within a simple hopping model for glassy materials. (The model is Bouchaud's model of glasses (J.-P. Bouchaud, J. Phys. I France 2, 1705 (1992)), as extended to describe rheological properties (P. Sollich, F. Lequeux, P. Hébraud, M.E. Cates, Phys. Rev. Lett. 78, 2020 (1997)).) We investigate the breakdown, near the glass transition, of the (generalized) Stokes-Einstein relation between self-diffusion of a tracer particle and the (frequency-dependent) viscosity of the system as a whole. This stems from the presence of a broad distribution of relaxation times of which different moments control diffusion and rheology. We also investigate the effect of flow (oscillatory shear) on self-diffusion and show that this causes a finite diffusivity in the temperature regime below the glass transition (where this was previously zero). At higher temperatures the diffusivity is enhanced by a power law frequency dependence that also characterises the rheological response. The relevance of these findings to soft glassy materials (foams, emulsions etc.) as well as to conventional glass-forming liquids is discussed. Received 31 August 1998 and Received in final form 25 January 1999     

Surface effects in field-induced smectic transitions

We study the surface behavior of a semi-infinite smectogenic sample bounded by a solid wall, in the presence of an external electric field. Our analysis is performed in the framework of a Landau-de Gennes theory. For the sake of simplicity, we consider only the case in which, in the absence of field and surfaces, a direct isotropic to smectic-A transition occurs, while in the presence of the electric field a nonspontaneous nematic phase appears. Two new surface phases are identified, namely a parasmectic and a surface-induced smectic phase. The shifts in the transition temperatures and the critical behavior of the surface states are analyzed. Received: 25 August 1997 / Accepted: 23 January 1998     

Non-linear evolution of step meander during growth of a vicinal surface with no desorption

Step meandering due to a deterministic morphological instability on vicinal surfaces during growth is studied. We investigate nonlinear dynamics of a step model with asymmetric step kinetics, terrace and line diffusion, by means of a multiscale analysis. We give the detailed derivation of the highly nonlinear evolution equation on which a brief account has been given [6]. Decomposing the model into driving and relaxational contributions, we give a profound explanation to the origin of the unusual divergent scaling of step meander (where F is the incoming atom flux). A careful numerical analysis indicates that a cellular structure arises where plateaus form, as opposed to spike-like structures reported erroneously in reference [6]. As a robust feature, the amplitude of these cells scales as t ^1/2, regardless of the strength of the Ehrlich-Schwoebel effect, or the presence of line diffusion. A simple ansatz allows to describe analytically the asymptotic regime quantitatively. We show also how sub-dominant terms from multiscale analysis account for the loss of up-down symmetry of the cellular structure. Received 4 May 2000 and Received in final form 8 September 2000     

Single-particle thermal diffusion of charged colloids: Double-layer theory in a temperature gradient

The double-layer contribution to the single-particle thermal diffusion coefficient of charged, spherical colloids with arbitrary double-layer thickness is calculated and compared to experiments. The calculation is based on an extension of the Debye-Hückel theory for the double-layer structure that includes a small temperature gradient. There are three forces that constitute the total thermophoretic force on a charged colloidal sphere due to the presence of its double layer: i) the force F _W that results from the temperature dependence of the internal electrostatic energy W of the double layer, ii) the electric force F _el with which the temperature-induced non-spherically symmetric double-layer potential acts on the surface charges of the colloidal sphere and iii) the solvent-friction force F _sol on the surface of the colloidal sphere due to the solvent flow that is induced in the double layer because of its asymmetry. The force F _W will be shown to reproduce predictions based on irreversible-thermodynamics considerations. The other two forces F _el and F _sol depend on the details of the temperature-gradient-induced asymmetry of the double-layer structure which cannot be included in an irreversible-thermodynamics treatment. Explicit expressions for the thermal diffusion coefficient are derived for arbitrary double-layer thickness, which complement the irreversible-thermodynamics result through the inclusion of the thermophoretic velocity resulting from the electric- and solvent-friction force.     

First-principles studies of kinetics in epitaxial growth of III–V semiconductors

We demonstrate how first-principles calculations using density-functional theory (DFT) can be applied to gain insight into the molecular processes that rule the physics of materials processing. Specifically, we study the molecular beam epitaxy (MBE) of arsenic compound semiconductors. For homoepitaxy of GaAs on GaAs (001), a growth model is presented that builds on results of DFT calculations for molecular processes on the β2-reconstructed GaAs (001) surface, including adsorption, desorption, surface diffusion, and nucleation. Kinetic Monte Carlo simulations on the basis of the calculated energetics enable us to model MBE growth of GaAs from beams of Ga and As₂ in atomistic detail. The simulations show that island nucleation is controlled by the reaction of As₂ molecules with Ga adatoms on the surface. The analysis reveals that the scaling laws of standard nucleation theory for the island density as a function of growth temperature are not applicable to GaAs epitaxy. We also discuss heteroepitaxy of InAs on GaAs (001), and report first-principles DFT calculations for In diffusion on the strained GaAs substrate. In particular, we address the effect of heteroepitaxial strain on the growth kinetics of coherently strained InAs islands. The strain field around an island is found to cause a slowing down of material transport from the substrate towards the island, and thus helps to achieve more homogeneous island sizes. Received: 2 May 2001 / Accepted: 23 July 2001 / Published online: 3 April 2002     

Simulation study of lateral diffusion in lipid-sterol bilayer mixtures

We employ off-lattice Monte Carlo simulations to study lateral diffusion in lipid-sterol bilayers using a two-dimensional model system which has been designed to simulate the experimental phase diagrams of both lipid-cholesterol and lipid-lanosterol systems. We focus on the effects of varying sterol concentration and temperature on the tracer diffusion coefficient, D, which characterizes the lateral motion of single tagged lipids in a bilayer. Generally, we find that increasing the cholesterol concentration suppresses D due to an increased conformational ordering of lipid chains. We argue that this effect competes with an increase in the average free area per lipid, which favours an increase in D. At temperatures close to the main transition temperature, the competition between the two effects leads to intriguing behavior of D. Overall, the model results are in excellent qualitative agreement with available experimental results for lipid-cholesterol mixtures. Additional studies of a model lipid-lanosterol system, for which experimental diffusion results are not available, predict that the presence of lanosterol has a smaller effect than cholesterol on reducing D relative to the pure lipid system. We conclude that the molecular model employed contains the essential features required to describe many of the qualitative features of the lateral diffusion behavior in lipid-sterol systems. Received 24 November 2000 and Received in final form 30 April 2001     

A DMRG study of the -symmetric Heisenberg chain

The spin one-half Heisenberg chain with U q [ SU (2)] symmetry is studied via density-matrix renormalization. Ground-state energy and q-symmetric correlation functions are calculated for the non-Hermitian case with integer r. This gives bulk and surface exponents for (para)fermionic correlations in the related Ising and Potts models. The case of real q corresponding to a diffusion problem is treated analytically. Received: 18 February 1998 / Accepted: 17 March 1998     

Disorder driven roughening transitions of elastic manifolds and periodic elastic media

The simultaneous effect of both disorder and crystal-lattice pinning on the equilibrium behavior of oriented elastic objects is studied using scaling arguments and a functional renormalization group technique. Our analysis applies to elastic manifolds, e.g., interfaces, as well as to periodic elastic media, e.g., charge-density waves or flux-line lattices. The competition between both pinning mechanisms leads to a continuous, disorder driven roughening transition between a flat state where the mean relative displacement saturates on large scales and a rough state with diverging relative displacement. The transition can be approached by changing the impurity concentration or, indirectly, by tuning the temperature since the pinning strengths of the random and crystal potential have in general a different temperature dependence. For D dimensional elastic manifolds interacting with either random-field or random-bond disorder a transition exists for 2<D<4, and the critical exponents are obtained to lowest order in . At the transition, the manifolds show a superuniversal logarithmic roughness. Dipolar interactions render lattice effects relevant also in the physical case of D=2. For periodic elastic media, a roughening transition exists only if the ratio p of the periodicities of the medium and the crystal lattice exceeds the critical value . For p<p _c the medium is always flat. Critical exponents are calculated in a double expansion in and and fulfill the scaling relations of random field models. Received 28 August 1998     

Growth of carbon nanotubes (CNTs) on metallic underlayers by diffusion plasma-enhanced chemical vapour deposition (DPECVD)

Here, we demonstrate the low-temperature (480–612 °C) synthesis of carbon nanotubes (CNTs) on different metallic underlayers (i.e., NiV, Ir, Ag, Pt, W, and Ta) using diffusion (dc) plasma-enhanced (~20 W, −600 V) chemical vapour deposition (DPECVD). The catalyst used is bi-layered Fe/Al and the feedstock used is a mixture of C₂H₂ and NH₃ (1:4). The crucial component is the diffusion of radical ions and hydrogen generated such as H₂/H⁺/H₂⁺/NH₃⁺/CH₂⁺/C₂H₂⁺ (which are confirmed by in-situ mass spectroscopy) from the nozzle, where it is inserted for most effective plasma diffusion between a substrate and a gas distributor.