Nacer
Boudjemia
*a,
Kari
Jänkälä
a,
Tatsuo
Gejo
bc,
Kiyonobu
Nagaya
bd,
Kenji
Tamasaku
b,
Marko
Huttula
a,
Maria Novella
Piancastelli
bef,
Marc
Simon
be and
Masaki
Oura
b
aNano and Molecular Systems Research Unit, University of Oulu, P.O. Box 3000, 90014 Oulu, Finland. E-mail: nacer.boudjemia@oulu.fi
bRIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5148, Japan
cGraduate School of Materials Science, University of Hyogo, Kamigori-cho, Ako-gun, Hyogo 678-1297, Japan
dDepartment of Physics, Kyoto University, Kyoto 606-8502, Japan
eSorbonne Université, CNRS, Laboratoire de Chimie Physique-Matière et Rayonnement, LCPMR, F-75005, Paris, France
fDepartment of Physics, and Astronomy, Uppsala University, SE-75120 Uppsala, Sweden
First published on 11th February 2019
Hard X-ray electron spectroscopic study of iodine 1s and 2s photoionization of iodomethane (CH3I) and trifluoroiodomethane (CF3I) molecules is presented. The experiment was carried out at the SPring-8 synchrotron radiation facility in Japan. The results are analyzed with the aid of relativistic molecular and atomic calculations. It is shown that charge redistribution within the molecule is experimentally observable even for very deep levels and is a function of the number of electron vacancies. We also show that the analysis of Auger spectra subsequent to hard X-ray photoionization can be used to provide insight into charge distribution in molecules and highlight the necessity of quantum electrodynamics corrections in the prediction of core shell binding energies in molecules that contain heavy atoms.
The main aim of the present work is to investigate whether the molecular environment plays a role in such deep ionization processes, and to assess whether widespread concepts such as electronegativity and charge distribution inside a molecule extend down to very deep levels, or should the levels be considered purely atomic in nature.
Photoelectron spectra of the 1s and 2s core levels of both molecules have been measured, and compared with Dirac–Fock theoretical calculations. Core-hole lifetimes have been extracted, and compared with previously reported values based on fluorescence measurements. To complement the information from photoelectron spectra, Auger LMX spectra have also been measured. As in the previous work on Kr5 and Xe,6 the comparison of the LMX decay above and below the K-edge is informative about the importance of radiative versus nonradiative decay.
A key point in the theoretical description we use has been to perform calculations for the molecular systems, for the iodine isolated atom, and for the iodine negative ion (isoelectronic with xenon). The rationale is to assess the importance of charge redistribution within the molecule even for deep core excitation and relaxation processes.
As results we show that there are measurable chemical shifts for both 1s and 2s ionization energies between the two molecules (0.7 eV for the 1s level and 1.3 eV for the 2s level). To our knowledge, this is the first time that chemical shift is observed in the ionization of such deep levels. An interesting result of the theoretical simulations including a neutral iodine atom and a Xe-like iodine negative ion is that the best agreement with the experimental ionization energies is obtained when the neutral atom is considered. A charge distribution analysis indicates that while the difference in charge withdrawing by the fluorine atoms compared to the hydrogen atoms is strong, this effect is compensated by charge redistribution on the carbon site, so that the iodine atom is neutral in both cases. At variance with this result, in the analysis of the Auger spectra we show that the negative charge on the iodine site increases significantly when the total charge is 2. In other words, charge redistribution within the molecule is a function of the number and of the depth of the electron vacancies.
The energy resolution of the analyzer was 781 meV and 313 meV for the two pass energies, 500 eV and 200 eV, used in the high and low photon energy measurements, respectively. The photon energy bandwidth was estimated to be 4.83 eV at 35.46 keV and 1.30 eV at 10.0 keV. The energy scale of the photon beam at 35.5 keV was calibrated by using the 1s photoelectron line of Xe6 and at 10.0 keV using the 2s and 2p photoelectron lines of Xe atoms.8 The kinetic energy scale of the spectrometer was calibrated with the aid of the well-known KLX Auger electron spectrum of Ar9 and LMM spectrum of Xe.10,11
The molecular states were further studied by calculating the partial charges at different sites. This was done by performing a natural population analysis (NPA)21 with the JANPA22 package using the CCSD molecular orbitals from the ORCA calculations. Partial charges were calculated using the NPA method due to its reliability in the case of bonds with high ionic character.21
The DIRAC code was used to obtain the absolute 1s and 2s ionization energies of I in CH3I, CF3I, I and I− forms within the Δ-SCF framework using the geometries obtained from ORCA calculations. For these calculations the Dyall's triple-ζ basis set23 was adopted. The calculations were carried out at two levels, using the bare Dirac–Coulomb Hamiltonian and including the Gaunt interaction.24,25
In addition to molecular calculations, single-atom configuration interaction Dirac–Fock simulations were performed. The results are used to check the quality of the molecular calculations, to obtain quantum electrodynamics (QED) corrections not present in the molecular codes and to analyze the experimental LMX Auger spectra subsequent to K and L photoionization, and KL fluorescence decay in the binding energy region of the iodine K-shell of CH3I and CF3I molecules. The atomic wave functions were calculated using the relativistic GRASP2K26 program. The configurations included for the calculation of the Auger decay spectra were selected by the same procedure as presented in ref. 10. The Auger spectra and photoionization cross sections were modeled by using the RATIP27 package, and fluorescence decay rates were obtained from the REOS28 program using a length gauge.
From the fits we extracted lifetime broadening of the iodine 2s core–hole states in CH3I and CF3I molecules. The values are 2.71 ± 0.13 eV, and 2.85 ± 0.14 eV, respectively, which corresponds to the lifetimes of 241.1 ± 11 as and 230.6 ± 10 as. The values differ considerably from the value of 3.46 eV29 reported for solid iodine, obtained using X-ray fluorescence measurements. The values are also notably smaller than those of the calculated widths using GRASP2K in the present work for atomic iodine (3.24 eV) and the Xe-like iodine anion (3.26 eV), which also indicates that the chemical environment should not affect the lifetime of the 2s hole of iodine considerably. It is emphasized that in the present measurement the total full width at half-maximum (FWHM) of the 2s lines were about 3.4 eV, which combined with the known photon bandpass of about 1.3 eV and analyzer broadening of about 0.8 eV means that the discrepancy cannot be explained solely by fitting uncertainties. The reason for this observed anomaly, however, remained unknown. We assume that the maximum of the peak corresponds to the GS → I 2s−1(ν0 = 0) transition and it is noted that the same discrepancy is also observed in the case of 2s ionization of atomic Xe, which means that the anomaly does not arise from molecular effects.
Fig. 2 presents the 1s photoelectron spectra of CH3I and CF3I measured at the photon energy of 35440 ± 5 eV. The spectra consist of a single line that corresponds to the direct 1s ionization of iodine and a satellite region at higher binding energy. In comparison to the 2s case, due to the larger lifetime broadening and photon bandpass, the satellite features are now seen as an asymmetric shoulder on the right-hand side of the main line. Due to the overlap and lack of precise knowledge about the shape of the satellite part, fitting of the 1s spectrum is somewhat complicated. In the present work the fitting was carried out by making an approximation that the relative position of the satellite spectrum with respect to the main line does not change between 1s and 2s ionization. The approximation can be justified by the fact that both 1s and 2s holes are very deep in the system and therefore relaxation of the valence orbitals in the final state is practically the same in both cases. Therefore the shoulders in the 1s spectra were fitted by fixing the relative positions of the four profiles obtained from the fit of the better resolved satellite part of the 2s spectrum. In addition a fifth profile centered at about 33260 eV was added to cover the far extending tail of the satellite part.
The shapes of the iodine 1s spectra of CH3I and CF3I are similar to the 1s spectrum of atomic Xe.5,6 The main line is accompanied by a satellite shoulder and a very broad tail that extends about 140 eV above the main line. Such similarity is expected since the iodine 1s ionized state of CH3I and CF3I is not dissociative and the large lifetime broadening smears out vibrational and rotational effects. However, after a careful analysis it is possible to gain information from the spectra, as will be discussed in the following.
The 1s main lines in Fig. 2 have FWHM of about 11.63 eV for CH3I and 11.69 eV for CF3I. Taking into account the experimental contribution of 4.84 eV, we are able to deduce the iodine 1s lifetime broadening to be 9.5 ± 0.6 eV for CH3I and 9.6 ± 0.6 eV for CF3I, corresponding to the lifetimes of 69 ± 4 as and 68 ± 4 as, respectively. The values are the same as those5 measured with a slightly different photon energy and well in line with the lifetime broadening of 9.6 ± 0.2 eV (68 ± 2 as) of the 1s ionized state of Xe.6 Interestingly, however, similar to the 2s case the lifetime broadening values obtained from the present experiment are smaller than the value of 10.6 eV obtained from fluorescence measurements of solid iodine,29 thus showing the same trend as seen in the case of 1s ionization of Xe.6 Comparison to the calculated broadening value of 10.3 eV using GRASP2K for atomic and Xe-like iodine also shows the same behavior as in the 2s case that the theory overestimates the width against the present measurement, but agrees well with the fluorescence data from ref. 29.
The results for the experimental 1s and 2s binding energies of iodine in CH3I and CF3I molecules are summarized in Table 1. The table also contains the results from theoretical calculations for I in the studied molecules, as well as in atomic neutral and Xe-like forms. The atomic calculations were carried out to monitor the quality of the chosen basis set for molecular calculations and to obtain estimates for QED corrections that, to our knowledge, are not currently implemented in commonly available molecular codes, namely Breit interaction, vacuum polarization and self energy, which make an important contribution to deep core binding energies of heavy atoms (see, e.g., ref. 30–32 and references therein).
Iodine 1s−1 | CH3I | CF3I | I | I− |
---|---|---|---|---|
a Corrections from atomic calculations. b Value for solid iodine from ref. 8. | ||||
DIRAC | ||||
D–C | 33293.7 | 33295.0 | 33295.2 | 33286.2 |
+G | 33210.1 | 33211.4 | 33211.6 | 33202.4 |
+B, VP, SEa | 33177.0 | 33178.3 | 33178.5 | 33169.5 |
GRASP2K | ||||
D–C | — | — | 33295.4 | 33286.2 |
+B | — | — | 33217.8 | 33208.5 |
+B, VP, SE | — | — | 33178.7 | 33169.5 |
Experiment | 33175.2 | 33175.9 | — | 33169.69(89)b |
Comparison of the atomic calculations with bare Dirac–Coulomb Hamiltonian (lines D–C) between DIRAC and GRASP2K codes shows a remarkably good agreement. The results for both I and I− are within a maximum of 0.2 eV from each other, which shows that the selected basis set (Dyall's triple-ζ) is of sufficient quality for the present task. Comparing the addition of Gaunt correction in DIRAC against the full Breit interaction (lines +G and +B) in GRASP2K shows that in the case of 2s ionization the two differ by 0.5 eV, but the difference increases drastically to 6.1 eV in the case of 1s ionization. We thus arrive at the conclusion, which is the same as the example in ref. 33, that even though the magnetic part (i.e. Gaunt) dominates the Breit interaction, the retardation term needs to be included in the calculation of core ionization energies of heavy atoms. The Dirac–Coulomb–Breit Hamiltonian is, however, still not sufficient, because the QED terms vacuum polarization and especially self energy become significant (see, e.g.ref. 6, 30 and 33). In the present case for 2s the combined effect of the two terms is about 5.1 eV and for 1s it is 39.1 eV, which are far above the experimental resolution. As already mentioned, the considered QED terms are not available even in the relativistic DIRAC code. Therefore we made an assumption that the terms do not change much when moving from an atomic to molecular environment, which is justified by the fact that the corrections are practically the same between I and I−. The results of this approximation are summarized in line “+B, VP, SE” of Table 1, where the total shift arising from full Breit interaction, vacuum polarization and self energy is taken from the atomic calculation to correct the bare Dirac–Coulomb result of the DIRAC calculations. This provides theoretical binding energies that are remarkably close to the experimental binding energies, in the case of 1s the difference is about 2 eV and in the case of 2s about 4.5 eV. These results are notably better than the ones obtained directly from DIRAC with a Gaunt correction that deviates only for 1s and 2s by about 35 eV and 9 eV from the experimental values, respectively. We may thus conclude that even though it seems that one may transfer the leading QED corrections from atomic to molecular calculations, there is a clear future need to implement them directly into the relativistic molecular codes.
The experiment shows that the chemical shift between CH3I and CF3I is from the 2s spectra −1.3 ± 0.05 eV and from the 1s spectra −0.7 ± 0.6 eV, when the calculation gives −1.4 eV for 2s and −1.3 eV for 1s. Deviation between the experiment and theory in the 1s case is in the region of error bars due to the large width of the photoelectron line and overlapping satellite region. Nevertheless, it is notable that the chemical shift is observed in the 1s case, which to the best of our knowledge makes the present data the highest energy photoelectron spectroscopic measurement where a molecular chemical shift is observed. We emphasize that the knowledge of the chemical shifts transferred all the way to the deepest orbitals of heavy atoms can prove to be useful in the site-specific study of complex systems that contain light and heavy atoms. This is because the photon energy required to ionize iodine 1s is so far above all thresholds of the second row elements that materials containing mainly these elements are almost transparent to hard X-rays.
Comparison of the chemical shifts between CH3I, CF3I, I and I− shows that the molecular values are surprisingly close to I rather than I− and that the chemical shift between CH3I and CF3I is fairly small, only about 1.3 eV, despite the fact that H atoms are electropositive and F atoms are highly electronegative. Also since iodine is electronegative, one might intuitively assume that in a molecular environment the atom would appear as a Xe-like anion. The reason for the behavior can be, however, understood from the distribution of partial charges in the molecule.
Table 2 shows the results of natural population analysis carried out for the molecules in the original state and by increasing the nuclear charge of iodine by one and two, which mimics the effect of removal of one or two core electrons, which will be discussed in the next chapter. In the present work, the neutral case, i.e. column “I(Z = 53)”, shows that, as expected, there is a huge difference in the partial charges in the H and F-sites. However, the charge is compensated by a change of partial charge in the C-site so that in CF3I the carbon loses almost a unit charge whereas in CH3I the carbon gains charge. The change is, however, not reflected in the I atoms because the CH3 and CF3 units as a whole are almost neutral. The small chemical shift between the two molecules can be understood by considering the C bridge acting as a charge buffer, which nullifies most of the effect arising from the change of H atoms to F atoms. This can also be seen from the calculated 1s binding energies of carbon which are 291.58 eV in CH3I and 299.46 eV in CF3I, thus giving a chemical shift of 7.9 eV. We notice that a comparably small chemical shift (1.05 eV) has been reported in the literature between the I 4d3/2 levels for the same two molecules,34 confirming our overall interpretation based on charge distribution. It also indicates that even though in the present case the iodine 1s and 2s chemical shifts are quite small, it should not be considered as a general rule.
I(Z = 53) | I(Z = 54) | I(Z = 55) | |
---|---|---|---|
![]() |
−0.67 | −0.50 | −0.36 |
C![]() |
0.20 | 0.25 | 0.30 |
CH3![]() |
0.08 | −0.24 | −0.54 |
![]() |
0.95 | 1.06 | 1.13 |
C![]() |
−0.36 | −0.29 | −0.20 |
CF3![]() |
0.13 | −0.18 | −0.52 |
The observed changes in the intensities between above and below the 1s threshold are highly similar to the same case in atomic xenon, discussed recently in ref. 6. The variations can be understood by the fact that below the 1s threshold the population of the initial states of the Auger decay is simply governed by the branching ratio of the 2s, 2p1/2 and 2p3/2 photoionization cross sections at the selected photon energy, but above the 1s threshold, population of 2p hole states increases dramatically because of dominant KL2,3 fluorescence decay. The consequence of this can be seen as an increase in the intensity of the L2,3MX Auger lines with respect to L1MX lines.
As expected, also the overall shapes of the measured Auger spectra are very similar to atomic Xe studied in ref. 6 and 10. In comparison to Xe, the whole spectrum is shifted by about 150 eV to smaller kinetic energy, but otherwise one can locate precisely the same spectral features and roughly the same relative intensities. The reason is that the decay takes place between very deep atomic-like orbitals, and photoionization does not directly affect the valence structure of the systems. Therefore the initial and final states of the decay remain as the bonding state. Moreover, due to the large lifetime broadening, all vibrational effects smear out. Since the LMX Auger spectrum of Xe has been recently analyzed in ref. 10, we do not repeat line identifications here. An interested reader can easily label the spectral features using the Xe case in Fig. 1 of ref. 10.
Comparison of the experimental Auger spectra of CH3I and CF3I molecules in Fig. 3 shows remarkable similarity. Line positions and intensities show no statistically significant differences in the present experimental resolution and noise level. The difficulty in finding differences is at least partly due to large natural width combined with considerable line overlap, which makes the study of individual Auger decay lines challenging. However, even the most isolated lines in the spectra appear to be almost identical in both molecules.
Comparison of the experimental spectra against the theory shows an interesting behavior. The intensities of both calculations, Xe-like and pure iodine, are in very good agreement with the experiment. However, the calculated line positions are found to be located in such a way that the pure iodine spectrum is at lower and Xe-like spectrum at higher kinetic energy, and the experimental spectrum is roughly in the center between the two cases. It should be noted that it is not likely that the behavior arises from computational uncertainties because in ref. 6 and 10 it was shown that in the case of atomic Xe the method predicts the positions of all line groups of the spectrum exceptionally well. The observation is in clear contrast to the photoionization results presented in Table 1, where in both molecules it is seen that the iodine site appears to be almost charge neutral “iodine-like” in the ground state. As in the case of photoionization, the behavior can be qualitatively understood by studying the partial charges in the molecules. The two rightmost columns of Table 2 show that in singly and doubly core ionized cases the iodine site gains a charge of about −0.2 e and −0.5 e, respectively. It thus means that when the core of the iodine site becomes more positive, it is able to attract more negative charge from the rest of the molecule and the chemical environment in the valence region of the site becomes more Xe-like, which in turn is reflected in the kinetic energies of the Auger lines. The values from Table 2 are well in line with the atomic calculations in Fig. 3 showing that core ionization changes the valence charge distribution in the molecule significantly, but even doubly core ionized core is not sufficient to make the site completely Xe-like in valence.
Theoretical calculations have been performed for the molecules, atomic iodine and Xe-like I− (negatively charged iodine). We show that for the 1s and 2s ionization energies the experimental values are in better agreement with theory when atomic iodine is considered, implying that the charge transfer from the CX3 moiety to the iodine center is negligible for deep core ionization. At variance with this result, the calculated energy values for the Auger transitions are in-between atomic iodine and Xe-like iodine. This in turn indicates that the charge distribution inside the molecules, and namely the capability of the electronegative iodine atom to draw charge from the rest of the molecule, is a detectable effect when the species are doubly charged (as in Auger decay) even for levels as deep as the ones involved in the LMX Auger transitions. This conclusion is supported by calculated charge distributions within the molecules.
This journal is © the Owner Societies 2019 |