Core-hole localization and ultra-fast dissociation in SF6

Resonant inelastic x-ray scattering spectra excited at the fluorine K resonances of SF6 have been recorded. While a small but significant propensity for electronically parity-allowed transitions is found, the observation of parity-forbidden electronic transitions is attributed to vibronic coupling that breaks the global inversion symmetry of the electronic wavefunction and localizes the core hole. The dependence of the scattering cross section on the polarization of the incident radiation and the scattering angle is interpreted in terms of local π/σ symmetry around the S–F bond. This symmetry selectivity prevails during the dissociation that occurs during the scattering process.


Introduction
Electronic states in a inversion-symmetric system have well-defined parity. Hence, a fluorine 1s core hole in the SF 6 molecule is typically described in terms of a multiconfigurational symmetry-adapted wavefunction, comprising core holes at various fluorine sites. If inversion symmetry is broken the sites become inequivalent and the core hole can localize Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. [1,2]. In resonant inelastic x-ray scattering (RIXS) symmetry breaking due to ultrafast electronic-vibronic coupling is manifested in population of electronically parity forbidden final states [3][4][5]. Even in such cases, however, RIXS spectra are sensitive to the symmetry of the wavefunction. There is still a propensity for populating 'allowed' relative to 'forbidden' [3] final states if the geometry is close to being inversion symmetric. Moreover, the local symmetry of molecular orbitals can be probed by varying the angle between the polarization direction of the incident radiation and the direction of the scattered radiation. Symmetry considerations are essential in many fields of physics, and the role of global and local symmetry has e.g., been investigated for transmission and reflection in scattering systems [6,7]. The global symmetry breaking discussed in this paper concerns the electronic part of the wavefunction, while the symmetry selectivity due to the dipole approximation (dipole selection rules) may very well be operational for the total wavefunction which includes both electronic and vibronic degrees of freedom [5,8].
In a crude ionic picture each F atom in SF 6 has a filled 2p shell, with atomic orbitals oriented either parallel or perpendicular to the S-F bond forming local σ or π type bonds, respectively. Electrostatic interactions and covalency lift the degeneracy so that molecular orbitals of various symmetry are energetically separated. Together with the relative simplicity of s ↔ p transitions and the initial inversion symmetry, this makes the SF 6 molecule a showcase for symmetry sensitivity in RIXS. x-ray emission spectra (XES) of SF 6 have earlier been measured with electron excitation [9,10], and RIXS have been measured at the K edge [11] and the L edge [12,13] of sulfur.
Here we present RIXS spectra of SF 6 excited in the vicinity of the F K edge. Indeed, electronically dipole forbidden transitions are observed, demonstrating dynamic breaking of the global electronic inversion symmetry and suggesting core-hole localization. Only a small but significant propensity for parity allowed transitions is found. A dependence of the polarization of the incident radiation is observed, which is interpreted in terms local dipole selection rules related to the local π/σ symmetry around the S-F bond. Analysis of the electron emission has earlier demonstrated that dissociation can occur prior to the core hole decay [14][15][16][17], and indeed a sharp peak suggesting atomic transitions is observed in the present RIXS spectra. The polarization dependence indicates that local π/σ symmetry is retained, and our analysis shows that the sharp peak primarily is associated with emission at large interatomic distances during dissociation.

Experimental
RIXS spectra were measured at the F K edge using the AERHA spectrometer [20] mounted at the SEXTANTS beamline [21] of the SOLEIL light source and a gas cell with a 100 nm thick Si 3 N 4 membrane, through which both primary and secondary radiation penetrated. The SF 6 sample pressure was 1 bar, and the angle between the incident radiation and the membrane surface was 25 • . At a scattering angle of 85 • , spectra were measured with both vertical and horizontal linear polarization of the incident radiation.
The absolute energy scale of the monochromator was calibrated by matching the fluorescence yield spectrum to the published x-ray absorption spectrum in reference [22]. The spectrometer energy scale was tied to the monochromator energy via elastic scattering. From the elastic peaks the overall energy resolution was estimated to be 0.3 eV.
The XES spectrum is best aligned to the valence-band photoemission spectrum (PES) [18] using 694.6 eV for the F 1s binding energy, which complies with the value given by Siegbahn et al [19]. A precise determination of the BE from XES is hampered as the ionisation is accompanied by a large number of unresolved vibrational excitations that contribute to the width and the general shape of the PES peak. Whether the adiabatic or the vertical binding energy is the most appropriate reference for XES is influenced by the excitation-emission dynamics, and our results do not rule out other determinations of the fluorine 1s binding energy [23,24].

Theory
The local symmetry and character of different transitions are analyzed in the molecular orbital frame work of the electronic ground state, augmented with classical simulations of dynamical effects in the core-ionized state, as employed in previous studies of complex systems [25,26]. The electronic structure of the SF 6 molecule was described on the level of density functional theory in the CP2K software suite (version 6.0) [27][28][29][30] using the Gaussian augmented plane wave (GAPW) method [31], the hybrid PBE0 functional [32,33] and augcc-VTZP basis sets [34]. A hybrid functional was required to obtain a good description of the chemical bonding and an accurate ground state geometry. Combined with the Gaussian basis, an auxiliary plane wave basis (with a cut-off of 600 Ry) was used to describe the electron density in a cubic simulation cell of dimension 20 Å, combined with a Poisson solver to model isolated systems [35]. The geometry optimization yields O h symmetry and a S-F bond distance of 1.575 Å, which compares well with the experimental value of 1.564 Å [36].
Simulated fluorine K-edge XES were derived from ground state Kohn-Sham orbitals as implemented in the GAPW method [37]. The discrete XES transitions and the density of states were convoluted with a Gaussian broadening with a fullwidth half-maximum of 0.3 eV and ad hoc energy shifts were added for direct comparison to the measured XES and PES spectra. However, the shifts were chosen such that the theoretical PES and XES spectra are aligned in figure 1, since in the absence of dynamical effects they are modelled on the basis of the same Kohn-Sham orbitals. The influence of nuclear dynamics during the core-hole lifetime was investigated classically by performing classical ab initio molecular dynamics simulations to follow the dynamics in the F 1s core ionized state, in analogy with previous XES simulations [25]. The core ionization was modelled in the Z + 1 approximation by replacing SF 6 with SF 5 Ne + . The core-excited state trajectory was followed for 20 fs with a time-step of 0.5 fs which is sufficient for capturing nuclear dynamics during an estimated F 1s lifetime with an exponential decay of 3.3 fs [17]. The XES spectra of SF 6 were calculated for configurations along the trajectory in the core-ionized state and used for deriving a life-time average according to the exponential decay.

Results and discussion
The XES spectrum excited at 751.4 eV, far above any sharp resonances, is dominated by a 3-4 eV broad structure peaking at around 678 eV (figure 1). This peak can be attributed to  [18]. The XES spectrum is aligned to the PES spectrum using 694.6 eV for the core level binding energy [19]. The theoretical spectra are shifted to fit the experimental PES.
transitions from the F 1s core hole states to states with a hole in F 2p-derived outer orbitals.
In a pure ionic-bond picture the SF 6 molecule binds the S 6+ and F − closed shell atoms. It has O h symmetry with the formal electronic configuration: 1g (S1s 2 )2a 2 1g 1e 4 g 1t 6 1u (F1s 12 )3a 2 1g (S2s 2 )2t 6 1u (S2p 6 )4a 2 1g 3t 6 1u 2e 4 g (F2s 12 )5a 2 1g 4t 6 1u 1t 6 2g 3e 4 g 5t 6 1u 1t 6 2u 1t 6 1g (F2p 36 ). Assuming that the symmetry is unbroken the fluorine core orbitals are of a 1g , e g and t 1u symmetry. States with core vacancy of corresponding symmetries can be built in linear combinations of local F 1s hole states, and vice versa. As emission associated with electron transfer from all valence orbitals to the core levels is allowed, and as there is no significant symmetry selectivity for excitation high above the ionization limit the discussion of local versus global symmetry selectivity becomes irrelevant for the high-energy excited spectrum. This will, however, become crucial in the discussion of the resonantly excited spectra, and for this purpose we show the Kohn-Sham valence orbitals in figure 2. S-F bonding/antibonding combinations involving sulfur 3s orbitals constitute orbitals of a 1g symmetry, which for local S-F bonds unambiguously give σ symmetry (figure 2). Bonding/antibonding combinations with S 3p are of t 1u symmetry, and as seen in figure 2 this global symmetry allows for both π and σ local S-F symmetry. A feature corresponding to the bonding 4t 1u orbital is found around 672 eV emission (figure 1). There are weaker S-F interactions allowed via S 3d combinations giving rise to 3e g with local σ symmetry and 1t 2g with local π character, with intensity found at 675 eV emission energy (figure 1). The other global symmetries primarily correspond to 'non-bonding' orbitals. The 1t 1g and 1t 2u orbitals correspond to f symmetry at the sulfur, and are nonbonding, whereas 5t 1u is simply strongly localized on fluorine.
In the calculations, we observe that the 4t 1u has strong σ character and that 5t 1u has predominantly π character. Intensity corresponding to the 'non-bonding' orbitals make up the main 3-4 eV broad structure in the emission spectrum peaking at around 678 eV (figure 1). Notice that the corresponding antibonding a * 1g , t * 1u and t * 2g and e * g orbitals are unoccupied, and available for excitation in x-ray absorption.
The assignment of the features in the XES spectrum (figure 1) is accomplished by comparing to the He II excited PES of reference [18], and our theoretical analysis based on DFT calculations. The simulated XES and total density of states (TDOS) are aligned to the PES spectrum while the alignment of the two experimental scales corresponds to a F 1s binding energy of 694.6 eV.
Apart from the main features a satellite assigned to multiple excitations in the excitations step is found at 681 eV. Such satellites are abundant in electron-excited XES [10], and the first such excitations are allowed at our highest excitation energy (751.4 eV). This satellite disappears as the excitation energy is lowered below the threshold for multiple excitations.
While the electronic final states in PES and XES are the same, relative shifts and different broadening of PES and XES features can be attributed to different vibrational dynamics. We observe that the XES features are broader than the corresponding PES features. As the S-F bond has ionic character, it is expected to weaken by the removal of a F 1s electron, an expectation that is supported by DFT, as seen in the molecular dynamics simulations presented below. Coupling to anti-symmetric vibrational modes may lower the electronic symmetry and localize the core hole to one of the fluorine atoms. In case of substantial nuclear rearrangement the scattering process may be regarded as incoherent in the sense that the original inversion symmetry becomes irrelevant for the interpretation. If the XES signal emanates from a fluorine atom with a prolonged S-F bond in a distorted molecule which has lost the O h symmetry, the final states would represent holes in valence orbitals that to some degree are localized at this specific fluorine atom, rather than electronically inversionsymmetric orbitals. Picturing molecular orbitals as a linear combination of atomic orbitals a prolongation of the S-F distance implies decreasing overlap between atomic orbitals and consequently a decreasing energy spacing. However, although our molecular dynamics simulation (see inset in figure 3) confirms that the bond associated with the ionized atom is prolonged (reaching 2 Å after only 10 fs) upon core ionization, the reduced spacing in XES as compared to PES is related to changes in intensity due to rehybridization rather than reduced covalent splitting of bonding and anti-bonding orbitals. The symmetry assignment of the calculated XES (t = 0) is presented at the bottom of figure 3. In terms of local symmetry, the 5a 1g , and 3e g are of pure σ symmetry, whereas 4t 1u and 5t 1u are dominated by σ and π, respectively. The other features 1t 2g , 1t 2u and 1t 1g have π symmetry. The change in spectral shape is followed along the core-ionized state dynamics. For the geometry at 10 fs, many features have shifted and lost most of their intensity due to the elongation (breaking) of the S-F bond to the core-ionized F atom and only one low energy σ and two degenerate high-energy π features remain. At full dissociation into SF 5 and F fragments and neglecting spin-orbit coupling, these will merge into a single feature due to the spherical symmetry to the F atom. At long time evolution 6-10 fs, we observe a narrowing of the spacing between the σ and π features as a consequence of the reduced overlap between S and F. In figure 3 we see that overall the spectral features both change in intensity and red-shift in emission energy as a function of distortion. In this classical analysis of dynamical effects in XES [25], it is assumed that the measured spectrum comprises features due to transitions at various times during the life of the core hole. As a consequence, the features associated with the σ bonding orbitals in the lowenergy part of the spectrum are due to early transitions at relatively short S-F distance, whereas the main feature has a large contribution from later transitions at relatively long S-F distance.
In figure 4, the x-ray absorption spectrum is presented in panel figure 4(a) together with the RIXS at different resonances in panels figures 4(b)-(f). To begin with we start by giving attention to the excitation-energy dependence when the energy is tuned to the shape resonances above the ionization limit, corresponding to unoccupied orbitals of t 1u , t 2g , and e g symmetry at 694.6 eV, 699.8 eV, and 713.1 eV, respectively (figures 4(c), (e), (f)). Apart from the expected elimination of the high-energy satellite the spectra are very similar to the high-energy excited spectrum, indicating that the excited electron couples only little to rest of the electron system. This is observed in spite of the fact that the final states assigned to the features at 672 eV and 675 eV have different parity, 4t 1u and 1t 2g , respectively. If inversion symmetry is preserved one expects a relative enhancement of the gerade final state upon resonant excitation of a shape resonance associated with at gerade orbital, i.e., excitation of e g (713.1 eV), and t 2g (699.8 eV) would open primarily ungerade core holes, and emission to the gerade 1t 2g final state would be enhanced, relative to the ungerade 4t 1u state. Conversely, excitation at t 1u (694.6 eV and 696.3 eV) enhances emission to 4t 1u relative to the gerade 1t 2g final state. We see that this indeed happens in figure 4(d). This implies that the original inversion symmetry does influence the electronic transitions in spite of the coupling to symmetrybreaking vibrational modes. Note, however, that this is a small effect, and one can conclude that the gerade/ungerade selection rule applies only in propensity form.
There is a significant polarization dependence of the relative spectral intensities (figures 4(b), (c), (e), (f)), which can be qualitatively described in both the local and the global symmetric pictures. In the symmetric picture, e.g. excitation of an electron to a gerade a 1g valence orbital can only open a t 1u core hole. With horizontal polarization this implies a combination of core holes primarily on the left and right fluorine atoms, whereas vertical polarization implies a combination of core holes on the up and down atoms. The same site selectivity is  [22]. Note that the spectrum deviates from an absorption spectrum due to saturation effects, and is used only for reference. Panel (b), (c), (e), (f): RIXS spectra excited with vertical (red) and horizontal (blue) polarization of the incident radiation, at the excitation energies marked with vertical lines in panel (a). A small but significant polarization dependence is observed in panel (e) and (f), while a dramatic effect is observed in panel (b). Panel (d) highlights the excitation-energy dependence of the two well-resolved features that correspond to orbitals of different parity. expected in a localized picture, considering that the a 1g orbital has local σ symmetry (see figure 2). Here, excitation of a F 1s electron to σ orbital opens a nodal plane perpendicular to the polarization direction, implying that primarily σ orbitals along the bonds in the same direction as the polarization are excited. Horizontal polarization thus primarily excites core holes on the left and right fluorine atoms, whereas vertical polarization excites up and down atoms.
Corresponding considerations apply for the emission. While σ type orbitals radiate primarily perpendicularly to the bond, transitions between orbitals of local π and 1s orbitals radiate also in the bond direction.
Therefore, we expect that in our experimental geometry, horizontal polarization and excitation to resonances of local σ symmetry relatively enhances emission from orbitals of local π symmetry. Analogously, vertical polarization and excitation to resonances of local σ symmetry relatively enhances emission from orbitals of local σ symmetry. The situation is illustrated in figure 5. Similar simplified consideration can be made for Figure 5. A sketch of the π/σ asymmetry during breakup. As the first resonance has local σ character horizontal polarization excites the blue atoms left and right (a), so that emission from π orbitals are emphasized in the emission (b). Vertical polarization excites the blue atoms up and down (c) so that also σ orbitals are also observed in the emission (d).
excitation to orbitals of local π symmetry. This primarily selects top and bottom fluorine atoms as well as the front and back atoms. The angular anisotropy is less in this case as both π and σ orbitals radiate perpendicularly to the bond. A rigorous treatment of polarization dependence in RIXS can be found in reference [5].
The observed trends in the polarization dependence (figures 4(b), (c), (e), (f)) are fully in line with this expectation. We recall that the main peak at 678 eV is attributed to final states of predominantly π symmetry, whereas the features below 675 eV are associated with σ bonding states. At 713.1 eV (figure 4(f)) the σ shape resonance is excited, and π emission is relatively enhanced at horizontal polarization. At 699.8 eV (figure 4(e)) the π shape resonance is excited, and π emission is relatively enhanced at vertical polarization. At 694.6 eV (figure 4(c)) the t 1u resonance comprises both π and σ symmetry, and there is no clear polarization dependence at the present level of accuracy.
We can conclude that, while the global inversion symmetry is broken during the scattering process, and gives a slight influence on the spectra, the polarization dependence can be understood in terms of local π/σ symmetry.
At 688.9 eV (figure 4(b)) the first unoccupied a 1g orbital is excited. This spectrum takes a form which is dramatically different from spectra excited at higher excitation energies. Both the main peak, corresponding to the nominally non-bonding orbitals of local π symmetry, and the feature corresponding to the σ bonding t 1u orbital becomes sharper, and especially a narrow peak is formed just above 676 eV emission energy. From the electronic decay channel it is well known that ultrafast dissociation occurs upon resonant excitation at the first resonance [14][15][16][17], and it is tempting to assign this peak as emission from a free fluorine atom. We note, however, that there is a strong polarization dependence. As a free fluorine atom is spherically symmetric, with a vacancy in the symmetric 1s orbital and a full n = 2 shell, we expect no angular anisotropy in the radiative decay. According to the discussion above, the polarization dependence of relative intensities is due to the selection of atomic centres; left/right for horizontal polarization and up/down for vertical polarization. Together with the symmetry of the free atom, this picture does not let us expect any polarization dependence. We therefore conclude that the core hole decays during the dynamic dissociation process, and primarily prior to complete separation of the fluorine atom from the molecule.
Then the polarization can be readily understood if the π/σ symmetry is considered to be retained. Horizontal polarization preferentially excites the left and right atoms (blue atoms in figure 5(a)). In the geometry of our experimental set-up the radiation of the primarily non-bonding π orbitals is enhanced in the bond direction, which coincides with the direction of the spectrometer. These orbitals simulate 2p orbitals of a free fluorine atom (red atom in figure 5(b)) also before complete dissociation, and one may expect a corresponding sharp atomic-like peak in the emission spectrum. This simple consideration is consistent with the notion that the potential surfaces of a local core hole state and a state with a hole in a local non-bonding (π) orbital are similar, implying a sharp emission peak and only little vibrational broadening associated with the transition between the two.
At vertical polarization primarily up and down atoms are excited (blue atoms in figure 5(c)). Then emission from bonding σ orbitals are relatively enhanced in our geometry, and their bonding character is retained a larger atomic distances than the non-bonding orbitals (red atom in figure 5(d)). A hole in a bonding (σ) orbital typically has a potential surface which is significantly different from the core hole state and consequently transitions between the two are accompanied by substantial vibrational excitations. Consequently, the emission spectrum would have a less atomic character than the spectrum excited with horizontal polarization. This simple consideration captures the overall phenomenology. It implies that the sharp line in the spectrum is only quasiatomic, and that the decay principally occurs before complete dissociation.
The theoretical simulations (figure 3) confirm that there is interaction between the receding atom and the rest of the molecule, also after times longer than the scattering duration, manifested in a σ/π splitting. Although this splitting is not observed in the experimental spectra, the interaction is reflected in the polarization anisotropy.

Conclusion
We have presented and analyzed RIXS spectra of SF 6 excited in the vicinity of the F K edge. Dynamic symmetry breaking is demonstrated in the observation of electronically dipole forbidden states, and the initial-state inversion symmetry influences the spectra only slightly. On the other hand, the local π/σ symmetry is crucial for the interpretation of the spectra. Small significant effects are observed for excitation on the sharp shape resonances above the ionization limit, and a dramatic polarization dependence is observed upon excitation to the dissociative first resonance. We attribute this to radiative decay during the dissociation process, where emission from non-bonding orbitals appears atomic while σ orbitals still interact with the molecular fragment.