Search

We study the no-slip boundary conditions for water at a hydrophobic (graphite) surface using non-equilibrium molecular-dynamics simulations. For the planar Couette flow, we find a slip length of 64 nm at 1 bar and 300 K, decreasing with increasing system pressure to a value of 31 nm at 1000 bar. Changing the properties of the interface to from hydrophobic to strongly hydrophilic reduces the slip to 14 nm. Finally, we study the flow of water past an array of carbon nanotubes mounted in an inline configuration with a spacing of 16.4 x 16.4 nm. For tube diameters of 1.25 and 2.50 nm we find drag coefficients in good agreement with the macroscopic, Navier-Stokes values. For carbon nanotubes, the no-slip condition is valid to within the definition of the position of the interface.

We study the hydrophobic/hydrophilic behavior of carbon nanotubes using molecular dynamics simulations. The energetics of the carbon-water interface are mainly dispersive but in the present study augmented with a carbon quadrupole term acting on the charge sites of the water. The simulations indicate that this contribution is negligible in terms of modifying the structural properties of water at the interface. Simulations of two carbon nanotubes in water display a wetting and drying of the interface between the nanotubes depending on their initial spacing. Thus, initial tube spacings of 7 and 8 A resulted in a drying of the interface whereas spacing of > 9 A remain wet during the course of the simulation. Finally, we present a novel particle-particle-particle-mesh algorithm for long range potentials which allows for general (curvilinear) meshes and "black-box" fast solvers by adopting an influence matrix technique.

Recently developed semiempirical potential energy functions for the Ga-As-Si and Ga-As-Al systems have been applied here to determine the excess formation energy for GaAs clusters on GaAs(00-1), AlAs(00-1), Si(001), and one atomic layer As-covered Si(001) substrates as a function of cluster size and cluster shape by the Monte Carlo technique. Pyramidal type ledges on the GaAs clusters are found to be the favored ledge for the first three layers while an inverted-pyramidal type ledge is also favored in certain cases for the As1/Si(001) substrate. Cluster formation at ledges is compared with cluster formation on a flat terrace for the Si(001) and the As1/Si(001) substrates.

Semiempirical potential energy functions have been utilized for a variety of calculations in the Ga-Al-As system. Surface energies have been calculated for several orientations of GaAs. Ledge energies for the GaAs (0 0 -1) (As-terminated) surface show long-range interaction effects with the ledge energy increasing with spacing. GaAs (001)/AlAs (001) superlattices have been simulated for a range of interlayer spacings with the excess interfacial energy per interlayer increasing from 5 erg/sq cm at an interlayer spacing of 1 molecular layer (5.8 A) to 50 erg/sq cm at an interlayer spacing of 18 molecular layers (103.8 A).

Both strip and triangular clusters, composed of 2 -1 -1 line ledges, have been simulated on the Si (111) surface. The long-range ledge-ledge interaction and the surface stress tensor distribution have been evaluated for these two pill-box geometries using a semiempirical potential-energy function that incorporates both two-body and three-body contributions. The consequences of the ledge-ledge interaction on two-dimensional nucleation for Si (111) has been evaluated as a function of Si adatom supersaturation and shown to differ significantly from conventional theory, where such interaction is neglected.

Using a computer simulation technique with a semiempirical potential, a Si crystal was cleaved along the (111) plane. The pi-bonded chain structural features of the Si(111) cleavage surface are observed and found to be a consequence of the dynamics of this cleavage process and seem not to be influenced by the final energetics.

During this reporting period three investigations were carried out. The first area of research concerned the analysis of the structure-energy relationship in small clusters. This study is very closely related to the improvement of the potential energy functions which are suitable and simple enough to be used in atomistic simulation studies. Parameters obtained from ab initio calculations for dimers and trimers of Al were used to estimate energetics and global minimum energy structures of clusters continuing up to 15 Al atoms. The second research topic addressed modeling of the collision process for atoms impinging on surfaces. In this simulation study qualitative aspects of the O atom collision with a graphite surface were analyzed. Four different O/graphite systems were considered and the aftermath of the impact was analyzed. The final area of investigation was related to the simulation of thin amorphous Si films on crystalline Si substrates. Parameters obtained in an earlier study were used to model an exposed amorphous Si surface and an a-Si/c-Si interface. Structural details for various film thicknesses were investigated at an atomistic level.

A direct and convenient method is presented for deriving expressions which equate any thermodynamic state function to averages of specific dynamical functions and their fluctuations over the classical microcanonical distribution. Specific expressions are obtained for a variety of thermodynamic quantities. The effect of various entropy definitions on the results are assessed, and the latter are compared to previous work in the literature. The derived formulas are applied to the analysis of molecular-dynamics computer simulations.

From intermolecular force studies, it is now known that the overall non-additive contribution to the lattice enegy is positive so that analysis based on only pairwise additivity suggests a shallower intermolecular potential than the true value. Two body contributions alone are also known to be categorically unable to even qualitatively describe some configurations of molecular clusters in the gas phase or the general relaxation and reconstruction of fcc crystal surfaces. In addition, the many-body contribution was shown to play a key role in the stability of certain crystal structures. In these recent analyses, a relatively simple potential energy function (PEF), comprising only a two-body Mie-type potential plus a three-body Axilrod-Teller-type potential, was found to be extremely effective. This same parametric PEF is applied to describe the bulk stability and surface energy for the diamond cubic structure. To test the stability condition, the FCC, BCC, diamond cubic, graphite and beta-tin structures were considered.

The results of an atomistic model for the simulation of the early stages of crack initiation in a two-dimensional triangular lattice are presented. In the current model, each particle in the system is treated discretely and assumed to be interacting with the surrounding particles via Lennard-Jones potentials. A uniaxial load (in incremental elongations) is applied to the rectangular two-dimensional slab in either the x or the y direction. After each incremental elongation the system is equilibrated using a static method. Initially, elastic behavior in the x and y directions is observed. Continued elongation results in plastic deformation. In lattices with point defects, the defects first move to the surface, creating vacancies which trigger plastic deformation.

Molecular dynamics was used to study the structure, dispersion and short-time behavior of ten-atom clusters adsorbed onto amorphous and crystalline substrates, in which the cluster atoms differed from the substrate atoms. Two adatom-substrate model systems were chosen; one, in which the interaction energy between adatom pairs was greater than that between substrate pairs, and the other, in which the reverse was true. At relatively low temperature ranges, increased dispersion of cluster atoms occurred: (a) on the amorphous substrate as compared to the FCC(100) surface, (b) with increasing reduced temperature, and (c) with adatom-substrate interaction energy stronger than adatom-adatom interaction. Two-dimensional clusters (rafts) on the FCC(100) surface displayed migration of edge atoms only, indicating a mechanism for the cluster rotation and shape changes found in experimental studies.

Minimum energy configurations of microclusters (up to six atoms) have been calculated using two- and three-body interactions. Structural changes were parametrically analyzed as a function of the intensity of three-body forces. The results are qualitative in nature; they indicate, however, that three-body interactions play an important role in the equilibrium structure of microclusters. The effect of the intensity of the three-body interactions on the structure of small clusters is not manifested in a continuous manner. Rather, changes in the energetically most stable structure occur abruptly. The results are in qualitative agreement with experimental observations as well as other calculations.

A rate equation describing the desorption of particles from substrates has been formulated using a simple classical model in terms of the particulate escape probabilities and escape velocities. Numerical results for desorption rates as a function of temperature using simple two-body potentials are calculated for various metallic systems. These results were obtained for the low coverage limits and are found to be in good agreement with experimental findings. The present formulation also provides an improved picture of the physical nature of the commonly discussed parameters, such as the frequency factor, desorption energy, and desorption temperature in relation to observed desorption phenomena.

Results of an atomistic calculation supporting two-dimensional diffusion of a Pt atom on a Pt(110) surface are presented. The calculations for diffusion energies indicate that interchannel diffusion by the replacement of a channel-wall atom by the adatom is energetically favorable.

The activation energies for diffusion were determined for gold, platinum and iridium adatoms on plane and plane PT surfaces and were found to be in good agreement with the measurements reported by Bassett and Webber. The Lennard-Jones pair potentials were used to model the interatomic forces, and relaxation of the substrate atoms in near proximity to the adatom was considered in detail. The present calculations clarify the mechanism of the observed two-dimensional diffusion of platinum and iridium atoms on a plane PT surface. The results are compared with those obtained using Morse potential functions and different relaxation techniques.