Search NASA⌕ Search

Engineering topics

Tibbits, P.

Publications and source records attributed to Tibbits, P..

Surface disordering of Pb(110)

A molecular dynamics simulation of Pb(110), using embedded atom method energy functional and two-body potential derived for the purpose, shows that near 400 K the three surface layers begin to disorder well before the bulk (interior) layers. Transfer of atoms from the first subsurface layer to the surface layer occurs. Disordering propagates beyond the top three layers above 550 K, accompanied by formation of an adlayer and vacancies in the top three layers. Behavior of the two-dimensional layer structure factors indicates that disordering is anisotropic. Simulation results are consistent with experimental observations of surface roughening near 400 K and more extensive surface ordering above 525 K. Results are consistent with simulations for Ni and Al.

Tibbits, P.↗

Molecular Dynamics study of Pb overlayer on Cu(100)

Isothermal-isobaric Molecular Dynamics (MD) simulation of a submonolayer Pb film in c(2x2) ordered structure adsorbed on a Cu(100) substrate showed retention of order to high T. The Embedded Atom Method (EAM) calculated the energy of atoms of overlayer and substrate. The time-averaged squared modulus of the two dimensional structure factor for the Pb overlayer measured the order of the overlayer. The results are for increasing T only, and require verification by simulated cooling.

Karimi, M.↗

Order/disorder and phase diagram of H on Pd(100)

A phase boundary for H-Pd(100) was calculated using the Metropolis (1953) algorithm and the embedded atom method (EAM) described by Daw and Foiles (1987). The calculated phase boundary agreed with an experimentally determined phase boundary in its curvature and the coverage at which maximum Tc appeared, but was about 125 K lower than the experimental phase boundary.

Tibbits, P.↗

Application of the embedded atom method to Pb and Be

We have derived the embedding energy functional and two-body potential of the embedded atom method using decreasing exponentials for both the electron density and the two body potential. The embedding function was obtained from the equation of state given by Rose et al. (1984). Because of the form of the embedding function, the equilibrium lattice constant, cohesive energy, and bulk modulus are automatically satisfied. The two parameters phi(e) and gamma of the two-body potential were determined by fitting to shear modulus and the single vacancy formation energy. Contributions of up to the third nearest neighbors were included in the evaluation of the charge density rho and the two-body potential phi. The stability and anisotropy of each structure were estimated and compared with the available experimental data.

Karimi, M.↗