“…Since experimental data for surface phonon frequencies [3] and another theoretical calculation [9] are available on Ni(9 7 7), in this paper we use the MD method and MAEAM potential to calculate the temperature dependence of the surface relaxation, phonon spectrum, mean square displacement of atoms, and the layer structure factor for the stepped (9 7 7) surface of Ni up to 1700 K, which is only 60 K below the bulk melting point. Present calculation is carried out for much higher temperature than real space GreenÕs function approach [9].…”
Section: Introductionmentioning
confidence: 99%
“…This method provides the information about the surface phonon density of states and permits a clear discrimination for the surface localized modes and surface resonances. Using the GreenÕs function method with the EAM potential, Rahman and coworkers [9,10] calculated the local vibrational densities of (9 7 7) for Ni and (2 1 1), (3 3 1), (5 1 1) for Cu. In addition to phonon frequencies, the GreenÕs function has also been used to calculate the mean square vibrational amplitude of atoms, the surface specific heat and other thermodynamic quantities.…”
“…Since experimental data for surface phonon frequencies [3] and another theoretical calculation [9] are available on Ni(9 7 7), in this paper we use the MD method and MAEAM potential to calculate the temperature dependence of the surface relaxation, phonon spectrum, mean square displacement of atoms, and the layer structure factor for the stepped (9 7 7) surface of Ni up to 1700 K, which is only 60 K below the bulk melting point. Present calculation is carried out for much higher temperature than real space GreenÕs function approach [9].…”
Section: Introductionmentioning
confidence: 99%
“…This method provides the information about the surface phonon density of states and permits a clear discrimination for the surface localized modes and surface resonances. Using the GreenÕs function method with the EAM potential, Rahman and coworkers [9,10] calculated the local vibrational densities of (9 7 7) for Ni and (2 1 1), (3 3 1), (5 1 1) for Cu. In addition to phonon frequencies, the GreenÕs function has also been used to calculate the mean square vibrational amplitude of atoms, the surface specific heat and other thermodynamic quantities.…”
“…The reduced coordination of the step atoms on a vicinal surface makes the average vibration amplitude of vicinal surface atoms significantly higher than the vibration amplitude on low-index surfaces [23][24][25]. As a consequence, the excess vibrational entropy of a vicinal surface is higher than that of the corresponding low-index facet combination.…”
mentioning
confidence: 96%
“…In the case of faceting, all steps are brought together to form a (111) terrace, and the step-free remainder of the surface has the (100) orientation. We assume for simplicity that the surface atoms are harmonic, isotropically vibrating Einstein oscillators, and that the step atoms on the (n11) surface have a vibration amplitude s step that is 20% higher than that on the (111) surface, s ͑111͒ [23][24][25]. The other atoms at the (n11) surface are assumed to have the same vibration amplitude as the terrace atoms on the (100) surface [25].…”
mentioning
confidence: 99%
“…As a consequence, the excess vibrational entropy of a vicinal surface is higher than that of the corresponding low-index facet combination. For a proper calculation of the differences between the excess entropies, the full vibrational energy spectrum has to be computed for each of the surface orientations involved [23,24]. Here, we give only a simple estimate of the effect for a surface with the (n11) orientation.…”
Molecular‐dynamics (MD) simulation is a well‐developed numerical technique that involves the use of a suitable algorithm to solve the classical equations of motion for atoms interacting with a known interatomic potential. This method has been used for several decades now to illustrate and understand the temperature and pressure dependencies of dynamical phenomena in liquids, solids, and liquid‐solid interfaces.
MD simulation techniques are also well suited for studying surface phenomena, as they provide a qualitative understanding of surface structure and dynamics. This article considers the use of MD techniques to better understand surface disorder and premelting. Specifically, it examines the temperature dependence of structure and vibrational dynamics at surfaces of face‐centered cubic (fcc) metals—mainly Ag, Cu, and Ni. It also makes contact with results from other theoretical and experimental methods.
While the emphasis in this article is on metal surfaces, the MD technique has been applied over the years to a wide variety of surfaces including those of semiconductors, insulators, alloys, glasses, and simple or binary liquids. A full review of the pros and cons of the method as applied to these very interesting systems is beyond the scope of this article. However, it is worth pointing out that the success of the classical MD simulation depends to a large extent on the exactness with which the forces acting on the ion cores can be determined. On semiconductor and insulator surfaces, the success of
ab initio
molecular dynamics simulations has made them more suitable for such calulations, rather than classical MD simulations.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.