Decoherent phonon effects in fast atom-surface scattering

L. Frisco,M. S. Gravielle
DOI: https://doi.org/10.1016/j.nimb.2023.03.039
2023-03-15
Abstract:Understanding the influence of phonon-mediated processes on grazing-incidence fast atom diffraction (GIFAD) patterns is relevant for its use as a surface analysis technique. In this work, we apply the Phonon-Surface Initial Value Representation (P0-SIVR) approximation to study lattice vibration effects on GIFAD patterns for the He-LiF(001) system at room temperature. The P0-SIVR approach is a semiquantum method that describes the zero-phonon atom-surface scattering including the contribution of intermediate phonon transitions. The main features introduced by thermal lattice vibrations in the angular distributions of the scattered helium atoms are investigated by considering normal energies in the 0.1 -3 eV range. In all the cases, thermal fluctuations introduce a wide polar spread that transforms the interference maxima into elongated strips. We found that the polar width of these fringes does not depend on the normal energy, as it was experimentally observed. In addition, when the normal energy increases, not only the relative intensities of interference peaks are affected by the crystal vibrations, but also the visibility of the interference structures, which disappear completely for normal energies approximately equal to or higher than 3 eV. These findings agree fairly well with the experimental data reported in Phys. Chem. Chem. Phys. 237615 (2021). However, it should be noticed that the simulated polar widths underestimate the experimentally-derived limit, suggesting that there are other mechanisms, such as inelastic phonon processes, that contribute to the polar dispersion of the GIFAD patterns.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the influence of lattice vibrations (i.e., phonon effects) on Grazing - Incidence Fast Atom Diffraction (GIFAD) patterns, especially under different normal - incidence energies. Specifically, the author uses the He/LiF(001) system model at room temperature and the Phonon - Surface Initial Value Representation (P - SIVR) approximation method to explore how thermal lattice vibrations change the scattering - angle distribution of helium atoms. ### Main problems and findings 1. **Broad polar - angle dispersion**: - The paper points out that when thermal lattice vibrations are considered, the interference maxima will become strips extending in the polar - angle direction, rather than the thin rings (Laue circles) produced by the ideally stationary crystal surface. This broad polar - angle dispersion phenomenon has also been observed in experiments. 2. **Dependence on normal energy**: - Research shows that the polar - angle width does not depend on the normal energy, which is consistent with experimental observations. However, as the normal energy increases, the visibility of the interference structure gradually decreases and completely disappears at about 3 eV. This indicates that at high normal energies, thermal vibrations have a significant decoherence effect on the scattering pattern. 3. **Comparison with experimental data**: - The simulated polar - angle width underestimates the experimental results, suggesting that there may be other mechanisms (such as inelastic phonon processes) causing polar - angle dispersion. The paper also mentions that the current simulation results differ from some experimental data, probably because the used atom - surface interaction potential model is not accurate enough. ### Formula summary The main formulas involved in the paper include: - **P0 - SIVR scattering probability**: \[ \frac{dP(P0 - SIVR)}{d\Omega_f}=K_f^2\left|A(P0 - SIVR)\right|^2 \] where \(K_f = K_i\), and \(\Omega_f = (\theta_f,\phi_f)\) represents the final scattering direction. - **Effective P0 - SIVR amplitude**: \[ A(P0 - SIVR)=\int dR_o f(R_o)\int dK_o g(K_o)\int du_o a_0(R_o,K_o,u_o) \] - **Log - normal function fitting**: \[ P(\theta_f)=\frac{A}{\omega\theta_f}\exp\left[-\frac{2[\ln(\theta_f / \theta_c)]^2}{\omega^2}\right] \] These formulas are used to describe the scattering probability and angle distribution and to verify the validity of the theoretical model through numerical calculations.