A computational approach is presented for prediction and interpretation of core-level spectra of complex molecules. Applications are presented for several isolated organic molecules, sampling a range of chemical bonding and structural motifs. Comparison with gas phase measurements indicate that spectral lineshapes are accurately reproduced both above and below the ionization potential, without resort to ad hoc broadening. Agreement with experiment is significantly improved upon inclusion of vibrations via molecular dynamics sampling. We isolate and characterize spectral features due to particular electronic transitions enabled by vibrations, noting that even zero-point motion is sufficient in some cases.
PACS numbers: UnknownWhen applied to molecular systems, core level spectroscopies are powerful probes of both occupied and unoccupied electronic states, uniquely revealing intimate details of both intra-and inter-molecular interactions [1]. Methods involving x-ray absorption (XAS, NEXAFS, XANES) or x-ray photo-electron spectroscopy (XPS) are increasingly being applied to complex molecular systems, including nucleotides, peptides and large organic molecules [2]. However, a major limitation of this technology is the fact that extraction of molecular information from these experiments often depends explicitly on comparisons with theoretical calculations, which are extremely challenging to perform at experimental accuracy. In this Letter, we describe the extension of a recently developed method for predicting core-level spectra of condensed phases [3] to isolated organic molecules -pyrrole, s-triazine, pyrrolidine and glycine -which demonstrates qualitative improvements over existing methods [4][5][6] in comparison with experiment and provides new insights into the origins of particular spectral features in terms of coupling of electronic and vibrational degrees of freedom.The challenges for simulating gas phase core-level spectra are maintaining accuracy in the following areas: (1) description of the core-hole excited state; (2) representation of both bound excitonic states below the ionization potential (IP) and resonance states in the continuum above the IP; and (3) inclusion of vibrational effects, either due to experiments being performed near room temperature, or from intrinsic zero-point motion.Density functional theory (DFT) [7,8] has proved accurate in reproducing the excitation energies associated with core-level spectra via total energy differences (socalled ∆SCF or ∆KS) [9]. Accordingly, we model the lowest energy core-level excited state self-consistently using a full core-hole and excited electron (XCH) [3]. This is particularly important for molecular systems, where screening of the core-hole excitation is greatly enhanced by the presence of the excited electron, which can be strongly bound to the core-hole in the lowest energy excited state. In contrast, for non-molecular condensed phases, such as covalent and ionic crystalline solids, the inherent dielectric screening of the valence charge density o...