Subatomic-scale resolution with SPM: Co adatom on p(2 × 1)Cu(110):O.

We study the limits of SPM subatomic resolution in imaging orbital magnetic features on a model system of a Co atom on a p(2 × 1)Cu(110):O surface. We show that scanning tunneling spectroscopies allow the determination of the occupation of the Co d shells and the value of the Hubbard U in the DFT + U modeling, and that standard near-contact AFM can in principle image the asymmetry due to partial filling of the d shells at close distances in the small-amplitude regime. Due to the partially ionic character of Co, a faint asymmetry is predicted to also arise in the electrostatic force. We anticipate these features to be even stronger for a transition metal adsorbate featuring larger departures from sphericality in charge/spin densities.

Optical Gaps in Pristine and Heavily Doped Silicon Nanocrystals: DFT versus Quantum Monte Carlo Benchmarks.

We present a time-dependent density functional theory (TDDFT) study of the optical gaps of light-emitting nanomaterials, namely, pristine and heavily B- and P-codoped silicon crystalline nanoparticles. Twenty DFT exchange-correlation functionals sampled from the best currently available inventory such as hybrids and range-separated hybrids are benchmarked against ultra-accurate quantum Monte Carlo results on small model Si nanocrystals. Overall, the range-separated hybrids are found to perform best. The quality of the DFT gaps is correlated with the deviation from Koopmans' theorem as a possible quality guide. In addition to providing a generic test of the ability of TDDFT to describe optical properties of silicon crystalline nanoparticles, the results also open up a route to benchmark-quality DFT studies of nanoparticle sizes approaching those studied experimentally.

Promoting Atoms into Delocalized Long-Living Magnetically Modified State Using Atomic Force Microscopy.

We report on a low-temperature atomic force microscropy manipulation of Co atoms in ultrahigh vacuum on an oxidized copper surface in which the manipulated atom is kept delocalized above several surface unit cells over macroscopic times. The manipulation employed, in addition to the ubiquitous short-range tip-generated chemical forces, also long-range forces generated via Friedel oscillations of the metal charge density due to Co nanostructures prearranged on the surface by lateral manipulation. We show that our manipulation protocol requires mechanical control of the spin state of the Co atom.

Charged vanadium-benzene multidecker clusters: DFT and quantum Monte Carlo study.

Using explicitly correlated fixed-node quantum Monte Carlo and density functional theory (DFT) methods, we study electronic properties, ground-state multiplets, ionization potentials, electron affinities, and low-energy fragmentation channels of charged half-sandwich and multidecker vanadium-benzene systems with up to 3 vanadium atoms, including both anions and cations. It is shown that, particularly in anions, electronic correlations play a crucial role; these effects are not systematically captured with any commonly used DFT functionals such as gradient corrected, hybrids, and range-separated hybrids. On the other hand, tightly bound cations can be described qualitatively by DFT. A comparison of DFT and quantum Monte Carlo provides an in-depth understanding of the electronic structure and properties of these correlated systems. The calculations also serve as a benchmark study of 3d molecular anions that require a balanced many-body description of correlations at both short- and long-range distances.

Atomic force microscopy identification of Al-sites on ultrathin aluminum oxide film on NiAl(110).

Ultrathin alumina film formed by oxidation of NiAl(110) was studied by non-contact atomic force microscopy in an ultra high vacuum at room temperature with the quest to provide the ultimate understanding of structure and bonding of this complicated interface. Using a very stiff Si cantilever with significantly improved resolution, we have obtained images of this system with unprecedented resolution, surpassing all the previous results. In particular, we were able to unambiguously resolve all the differently coordinated aluminum atoms. This is of importance as the previous images provide very different image patterns, which cannot easily be reconciled with the existing structural models. Experiments are supported by extensive density functional theory modeling. We find that the system is strongly ionic and the atomic force microscopy images can reliably be understood from the electrostatic potential which provides an image model in excellent agreement with the experiments. However, in order to resolve the finer contrast features we have proposed a more sophisticated model based on more realistic approximants to the incommensurable alumina interface.

Vertical atomic manipulation with dynamic atomic-force microscopy without tip change via a multi-step mechanism.

Manipulation is the most exciting feature of the non-contact atomic force microscopy technique as it allows building nanostructures on surfaces. Usually vertical manipulations are accompanied by an abrupt tip modification leading to a change of contrast. Here we report on low-temperature experiments demonstrating vertical manipulations of 'super'-Cu atoms on the p(2 × 1) Cu(110):O surface, both extractions to and depositions from the tip, when the imaging contrast remains the same. These results are rationalized employing a novel and completely general method that combines density functional theory calculations for obtaining energy barriers as a function of tip height and a Kinetic Monte Carlo algorithm for studying the tip dynamics and extraction of manipulation statistics. The model reveals a novel multi-step manipulation mechanism combining activated jumps of 'super'-Cu atoms to/from the tip with their drag by and diffusion on the tip.

Quantum Monte Carlo Study of π-Bonded Transition Metal Organometallics: Neutral and Cationic Vanadium-Benzene and Cobalt-Benzene Half Sandwiches.

We present accurate quantum Monte Carlo (QMC) calculations that enabled us to determine the structure, spin multiplicity, ionization energy, dissociation energy, and spin-dependent electronic gaps of neutral and positively charged vanadium-benzene and cobalt-benzene systems. From total/ionization energy, we deduce a sextet (quintet) state of neutral (cationic) vanadium-benzene systems and quartet (triplet) state of the neutral (cationic) cobalt-benzene systems. Vastly different energy gaps for the two spin channels are predicted for the vanadium-benzene system and broadly similar energy gaps for the cobalt-benzene system. For this purpose, we have used a multistage combination of techniques with consecutive elimination of systematic biases except for the fixed-node approximation in QMC. Our results significantly differ from the established picture based on previous less accurate calculations and point out the importance of high-level many-body methods for predictive calculations of similar transition metal-based organometallic systems.

Critical Importance of van der Waals Stabilization in Strongly Chemically Bonded Surfaces: Cu(110):O.

We provide strong evidence that different reconstructed phases of the oxidized Cu(110) surface are stabilized by the van der Waals (vdW) interactions. These covalently bonded reconstructed surfaces feature templates that are an integral part of the surfaces and are bonded on the bare metal surface by a combination of chemical and physical bonding. The vdW stabilization in this class of systems affects predominantly the intertemplate Cu-O interactions in structures sparsely populated by these templates. The conventional dispersionless density functional theory (DFT) methods fail to model such systems. We find a failure to describe the thermodynamics of the different phases that are formed at different oxygen exposures and spurious minima on the potential energy surface of a diffusing surface adatom. To overcome these issues, we employ a range of different DFT methods that account for the missing vdW correlations. Surprisingly, despite vast conceptual differences in the different formulations of these methods, they yield physically identical results for the Cu(110):O surface phases, provided the massive screening effects in the metal are taken into account. Contrary, the vibrational contribution does not consistently stabilize the experimentally observed surface structures. The van der Waals surface stabilization, so far deemed to play only a minor role in hard-bonded surfaces, is suggested to be a more general key feature for this and other related surfaces.

Spin multiplicity and symmetry breaking in vanadium-benzene complexes.

We present accurate quantum Monte Carlo (QMC) calculations which enabled us to determine the structure, spin multiplicity, ionization energy, dissociation energy, and spin-dependent electronic gaps of the vanadium-benzene system. From total and ionization energy we deduce a high-spin state with vastly different energy gaps for the two spin channels. For this purpose we have used a multistage combination of techniques with consecutive elimination of systematic biases except for the fixed-node approximation in QMC calculations. Our results significantly differ from the established picture based on previous less accurate calculations and point out the importance of high-level many-body methods for predictive calculations of similar transition metal-based organometallic systems.

Disentanglement of triplet and singlet states of azobenzene: direct EELS detection and QMC modeling.

Singlet and triplet excited states of trans-azobenzene have been measured in the gas phase by electron energy loss spectroscopy (EELS). In order to interpret the strongly overlapping singlet and triplet bands in the spectra a set of large-scale correlated quantum Monte-Carlo (QMC) simulations was performed. The EELS/QMC combination of methods yields an excellent agreement between theory and experiment and for the two low-lying excited singlet and two low-lying triplet states permitted their unambiguous assignment. In addition, EELS revealed two overlapping electronic states in the band commonly assigned as S(2), the lower one with a pronounced vibrational structure, the upper one structureless. Finally, the agreement between theory and experiment was shown to further increase by taking computationally into account the finite temperature effects.

Ground and excited electronic states of azobenzene: a quantum Monte Carlo study.

Large-scale quantum Monte Carlo (QMC) calculations of ground and excited singlet states of both conformers of azobenzene are presented. Remarkable accuracy is achieved by combining medium accuracy quantum chemistry methods with QMC. The results not only reproduce measured values with chemical accuracy but the accuracy is sufficient to identify part of experimental results which appear to be biased. Novel analysis of nodal surface structure yields new insights and control over their convergence, providing boost to the chemical accuracy electronic structure methods of large molecular systems.

Thermodynamics of catalytic formation of dimethyl ether from methanol in acidic zeolites.

We present a theoretical study of the formation of the first intermediate, dimethyl ether, in the methanol to gasoline conversion within the framework of an ab initio molecular dynamics approach. The study is performed under conditions that closely resemble the reaction conditions in the zeolite catalyst including the full topology of the framework. The use of the method of thermodynamic integration allows us to extract the free-energy profile along the reaction coordinate. We find that the entropic contribution qualitatively alters the free-energy profile relative to the total energy profile. Different transition states are found from the internal and free energy profiles. The entropy contribution varies significantly along the reaction coordinate and is responsible for stabilizing the products and for lowering the energy barrier. The hugely inhomogeneous variation of the entropy can be understood in terms of elementary processes that take place during the chemical reaction. Our simulations provide new insights into the complex nature of this chemical reaction.

Silicon clusters of intermediate size: energetics, dynamics, and thermal effects

A combination of the best available theoretical techniques for energetics, dynamics, and thermodynamics is employed in an extensive study of Si(n) ( n = 20,25) clusters. For T = 0 we solve the electronic structure by the density functional and the highly accurate quantum Monte Carlo approaches. Finite temperature and dynamical effects are investigated by the ab initio molecular dynamics method. This combination of methods enables us to find several new low-energy isomers and to explain the differences in properties, behavior, and stability of elongated versus compact types of structures and to elucidate the origin of the existing discrepancies between theory and experiments.