Reply to the 'Comment on "Perturbation theory of scattering for grazing-incidence fast-atom diffraction"' by G. A. Bocan, H. Breiss, S. Szilasi, A. Momeni, E. M. S. Casagrande, E. A. Sánchez, M. S. Gravielle and H. Khemliche, Phys. Chem. Chem. Phys., 2023, , DOI: 10.1039/D3CP02486E

In this Reply, we show that criticisms of perturbation theory for grazing-incidence fast-atom diffraction (GIFAD) are ill-founded. We show explicitly that our formulation (W. Allison, S. Miret-Artés and E. Pollak, Phys. Chem. Chem. Phys., 2022, 24, 15851) provides a similar precision in describing the observed phenomena as ab initio potentials. Since that is the main criterion to distinguish between methods, it seems reasonable to conclude that the perturbation approach using a Morse-type potential reproduces the essential aspects of the dynamics correctly. In addition we expand on the historical context and summarize the physical insights provided by our methods.

The stochastic wave function method for diffusion of alkali atoms on metallic surfaces

The stochastic wave function method is proposed to study the diffusion regimes of alkali atoms on metallic surfaces. The Lindblad approach, based on the microscopic Hamiltonian information in the Caldeira-Leggett model, is presented and numerical calculations of the dynamics are carried out to characterize surface diffusion for two different systems: Na-Cu(111) and Li-Cu(111). Calculations of the intermediate scattering function for an isolated adsorbate are compared, in the Brownian limit, with results deduced from helium spin-echo (HeSE) experiments after reducing them to single adsorbate dynamics. To illustrate the method we present the dependence on momentum transfer and the temperature dependency. Results show that the experiment can be described at a quantitative level by the 1-D quantum model (reduced dimensionality).

Atom-surface scattering in the classical multiphonon regime

Many experiments that utilize beams of incident atoms colliding with surfaces as a probe of surface properties are carried out at large energies, high temperatures and with large mass atoms. Under these conditions the scattering process does not exhibit quantum mechanical properties such as diffraction or single-phonon excitation, but rather can be treated with classical physics. This is a review of work carried out by the authors over a span of several years to develop theoretical frameworks using classical physics for describing the scattering interactions of atom with surfaces.

Surface diffusion within the Caldeira-Leggett formalism

Surface diffusion is described in terms of the intermediate scattering function in the time domain and reciprocal space. Two extreme time regimes are analyzed, ballistic (very short times) and Brownian or diffusive (very long times). This open dynamics is studied from the master equation for the reduced density matrix within the Caldeira-Leggett formalism. Several characteristic magnitudes in this decoherence process such as the coherence length, ensemble width and purity of the density matrix are analyzed. Furthermore, for flat surfaces, the surface diffusion is considered for the Schrödinger cat states and identical adsorbates or adparticles, bosons and fermions. The analytical results are compared with those issued from solving the Lindblad master equation through the stochastic wave function method. This numerical analysis is extended to be applied to corrugated surfaces.

Momentum-Space Decoherence of Distinguishable and Identical Particles in the Caldeira-Leggett Formalism

In this work, momentum-space decoherence using minimum and nonminimum-uncertainty-product (stretched) Gaussian wave packets in the framework of Caldeira-Leggett formalism and under the presence of a linear potential is studied. As a dimensionless measure of decoherence, purity, a quantity appearing in the definition of the linear entropy, is studied taking into account the role of the stretching parameter. Special emphasis is on the open dynamics of the well-known cat states and bosons and fermions compared to distinguishable particles. For the cat state, while the stretching parameter speeds up the decoherence, the external linear potential strength does not affect the decoherence time; only the interference pattern is shifted. Furthermore, the interference pattern is not observed for minimum-uncertainty-product-Gaussian wave packets in the momentum space. Concerning bosons and fermions, the question we have addressed is how the symmetry of the wave functions of indistinguishable particles is manifested in the decoherence process, which is understood here as the loss of being indistinguishable due to the gradual emergence of classical statistics with time. We have observed that the initial bunching and anti-bunching character of bosons and fermions, respectively, in the momentum space are not preserved as a function of the environmental parameters, temperature, and damping constant. However, fermionic distributions are slightly broader than the distinguishable ones and these similar to the bosonic distributions. This general behavior could be interpreted as a residual reminder of the symmetry of the wave functions in the momentum space for this open dynamics.

Ag5-induced stabilization of multiple surface polarons on perfect and reduced TiO2 rutile (110)

The recent advent of cutting-edge experimental techniques allows for a precise synthesis of subnanometer metal clusters composed of just a few atoms, opening new possibilities for subnanometer science. In this work, via first-principles modeling, we show how the decoration of perfect and reduced TiO₂ surfaces with Ag₅ atomic clusters enables the stabilization of multiple surface polarons. Moreover, we predict that Ag₅ clusters are capable of promoting defect-induced polarons transfer from the subsurface to the surface sites of reduced TiO₂ samples. For both planar and pyramidal Ag₅ clusters, and considering four different positions of bridging oxygen vacancies, we model up to 14 polaronic structures, leading to 134 polaronic states. About 71% of these configurations encompass coexisting surface polarons. The most stable states are associated with large inter-polaron distances (>7.5 Å on average), not only due to the repulsive interaction between trapped Ti³⁺ 3d¹ electrons, but also due to the interference between their corresponding electronic polarization clouds [P. López-Caballero et al., J. Mater. Chem. A 8, 6842-6853 (2020)]. As a result, the most stable ferromagnetic and anti-ferromagnetic arrangements are energetically quasi-degenerate. However, as the average inter-polarons distance decreases, most (≥70%) of the polaronic configurations become ferromagnetic. The optical excitation of the midgap polaronic states with photon energy at the end of the visible region causes the enlargement of the polaronic wave function over the surface layer. The ability of Ag₅ atomic clusters to stabilize multiple surface polarons and extend the optical response of TiO₂ surfaces toward the visible region bears importance in improving their (photo-)catalytic properties and illustrates the potential of this new generation of subnanometer-sized materials.

Surface Structure of Bi(111) from Helium Atom Scattering Measurements. Inelastic Close-Coupling Formalism

Elastic and inelastic close-coupling (CC) calculations have been used to extract information about the corrugation amplitude and the surface vibrational atomic displacement by fitting to several experimental diffraction patterns. To model the three-dimensional interaction between the He atom and the Bi(111) surface under investigation, a corrugated Morse potential has been assumed. Two different types of calculations are used to obtain theoretical diffraction intensities at three surface temperatures along the two symmetry directions. Type one consists of solving the elastic CC (eCC) and attenuating the corresponding diffraction intensities by a global Debye-Waller (DW) factor. The second one, within a unitary theory, is derived from merely solving the inelastic CC (iCC) equations, where no DW factor is necessary to include. While both methods arrive at similar predictions for the peak-to-peak corrugation value, the variance of the value obtained by the iCC method is much better. Furthermore, the more extensive calculation is better suited to model the temperature induced signal asymmetries and renders the inclusion for a second Debye temperature for the diffraction peaks futile.

Communication: Semiclassical perturbation theory for the quantum diffractive scattering of atoms on thermal surfaces

Inspired by the semiclassical perturbation theory of Hubbard and Miller [J. Chem. Phys. 80, 5827 (1984)], we derive explicit expressions for the angular distribution of particles scattered from thermal surfaces. At very low surface temperature, the observed experimental background scattering is proportional to the spectral density of the phonons. The angular distribution is a sum of diffraction peaks and a broad background reflecting the spectral density. The theory is applied to measured angular distributions of Ne, Ar, and Kr scattered from a Cu(111) surface.

Dissipative geometric phase and decoherence in parity-violating chiral molecules

Within a generalized Langevin framework for open quantum systems, the cyclic evolution of a two-level system is analyzed in terms of the geometric phase extended to dissipative systems for Ohmic friction. This proposal is applied to the dynamics of chiral molecules where the tunneling and parity violating effects are competing. The effect of different system-bath coupling functions in the dissipated energy is shown to be crucial to understand the behavior of the geometric phase as well as the decoherence displayed by the corresponding interference patterns.

Communication: Quantum Zeno-based control mechanism for molecular fragmentation

A quantum control mechanism is proposed for molecular fragmentation processes within a scenario grounded on the quantum Zeno effect. In particular, we focus on the van der Waals Ne-Br(2) complex, which displays two competing dissociation channels via vibrational and electronic predissociation. Accordingly, realistic three-dimensional wave packet simulations are carried out by using ab initio interaction potentials recently obtained to reproduce available experimental data. Two numerical models to simulate the repeated measurements are reported and analyzed. It is found that the otherwise fast vibrational predissociation is slowed down in favor of the slow electronic (double fragmentation) predissociation, which is enhanced by several orders of magnitude. Based on these theoretical predictions, some hints to experimentalists to confirm their validity are also proposed.

Quantum Zeno and anti-Zeno effects in surface diffusion of interacting adsorbates

Surface diffusion of interacting adsorbates is here analyzed within the context of two fundamental phenomena of quantum dynamics, namely the quantum Zeno effect and the anti-Zeno effect. The physical implications of these effects are introduced here in a rather simple and general manner within the framework of non-selective measurements and for two (surface) temperature regimes: high and very low (including zero temperature). The quantum intermediate scattering function describing the adsorbate diffusion process is then evaluated for flat surfaces, since it is fully analytical in this case. Finally, a generalization to corrugated surfaces is also discussed. In this regard, it is found that, considering a Markovian framework and high surface temperatures, the anti-Zeno effect has already been observed, though not recognized as such.

A theoretical study of rotational and translational diffusion dynamics of molecules with a six-fold point symmetry adsorbed on a hexagonal lattice by neutron scattering

A complete analytical model for the rotational and translational diffusion of molecules with a six-fold point symmetry on a hexagonal lattice is presented. It can be applied, in particular, to the diffusion of benzene molecules adsorbed flat on the basal plane of graphite in the case of incoherent scattering. Under the weak hindered approximation, the classical mechanics framework and making use of the van Hove formalism of correlation functions, the intermediate scattering function and its Fourier transform, the scattering law, are both obtained. They can be expressed as sums of exponential decays or Lorentzian functions, respectively, containing the contribution of each of the dynamical processes taking place. In the case of benzene lying flat on the substrate we expect translational diffusion, continuous rotations of isolated molecules and hindered rotations of molecules within clusters. Each particular diffusive mechanism can be recognized owing to its particular signature in the dependence of the quasi-elastic broadening on the momentum transfer.

He atom scattering from ZnO surfaces: calculation of diffraction peak intensities using the close-coupling approach

Diffraction intensities of a molecular He beam scattered off the clean and water-covered ZnO(1010) surface have been simulated using a new potential model in conjunction with the close-coupling formalism. The effective corrugation functions for the systems He-H2O/ZnO(1010) and He-H2O/ZnO(1010) have been obtained from density functional theory calculations within the Esbjerg-Nørskov approximation. Using these data a potential model is constructed consisting of a corrugated Morse potential at small He-surface distances and a semiempiric attractive part at larger distances. The diffraction patterns obtained from close-coupling calculations agree with the experimental data within about 10%, which opens the possibility to simulate He diffraction from surfaces of any structural complexity and to verify surface and adsorbate structures proposed theoretically by employing this kind of analysis.

Phonon lineshapes in atom-surface scattering

Phonon lineshapes in atom-surface scattering are obtained from a simple stochastic model based on the so-called Caldeira-Leggett Hamiltonian. In this single-bath model, the excited phonon resulting from a creation or annihilation event is coupled to a thermal bath consisting of an infinite number of harmonic oscillators, namely the bath phonons. The diagonalization of the corresponding Hamiltonian leads to a renormalization of the phonon frequencies in terms of the phonon friction or damping coefficient. Moreover, when there are adsorbates on the surface, this single-bath model can be extended to a two-bath model accounting for the effect induced by the adsorbates on the phonon lineshapes as well as their corresponding lineshapes.

Two-bath model for activated surface diffusion of interacting adsorbates

The diffusion and low vibrational motions of adsorbates on surfaces can be well described by a purely stochastic model, the so-called interacting single adsorbate model, for low-moderate coverages (theta < or approximately equal to 0.12). Within this model, the effects of thermal surface phonons and adsorbate-adsorbate collisions are accounted for by two uncorrelated noise functions, which arise in a natural way from a two-bath model based on a generalization of the one-bath Caldeira-Leggett Hamiltonian. As an illustration, the model is applied to the diffusion of Na atoms on a Cu(001) surface with different coverages.

A causal look into the quantum Talbot effect

A well-known phenomenon in both optics and quantum mechanics is the so-called Talbot effect. This near field interference effect arises when infinitely periodic diffracting structures or gratings are illuminated by highly coherent light or particle beams. Typical diffraction patterns known as quantum carpets are then observed. Here the authors provide an insightful picture of this nonlocal phenomenon as well as its classical limit in terms of Bohmian mechanics, also showing the causal reasons and conditions that explain its appearance. As an illustration, theoretical results obtained from diffraction of thermal He atoms by both N-slit arrays and weak corrugated surfaces are analyzed and discussed. Moreover, the authors also explain in terms of what they call the Talbot-Beeby effect how realistic interaction potentials induce shifts and distortions in the corresponding quantum carpets.

Line shape broadening in surface diffusion of interacting adsorbates with quasielastic he atom scattering

The experimental line shape broadening observed in adsorbate diffusion on metal surfaces with increasing coverage is usually related to the nature of the adsorbate-adsorbate interaction. Here we show that this broadening can also be understood in terms of a fully stochastic model just considering two noise sources: (i) a Gaussian white noise accounting for the surface friction, and (ii) a shot noise replacing the physical adsorbate-adsorbate interaction potential. Furthermore, contrary to what could be expected, for relatively weak adsorbate-substrate interactions the opposite effect is predicted: line shapes get narrower with increasing coverage.

Generalized Chudley-Elliott vibration-jump model in activated atom surface diffusion

Here the authors provide a generalized Chudley-Elliott expression for the activated atom surface diffusion which takes into account the coupling between both low-frequency vibrational motion (namely, the frustrated translational modes) and diffusion. This expression is derived within the Gaussian approximation framework for the intermediate scattering function at low coverage. Moreover, inelastic contributions (arising from creation and annihilation processes) to the full width at half maximum of the quasielastic peak are also obtained.

Surface diffusion and low vibrational motion with interacting adsorbates: a shot noise description

Here, an approach in terms of shot noise is proposed to study and characterize surface diffusion and low vibrational motion when having interacting adsorbates on surfaces. In what we call statistical limit, that is, at long times and high number of collisions, one expects that diffusing particles display an essential Markovian behavior. Accordingly, the action of the pairwise potentials accounting for particle-particle collisions is equivalent to considering a shot noise acting on a single particle. We call this approach the interacting single adsorbate approximation, which gathers three important advantages: (i) the dynamics underlying surface diffusion and low vibrational motion can be easily understood in terms of relatively simple stochastic processes; (ii) from our model, appropriate (and well justified) working formulas are easily obtained, which explain the results arising from more complicated (but commonly used) molecular dynamics simulations within the Langevin formulation; and (iii), at the same time, it is less demanding computationally than the latter type of calculations. In order to illustrate the application of this model, numerical results are presented. Specially, our model reproduces the experimental observation regarding the broadening of the quasielastic peak ruling surface diffusion.

Quasi-elastic peak lineshapes in adsorbate diffusion on nearly flat surfaces at low coverages: the motional narrowing effect in Xe on Pt(111)

Quasi-elastic helium atom scattering measurements have provided clear evidence for a two-dimensional free gas of Xe atoms on Pt(111) at low coverages. Increasing the friction due to the surface, a gradual change of the shape of the quasi-elastic peak is predicted and analysed for this system in terms of the so-called motional narrowing effect. The type of analysis presented here for the quasi-elastic peak should be prior to any deconvolution procedure carried out in order to better extract information from the process, e.g. diffusion coefficients and jump distributions. Moreover, this analysis also provides conditions for the free gas regime different to those reported earlier.

Quantum trajectories in elastic atom-surface scattering: threshold and selective adsorption resonances

The elastic resonant scattering of He atoms off the Cu(117) surface is fully described with the formalism of quantum trajectories provided by Bohmian mechanics. Within this theory of quantum motion, the concept of trapping is widely studied and discussed. Classically, atoms undergo impulsive collisions with the surface, and then the trapped motion takes place covering at least two consecutive unit cells. However, from a Bohmian viewpoint, atom trajectories can smoothly adjust to the equipotential energy surface profile in a sort of sliding motion; thus the trapping process could eventually occur within one single unit cell. In particular, both threshold and selective adsorption resonances are explained by means of this quantum trapping considering different space and time scales. Furthermore, a mapping between each region of the (initial) incoming plane wave and the different parts of the diffraction and resonance patterns can be easily established, an important issue only provided by a quantum trajectory formalism.

Quantum trajectories in atom-surface scattering with single adsorbates: the role of quantum vortices

In this work, a full quantum study of the scattering of He atoms off single CO molecules, adsorbed onto the Pt(111) surface, is presented within the formalism of quantum trajectories provided by Bohmian mechanics. By means of this theory, it is shown that the underlying dynamics is strongly dominated by the existence of a transient vortitial trapping with measurable effects on the whole diffraction pattern. This kind of trapping emphasizes the key role played by quantum vortices in this scattering. Moreover, an analysis of the surface rainbow effect caused by the local corrugation that the CO molecule induces on the surface, and its manifestation in the corresponding intensity pattern, is also presented and discussed.

Bound-state energies in argon trimers via a variational expansion: The effects from many-body corrections

In this paper we study the bound-state energies and geometries of Ar(3) for J=0, using the distributed Gaussian functions method that provides a configurational description of the different structures contributing to these states. Atom-atom potentials are employed and three-body long-range effects are also included in the computational treatment by adding to the sum of potentials the Axilrod-Teller triple-dipole correction for the whole rotationless energy spectrum. An estimate of the total number of bound states for the Ar trimer is given. With respect to previous calculations, limited to the lower-lying states, our results show slightly larger nonadditive effects and are further able to predict the full range of the bound spectrum. Changes on the geometries of a large part of the vibrationally excited states of Ar(3) when the Axilrod-Teller term is included in the molecular potential are found by the present study.

A complete configurational study for the bound states of Ne trimers

The structural properties and the energetics of the ground and the excited bound states of Ne(3) for zero total angular momentum are examined using different modelings for the two-body interactions. We employ a method consisting of a variational approach with a distributed Gaussian functions (DGF) basis set expansion. We discuss at length the advantages and possible limitations of such an approach, comparing it to other methods which have been applied in the literature to the same system. The DGF method turns out to be very accurate in giving us the bound states energetics and also provides in a natural way a convincing pictorial description of all the states, including those with dominant linear configurations. Additional bound states are found for the Ne(3) system with respect to those indicated in previous works and we suggest a "stabilization" procedure that can be used to assess the truly bound nature of a state. Some considerations on the relative reliability of the examined two-body interactions are also reported.

Collisional line shapes for low frequency vibrations of adsorbates on a metal surface

The dynamics of atoms or molecules adsorbed on a metal surface, and excited by collisions with an atomic beam, are treated within a theory that includes energy dissipation into lattice vibrations by means of a frequency and temperature dependent friction function. The theory provides dynamic structure factors for energy transfer derived from collisional time correlation functions. It describes the relaxation of a vibrationally excited atom or molecule within a model of a damped quantum harmonic oscillator bilinearly coupled to a bath of lattice oscillators. The collisional time correlation function is generalized to include friction effects and is applied to the vibrational relaxation of the frustrated translation mode of Na adsorbed on a Cu(001) surface, CO on Cu(001), and CO on Pt(111), following excitation by collisions with He atoms. Results for the frequency shift and width of line shapes versus surface temperature are in very good agreement with experimental measurements of inelastic He atom scattering. Our interpretation of the experimental results provides insight on the relative role of phonon versus electron-hole relaxation.

Chaotic transport of particles in two-dimensional periodic potentials driven by ac forces

The diffusive and directed transport of particles in a two-dimensional periodic potential subjected to frictional and time-periodic forces is analyzed in detail. The model represents diffusion of atoms adsorbed on metal surfaces under an applied ac electric field (surface electromigration) in the low-temperature limit. The second dimension and the potential energy coupling are shown to play an important role on both diffusion and net currents, depending on the direction of the drive. A properly chosen biharmonic field is able to control the directed ratchetlike dynamics of atoms on symmetric surfaces, since current reversals take place by different stabilization of attractors. Reversals identified with hysteresis loops between periodic running attractors are robust against an increase of the second harmonic amplitude, and against temperature effects inside the experimental range for measurements of surface diffusion.