D.
Koulentianos‡
*abc,
S.
Carniato‡
a,
R.
Püttner‡
d,
J. B.
Martins
a,
O.
Travnikova
ae,
T.
Marchenko
ae,
L.
Journel
ae,
R.
Guillemin
ae,
I.
Ismail
ae,
D.
Céolin
e,
M. N.
Piancastelli
af,
R.
Feifel
b and
M.
Simon
ae
aSorbonne Université, CNRS, Laboratoire de Chimie Physique-Matière et Rayonnement, F-75005 Paris Cedex 05, France
bDepartment of Physics, University of Gothenburg, Origovägen 6B, SE-412 96 Gothenburg, Sweden
cChemical Sciences and Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, IL 60439, USA. E-mail: koulentianos@anl.gov
dFachbereich Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany
eSynchrotron SOLEIL, L’Orme des Merisiers, Saint-Aubin, BP 48, F-91192 Gif-surYvette Cedex, France
fDepartment of Physics and Astronomy, Uppsala University, Box 516, SE-751 20 Uppsala, Sweden
First published on 22nd April 2021
Using synchrotron radiation in the tender X-ray regime, a photoelectron spectrum showing the formation of single site double-core-hole pre-edge states, involving the K shell of the O atom in CO, has been recorded by means of high-resolution electron spectroscopy. The experimentally observed structures have been simulated, interpreted and assigned, employing state-of-the-art ab initio quantum chemical calculations, on the basis of a theoretical model, accounting for their so-called direct or conjugate character. Features appearing above the double ionization threshold have been reproduced by taking into account the strong mixing between multi-excited and continuum states. The shift of the σ* resonance below the double ionization threshold, in combination with the non-negligible contributions of multi-excited configurations in the final states reached, gives rise to a series of avoided crossings between the different potential energy curves.
In what concerns molecular DCH continuum states, the creation of the two vacancies may occur either in the core shell of one atom (denoted as K−2, for the case of K-shell vacancies), or in the core shells of two different atoms (denoted as K−1K−1 when different K shells are involved). The first case is referred to as single site (ss) DCH states, while the second as two site (ts) DCH states. As it was predicted in ref. 23, ts-DCH states will be characterized by enhanced chemical-shifts, allowing for a more sensitive probe of the chemical environment of an atom in a molecule, compared to the traditional electron spectroscopy for chemical analysis (ESCA) technique, developed by K. Siegbahn.25 A theoretical investigation of this feature of ts-DCH states has been reported in ref. 26 for the formamide molecule and the nucleobases. On the other hand, ss-DCH states will be characterized by large orbital relaxation effects, as it was again theoretically predicted in ref. 23, and was further validated afterwords.4,16
Concerning the experimental observation of such states, upon single photon absorption using SR, the works of Eland et al.1 and Lablanquie et al.2 reported for the first time on the formation of ss-DCH continuum states in different small light-atom (C, N, O) containing molecules. In both ref. 1 and 2, the DCH continua were measured by means of time-of-flight (TOF) multi-electron coincidence spectroscopy, using a magnetic bottle spectrometer. With the same experimental setup, the ts-DCH continuum of the triply bonded C2H2 molecule was recorded shortly later,3 and subsequently, Nakano et al.5 measured the ss and ts-DCH continua of CO, N2 and the C2H2n (n = 1,2,3) series. Finally, this technique was used as well in order to perform DCH measurements in bigger molecules like propylene6 (C3H6) and benzene7 (C6H6). The limitations on the observation of the elusive ts-DCH states, using SR, have been demonstrated and discussed in ref. 6, for the case of ethyl trifluoroacetate (C4H5F3O2). Such limitations have been shown to be overcome, using instead of SR, the intense femtosecond pulses of an XFEL.18–22
Beyond DCH continuum states, the formation of a fully depleted K-shell through a simultaneous core-ionization core-excitation mechanism was observed in ref. 1, 8 and 10. Here, like in the case of DCH continuum states, a significant contraction of the valence orbitals takes place upon K-shell photoionization, which in turn leads to the shake-up of the second K-shell electron, instead of its shake-off like in the formation of a DCH continuum state. This type of DCH states is usually denoted as K−2V, with V referring to the unoccupied valence or Rydberg orbital that the core-excited electron occupies. The interest for molecular states of such form stems from the fact that they access, in highly-symmetric molecules, final-state orbitals which are non-reachable through conventional near edge X-ray absorption fine structure (NEXAFS) spectroscopy. A detailed theoretical description of DCH pre-edge states has been given by Carniato et al.,9 interpreting the observed spectral features in terms of their direct (dipolar ionization-monopolar excitation) and conjugate (monopolar ionization-dipolar excitation) nature, upon single photon absorption, similarly to the theoretical description of valence shake-ups initially given by Martin and Shirley.27 In general, the K−2V signatures have been treated so far like NEXAFS features of the K−1 cation.
For the experimental observation of DCH pre-edge states, beyond TOF electron spectroscopy, high-resolution single electron spectroscopy, using a hemispherical electron analyser, has been proven to be a very efficient tool. The first experimental studies, using this technique, reported on the formation of K−1L−1V states in argon11 and K−2V in neon.12 Concerning the observation of molecular K−2V states, Carniato and co-workers13 recorded the core-ionized core-excited states involving the K shell of C atom in CO2, while Feifel et al.14 reported on the K−2V states of CS2 and SF6, involving the K shells of each and every atom. Both ref. 13 and 14, are two showcases for the advantages of K−2V spectroscopy, as for the centro-symmetric molecules studied there, the forbidden in NEXAFS spectroscopy g → g transitions were clearly observed. Accordingly, following the work done in ref. 11, the K−1L−1V DCH states of the iso-electronic to argon HCl molecule15 were investigated. Furthermore, a study of the K−2V states in CH3CN, a molecule possessing two non-equivalent carbon atoms, along with a discussion of the factors contributing to the chemical-shift between them, can be found in ref. 16. Finally, a high-resolution electron spectrum, recording the K−2V states of water can be found in ref. 17. More recently, it was demonstrated from Mazza and coworkers,28 for the case of Ne, how such K−2V resonances can be reached using instead of SR the light produced at an XFEL. XFEL studies are expected to shed light on the formation of DCH pre-edge states on third and higher row elements, which are characterized by very small cross-sections,29 making their experimental observation using SR impracticable. Regarding the neutrally excited DCH states of the form K−2VV′, where both core electrons have been excited to the V and V′ unoccupied molecular/Rydberg orbitals, though such states are known to be accessible through X-ray absorption spectroscopy (XAS),30,31 yet no experimental evidence on their Auger decay is available to date. Finally, another interesting feature of DCH spectroscopy can be found in the probe of the induced nuclear dynamics, upon the formation of a double inner-shell vacancy. On one hand, the (sub)femto-second lifetimes of DCH states would not allow for a significant change of the molecular geometry, but on the other “hand”, their highly repulsive potential energy surfaces, could even result in ultra-fast dissociation.32,33
Within the frame of the current work, the formation of K−2V DCH states, involving the K shell of O atom in CO has been recorded by means of single electron high-resolution spectroscopy. The experimental photoelectron spectrum has been reproduced by theoretical simulations, taking into account the direct and conjugate character of each observed transition, with experiment and theory being in very good agreement. Previously known aspects of K−2V spectroscopy, like the increased term values (TV) of the different pre-edge states and the shift of the σ* transition below threshold, are present here. Nevertheless, as several avoided crossings appear in the potential energy curves (PEC) of the different final states, an interpretation of the observed spectral features as the K−1 analogue of the regular CO NEXAFS spectrum34 should be treated with caution. A behaviour similar to the present K−2V spectrum has been observed in the K−1V absorption spectrum of O2 molecule.35–37 In the latter case, the below-threshold position of the σ* resonance, within the Rydberg region, causes a series of avoided crossings, making the assignment of the different resonances non-trivial.
The energy calibration of the acquired spectra was performed in two steps. First the LMM Auger spectrum of argon was recorded, and the L3M2,3M2,3(1S0) transition at 201.09 eV40 kinetic energy was used to calibrate the spectrometer. Subsequently, the photon energy was calibrated based on the binding energy of 248.63 eV for argon 2p3/2−1.41
The initial and final states have been described with a post Hartree–Fock (HF) configuration interaction (CI) procedure, using for oxygen an aug-cc-PCV5Z basis set42 augmented by additional diffuse (s,p,d) orbitals, in order to better take into account Rydberg states and the relaxation of the valence orbitals in the presence of double core vacancies. A lighter cc-pVTZ basis set42 was used for the neighboring carbon atom.
Our approach relies on the same set of orthogonal molecular orbitals (MO) for both the ground state of the neutral molecule, and the DCH states of the singly-ionized molecular ions. We made the choice to consider a set consisting of HF MOs optimized for the DCH K−2 ion. The final-state CI wave functions of the molecular ions are expanded over DCH (2h-1e) configurations, DCH/valence-hole (3h-2e) configurations and DCH/double-valence-hole (4h-3e) configurations.
The initial ground state has also been described at a CI level of theory including the main K2 closed-shell configuration, 1σ2(O1s)2σ2(C1s)3σ24σ21π45σ2, as well as the single core-hole, single core-hole and single valence-hole and double valence-hole excited configurations. In principle, in the HF Hamiltonian of an N-particle system, single-core excitations will not contribute because of the Brillouin theorem.43 However, here the initial neutral ground state is built from the same set of the doubly core-ionized K−2 MOs as for K−2V, i.e. quite far from the expected MOs corresponding to the neutral atom. Under these peculiar conditions, single-core and single core-valence excitations in the CI expansion of the initial state give tiny contributions, but are essential for the description of the direct process of the basic K−2V and K−2Vv′ → v′′ multi-excited states, (see previous works on this subject8–10,14). The large expansion used in this description allows compensating the relaxation effects induced by the DCH in the reference set. The choice of a unique reference set of molecular orbitals has also the advantage of providing a simple identification of contributing terms in, respectively, direct and conjugate cross sections. First, this unique set gives the dichotomous (0,1) values for the overall integrals taking part in the direct amplitudes. Second, it reduces the (N − 1)-electron dipolar matrix elements taking part in the conjugate amplitudes, to single-electron matrix elements. The latter confirms that the CI expansions of the initial and final states are well balanced. All the configurations present in the description of the final state can be reached both through the direct and the conjugate paths, except for the (4h-3e) ones that would need the presence of (3h-3e) configurations in the ground state to be reached through the direct process.
While the peaks below the double ionization potential (DIP) can be essentially described as an ejection of one core electron accompanied by a single discrete transition, the situation is less straightforward in the region above the DIP, where one can expect the effect of coupling between the discrete multi-excited states and the double-ionization continuum. The simulation of this region taking into account the coupling effect is out of reach due to a huge number of CI configurations, but in order to adequately reproduce this region we adopted the following strategy.
In order to calculate the continuum profile of the spectrum a static-exchange (STEX-HF) approximation technique44–49 has been adapted and combined with the Stieltjes imaging technique.50,51 Our analysis reveals three main families of the CI states contributing to the spectrum above the DIP. In the first case, more than 90% of a CI state consists of a single (2h-1e) configuration. Since the corresponding core-excitation in the presence of a core-hole involves a vacant orbital located above the DIP, this state was considered as a pure continuum state. In the second case, more than 90% of a CI state consists of (3h-2e) and (4h-3e) configurations corresponding to discrete multi-excitations where core-excitations involve vacant orbitals located below the DIP. In this situation, similarly to the one encountered for the region below the DIP, the CI state was treated as a pure discrete state. In the third case, a CI state exhibits strong mixing between a single “continuum” state and various discrete “multiple-excitations” states. Therefore, for the sake of simplicity, these states and their associated cross sections were considered in the pool of the continuum states. The contribution of the discrete multi-excitations is superimposed to the simulated double continuum. Here we have assumed that there is no coupling/interference between the discrete and the double continuum parts.
Apart from this consideration, there is also a certain imbalance in the description of these core/valence excited electronic states compared to the final states in the region below the DIP. Indeed, if a correlation treatment at a CI level of theory, including single and double excitations CI(SD), is acceptable to treat formally K−2V single inner-shell excitations like in DCH (2h-1e) final states, namely the case for most of the peaks under the DIP, a larger degree of correlation up to CI(SDTQ), i.e. taking into account also triple and quadruple excitations, is in principle required to be taken into account for the description of K−2Vv′ → v′′ states in the region above the DIP. This imbalance leads to calculated transition energies which are potentially too large. However, the calculation of electronic transitions becomes computationally very demanding with the inclusion of quadruple excitations leading to a very large configuration space, not possible to be treated with our home code. Thus, we can consider that the calculated K−2Vv′ → v′′ photoelectron satellites have too large energies (by 1 eV to 2 eV) compared to the experimental data. For better comparison between experimental and theoretical results, we have shifted the calculated K−2Vv′ → v′′ structures by 1 eV.
In the simulations, we make an assumption that the angular distribution parameter β associated with the direct shake-up is equal to 2 (p-wave) for the O 1s photoelectron, and equal to 0 for the conjugate shake-up (s-wave). Considering a half acceptance angle of 20° ± 10° for the entrance lens of the analyzer, a factor of 2.8 ± 0.2 between the direct and the conjugate shake-up pathways is taken into account.12
The potential curve of the 1Σ+ ground state of neutral CO was simulated with a Morse potential:
V(r) = D·(1 − e−ρ(r−r0))2 | (1) |
In order to analyze the nature and properties of the most important transitions in the O K−2V spectrum of CO, a series of single point energy calculations at the DFT level of theory have been performed. In these calculations only one configuration within each state has been considered, namely the one with the highest coefficient, while all the other configurations have been neglected. Single point of energy with mono-determinantal technique must be taken with caution when strong interaction of few configurations occurs. However, this approach is useful to analyze various properties associated with each main configuration interacting at play. For that we have estimated a few different parameters, including (i) the mean value of the square root of the squared radius of the molecular orbital, (ii) the gradient (in eV Å−1) at the initial equilibrium distance (1.128 Å); (iii) the width of the Franck–Condon envelope and (iv) the effective quantum number n* deduced from the hydrogen-like formula , (α0 being the bohr radius and
= 0,1 for s and p states) and the corresponding term values compared to ab initio DFT/BL3LYP values.
For the 2π* orbital, the quantity was calculated at 1.44 Å at the HF level of theory and 1.41 Å at the DFT/B3LYP level of theory (see Tables S2 and S3 in the ESI†), thus close to that of the highest occupied molecular orbital (HOMO) 5σ (1.61 Å), indicating in this way a valence rather than a Rydberg character. The theoretical position of this transition was calculated at the DFT/B3LYP level of theory at 1153.8 eV, a value differing from the experimental value by 2.2 eV. It should be noted that such mismatches are systematic within a DFT treatment of DCH processes.5–7 For that reason, the theoretical spectrum was blue-shifted to 1156 eV, in order to match the experimental excitation energy to the 1(2Π) state.
The same shift has been applied to the value of the K−2 threshold (indicated with an orange bar in Fig. 1), resulting in a value of 1177 eV binding energy, in fairly good agreement with the experimental value of 1178.0 ± 0.8 eV, measured by Lablanquie et al.2 From the aforementioned numbers, a TV of ≈22 eV is estimated for the transition to the 2π* MO, a value which is ≈2.5 times higher than the 8.97 eV estimated for the singly core-excited K−12π* state in ref. 34, where the vibrational sub-structure of the transition was also recorded, by means of photo-absorption spectroscopy. With the latter value being in excellent agreement with the TV of of 9.0 eV estimated in Sham et al.56 by means of NEXAFS spectroscopy, and with the one of 8.3 eV obtained by Hitchcock and Brion57 through electron energy loss spectroscopy (EELS), it is evident that the fully depleted K shell causes a strong attractive field, which shifts these core-ionized core-excited states away from threshold, as has been already found in several previous works focusing on the observation of K−2V states.13,14,16
A final point to be discussed, before moving on to region II, is the slope of the PEC of the K−22π* state, shown as a circle-marked light blue line in the middle panel of Fig. 2. As it was mentioned above, the PEC is strongly dissociative, which is also expected from the Z + 2 approximation as well. More specifically, it is known that the K−2π* DCH state, for the iso-electronic case of N2 is stable10 and its equivalent Z + 2 approximation, NF+ in the ground state, is also known to be stable.58 In contrast, the equivalent of O K−22π* in CO is given by the noble-gas containing cation CNe+, which is expected to be dissociative in the Franck–Condon region of CO in the ground state. As can be seen from Fig. 2, the extent of the Franck–Condon region is from ≈1.08 Å to ≈1.18 Å. We note that Frenking et al.59 calculated the ground state of CNe+ to be weakly bound with a dissociation energy De ≈ 130 meV and a bond length of 2.077 Å.
According to the Condon reflection approximation,60 the slope of the PEC will depend on the full width at half maximum (FWHM) of the Gaussian profile describing the line-shape of the recorded transition, through S = aσ, as has been demonstrated in ref. 61, where S is the slope of the potential energy curve, σ the FWHM of the Gaussian and a2 = μω/ħ, with μ being the reduced mass of the molecule and ω its vibrational frequency in the electronic ground state. The line-shape of the transition will formally be described by a Voigt profile, being the result of the convolution of a Lorentzian profile accounting for the lifetime broadening, with a Gaussian profile having two contributions. The first is due to the experimental resolution, while the second is the result of the dissociative broadening. Yet, as it is explained below by giving the relative widths of each contribution, a Gaussian profile can adequately describe the above mentioned resonance, thus the value of the slope can be estimated from the width of the peak. From a fit analysis using a Gaussian profile, we obtained a value of σ = 1.90 ± 0.15 eV. Beyond the dissociative broadening, a Lorentzian profile of ΓDCH = 513 meV, should also be used. The value of ΓDCH above has been set equal to 3 × ΓSCH,12,62 with ΓSCH the lifetime of the O K−1 SCH state of CO taken from ref. 63. Finally, convolution with a Gaussian having a FWHM of 1 eV should be included, in order to properly account for the experimental resolution. From the values of the FWHM of the above mentioned profiles, it can be seen that the Gaussian part has a major contribution in the line-shape of the peak corresponding to the K−22π* state, therefore it was chosen for our fit analysis, instead of the more complex Voigt function. From the obtained value of σ, the value of the slope is estimated at −39.9 ± 3.2 eV Å−1. This value is in relatively good agreement with the theoretical value of −33.2 eV Å−1 obtained at the DFT/B3LYP level of theory (Table S3 in the ESI†), yet it is significantly different from the CI numerical gradient of −24.9 eV Å−1. The reason for the latter mismatch can be found in correlation effects, which are not treated properly within the CI method framework.
For the first peak at 1165.80 eV, the K−26σ configuration was found from our theoretical investigations to be ≈77% of the CI state, at the equilibrium distance of the electronic ground state (see Table S1, ESI†), with the final state MO having a strong oxygen 3s character, with a less intense, yet non-negligible, carbon 3p (lying along the inter-nuclear axis) contribution (see Fig. S2 and Tables S2 and S3 of the ESI†). As shown in Fig. 1, this structure has a dominant direct component, though its conjugate (Σ+) character is not negligible. In addition, from our theoretical results it was found that a strong mixing with the K−27σ* and the K−25σ12π*2 configurations takes place. The 7σ* orbital of the former configuration shows a nodal structure typical for the unoccupied σ* orbital obtained in the LCAO approach from the 2pσ orbitals of carbon and oxygen. In the latter configuration an electron from the 5σ MO has been promoted to the 2π*, being now doubly occupied. We note that the latter configuration is allowed under the ΔΛ = 0, ±1 selection rule.
For the 6σ orbital, the was estimated at 1.85 Å at the DFT/B3LYP level of theory (Table S3, ESI†), concomitant with a quantum defect of δ = 0.8, or equivalently a small effective quantum number n* = n − δ = 2.2. The values of
and n*, in the HF level of theory, are predicted at 2.07 Å and 2.20, respectively (Table S2, ESI†). Note that the size of the 6σ orbital is comparable to the size of the valence orbitals 2π* and 7σ*, although from its nodal structure it has a strong O 3s and C 3p character (see Fig. S2, ESI†). Obviously, the high positive charge due to the DCH, seen by this orbital, causes a significant contraction. As a result, it does not experience the molecular dication as almost point-like and cannot be considered a pure Rydberg orbital, but has to be considered as an orbital with significant valence character, i.e. it has an influence on the molecular bond. Despite its size, the shape of the 6σ orbital can be considered mainly as 3s-like (see Fig. S2 in the ESI†). In the SCH studies of ref. 34, 56 and 57, a similar K−13sσ transition to a Rydberg orbital was identified, with its TV estimated at ≈3.6 eV in ref. 34 and 56. From the recorded peak position and the calculated K−2 threshold, the TV of the 1(2Σ+) state is estimated at ≈11 eV, thus increased by a factor of 3 compared to the SCH case, for the same reasons as explained above.
The PEC of the 1(2Σ+) state is shown with triangle-marked red solid line in the upper panel of Fig. 2. At the internuclear distance of ≈1.2 Å the slope becomes steeper due to an avoided crossing with the 2(2Σ+) state (star-marked dashed black line). The energetic limit at large internuclear distances is identical to that of the K−22π* state, because the same atomic states are involved. Thus we can conclude that at larger bond distances the K−27σ* configuration, including the anti-bonding σ* orbital, becomes dominant.
The CI gradient at the equilibrium distance associated to the 1(2Σ+) state was found at −27.2 eV Å−1, and the deduced CI Franck–Condon FWHM at 1.5 eV. Additionally, if the O K−2 core-hole lifetime and the experimental resolution are taken into account, this value becomes 2.15 eV, which is in reasonable agreement with the experimental value of ≈2.5 eV. Considering on the other hand a single-point energy for the K−26σ configuration, the DFT gradient was estimated at −40.0 eV Å−1 (Table S3 in the ESI†), with the deduced numerical FWHM of the corresponding Franck–Condon envelope at 2.25 eV. However, it should be pointed out here that there are two possible explanations for the differences between experiment and theory. First, the Franck–Condon approximation is not strictly valid due to the avoided crossing at the upper limit of the Franck–Condon region. Second, contributions of transitions to the 3(2Σ+) state also show up at higher binding energies of the peak in region II.
The PEC shown by the star-marked black dashed line corresponds to a 2(2Σ+) state, which mainly consists of the configurations K−25σ12π*2 as well as K−25σ12π′*2 at the equilibrium distance of the ground state, i.e. it is reached after a K → 2π* transition accompanied by an additional 5σ → 2π* excitation. Within the framework of our model, the intensity of this line contributes only through the direct path from the initial O K15σ12π*2 configuration of the ground state, which leads to K−25σ12π*2 configuration after photon absorption and K-shell ionization. However, such core-valence excited configurations play a minor role in the neutral ground state and the final cross-section is very weak. We note here that the K−25σ12π*2 state can be formed as well through the doubly core-excited K−25σ22π*2 configuration, which contributes to the ground state, after valence ionization of the 5σ MO. This valence ionization channel could interfere with the core-valence direct path described above, but has not been considered in the present study. The main relevance of the K−25σ12π*2 configuration for the present study is probably that it leads to avoided crossings.
The weak structure, theoretically predicted at 1167.14 eV, is caused by a transition to the 3(2Σ+) state. Its PEC is shown with the square-marked green solid line in Fig. 2. For this state the same configurations as for the 1(2Σ+) state should interact, but in addition, a K−28σ configuration is also present there. Nonetheless, as the two lines caused by the transitions to the 1(2Σ+) and 3(2Σ+) state are very close in energy, they can not be resolved experimentally.
As shown in the upper panel of Fig. 2, a second avoided crossing within the Franck–Condon region, is present between the 4(2Σ+) and the 5(2Σ+) PECs, while in the middle panel of Fig. 2, an avoided crossing in the edges of the Franck–Condon region between the 3(2Π) and 4(2Π) PECs is also evident. Because of their presence, the given assignments are valid only within the Franck–Condon region and a one-to-one correspondence with the NEXAFS study of the neutral species in ref. 34 should not be attempted. This can be explained by the below-threshold shift of the 7σ* resonance and states of the type K−25σ−12π*(3Π)nσ(2Π) which are both strongly dissociative. This results in the diabatic picture in a crossing with the K−2Ryd states which are metastable with a potential-well depth of about 1 eV and an equilibrium distance of ≈1.55 Å found when the present PECs are transfered to a diabatic picture. Similar values, namely a dissociation energy of 1.6 eV and an equilibrium distance of 1.67 Å have been found in calculations performed for the Z + 2 molecule CNe2 +, see ref. 66. In case of the K−1 resonant excitations, these states with a different bond character are well separated in energy. The significantly different behaviour of the cationic and neutral CO, found here, is in line with the results of a recent theoretical and experimental study by Couto and co-workers,67 where the NEXAFS spectra, involving both the O and C K edges were recorded for CO+, after the creation of an initial 5σ−1 vacancy.
Before we discuss the PECs of the states 6(2Σ+) to 8(2Σ+) and in particular 3(2Π) to 5(2Π), we want to point out that the state K−25σ−11π−1(3Π)6σ has also been calculated at a binding energy of 1174.6 eV and could contribute to the peak IV at 1172 as a conjugate channel. This finding is in agreement with calculations for the dication CNe2+, see ref. 66, which represents the Z + 2 approximation for the K−2 states of CO. In CNe2+ the states σ−1π*1(3Π), σ−1π*1(1Π), and σ−1π*1(3Σ) are 4.6 eV, 8.7 eV, and 10.8 eV, respectively, above the K−2 ground state. These values are in reasonable agreement with the splittings calculated for the double-shake states. Based on these considerations we can expect dissociative K−25σ−11π−1(3Π)nλ states below the K−2 threshold, which have a significant influence on the PECs.
The dissociative PECs of the states 6(2Σ+) to 8(2Σ+) are shown in the upper panel of Fig. 2 with the dashed purple, dashed-dotted cyan and dotted black lines, respectively. The PECs of the 3(2Π) to 5(2Π) states are shown in the middle panel of Fig. 2, with the x-marked dotted brown, black dashed and dotted green lines, respectively. The strong interactions of the final states corresponding to the above mentioned configurations are reflected again in the avoided crossings between the different PECs. The most prominent level crossings are at ≈1.33 Å between the 2(2Π) and 3(2Π) states, and at ≈1.20 Å between the 3(2Π) and 4(2Π) states. By converting the adiabatic picture to a diabatic picture, i.e. by eliminating the avoided level crossings, one obtains a strongly dissociative state, which is described in the limit of long internuclear distances by the cross-marked pink dotted curve labeled 2(2Π) and has at the equilibrium distance of the ground state an energy of ≈1168 eV. The latter value suggests in agreement with our calculations that the dissociative configuration is given by K−25σ−11π−1(3Π)6σ, see above. This configuration is expected to be dissociative and crosses the K−2Ryd states, which are expected to be weakly bound in the diabatic picture, see above.
In the PECs a number of additional avoided level crossings can be observed, so for example at ≈1.60 Å between 3(2Π) and 4(2Π) as well as at ≈1.35 Å between 4(2Π) and 5(2Π) states. These avoided crossings are probably due to other dissociative K−25σ−11π−1(3Π)nλ configurations. In particular, at the equilibrium distance of the ground state, a double shake-up K−25σ−11π−1(3Π)7σ* state is theoretically predicted at ≈ 1177 eV binding energy. The corresponding peak can be seen in the calculations displayed in Fig. 1, together with two more peaks at ≅1178 and ≅1179 eV, which probably belong to transitions to the same kind of states.
The slight mismatch between the experimental and the theoretical intensities can be explained by the formation of a well known background, which has been observed in several K−2V studies.11,15,16 This background is caused by the higher n unresolved Rydberg states below threshold, as well as the slow electron of the double ionization process above threshold, and can be modeled by an arctan function,68 being the result of the convolution of a step function with the Lorentzian profile describing the lifetime broadening of the core-excited state.
From a more general point of view we found that the O K−2V spectrum of CO is more complex than the O K−1V of the same molecule. In the latter spectrum the O 1s → 2π* and the O 1s → Ryd resonances below the ionization threshold are well separated from the O 1s → σ* shape resonance and the double excitations. This separation does not exist for the DCH states and leads to a number of avoided level crossings and consequently to significant mixing of different configurations. It can be explained by three findings. First, the unoccupied σ* valence orbital, which is normally present above the ionization threshold in SCH studies and is known as shape resonance, is in the present study significantly below threshold. Second, in the present case the K−2v′ → v′′ satellites are energetically closer to the K−2 state than the K−1v′ → v′′ satellites to the K−1 threshold. Third, the term value of the excited state orbital V in the K−2Vv′ → v′′ double shake-up states is larger than in the K−1Vv′ → v′′ doubly exited states due to the higher charge seen by this orbital in the DCH states. Since the last argument is expected to hold in general and since the shift of the unoccupied σ* valence orbital from above threshold in SCH states to below threshold in DCH states has been observed already for a number of molecules, we expect that the higher complexity of the K−2V spectra is a general phenomenon. This result suggests that an interpretation of K−2V structures as cationic K−1 absorption lines should be treated with caution, as the assignment of the different features might become non-trivial due to the avoided crossings, especially if a significant change in the equilibrium geometry takes place. This is not the case in the present work, as the equivalent bond length of the molecule in its electronic ground state is at 1.128 Å,69 while the bond length of the O doubly core-ionized CO has been calculated at 1.160 Å.70 However, further investigations on similar systems like N2, NO, and O2 are necessary to shed more light on the matter. Finally, multi-shake-up excitations involving the occupied MOs of the CO molecule have been observed and give rise to an intense above-threshold spectral feature.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cp00607j |
‡ These authors contributed equally to this work. |
This journal is © the Owner Societies 2021 |