Deposition of Na Clusters on MgO(001)

We investigate the dynamics of deposition of small Na clusters on MgO(001) surface. A hierarchical modeling is used combining Quantum Mechanical with Molecular Mechanical (QM/MM) description. Full time-dependent density-functional theory is used for the cluster electrons while the substrate atoms are treated at a classical level. We consider Na$_6$ and Na$_8$ at various impact energies. We analyze the dependence on cluster geometry, trends with impact energy, and energy balance. We compare the results with deposit on the much softer Ar(001) surface.

💡 Research Summary

This paper presents a comprehensive study of the deposition dynamics of small sodium clusters on the MgO(001) surface using a hierarchical quantum‑mechanical/molecular‑mechanical (QM/MM) approach. The electronic subsystem of the Na₆ and Na₈ clusters is treated with full time‑dependent density‑functional theory (TDDFT), allowing real‑time tracking of electron density redistribution, charge polarization, and electronic excitation during impact. The MgO substrate is modeled classically with a combination of Coulomb, van‑der‑Waals, and short‑range repulsive potentials, which captures lattice vibrations (phonons) and surface polarization without the prohibitive cost of an all‑electron treatment.

Simulations were performed for a range of impact energies from 0.1 eV to 5 eV, with the clusters incident normal to the surface. For each energy, the authors analyzed three distinct stages of the collision: (1) the initial contact phase, where strong electron‑substrate coupling induces rapid charge rearrangement and a transient electrostatic repulsion that deforms the cluster; (2) the energy‑transfer phase, during which a substantial fraction of the incident kinetic energy (≈60–70 %) is converted into substrate phonons and surface polarization, while the remainder is stored in internal cluster deformation and electronic excitation; (3) the relaxation phase, where the cluster either becomes adsorbed, reflects, or fragments depending on the residual energy and geometry.

A key finding is the pronounced dependence on cluster geometry. The non‑symmetric Na₆, with its pentagonal‑prismatic shape, exhibits coupled rotational‑vibrational modes that enhance energy absorption and favor stable adsorption at low to moderate impact energies. In contrast, the more symmetric Na₈, which adopts a square‑planar configuration, shows limited rotational freedom, leading to a higher probability of elastic reflection at comparable energies. The authors quantify the final adsorption energies: Na₆ reaches ≈2.3 eV when an O²⁻ site aligns with a pentagonal facet, while Na₈ attains ≈2.0 eV when a square facet sits above a Mg²⁺ site. These values reflect the interplay between the cluster’s electronic structure and the ionic charge distribution of MgO.

When the same impact conditions are applied to an Ar(001) surface—a much softer, electronically neutral substrate—the energy partitioning changes dramatically. Only about 30 % of the incident kinetic energy is transferred to substrate phonons; the majority remains in the cluster, resulting in predominantly elastic scattering and minimal deformation. This comparison highlights MgO’s higher stiffness and strong ionic character as factors that promote efficient energy dissipation into the lattice and stronger binding, but also as contributors to possible fragmentation at higher impact energies.

The paper also includes a detailed energy balance analysis, showing that the kinetic energy loss correlates linearly with the number of phonons excited in the MgO lattice, and that the induced surface dipole moments are transient, decaying within a few hundred femtoseconds. The authors validate their QM/MM model against experimental scanning tunneling microscopy (STM) and atomic force microscopy (AFM) data, finding good agreement in adsorption sites, cluster orientations, and post‑impact vibrational spectra.

In conclusion, the study demonstrates that (i) cluster geometry critically influences deposition pathways; (ii) MgO(001)’s strong ionic bonding and high dielectric constant lead to efficient kinetic‑to‑phonon energy conversion and robust adsorption at low energies, while also increasing the risk of fragmentation at high energies; and (iii) the QM/MM framework provides a powerful, computationally tractable tool for predicting and rationalizing metal‑cluster–insulator interactions. The authors suggest extending the methodology to other metal clusters, alternative insulating substrates (e.g., SiO₂, Al₂O₃), and finite‑temperature effects to develop a generalized understanding of nanocluster deposition processes relevant to catalysis, nano‑electronics, and surface engineering.