Florian Libisch

Hot Electrons Do the Impossible: Plasmon-Induced Dissociation of H2 on Au

Heterogeneous catalysis is of paramount importance in chemistry and energy applications. Catalysts that couple light energy into chemical reactions in a directed, orbital-specific manner would greatly reduce the energy input requirements of chemical transformations, revolutionizing catalysis-driven chemistry. Here we report the room temperature dissociation of H(2) on gold nanoparticles using visible light. Surface plasmons excited in the Au nanoparticle decay into hot electrons with energies between the vacuum level and the work function of the metal. In this transient state, hot electrons can transfer into a Feshbach resonance of an H(2) molecule adsorbed on the Au nanoparticle surface, triggering dissociation. We probe this process by detecting the formation of HD molecules from the dissociations of H(2) and D(2) and investigate the effect of Au nanoparticle size and wavelength of incident light on the rate of HD formation. This work opens a new pathway for controlling chemical reactions on metallic catalysts.

Photovoltaic effect in an electrically tunable van der Waals heterojunction.

Semiconductor heterostructures form the cornerstone of many electronic and optoelectronic devices and are traditionally fabricated using epitaxial growth techniques. More recently, heterostructures have also been obtained by vertical stacking of two-dimensional crystals, such as graphene and related two-dimensional materials. These layered designer materials are held together by van der Waals forces and contain atomically sharp interfaces. Here, we report on a type-II van der Waals heterojunction made of molybdenum disulfide and tungsten diselenide monolayers. The junction is electrically tunable, and under appropriate gate bias an atomically thin diode is realized. Upon optical illumination, charge transfer occurs across the planar interface and the device exhibits a photovoltaic effect. Advances in large-scale production of two-dimensional crystals could thus lead to a new photovoltaic solar technology.

The Multiradical Character of One- and Two-Dimensional Graphene Nanoribbons

When is an acene stable? The pronounced multiradical character of graphene nanoribbons of different size and shape was investigated with high-level multireference methods. Quantitative information based on the number of effectively unpaired electrons leads to specific estimates of the chemical stability of graphene nanostructures.

Wave-function mapping of graphene quantum dots with soft confinement.

Using low-temperature scanning tunneling spectroscopy, we map the local density of states of graphene quantum dots supported on Ir(111). Because of a band gap in the projected Ir band structure around the graphene K point, the electronic properties of the QDs are dominantly graphenelike. Indeed, we compare the results favorably with tight binding calculations on the honeycomb lattice based on parameters derived from density functional theory. We find that the interaction with the substrate near the edge of the island gradually opens a gap in the Dirac cone, which implies soft-wall confinement. Interestingly, this confinement results in highly symmetric wave functions. Further influences of the substrate are given by the known moiré potential and a 10% penetration of an Ir surface resonance into the graphene layer.

Embedded correlated wavefunction schemes: theory and applications.

Conspectus Ab initio modeling of matter has become a pillar of chemical research: with ever-increasing computational power, simulations can be used to accurately predict, for example, chemical reaction rates, electronic and mechanical properties of materials, and dynamical properties of liquids. Many competing quantum mechanical methods have been developed over the years that vary in computational cost, accuracy, and scalability: density functional theory (DFT), the workhorse of solid-state electronic structure calculations, features a good compromise between accuracy and speed. However, approximate exchange-correlation functionals limit DFTs ability to treat certain phenomena or states of matter, such as charge-transfer processes or strongly correlated materials. Furthermore, conventional DFT is purely a ground-state theory: electronic excitations are beyond its scope. Excitations in molecules are routinely calculated using time-dependent DFT linear response; however applications to condensed matter are still limited. By contrast, many-electron wavefunction methods aim for a very accurate treatment of electronic exchange and correlation. Unfortunately, the associated computational cost renders treatment of more than a handful of heavy atoms challenging. On the other side of the accuracy spectrum, parametrized approaches like tight-binding can treat millions of atoms. In view of the different (dis-)advantages of each method, the simulation of complex systems seems to force a compromise: one is limited to the most accurate method that can still handle the problem size. For many interesting problems, however, compromise proves insufficient. A possible solution is to break up the system into manageable subsystems that may be treated by different computational methods. The interaction between subsystems may be handled by an embedding formalism. In this Account, we review embedded correlated wavefunction (CW) approaches and some applications. We first discuss our density functional embedding theory, which is formally exact. We show how to determine the embedding potential, which replaces the interaction between subsystems, at the DFT level. CW calculations are performed using a fixed embedding potential, that is, a non-self-consistent embedding scheme. We demonstrate this embedding theory for two challenging electron transfer phenomena: (1) initial oxidation of an aluminum surface and (2) hot-electron-mediated dissociation of hydrogen molecules on a gold surface. In both cases, the interaction between gas molecules and metal surfaces were treated by sophisticated CW techniques, with the remainder of the extended metal surface being treated by DFT. Our embedding approach overcomes the limitations of conventional Kohn-Sham DFT in describing charge transfer, multiconfigurational character, and excited states. From these embedding simulations, we gained important insights into fundamental processes that are crucial aspects of fuel cell catalysis (i.e., O2 reduction at metal surfaces) and plasmon-mediated photocatalysis by metal nanoparticles. Moreover, our findings agree very well with experimental observations, while offering new views into the chemistry. We finally discuss our recently formulated potential-functional embedding theory that provides a seamless, first-principles way to include back-action onto the environment from the embedded region.

Electron-hole crossover in graphene quantum dots.

We investigate the addition spectrum of a graphene quantum-dot in the vicinity of the electron-hole crossover as a function of perpendicular magnetic field. Coulomb-blockade resonances of the 50 nm wide dot are visible at all gate voltages across the transport gap ranging from hole to electron transport. The magnetic field dependence of more than 50 states displays the unique complex evolution of the diamagnetic spectrum of a graphene dot from the low-field regime to the Landau regime with the n=0 Landau level situated in the center of the transport gap marking the electron-hole crossover. The average peak spacing in the energy region around the crossover decreases with increasing magnetic field. In the vicinity of the charge neutrality point we observe a well resolved and rich excited state spectrum.

Graphene quantum dots: Beyond a Dirac billiard

We present realistic simulations of quantum confinement effects in phase-coherent graphene quantum dots with linear dimensions of 10\char21{}40 nm. We determine wave functions and energy-level statistics in the presence of disorder resulting from edge roughness, charge impurities, or short-ranged scatterers. Marked deviations from a simple Dirac billiard for massless fermions are found. We find a remarkably stable dependence of the nearest-neighbor level spacing on edge roughness suggesting that the roughness of fabricated devices can be possibly characterized by the distribution of measured Coulomb blockade peaks.

Dynamically encircling an exceptional point for asymmetric mode switching.

Physical systems with loss or gain have resonant modes that decay or grow exponentially with time. Whenever two such modes coalesce both in their resonant frequency and their rate of decay or growth, an ‘exceptional point’ occurs, giving rise to fascinating phenomena that defy our physical intuition. Particularly intriguing behaviour is predicted to appear when an exceptional point is encircled sufficiently slowly, such as a state-flip or the accumulation of a geometric phase. The topological structure of exceptional points has been experimentally explored, but a full dynamical encircling of such a point and the associated breakdown of adiabaticity have remained out of reach of measurement. Here we demonstrate that a dynamical encircling of an exceptional point is analogous to the scattering through a two-mode waveguide with suitably designed boundaries and losses. We present experimental results from a corresponding waveguide structure that steers incoming waves around an exceptional point during the transmission process. In this way, mode transitions are induced that transform this device into a robust and asymmetric switch between different waveguide modes. This work will enable the exploration of exceptional point physics in system control and state transfer schemes at the crossroads between fundamental research and practical applications.

Size-extensivity-corrected multireference configuration interaction schemes to accurately predict bond dissociation energies of oxygenated hydrocarbons

Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this schemes use on larger molecules and for more complex regions of combustion PESs.

A comparison of singlet and triplet states for one- and two-dimensional graphene nanoribbons using multireference theory

This study examines the radical nature and spin symmetry of the ground state of the quasi-linear acene and two-dimensional periacene series. For this purpose, high-level ab initio calculations have been performed using the multireference averaged quadratic coupled cluster theory and the COLUMBUS program package. A reference space consisting of restricted and complete active spaces is taken for the π-conjugated space, correlating 16 electrons with 16 orbitals with the most pronounced open-shell character for the acenes and a complete active-space reference approach with eight electrons in eight orbitals for the periacenes. This reference space is used to construct the total configuration space by means of single and double excitations. By comparison with more extended calculations, it is shown that a focus on the π space with a 6-31G basis set is sufficient to describe the major features of the electronic character of these compounds. The present findings suggest that the ground state is a singlet for the smaller members of these series, but that for the larger ones, singlet and triplet states are quasi-degenerate. Both the acenes and periacenes exhibit significant polyradical character beyond the traditional diradical.