Reaction dynamics of atomic hydrogen with the hydrogenated Si(001) (2×1) surface

A hierarchical effective mode approach to phonon-driven multilevel vibrational relaxation dynamics at surfaces

We discuss an efficient Hierarchical Effective Mode (HEM) representation of a high-dimensional harmonic oscillator bath, which describes phonon-driven vibrational relaxation of an adsorbate-surface system, namely, deuterium adsorbed on Si(100). Starting from the original Hamiltonian of the adsorbate-surface system, the HEM representation is constructed via iterative orthogonal transformations, which are efficiently implemented with Householder matrices. The detailed description of the HEM representation and its construction are given in the second quantization representation. The hierarchical nature of this representation allows access to the exact quantum dynamics of the adsorbate-surface system over finite time intervals, controllable via the truncation order of the hierarchy. To study the convergence properties of the effective mode representation, we solve the time-dependent Schrödinger equation of the truncated system-bath HEM Hamiltonian, with the help of the multilayer extension of the Multiconfigurational Time-Dependent Hartree (ML-MCTDH) method. The results of the HEM representation are compared with those obtained with a quantum-mechanical tier-model. The convergence of the HEM representation with respect to the truncation order of the hierarchy is discussed for different initial conditions of the adsorbate-surface system. The combination of the HEM representation with the ML-MCTDH method provides information on the time evolution of the system (adsorbate) and multiple effective modes of the bath (surface). This permits insight into mechanisms of vibration-phonon coupling of the adsorbate-surface system, as well as inter-mode couplings of the effective bath.

Angular distributions of H-induced HD and D2 desorptions from the Si(100) surfaces

We measured angular distributions of HD and D2 molecules desorbed via the reactions H+DSi(100)-->HD [abstraction (ABS)] and H+DSi(100)-->D2 [adsorption-induced-desorption (AID)], respectively. It was found that the angular distribution of HD molecules desorbed along ABS is broader than that of D2 molecules desorbed along AID, i.e., the former could be fit with cos(2.0+/-0.2) theta, while the latter with cos(5.0+/-0.5) theta. This difference of the angular distributions between the two reaction paths suggests that their dynamic mechanisms are different. The observed cos2 theta distribution for the ABS reaction was reproduced by the classical trajectory calculations over the London-Eyring-Polanyi-Sato potential-energy surfaces. The simulation suggests that the HD desorption along the ABS path takes place along the direction of Si-D bonds, but the apparent angular distribution is comprised of multiple components reflecting the different orientations of D-occupied Si dimers in the (2 x 1) and (1 x 2) double domain structures.

Internal pressure and surface tension of bare and hydrogen coated silicon nanoparticles

We present a study of internal pressure and surface tension of bare and hydrogen coated silicon nanoparticles of 2-10 nm diameter as a function of temperature, using molecular dynamics simulations employing a reparametrized Kohen-Tully-Stillinger interatomic potential. The internal pressure was found to increase with decreasing particle size but the density was found to be independent of the particle size. We showed that for covalent bond structures, changes in surface curvature and the associated surface forces were not sufficient to significantly change bond lengths and angles. Thus, the surface tension was also found to be independent of the particle size. Surface tension was found to decrease with increasing particle temperature while the internal pressure did not vary with temperature. The presence of hydrogen on the surface of a particle significantly reduces surface tension (e.g., drops from 0.83 J/m(2) to 0.42 J/m(2) at 1500 K). The computed pressure of bare and coated particles was found to follow the classical Laplace-Young equation.