Modeling of Surface Processes as Exemplified by Hydrocarbon Reactions

Temperature control in molecular dynamic simulations of non-equilibrium processes

Thermostats are often used in various condensed matter problems, e.g. when a biological molecule undergoes a transformation in a solution, a crystal surface is irradiated with energetic particles, a crack propagates in a solid upon applied stress, two surfaces slide with respect to each other, an excited local phonon dissipates its energy into a crystal bulk, and so on. In all of these problems, as well as in many others, there is an energy transfer between different local parts of the entire system kept at a constant temperature. Very often, when modelling such processes using molecular dynamics simulations, thermostatting is done using strictly equilibrium approaches serving to describe the NVT ensemble. In this paper we critically discuss the applicability of such approaches to non-equilibrium problems, including those mentioned above, and stress that the correct temperature control can only be achieved if the method is based on the generalized Langevin equation (GLE). Specifically, we emphasize that a meaningful compromise between computational efficiency and a physically appropriate implementation of the NVT thermostat can be achieved, at least for solid state and surface problems, if the so-called stochastic boundary conditions (SBC), recently derived from the GLE (Kantorovich and Rompotis 2008 Phys. Rev. B 78 094305), are used. For SBC, the Langevin thermostat is only applied to the outer part of the simulated fragment of the entire system which borders the surrounding environment (not considered explicitly) serving as a heat bath. This point is illustrated by comparing the performance of the SBC and some of the equilibrium thermostats in two problems: (i) irradiation of the Si(001) surface with an energetic CaF(2) molecule using an ab initio density functional theory based method, and (ii) the tribology of two amorphous SiO(2) surfaces coated with self-assembled monolayers of methyl-terminated hydrocarbon alkoxylsilane molecules using a classical atomistic force field. We discuss the differences in behaviour of these systems due to applied thermostatting, and show that in some cases a qualitatively different physical behaviour of the simulated system can be obtained if an equilibrium thermostat is used.

Computational view of surface based organic mass spectrometry

Surface based mass spectrometric approaches fill an important niche in the mass analysis portfolio of tools. The particular niche depends on both the underlying physics and chemistry of molecule ejection as well as experimental characteristics. In this article, we use molecular dynamics computer simulations to elucidate the fundamental processes giving rise to ejection of organic molecules in atomic and cluster secondary ion mass spectrometry (SIMS), massive cluster impact (MCI) mass spectrometry, and matrix-assisted laser desorption ionization (MALDI) mass spectrometry. This review is aimed at graduate students and experimental researchers.

Effect of cluster size in kiloelectronvolt cluster bombardment of solid benzene

Emission of benzene molecules by 5-keV cluster bombardment of a range of carbon projectiles from C6H6 to C180 is studied by a coarse-grained molecular dynamics (MD) technique. This approach permits calculations that are not feasible using more complicated potential energy functions, particularly as the interesting physics associated with the ion impact event approaches the mesoscale. These calculations show that the highest ejection yields are associated with clusters that deposit their incident energy 15-20 A below the surface. The highest yield for the projectiles is produced by the C20 and C60 projectiles. The results from the MD simulations are also compared favorably to an analytical model based on fluid dynamics to describe the energy deposition. The analytical model is then utilized to extend the range of the calculations to higher incident energies. The issue of the relative amount of chemical fragmentation and intact molecular desorption is also examined for the benzene crystal. These results show that damage accumulation at high-incident fluence should not be problematic and that it should be possible to perform molecular depth profiling via secondary ion mass spectrometry experiments. In general, the approach presented here illustrates the power of combining a simplified MD method with analytical strategies for describing a length scale that is difficult to achieve with traditional MD calculations.

IDEA: Interface dynamics and energetics algorithm

IDEA, interface dynamics and energetics algorithm, was implemented, in FORTRAN, under different operating systems to mimic dynamics and energetics of elementary events involved in interfacial processes. The code included a parallel elaboration scheme in which both the stochastic and the deterministic components, involved in the developed physical model, worked simultaneously. IDEA also embodied an optionally running VISUAL subroutine, showing the dynamic energy changes caused by the surface events, e.g., occurring at the gas-solid interface. Monte Carlo and ordinary differential equation system subroutines were employed in a synergistic way to drive the occurrence of the elementary events and to manage the implied energy flows, respectively. Biphase processes, namely isothermal and isobaric adsorption of carbon monoxide on nickel, palladium, and platinum surfaces, were first studied to test the capability of the code in modeling real frames. On the whole, the simulated results showed that IDEA could reproduce the inner characteristics of the studied systems and predict properties not yet experimentally investigated.

Sputtering of Water Ice Induced by C60 Bombardment: Onset of Plume Formation

The interaction of a 5 keV C(60) projectile with amorphous water ice is studied using molecular dynamics computer simulations. The energetic C(60) molecule causes large-scale collisional events in the subsurface region, involving more than 10(4) water molecules in a time of less than 3 ps. The energy deposited in the sample is sufficiently large to turn the ice into a superheated and superdense gas. The gas is expelled into the vacuum, leading to the formation of a flow that manifests itself in the angular and velocity distributions of emitted water molecules.

Microscopic Insights into the Sputtering of Thin Organic Films on Ag{111} Induced by C60 and Ga Bombardment

Molecular dynamics computer simulations have been employed to model the bombardment of Ag{111} covered with three layers of C6H6 by 15 keV Ga and C60 projectiles. The study is aimed toward examining the mechanism by which molecules are desorbed from surfaces by energetic cluster ion beams and toward elucidating the differences between cluster bombardment and atom bombardment. The results show that the impact of the cluster on the benzene-covered surface leads to molecular desorption during the formation of a mesoscopic scale impact crater via a catapulting mechanism. Because of the high yield of C6H6 with both Ga and C60, the yield enhancement is observed to be consistent with related experimental observations. Specific energy and angle distributions are shown to be associated with the catapult mechanism.

Enhancement of Sputtering Yields Due to C60versus Ga Bombardment of Ag{111} As Explored by Molecular Dynamics Simulations

The mechanism of enhanced desorption initiated by 15-keV C60 cluster ion bombardment of a Ag single crystal surface is examined using molecular dynamics computer simulations. The size of the model microcrystallite of 165,000 atoms and the sophistication of the interaction potential function yields data that should be directly comparable with experiment. The C60 model was chosen since this source is now being used in secondary ion mass spectrometry experiments in many laboratories. The results show that a crater is formed on the Ag surface that is approximately 10 nm in diameter, a result very similar to that found for Au3 bombardment of Au. The yield of Ag atoms is approximately 16 times larger than for corresponding atomic bombardment with 15-keV Ga atoms, and the yield of Ag3 is enhanced by a factor of 35. The essential mechanistic reasons for these differences is that the C60 kinetic energy is deposited closer to the surface, with the deeply penetrating energy propagation occurring via a nondestructive pressure wave. The numbers predicted by the model are testable by experiment, and the approach is extendable to include the study of organic overlayers on metals, a situation of growing importance to the SIMS community.

Molecule liftoff from surfaces

Molecular dynamics simulations have been used to model the kiloelectronvolt particle bombardment of organic layers on metal substrates such as occurs in the analytical techniques of secondary ion mass spectrometry and fast atom bombardment mass spectrometry. Vignettes of insights gained from the simulations along with comparisons to experimental data are presented in this Account. Topics include intact molecular ejection vs fragmentation, prediction of reaction pathways, influence of the substrate, and quantitative predictions of energy and angular distributions.

Surface mass spectrometry of molecular species

This tutorial discusses the predominant methods available for surface mass spectrometry (MS) of molecular species: thermal desorption spectroscopy, laser desorption MS, secondary ion MS, post-ionization of desorbed neutrals and surface matrix-assisted laser desorption/ionization. Each of these has the capability to analyze molecular species that are chemisorbed, physisorbed, covalently bound to or the predominant component of a solid surface. These surface MS methods are briefly described, then their capabilities demonstrated using data predominantly from the authors' work. Comparisons are made with related methods in conventional MS. A very brief discussion is provided on the importance of complementing surface MS data with data from x-ray photoelectron spectroscopy and other surface analysis tools.