Performance of correlation functionals in ab initio chemisorption cluster-model calculations: Alkali metals on Si(111)

The performance of different correlation functionals has been tested for alkali metals, Li to Cs, interacting with cluster models simulating different active sites of the Si(111) surface. In all cases, the ab initio Hartree-Fock density has been obtained and used as a starting point. The electronic correlation energy is then introduced as an a posteriori correction to the Hartree-Fock energy using different correlation functionals. By making use of the ionic nature of the interaction and of different dissociation limits we have been able to prove that all functionals tested introduce the right correlation energy, although to a different extent. Hence, correlation functionals appear as an effective and easy way to introduce electronic correlation in the ab initio Hartree-Fock description of the chemisorption bond in complex systems where conventional configuration interaction techniques cannot be used. However, the calculated energies may differ by some tens of eV. Therefore, these methods can be employed to get a qualitative idea of how important correlation effects are, but they have some limitations if accurate binding energies are to be obtained.

From nonwetting to prewetting: The asymptotic behavior of 4He drops on alkali substrates

We investigate the spreading of 4He droplets on alkali-metal surfaces at zero temperature, within the frame of finite range density-functional theory. The equilibrium configurations of several 4HeN clusters and their asymptotic trend with increasing particle number N, which can be traced to the wetting behavior of the quantum fluid, are examined for nanoscopic droplets. We discuss the size effects inferring that the asymptotic properties of large droplets correspond to those of the prewetting film.

Condensation of helium in nanoscopic alkali wedges at zero temperature

We present a complete calculation of the structure of liquid 4He confined to a concave nanoscopic wedge, as a function of the opening angle of the walls. This is achieved within a finite-range density functional formalism. The results here presented, restricted to alkali metal substrates, illustrate the change in meniscus shape from rather broad to narrow wedges on weak and strong alkali adsorbers, and we relate this change to the wetting behavior of helium on the corresponding planar substrate. As the wedge angle is varied, we find a sequence of stable states that, in the case of cesium, undergo one filling and one emptying transition at large and small openings, respectively. A computationally unambiguous criterion to determine the contact angle of 4He on cesium is also proposed.

Comment on "Reappraisal of experimental values of third-order elastic constants of some cubic semiconductors and metals"

In a recent paper A. S. Johal and D. J. Dunstan [Phys. Rev. B 73, 024106 (2006)] have applied multivariate linear regression analysis to the published data of the change in ultrasonic velocity with applied stress. The aim is to obtain the best estimates for the third-order elastic constants in cubic materials. From such an analysis they conclude that uniaxial stress data on metals turns out to be nearly useless by itself. The purpose of this comment is to point out that by a proper analysis of uniaxial stress data it is possible to obtain reliable values of third-order elastic constants in cubic metals and alloys. Cu-based shape memory alloys are used as an illustrative example.

Adsorption energy and spin state of first-row transition metals adsorbed on MgO(100)

Slab and cluster model spin-polarized calculations have been carried out to study various properties of isolated first-row transition metal atoms adsorbed on the anionic sites of the regular MgO(100) surface. The calculated adsorption energies follow the trend of the metal cohesive energies, indicating that the changes in the metal-support and metal-metal interactions along the series are dominated by atomic properties. In all cases, except for Ni at the generalized gradient approximation level, the number of unpaired electron is maintained as in the isolated metal atom. The energy required to change the atomic state from high to low spin has been computed using the PW91 and B3LYP density-functional-theory-based methods. PW91 fails to predict the proper ground state of V and Ni, but the results for the isolated and adsorbed atom are consistent within the method. B3LYP properly predicts the ground state of all first-row transition atom the high- to low-spin transition considered is comparable to experiment. In all cases, the interaction with the surface results in a reduced high- to low-spin transition energy.