noloco: An efficient implementation of van der Waals density functionals based on a Monte-Carlo integration technique.

The treatment of van der Waals interactions in density functional theory is an important field of ongoing research. Among different approaches developed recently to capture these non-local interactions, the van der Waals density functional (vdW-DF) developed in the groups of Langreth and Lundqvist is becoming increasingly popular. It does not rely on empirical parameters, and has been successfully applied to molecules, surface systems, and weakly-bound solids. As the vdW-DF requires the evaluation of a six-dimensional integral, it scales, however, unfavorably with system size. In this work, we present a numerically efficient implementation based on the Monte-Carlo technique for multi-dimensional integration. It can handle different versions of vdW-DF. Applications range from simple dimers to complex structures such as molecular crystals and organic molecules physisorbed on metal surfaces.

Electronic interactions between "pea" and "pod": the case of oligothiophenes encapsulated in carbon nanotubes.

One of the most challenging strategies to achieve tunable nanophotonic devices is to build robust nanohybrids with variable emission in the visible spectral range, while keeping the merits of pristine single-walled carbon nanotubes (SWNTs). This goal is realized by filling SWNTs ("pods") with a series of oligothiophene molecules ("peas"). The physical properties of these peapods are depicted by using aberration-corrected high-resolution transmission electron microscopy, Raman spectroscopy, and other optical methods including steady-state and time-resolved measurements. Visible photoluminescence with quantum yields up to 30% is observed for all the hybrids. The underlying electronic structure is investigated by density functional theory calculations for a series of peapods with different molecular lengths and tube diameters, which demonstrate that van der Waals interactions are the bonding mechanism between the encapsulated molecule and the tube.

Epitaxy of rodlike organic molecules on sheet silicates--a growth model based on experiments and simulations.

During the last years, self-assembled organic nanostructures have been recognized as a proper fundament for several electrical and optical applications. In particular, phenylenes deposited on muscovite mica have turned out to be an outstanding material combination. They tend to align parallel to each other forming needlelike structures. In that way, they provide the key for macroscopic highly polarized emission, waveguiding, and lasing. The resulting anisotropy has been interpreted so far by an induced dipole originating from the muscovite mica substrate. Based on a combined experimental and theoretical approach, we present an alternative growth model being able to explain molecular adsorption on sheet silicates in terms of molecule-surface interactions only. By a comprehensive comparison between experiments and simulations, we demonstrate that geometrical changes in the substrate surface or molecule lead to different molecular adsorption geometries and needle directions which can be predicted by our growth model.

A momentum space view of the surface chemical bond.

Well-ordered and oriented monolayers of conjugated organic molecules can offer new perspectives on surface bonding. We will demonstrate the importance of the momentum distribution, or symmetry, of the adsorbate molecules' π orbitals in relation to the states available for hybridization at the metal surface. Here, the electronic band structure of the first monolayer of sexiphenyl on Cu(110) has been examined in detail with angle-resolved ultraviolet photoemission spectroscopy over a large momentum range and will be compared to measurements of a multilayer thin film and to density functional calculations. In the monolayer, the one-dimensional intramolecular band structure can still be recognized, allowing an accurate determination of orbital modification upon bonding and the relative energetic positions of the electronic levels. It is seen that the character of the molecular π orbitals is largely maintained despite strong mixing between Cu and molecular states and that the lowest unoccupied molecular orbital (LUMO) is filled by hybridization with Cu s,p states rather than through a charge transfer process. It is also shown that the momentum distribution of the substrate states involved and the periodicity of the molecular overlayer play a large role in the final E(k) distribution of the hybrid states. The distinct momentum distribution of the LUMO, interacting with the Cu substrate s,p valence bands around the gap in the surface projection of the bulk band structure, make this system a particularly illustrative example of momentum resolved hybridization. This system demonstrates that, for hybridization to occur, not only do states require overlap in energy and space, but also in momentum.

Effect of rhenium on the dislocation core structure in tungsten.

Despite exhibiting the highest melting point of all metals, the technological use of tungsten is hampered by its room-temperature brittleness. Alloying with Re significantly ductilizes the material which has been assigned to modified properties of the 1/2(111) screw dislocation. Using density functional theory, we show that alloying induces a transition from a symmetric to an asymmetric core and a reduction in Peierls stress. This combination ductilizes the alloy as the number of available slip planes is increased and the critical stress needed to start plastic deformation is lowered.

Reconstruction of molecular orbital densities from photoemission data.

Photoemission spectroscopy is commonly applied to study the band structure of solids by measuring the kinetic energy versus angular distribution of the photoemitted electrons. Here, we apply this experimental technique to characterize discrete orbitals of large pi-conjugated molecules. By measuring the photoemission intensity from a constant initial-state energy over a hemispherical region, we generate reciprocal space maps of the emitting orbital density. We demonstrate that the real-space electron distribution of molecular orbitals in both a crystalline pentacene film and a chemisorbed p-sexiphenyl monolayer can be obtained from a simple Fourier transform of the measurement data. The results are in good agreement with density functional calculations.

Time-dependent density functional theory versus Bethe-Salpeter equation: an all-electron study.

We perform first-principle calculations for the macroscopic dielectric function within an all-electron full-potential framework, i.e., the linearized augmented planewave (LAPW) method. To this extent we pursue two different routes, which are, on the one hand, many-body perturbation theory (MBPT) by solving the Bethe-Salpeter equation (BSE), and time-dependent density functional theory (TDDFT), on the other hand. The implementation into one program package, i.e., the EXC!TiNG code, allows a direct comparison of these approaches, in terms of accuracy as well as efficiency. We briefly review the theory, and describe in detail the quantities specific to the LAPW method. As an example we show the results for GaAs which is well-investigated in the literature, and hence serves as a stringent test case.

Influence of the core-valence interaction and of the pseudopotential approximation on the electron self-energy in semiconductors.

State-of-the-art theory addresses single-electron excitations in condensed matter by linking density-functional theory (DFT) with many-body perturbation theory. In actual calculations it is common to employ the pseudopotential (PP) approach, where pseudo-wave-functions enter the calculation of the self-energy, and the core-valence interaction is treated at the DFT level. In this Letter we present accurate all-electron calculations of the self-energy and systematically compare the results to those of PP calculations. The analysis for a range of different materials reveals that both above mentioned approximations are indeed problematic.

Complex quasiparticle band structure induced by electron-phonon interaction: band splitting in the 1x1H/W110 surface.

We show that the self-consistent solution of the complex Dyson equation for the electron-phonon (EP) problem introduces many body effects which are often observed in photoemission experiments. The formalism is applied to the H covered W(110) surface, using first-principles results for the electronic and vibrational structure. We demonstrate that the measured spin-polarized surface band splitting [Phys. Rev. Lett. 84, 2925 (2000)10.1103/PhysRevLett.84.2925; 89, 216802 (2002)] can be traced back to different quasiparticle (QP) states induced by EP coupling. Despite the breakdown of the single QP picture, the spectral functions are very well represented by the predicted multiple QP structure.

Characterization of step-edge barriers in organic thin-film growth.

Detailed understanding of growth mechanisms in organic thin-film deposition is crucial for tailoring growth morphologies, which in turn determine the physical properties of the resulting films. For growth of the rodlike molecule para-sexiphenyl, the evolution of terraced mounds is observed by atomic force microscopy. Using methods established in inorganic epitaxy, we demonstrate the existence of an additional barrier (0.67 electron volt) for step-edge crossing-the Ehrlich-Schwoebel barrier. This result was confirmed by transition state theory, which revealed a bending of the molecule at the step edge. A gradual reduction of this barrier in the first layers led to an almost layer-by-layer growth during early deposition stage. The reported phenomena are a direct consequence of the complexity of the molecular building blocks versus atomic systems.

Importance of van der Waals interaction for organic molecule-metal junctions: adsorption of thiophene on Cu(110) as a prototype.

We report ab initio calculations for the interface energetics of a weakly adsorbed organic molecule on a metal surface, which serves as a model interface relevant for organic electronics. The studied thiophene ring is found to be physisorbed on the Cu(110) surface with an adsorption energy of -0.50 eV. Nonlocal correlations, i.e., van der Waals interactions, are solely responsible for the binding in this weakly interacting system, and the choice of the proper exchange-correlation function is crucially important. The adsorption of thiophene lowers the metal work function due to the formation of surface dipoles while no sizable charge transfer is found.

Lowest optical excitations in molecular crystals: bound excitons versus free electron-hole pairs in anthracene.

By solving the Bethe-Salpeter equation for the electron-hole Green function for crystalline anthracene we find the lowest absorption peak generated by strongly bound excitons or by a free electron-hole pair, depending on the polarization direction being parallel to the short or the long molecular axis, respectively. Both excitations are shifted to lower energies by pressure. The physical difference of these excitations is apparent from the electron-hole wave functions. Our findings are a major contribution to solve the long-standing puzzle about the nature of the lowest optical excitations in organic materials.

Suppression of electron-hole correlations in 3D polymer materials.

We present an ab initio study of the optical absorption spectra of isolated as well as crystalline trans-polyacetylene. We include excitonic effects by solving the Bethe-Salpeter equation for the electron-hole two-particle correlation function. We observe that the strength of the electron-hole interactions drastically reduces when going from an isolated polymer chain to a crystalline arrangement. This is not only a result of enhanced screening in the 3D material, but also of the increased spatial extent of the exciton perpendicular to the polymer chains. We point out that these findings apply to crystalline phases of conjugated polymers and molecular crystals in general.