Energy Level Alignment at Molecule-Metal Interfaces from an Optimally-Tuned Range-Separated Hybrid Functional

Zhenfei Liu,David A. Egger,Sivan Refaely-Abramson,Leeor Kronik,Jeffrey B. Neaton

DOI: https://doi.org/10.1063/1.4975321

2017-03-07

Abstract:The alignment of the frontier orbital energies of an adsorbed molecule with the substrate Fermi level at metal-organic interfaces is a fundamental observable of significant practical importance in nanoscience and beyond. Typical density functional theory calculations, especially those using local and semi-local functionals, often underestimate level alignment leading to inaccurate electronic structure and charge transport properties. In this work, we develop a new fully self-consistent predictive scheme to accurately compute level alignment at certain classes of complex heterogeneous molecule-metal interfaces based on optimally-tuned range-separated hybrid functionals. Starting from a highly accurate description of the gas-phase electronic structure, our method by construction captures important nonlocal surface polarization effects via tuning of the long-range screened exchange in a range-separated hybrid in a non-empirical and system-specific manner. We implement this functional in a plane-wave code and apply it to several physisorbed and chemisorbed molecule-metal interface systems. Our results are in quantitative agreement with experiments, both the level alignment and work function changes. Our approach constitutes a new practical scheme for accurate and efficient calculations of the electronic structure of molecule-metal interfaces.

Mesoscale and Nanoscale Physics,Materials Science

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the problem of energy level alignment at the molecule - metal interface. Specifically, it addresses the following challenges:

1. **Accuracy of energy level alignment**: Traditional density functional theory (DFT), especially methods using local and semi - local functionals, often underestimates this alignment when calculating the energy level alignment at the molecule - metal interface, leading to inaccuracies in the electronic structure and charge transport properties.

2. **Handling of complex heterogeneous interfaces**: For some complex molecule - metal interfaces, especially those involving strong chemisorption or significant charge transfer, traditional methods are difficult to provide accurate results.

3. **Capturing of non - local effects**: The energy level alignment at the molecule - metal interface is affected by important non - local surface polarization effects, and these effects are difficult to accurately describe in traditional DFT.

To address these problems, the authors developed a new fully self - consistent prediction scheme based on the optimally - tuned range - separated hybrid functional (OT - RSH). This scheme captures non - local effects by adjusting the long - range screened exchange term and can be applied to various types of molecule - metal interface systems, including physisorption and chemisorption systems.

### Main objectives

- **Improve the accuracy of energy level alignment**: Ensure that the calculation results are quantitatively consistent with experimental data, especially in terms of energy level alignment and work function changes.
- **Expand the scope of application**: Be applicable not only to weakly - coupled physisorption systems but also to strongly - coupled chemisorption systems.
- **Provide an efficient calculation method**: While ensuring accuracy, reduce the calculation cost, making the electronic structure calculation of complex interface systems more feasible.

Through this method, researchers hope to more accurately understand and predict the functional characteristics of the molecule - metal interface, thereby providing theoretical support for applications in fields such as nanoscience, energy conversion, and storage.

### Formula summary

The main formulas involved in the paper are as follows:

1. **Decomposition of the Coulomb operator**:
\[
\frac{1}{|r - r'|}=\alpha+\beta\frac{\text{erf}(\gamma|r - r'|)}{|r - r'|}+\left(1-\left[\alpha+\beta\frac{\text{erf}(\gamma|r - r'|)}{|r - r'|}\right]\right)
\]

2. **Expression for the exchange - correlation energy**:
\[
E_{XC}=\alpha E_{X}^{\text{SR}}+(1 - \alpha)E_{\text{GGA}}^{\text{SR}}+(\alpha+\beta)E_{X}^{\text{LR}}+(1 - \alpha-\beta)E_{\text{GGA}}^{\text{LR}}+E_{\text{GGA}}^{C}
\]

3. **Objective function for the optimization parameter γ**:
\[
J(\gamma)=|\varepsilon_{\gamma}(N)+IP_{\gamma}(N)|^{2}+|\varepsilon_{\gamma}(N + 1)+IP_{\gamma}(N + 1)|^{2}
\]

4. **Orbital energy renormalization condition**:
\[
\varepsilon_{N}^{\gamma}(\beta)-\varepsilon_{N}^{\gamma}(\beta_{0})=P=\frac{1}{4(z - z_{0})}
\]

Through these formulas, the authors achieve accurate calculations of the energy level alignment at the molecule - metal interface.