Néstor F. Aguirre

Communication: Unraveling the 4He droplet-mediated soft-landing from ab initio-assisted and time-resolved density functional simulations: Au@4He300/TiO2(110)

An ab-initio-based methodological scheme for He-surface interactions and zero-temperature time-dependent density functional theory for superfluid (4)He droplets motion are combined to follow the short-time collision dynamics of the Au@(4)He300 system with the TiO2(110) surface. This composite approach demonstrates the (4)He droplet-assisted sticking of the metal species to the surface at low landing energy (below 0.15 eV/atom), thus providing the first theoretical evidence of the experimentally observed (4)He droplet-mediated soft-landing deposition of metal nanoparticles on solid surfaces [Mozhayskiy et al., J. Chem. Phys. 127, 094701 (2007) and Loginov et al., J. Phys. Chem. A 115, 7199 (2011)].

Helium mediated deposition: Modeling the He−TiO2(110)-(1×1) interaction potential and application to the collision of a helium droplet from density functional calculations

This paper is the first of a two-part series dealing with quantum-mechanical (density-functional-based) studies of helium-mediated deposition of catalytic species on the rutile TiO(2)(110)-(1×1) surface. The interaction of helium with the TiO(2)(110)-(1×1) surface is first evaluated using the Perdew-Burke-Ernzerhof functional at a numerical grid dense enough to build an analytical three-dimensional potential energy surface. Three (two prototype) potential models for the He-surface interaction in helium scattering calculations are analyzed to build the analytical potential energy surface: (1) the hard-corrugated-wall potential model; (2) the corrugated-Morse potential model; and (3) the three-dimensional Morse potential model. Different model potentials are then used to study the dynamics upon collision of a (4)He(300) cluster with the TiO(2)(110) surface at zero temperature within the framework of a time-dependent density-functional approach for the quantum fluid [D. Mateo, D. Jin, M. Barranco, and M. Pi, J. Chem. Phys. 134, 044507 (2011)] and classical dynamics calculations. The laterally averaged density functional theory-based potential with an added long-range dispersion interaction term is further applied. At variance with classical dynamics calculations, showing helium droplet splashing out of the surface at impact, the time evolution of the macroscopic helium wave-function predicts that the helium droplet spreads on the rutile surface and leads to the formation of a thin film above the substrate. This work thus provides a basis for simulating helium mediated deposition of metallic clusters embedded within helium nanodroplets.

Including nuclear quantum effects into highly correlated electronic structure calculations of weakly bound systems

An interface between the APMO code and the electronic structure package MOLPRO is presented. The any particle molecular orbital APMO code [González et al., Int. J. Quantum Chem. 108, 1742 (2008)] implements the model where electrons and light nuclei are treated simultaneously at Hartree-Fock or second-order Möller-Plesset levels of theory. The APMO-MOLPRO interface allows to include high-level electronic correlation as implemented in the MOLPRO package and to describe nuclear quantum effects at Hartree-Fock level of theory with the APMO code. Different model systems illustrate the implementation: (4)He2 dimer as a protype of a weakly bound van der Waals system; isotopomers of [He-H-He](+) molecule as an example of a hydrogen bonded system; and molecular hydrogen to compare with very accurate non-Born-Oppenheimer calculations. The possible improvements and future developments are outlined.

Quantum solvent states and rovibrational spectra of small doped H3e clusters through the full-configuration-interaction nuclear orbital approach: The (H3e)N–Cl2(X) case (N≤4)

A full-configuration-interaction nuclear orbital treatment has been recently developed as a benchmark quantum-chemistry-like method to study small doped (3)He clusters [M. P. de Lara-Castells et al., J. Chem. Phys. 125, 221101 (2006)]. Our objective in this paper is to extend our previous study on ((3)He)(N)-Cl(2)(B) clusters, using an enhanced implementation that allows employing very large one-particle basis sets [M. P. de Lara-Castells et al., J. Chem. Phys. 131, 194101 (2009)], and apply the method to the ((3)He)(N)-Cl(2)(X) case, using both a semiempirical T-shaped and an ab initio He-dopant potential with minima at both T-shaped and linear conformations. Calculations of the ground and low-lying excited solvent states stress the key role played by the anisotropy of the He-dopant interaction in determining the global energies and the structuring of the (3)He atoms around the dopant. Whereas (3)He atoms are localized in a broad belt around the molecular axis in ground-state N-sized complexes with N=1-3, irrespective of using the T-shaped or the ab initio He-dopant potential function, the dopant species becomes fully coated by just four (3)He atoms when the He-dopant potential also has a minimum at linear configurations. However, excited solvent states with a central ring-type clustering of the host molecule are found to be very close in energy with the ground state by using the ab initio potential function. A microscopic analysis of this behavior is provided. Additional simulations of the molecular rovibrational Raman spectra, also including excited solvent states, provide further insights into the importance of proper modeling the anisotropy of the He-dopant interaction in these weakly bound systems and of taking into account the low-lying excitations.