Abstract

The magnetic dilution of vanadyl phthalocyanine (VOPc) within the isostructural diamagnetic titanyl phthalocyanine (TiOPc) affords promising molecular spin qubit platforms for solid-state quantum computing. The development of quantitative methods for determining how the interactions with a supporting substrate impact the electronic structure of the system are fundamental to determine their potential integration in physical devices. In this work we propose a combined approach based on x-ray absorption spectroscopy (XAS), atomic multiplet calculations, and density functional theory (DFT) to investigate the 3⁢𝑑 orbital level structure of VOPc on TiOPc/Ag(100). We characterize VOPc in different molecular environments realized by changing the thickness of the TiOPc interlayer and adsorption configuration on Ag(100). Depending on the molecular film structure, we find characteristic XAS features that we analyze using atomic multiplet calculations. We use a Bayesian optimization algorithm to accelerate the parameter search process in the multiplet calculations and identify the ground state properties, such as the 3⁢𝑑 orbital occupancy and splitting, as well as intra-atomic interactions. Our analysis indicates that VOPc retains its spin 𝑆= 1 2 character in all configurations. Conversely, the energy separation and sequence of the unoccupied V 3⁢𝑑 orbitals sensitively depend on the interaction with the surface and TiOPc interlayer. We validate the atomic orbital picture obtained from the multiplet model by comparison with DFT, which further allows us to understand the VOPc electronic properties using a molecular orbital description.