Abstract:The possibility to perform low-pressure gas-phase photoemission experiments has allowed the exploration at molecular scale of the complex correlation between electronic and structural properties of the matter. A significant advantage in the removal of the interference of the surrounding at low-pressure is also the easier correlation of the experimental results with theoretical ab-initio simulations. In the context of core ionization, the development of accurate methods for calculating inner-shell binding energies (BEs) offers a dual benefit. Firstly, it aids significantly in interpreting experimental data, helping to assign the contributions of all non-equivalent atoms, even in unresolved features arising from a molecular structure. Secondly, it allows for the anticipation of experimental results by accurately predicting spectral lines. In this context, we have developed and extensively tested a computational protocol based on plane-wave/pseudopotential density functional theory (PW-DFT) to predict XPS BEs in molecules and molecular aggregates. This protocol is based on a ∆SCF approach and has been tested comparing the present theoretical results with experimental XPS spectra collected from several molecular systems in gas phase, ranging from benzene derivates to biomolecules. Our calculations have been performed using semilocal and global/range-separated hybrid density functionals, containing increasing fractions of Hartree-Fock exact exchange (EXX). Specifically, PBE, B3LYP(20 % EXX), HSE(range separated with 25 % EXX at short range ) and BH&HLYP (50 % EXX) have been used for the assessment of the computational protocol. EOM-CCSD has been employed as reference method for comparison. Our PW-DFT approach demonstrated to be generally accurate and robust even using semilocal DFT, and also suitable for application to very large molecular systems and to organic thin films deposited on inorganic surfaces.

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is to develop and validate a computational protocol based on plane - wave / pseudopotential density - functional theory (PW - DFT) to accurately predict X - ray photoelectron spectroscopy (XPS) binding energies (BEs) in isolated molecules, molecular aggregates and periodic systems. Specifically, the research aims to: 1. **Improve the accuracy of experimental data interpretation**: By precisely calculating the inner - shell binding energies (BEs), help to assign the contributions of all non - equivalent atoms, even when unresolved features appear in the molecular structure. 2. **Predict experimental results**: By accurately predicting spectral lines, predict experimental results in advance, thus providing guidance for experimental design. 3. **Evaluate the effects of different exchange - correlation functionals**: Use functionals with different proportions of Hartree - Fock exact exchange (EXX) (such as PBE, B3LYP, HSE, BH&HLYP) to evaluate their contributions to computational accuracy and stability. 4. **Expand the application range**: Ensure that this method is applicable not only to isolated molecules, but also to very large molecular systems and cases where organic thin films are deposited on inorganic surfaces. ### Main research contents - **Theoretical framework**: Adopt the PW - DFT computational framework based on the ΔSCF method. - **Reference method**: Use EOM - CCSD as a reference method for comparison. - **Test systems**: Include various gas - phase molecular systems from benzene derivatives to biomolecules. - **Functional selection**: Use semi - local functionals (such as PBE) and global / range - separated hybrid functionals (such as B3LYP, HSE, BH&HLYP), and evaluate the influence of different EXX proportions. - **Treatment of periodic systems**: For periodic systems (such as doped surfaces or molecules adsorbed on metal surfaces), introduce the van der Waals density - functional (vdW - DF) to consider non - local long - range interactions. ### Key formulas - **Binding energy calculation formula**: \[ BE_{\text{C(1s)}} = 290.0+(E_{C1s}^{\text{(target)}} - E_{C1s}^{\text{(reference)}}) \] \[ BE_{\text{N(1s)}} = 405.6+(E_{N1s}^{\text{(target)}} - E_{N1s}^{\text{(reference)}}) \] \[ BE_{\text{O(1s)}} = 541.25+(E_{O1s}^{\text{(target)}} - E_{O1s}^{\text{(reference)}}) \] where \( E_{C1s}^{\text{(target)}} \), \( E_{N1s}^{\text{(target)}} \) and \( E_{O1s}^{\text{(target)}} \) are the electron energies of C(1s), N(1s) and O(1s) core ionization in the target molecule respectively; \( E_{C1s}^{\text{(reference)}} \), \( E_{N1s}^{\text{(reference)}} \) and \( E_{O1s}^{\text{(reference)}} \) are the electron energies of the corresponding core ionization in the reference molecules CH₄, NH₃ and CO₂ respectively. Through these methods and formulas, this research demonstrates the effectiveness and computational cost - effectiveness of the PW - DFT method in simulating the XPS binding energies of organic molecules and molecular aggregates.

Photoemission Spectroscopy of Organic Molecules Using Plane-Wave/Pseudopoential Density Functional Theory: A Comprehensive Computational Protocol for Isolated Molecules, Molecular Aggregates and Periodic Systems

Photoemission Spectroscopy and Orbital Imaging from Koopmans-Compliant Functionals

The Nuts and Bolts of Ab-Initio Core-Hole Simulations for K-shell X-Ray Photoemission and Absorption Spectra

Quantitative Predictions of Photoelectron Spectra in Amorphous Molecular Solids from Multiscale Quasiparticle Embedding

Molecular structure retrieval directly from laboratory-frame photoelectron spectra in laser-induced electron diffraction

Photoemission orbital tomography based on tight-binding approach: method and application to $π$-conjugated molecules

XPS core-level chemical shift by ab initio many-body theory

Modeling Liquid Photoemission Spectra: Path-Integral Molecular Dynamics Combined with Tuned Range-Separated Hybrid Functionals

Ab-initio angle and energy resolved photoelectron spectroscopy with time-dependent density-functional theory

Enhancing the Accuracy of XPS Calculations: Exploring Hybrid Basis Set Schemes for CVS-EOMIP-CCSD Calculations

Combination of the TD-DFT and the pixel counting method for determining electron capture cross sections for protons impacting on organic molecules

Photoemission orbital tomography for excitons in organic molecules

Combined Self-Consistent-Field and Spin-Flip Tamm-Dancoff Density Functional Approach to Potential Energy Surfaces for Photochemistry

Photoelectron angular distributions from strong-field ionization of oriented molecules

Molecular inner-shell photoabsorption/photoionization cross sections at core-valence-separated coupled cluster level: Theory and examples

Simulation of XPS C1s Spectra of Organic Monolayers by Quantum Chemical Methods

Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional

Excited-state potential-energy surfaces of metal-adsorbed organic molecules from Linear Expansion Δ-Self-Consistent Field Density-Functional Theory (ΔSCF-DFT)

Theoretical assessment of vibrationally resolved C1s X-ray photoelectron spectra of simple cyclic molecules

Photoionization and core resonances from range-separated density-functional theory: General formalism and example of the beryllium atom

A minimalist approach to 3D photoemission orbital tomography: algorithms and data requirements