Tomographic identification of all molecular orbitals in a wide binding energy range

arXiv:2501.05287 · physics.chem-ph, cond-mat.mes-hall, cond-mat.mtrl-sci · Submitted 2025-01-09 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: Today's paper: "Tomographic identification of all molecular orbitals in a wide binding energy range".

Mira: This study presents a comprehensive experimental and theoretical investigation using photoemission orbital tomography (POT) to identify all molecular orbitals in bisanthene adsorbed on a Cu(110) surface across a substantial binding…

Kai: First, who's behind it and why it matters.

Title and authors: Kai: We’ve just talked about the scope of this paper, specifically how they went about identifying all molecular orbitals on bisanthene across a broad binding energy range using photoemission orbital tomography. Now let’s talk about who put this work out there.

Mira: The authors are a solid team including Anja Haags, Dominik Brandstetter, Xiaosheng Yang, Larissa Egger, Hans Kirschner, Alexander Gottwald, Mathias Richter, Georg Koller, Francois C. Bocquet—and F. Stefan Tautz—all contributing to this experimental and theoretical effort.

Lev: It’s interesting to see a collaboration spanning both experimentalphysik IV A at RWTH Aachen University, NAWI Graz at the University of Graz, and the PTB in Berlin; that breadth suggests a very thorough approach to validating the results.

Kai: That diversity of institutions definitely lends credibility when you're talking about fundamental electronic structure investigations like this.

Mira: The title itself, "Tomographic identification of all molecular orbitals in a wide binding energy range," really tells you the ambition here is not just to find one or two important states, but to get a comprehensive picture of the entire electronic structure across a significant energy window.

Lev: A comprehensive identification across that whole range means they’re not just looking at the frontier orbitals near zero energy; they’re mapping out the entire potential landscape, which is important for understanding stability and reaction pathways.

Kai: Exactly, and it sets a high bar for what POT can achieve in terms of orbital resolution on these types of adsorbed systems.

Mira: And when you combine that comprehensive experimental reach with the subsequent functional benchmarking, the implication is that they’re providing a really strong foundation for how we should interpret DFT results in this specific area.

Lev: That foundation is what we need; if the benchmark itself is sound, then any quantum simulation built on top of it has a much higher chance of being meaningful.

Kai: So, we’ve seen the title and the authors, and they’re clearly setting out to deliver a detailed map of how these systems interact across their entire electronic spectrum.

Mira: And that ambition is what drives the rest of the paper, moving from experimental measurement to theoretical scrutiny.

The paper's summary: Kai: Now we’re getting into the actual summary of what this study actually achieved, which is quite detailed. They used photoemission orbital tomography to experimentally determine the orbital-resolved partial densities of states by analyzing the angular distribution of photoelectrons as a function of binding energy.

Mira: Specifically, they employed ARPES to measure momentum maps, which allowed them to extract these "orbital fingerprints" from the electron emission patterns.

Lev: From a hardware standpoint, extracting that information means they’ve successfully bypassed some of the inherent limitations of just measuring bulk properties and are getting direct access to the quantum states involved in the process.

Kai: They then took that raw data cube, Iexp(kx, ky; Eb), and used various numerical techniques—linear fitting, Monte Carlo fitting with Gaussian weight functions, and multilayer averaging—to deconvolute it into contributions from individual molecular Kohn-Sham orbitals.

Mira: That process of deconvolution is the heart of their experimental methodology; they are essentially solving a complex pattern recognition problem to separate the signal originating from different energy levels.

Lev: I’m thinking about those fitting methods—the Monte Carlo approach, for example—that require a lot of careful regularization to ensure you don't introduce artifacts when trying to fit something so complex.

Kai: The main result they present is that this process successfully yielded the experimental identification of fifteen π and twenty-three σ Kohn-Sham orbitals covering a binding energy range exceeding ten eV below the Fermi level.

Mira: That successful identification of all those specific molecular orbitals across that wide energy span is a significant accomplishment because it validates the method for comprehensive mapping, not just a spot check.

Lev: Getting that complete set of assignments means they’ve established a robust dataset for validating theoretical models, which is the crucial step before you can even think about running any complex simulations on real quantum hardware.

Kai: The paper also highlights how the choice of functional impacts the calculated charge transfer into states like pi(forty) and influences hybridization between molecular and metallic states.

Mira: That’s a vital finding because it shows that the theoretical description of charge movement is not static; it depends on whether you use PBE or HSE for the exchange-correlation term.

Lev: If we can tune those theoretical descriptions systematically, it gives us levers to control the physics in our simulations, which is exactly what error correction relies on.

Kai: The analysis of orbital characteristics further shows how orbital energies shift during adsorption—specifically noting an overall shift of roughly zero point five eV or lower for sigma orbitals upon adsorption.

Mira: That energy alignment detail is powerful because it gives us a physical picture of the bonding at the interface, showing how the molecular orbitals reorganize when they are actually stuck onto a substrate.

Lev: Understanding those shifts is necessary when we try to design materials for quantum systems where subtle changes in orbital energy can lead to entirely different ground states or transport properties.

The paper's improvements: Kai: So, the paper isn't just reporting the results; they suggest ways to improve this methodology, and those suggestions are pretty interesting for pushing the boundaries of what’s possible in photoemission tomography.

Mira: They point out that improved analyzer performance could sharpen the pDOS curves, which could potentially reveal systematic differences between Kohn-Sham and Dyson orbitals, not just in energy but also in k space.

Lev: If you can see differences in k space, that means the error isn't just an energy shift; it’s a failure of the orbital description itself, which is a much more fundamental problem when we think about simulating quantum phenomena.

Kai: They also suggest extending the methodology to unoccupied states using pump-probe ARPES experiments as a way to gain even more information.

Mira: Extending it to unoccupied states would be a huge step because right now they are focused heavily on occupied states, so seeing the unoccupied picture is essential for a complete picture of the electronic structure.

Lev: For error correction, we need full knowledge of the excited states too; if we only know about the occupied ones, we can’t fully characterize the Hamiltonian we are trying to simulate.

Kai: It seems like these improvements focus on increasing resolution and scope—better measurement sensitivity and looking at unoccupied states—to get that next level of detail.

Mira: And they are also pushing the theoretical side by suggesting we need better ways to rigorously compare KS orbitals with the experimentally measured Dyson orbitals, which is where their benchmark work really shines.

Lev: I think that push towards rigorous comparison between KS and Dyson orbitals is exactly what we need; it’s about building a better theoretical bridge between theory and experiment.

Conclusion: Kai: We’ve covered a lot today regarding the identification of all molecular orbitals in this paper, including the experimental setup, the functional benchmarking, and those suggested improvements for future work. In short, this paper successfully identified fifteen pi and twenty-three sigma Kohn-Sham orbitals for bisanthene on Cu(one hundred ten) over a binding energy range exceeding ten eV below the Fermi level.

Mira: To summarize, the main implication is that this work provides an essential benchmark for electronic structure methods, allowing us to quantitatively compare DFT calculations with experimental data and revealing how functional choice impacts charge transfer and hybridization at these interfaces.

Lev: For running on real hardware, having a reliable functional choice that minimizes systematic errors is key because it gives us confidence when we try to translate those models into actual quantum simulations.

Kai: This study successfully achieved its goal of retrieving the orbital structures and validating the use of DFT Kohn-Sham orbitals as reliable approximations for the experimentally measured Dyson orbitals.

Mira: The future direction suggested is that sharper analyzer performance could reveal systematic differences between Kohn-Sham and Dyson orbitals in k space, which would be incredibly valuable for refining our understanding of these systems.

Lev: Having that level of detail on the orbital structure is a necessary step toward building models that can accurately describe complex quantum physics.

Peter Grünberg Institut (PGI-3) · J¨ulich Aachen Research Alliance (JARA) · Experimentalphysik IV A, RWTH Aachen University · Institute of Physics, NAWI Graz, University of Graz · Physikalisch-Technische Bundesanstalt (PTB)

physics.chem-ph, cond-mat.mes-hall, cond-mat.mtrl-sci

Submitted: 2025-01-09

Updated: 2025-01-09

Comments: arXiv admin note: text overlap with arXiv:2209.11516

Journal ref: Phys. Rev. B 111, 165402 (2025)

DOI: 10.1103/PhysRevB.111.165402

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 83/100

The gist: This study presents a comprehensive experimental and theoretical investigation using photoemission orbital tomography (POT) to identify all molecular orbitals in bisanthene adsorbed on a Cu(110)

Key concepts

Photoemission Orbital Tomography (POT): This technique uses the angular distribution of emitted photoelectrons as a function of binding energy to extract orbital-resolved partial densities of states. It allows researchers to map out the electronic structure, identifying which specific molecular orbitals are responsible for which features in the experimental data.
Kohn-Sham Orbitals vs. Dyson Orbitals: Kohn-Sham orbitals are calculated within DFT, while Dyson orbitals are what is directly measured experimentally via photoemission. The study shows that Kohn-Sham orbitals serve as good approximations of the measured Dyson orbitals across a broad energy range, offering insight into their relationship.
Range-Separated Hybrid Functionals (HSE): This specific type of DFT functional was found to perform best for this organic/metal interface. It accurately predicts the electronic structure by combining different levels of exchange, significantly reducing errors compared to simpler functionals like PBE and B3LYP.
Intramolecular Dispersion: This refers to the energy contribution arising from electron correlation within the molecule itself. The study confirmed that this dispersion structure in bisanthene/Cu(110) is fundamentally derived from the known band structure properties of graphene.

Terminology

Monte Carlo — 19× in this paper
Benchmark — 13× in this paper · explained in 2 episodes
Analyzer — 12× in this paper · explained in 1 episode
Benchmarking — 11× in this paper · explained in 3 episodes
DFT Calculations — 11× in this paper · explained in 2 episodes
Alignment — 9× in this paper · explained in 1 episode
Photoemission Orbital Tomography — 9× in this paper
Density Functional Theory — 5× in this paper
Density Functional Theory (DFT) — 4× in this paper · explained in 2 episodes
Regularization — 4× in this paper · explained in 3 episodes
Consistency — 3× in this paper · explained in 3 episodes
Accuracy — 2× in this paper · explained in 1 episode
Core — 2× in this paper · explained in 2 episodes
European Research Council — 2× in this paper
Physikalisch-Technische Bundesanstalt — 2× in this paper
Quantum Physics — 1× in this paper · explained in 10 episodes
Spectral Signature — 1× in this paper · explained in 2 episodes

Summary

This study presents a comprehensive experimental and theoretical investigation using photoemission orbital tomography (POT) to identify all molecular orbitals in bisanthene adsorbed on a Cu(110) surface across a substantial binding energy range exceeding 10 eV. This work is significant because it provides an essential benchmark for electronic structure methods, allowing for the quantitative comparison of Density Functional Theory (DFT) calculations with experimental data, thereby revealing crucial insights into the nature of Kohn-Sham orbitals and their approximation to Dyson orbitals in organic/metal interfaces.

Experimental Identification via Photoemission Orbital Tomography (POT)

The core methodology involves using the angular distribution of photoelectrons as a function of binding energy to extract orbital-resolved partial densities of states (pDOS). The researchers employed ARPES, which allows for the measurement of momentum maps, to reveal orbital fingerprints. Key experimental steps included:

: Measuring photoemission intensity over a large binding energy window (e.g., 13 eV) and recording momentum maps at fixed binding energies or varying photon energies (hν = 35 eV and 57 eV). The angular distribution of the photoelectrons is related to the squared modulus of the Fourier transform of the initial-state Dyson orbital.

The resulting data cube, Iexp(kx, ky; Eb), is then deconvolved into contributions from individual molecular Kohn-Sham orbitals using various numerical approaches:

Linear fitting (minimizing Eq. 5 for each energy).
Monte Carlo fitting (using Gaussian weight functions with regularization constraints).
Multilayer averaging (MLA), which leverages established pattern recognition methods, such as a more general pattern recognition problem.

Theoretical Benchmarking and Functional Comparison

The experimentally determined orbital-resolved pDOS serves as a stringent benchmark for DFT calculations employing four different exchange-correlation functionals: PBE, HSE, PBE0, and B3LYP. The study systematically compares the calculated molecular-orbital-projected densities of states (MOPDOS) against experimental data. Key findings regarding functional performance include:

: The range-separated hybrid functional HSE performs best for the investigated organic/metal interface.

The comparison reveals that Kohn-Sham orbitals are good approximators of the Dyson orbitals measured in photoemission, over a wide range of binding energies and not just near the Fermi level. Furthermore, the analysis shows how different functionals affect orbital characteristics:

: The choice of functional influences the calculated amount of charge transfer into the LUMO (π(4,0)) and influences the degree of hybridization between molecular and metallic states.

Analysis of Orbital Characteristics and Dispersion

The study provides detailed insights into the electronic structure, including orbital energy alignments and intramolecular dispersion. The analysis focuses on:

Energy Level Alignment: Examining how orbital energies shift during adsorption by considering four steps: isolated molecule, adsorbed geometry (with structural distortion), free-standing monolayer, and finally the bisanthene/Cu(110) interface. The results show an overall shift of roughly 0.5 eV to lower energies for σ orbitals upon adsorption.
Intramolecular Dispersion: The reconstruction of the data cube from theoretical momentum maps and fitted weights confirms that the experimental data cube implicitly contains the intrinsic dispersion structure, which is derived from graphene's band structure. This demonstrates that the intramolecular dispersion in bisanthene/Cu(110) derives from the band structure of graphene.

Benchmarking DFT Functionals Against Experiment

A quantitative benchmarking exercise was conducted using mean errors (ME) and mean absolute errors (MAE) to compare DFT predictions with experimental binding energies. The results highlight functional strengths:

**: HSE is found to be the most accurate functional of the four tested ones for the present purpose, as it reduces MAE for both π and σ orbitals compared to B3LYP. **

The analysis also addresses systematic errors, noting that PBE exhibits significant self-interaction errors, leading to erroneously small binding energies (negative ME). The study concludes that while the overall agreement is strong, for individual orbitals, there may be systematic differences not only between the energies but also between the momentum distributions of Dyson and Kohn-Sham orbitals.

Conclusion and Future Outlook

The research successfully achieved its goal: the experimental identification of 15 π and 23 σ Kohn-Sham orbitals of bisanthene on Cu(110), covering a binding energy range of more than 10 eV below the Fermi level. The successful retrieval of orbital structures validates the use of DFT Kohn-Sham orbitals as reliable approximations for experimentally measured Dyson orbitals. Future directions suggest that improved analyzer performance could sharpen pDOS curves, potentially revealing systematic differences between Kohn-Sham and Dyson orbitals, not only in energy but also k space. Furthermore, the methodology can be extended to unoccupied states via pump-probe ARPES experiments.

Improvements for AI systems

Here are specific improvements that an AI system could implement based on the findings of this scientific paper, focusing on enhancing its capabilities in materials science and electronic structure prediction:

The following improvements leverage the core methodology of Photoemission Orbital Tomography (POT) combined with Density Functional Theory (DFT) benchmarking:

Incorporate a Kohn-Sham to Dyson Orbital Approximation module.
Implement an Orbital-Resolved Partial Density of States (pDOS) Extraction engine.
Develop a Functional Benchmarking and Error Quantification pipeline for DFT calculations, specifically targeting organic/metal interfaces.

The improved AI system can perform the following specific tasks:

Incorporate a Kohn-Sham to Dyson Orbital Approximation module:
Perform orbital energy predictions that are more rigorous than standard Kohn-Sham (KS) energies for wide binding energy ranges (> 10 eV), specifically by using the KS orbitals as approximations for Dyson orbitals, which rigorously account for the electron extraction process in photoemission spectroscopy. This allows the system to predict and evaluate electronic structure properties over a much broader energy spectrum.
Implement an Orbital-Resolved Partial Density of States (pDOS) Extraction engine:
Extract experimentally relevant pDOS data from raw photoemission intensity maps by deconvolving the angular distributions (using ARPES data as input), thereby providing orbital-resolved spectral information for molecular orbitals that were previously difficult to measure or identify.
Develop a Functional Benchmarking and Error Quantification pipeline:
Benchmark the performance of various DFT exchange-correlation functionals (PBE, PBE0, HSE, B3LYP) against experimental binding energies for organic/metal interfaces on a per-orbital basis (15 π and 23 σ orbitals).
Quantify systematic errors: The system can explicitly calculate Mean Errors (ME) and Mean Absolute Errors (MAE) for the predicted orbital energies across different functionals to determine which functional is most accurate for specific molecular orbital types (e.g., identifying HSE as superior for certain π and σ orbitals in this system).
Analyze hybridization effects: The AI can distinguish between band broadening due to hybridization with the substrate versus intrinsic orbital width, providing a more nuanced understanding of electronic structure formation at interfaces.
Predict structural/electronic shifts upon adsorption: By simulating the adsorption process (gas-phase, distorted geometry, free-standing layer, adsorbed geometry), the system can predict how molecular orbital energies shift due to bonding dipoles and substrate interactions for specific orbitals (e.g., predicting the energy shift of σ(0,8) upon adsorption).
Identify Unreliable Orbital Assignments: The system can flag assignments where experimental data falls outside the reliable fitting range (e.g., category 4 in Table IV), indicating that the experimental resolution or data cube size was insufficient for a particular orbital state.

Abstract

In the past decade, photoemission orbital tomography (POT) has evolved into a powerful tool to investigate the electronic structure of organic molecules adsorbed on surfaces. Here we show that POT allows for the comprehensive experimental identification of all molecular orbitals in a substantial binding energy range, in the present case more than 10 eV. Making use of the angular distribution of photoelectrons as a function of binding energy, we exemplify this by extracting orbital-resolved partial densities of states (pDOS) for 15 π and 23 σ orbitals from the experimental photoemission intensities of the prototypical organic molecule bisanthene (C 28 H 14) on a Cu(110) surface. In their entirety, these experimentally measured orbital-resolved pDOS for an essentially complete set of orbitals serve as a stringent benchmark for electronic structure methods, which we illustrate by performing density functional theory (DFT) calculations employing four frequently-used exchange-correlation functionals. By computing the respective molecular-orbital-projected densities of states of the bisanthene/Cu(110) interface, a one-to-one comparison with experimental data for an unprecedented number of 38 orbital energies becomes possible. The quantitative analysis of our data reveals that the range-separated hybrid functional HSE performs best for the investigated organic/metal interface. At a more fundamental level, the remarkable agreement between the experimental and the Kohn-Sham orbital energies over a binding energy range larger than 10,eV suggests that -- perhaps unexpectedly -- Kohn-Sham orbitals approximate Dyson orbitals, which would rigorously account for the electron extraction process in photoemission spectroscopy but are notoriously difficult to compute, in a much better way than previously thought.