invited papers
Theoretical studies of core excitation and ionization in molecular systems
aDepartment of Chemistry, Texas A&M University, PO Box 30012, College Station, TX 77842-3012, USA
*Correspondence e-mail: lucchese@mail.chem.tamu.edu
ab initio techniques and using a numerical representation of the continuum scattering orbitals. Here, core ionization of small molecular systems will be focused on. A brief review of previous results obtained using these methods is given. New results are then given for the N 1s ionization of NO and the N 1s ionization of N2. In the study of NO, the spin-coupling scheme which best describes the shape-resonant state and the energy dependence of the branching ratios for the singlet and triplet final ion states will be considered. In the study of N2 core ionization, the effects of the coupling of nearly degenerate core levels and the effects of correlation on the importance of interchannel coupling will be considered.
cross sections have been computed using correlated initial and final target states obtained with standardKeywords: molecules; photoionization; core ionization.
1. Introduction
Much progress has been made in the ability to compute the ) which we have used to study a number of core ionization processes. To consider larger systems, simpler treatments of the correlation effects are usually employed.
cross sections of molecules. For linear molecules the azimuthal symmetry reduces the difficulty of the calculation. Thus for these systems it is possible to perform very accurate calculations including both correlation in the initial and ionized target states and correlation of the photoelectron with the target electrons. In this paper we will discuss results obtained with the multichannel Schwinger method (MCSCI) (Stratmann & Lucchese, 1995There are a number of interesting resonance effects which occur at energies near the threshold for core ionization of a molecule (Chen et al., 1989). Just below the threshold for ionization there are transitions from the core level to valence and Rydberg levels of the molecule. These states then have a number of modes of decay. The excited states can decay by emission of light, by an Auger-type decay process, and by fragmentation. The Auger process is known as a resonant Auger decay. Furthermore, there are two modes of in the resonant Auger process. The first mode is a participator-type decay where the electron in the valence or Rydberg level participates in the decay process; the second mode is a spectator-type decay where the initially excited electron remains in its and the other electrons take part in the decay process. In our study on the resonant Auger ionization of CO due to the C 1s → 2π* excitation, we found that the MCSCI approach can obtain the width of the resonant state and can predict the degree to which the decay process is of the participator or spectator type (Botting & Lucchese, 1997).
When the photon energy is greater than the threshold for the core ionization process, the resonances that occur are qualitatively very similar to those found in ) or N2 (Lucchese et al., 1982) in which there are prominent σ*-type resonances which are trapped by an l = 3 angular momentum barrier. However, we have found in the of S 1s in SF6 that there is a resonance of t1u symmetry which is not similar to any of the unoccupied valence orbitals of that system and is characterized by an l = 9 angular momentum barrier (Natalense & Lucchese, 1999).
from the valence shell. One finds one-electron resonances that are due to angular momentum barriers. These resonances are usually fairly broad when compared with two-electron resonances. Often these shape resonances resemble unoccupied valence antibonding orbitals as in the of CO (Lucchese & McKoy, 1983One difference between valence ionization and core ionization is that generally in the valence region there are several valence ion states that are well described by single-hole states. At higher energies in the valence region there are usually shake-up states that can be characterized by two vacancies in the bound orbitals of the initial state and one virtual orbital occupied. These states often strongly mix with the inner valence single-hole states, leading to strong configuration mixing effects that can complicate the interpretation of the ). If the atoms are not equivalent, then there will be a more substantial splitting between the core levels.
cross sections. In the region of core ionization there are three possible cases which can occur. If there is only one atom in the molecule of a particular type, then in the region of the core ionization of that atom there will be only one hole state. If there are two or more atoms of a particular type, and if the atoms are equivalent, there will be a number of degenerate or nearly degenerate hole states. In these cases, there can be interesting coupling effects between the hole states and the vibrational motion of the molecule. One such effect is dynamic hole localization (Domcke & Cederbaum, 1977In a previous study of the C 1s of acetylene (Lin & Lucchese, 2000a), we have considered the case of chemically equivalent atoms. In that study we considered the effects of electron–electron correlation on the description of the σ* shape resonance which occurs at about 20 eV above the C 1s−1 threshold. It was found that, without the inclusion of correlation in the description of the target states, the effects of interchannel coupling were exaggerated. We further found quantitative agreement with experimental cross sections (Thomas et al., 1999) for the resolved 1σg-1 and 1σu-1 partial cross sections with a calculation including the two hole states and the four strongest shake-up states.
2. Computational methods
The method we have used to study the various et al., 1986; Bandarage & Lucchese, 1993; Stratmann & Lucchese, 1995; Stratmann et al., 1996) so we will only give a brief outline of the method here. In the MCSCI method, the photoionized state of a molecule is described by a wave function of the form
processes is the MCSCI method. In this method, the close-coupling approximation is used with the target states represented by (CI) wave functions. This leads to a set of coupled integro-differential equations that are solved using the Schwinger variational method. This method has been described in detail elsewhere (Lucchesewhere is the one-electron scattering wave function of channel p. represents the CI target ion wave function of channel p, which is a linear combination of spin-adapted N − 1 state functions (CSF) . and are antisymmetrized spin- and symmetry-adapted products of the corresponding N − 1 electron functions and the continuum one-electron function .
The wave function satisfies the time-independent Schrödinger equation
with the appropriate asymptotic boundary condition and where HN is the N-electron Hamiltonian and E is the total electronic energy. Given the molecular orbitals {} and the CI wave functions of the residual ion of different channels, we are able to obtain the coupled variational equations for the scattering wave functions,
where p, q are the channel indexes,
∊p is the asymptotic of the photoelectron in channel p, with = E-Ep, where Ep is the electronic energy of the residual ionic state in channel p, and Upq is the non-local optical potential (Bandarage & Lucchese, 1993).
The Lippmann–Schwinger equation is then
where and are the vectors of the channel-scattering wave function and Coulomb wave function, respectively. The potential in equation (5) is a Phillips–Kleinman pseudopotential (Lucchese et al., 1982; Bandarage & Lucchese, 1993), which enforces the orthogonality between a given scattering wave function and selected molecular orbitals. In (5), is the multichannel Coulomb Green's function matrix given by
The variational functional for the p can be constructed as
transition moment leading to the target state of channelwhere the quantity RiL(V) is a function of the electronic dipole operator, and {} forms the multichannel scattering basis set. The superscript L(V) denotes that this equation may be solved using either the length or the velocity approximation. This multichannel Schwinger result is further refined using Padé approximates to IpL(V) (Lucchese & McKoy, 1983).
The double-differential
of a given channel is expressed in length (or velocity) and the mixed form asand
The final cross sections and the asymmetry parameters are obtained by integrating over all orientations of the molecule in the laboratory frame.
The calculation of the MOLPRO program (Werner & Knowles, 1996). However, for valence ionization we have also found that a state-averaged multi-configuration self-consistent-field (MCSCF) calculation on a selection of important ion states has yielded very good results (Wells & Lucchese, 1999). Molecular natural orbitals generated from this step are used by the MCSCI code to construct the CSFs and total electron wave functions.
begins with the calculations of the initial and target states. In order to simplify the calculation, we use the same set of molecular orbitals in the CI calculations of the initial and all final states. For the ionization of the core levels, these common orbitals are usually obtained by a complete-active-space self-consistent-field (CASSCF) calculation on the ground state of the molecule of interest using theIn the following steps, the single center expansion (SCE) technique and numerical grid method have been employed in all integrations.
Different levels of the MCSCI calculations are performed to investigate how the different types of correlation contribute to the ); however, two-electron resonance processes such as autoionization cannot be studied.
process. Calculations based on the single-channel Schwinger Hartree–Fock (SCSHF) approximation contain no correlation. Multichannel Schwinger Hartree–Fock (MCSHF) approximation calculations show the results of target response to the photoelectron without including target correlation. Single-channel Schwinger (SCSCI) calculations show the effects of target correlation. The full MCSCI calculation has all correlation effects included. In addition, by selectively including and removing channels it is possible to consider which channels are important in various resonant features. In core ionization we are particularly interested in the effects of interchannel coupling between core-hole channels and in the effects of interchannel coupling between core-hole and shake-up channels. In the study of linear systems, the full MCSCI calculation can be performed with up to ∼20 ion states included in the close coupling expansion. For non-linear systems, we are only able to treat the problem at the SCSHF level. At this level, interesting shape resonance effects can be considered (Natalense & Lucchese, 19993. Results
3.1. N 1s ionization of NO
In our study of the RN–O = 1.15077 Å. In the first step, a restricted Hartree–Fock (RHF) calculation is performed on the ion ground state using the MOLPRO program (Werner & Knowles, 1996), followed by an MCSCF calculation including all core and valence orbitals. Dunning's aug-cc-pVTZ basis set with augmented diffuse functions was used in these calculations (Dunning, 1989). The natural orbitals are constructed from the MCSCF calculation and serve as the molecular basis set in the following MCSCI calculation.
of NO, the equilibrium structure of molecular NO is used in the calculations, in whichIn the next step, initial molecular states and target ion states are obtained by full CI (FCI) calculation with the graphical unitary group approach (GUGA) (Hinze, 1981) using the natural orbital set. SCSCI and two-channel Schwinger (2CSCI) calculations were performed using the initial and target states calculated above. The continuum scattering orbitals are obtained by performing a scattering calculation on a model potential using a set of general Gaussian-type orbitals and selecting a compact basis set from the final eigenfunctions (Stratmann et al., 1996).
The ground state of NO has the configuration of 1σO22σN23σ24σ25σ21π42π1. The N 1s ionization of NO leads to two Π states with different multiplicity, singlet or triplet states. The calculated (IP) is 409.59 eV for the 3Π state, and 410.44 eV for the 1Π state. The experimentally determined IPs are 410.3 and 411.8 eV, respectively (Siegbahn et al., 1969). The SCSCI and 2CSCI calculations produce very similar shapes of the total and partial cross sections. The partial cross sections from the 2CSCI calculation are plotted in Fig. 1(a). Sharp resonances occur in both channels but at slightly different positions. The maximum of the is 413.5 eV in the channel leading to the 1Π state and 414.5 eV in the channel leading to the 3Π state. These two shape resonances are also observed in the measured cross sections for the two channels, but are broader than those we obtained (Ruedel & Hergenhahn, 2000). By comparison with the single-channel scattering calculations, we observed little interchannel coupling between these two channels. The intensity ratio between these two channels is plotted as a function of photon energy in Fig. 1(b). In the high-energy limit, the intensity ratio 3Π/1Π is slightly higher than 3, which matches the spin statistic ratio. Near the threshold, however, the intensity ratio changes dramatically. The minimum of the ratio is ∼0.88, located at 413 eV, just around the first shape resonance. The maximum of the ratio is ∼5.24, located at 415 eV, close to the shape resonance in the channel leading to the 3Π state. Experiments have verified the existence of the structure in the measured intensity ratio, but the oscillatory structure near the threshold is much less significant.
The shape resonances have been verified as σ*-type resonances. Compared with previous calculations by Wallace et al. (1982), the positions of the shape resonances from our calculation are much closer to the threshold, and the intensity ratio is also different from the spin statistics that they used. For the N 1s core ionization of NO, the couplings between core and valence orbitals are strong enough that a multiplet-averaged approach is not feasible. Our results also indicate that the couplings between continuum scattering wave function and valence orbitals are much less significant than the core-valence coupling. Thus, the shape resonances in the photoelectron scattering from the singlet or triplet ion potentials are observed in different positions. This is in contrast to the behavior of the corresponding resonance in the 5σ−1 valence ionization of NO (Stratmann et al., 1996). In the valence region, the 5σ–kσ coupling is much stronger than the N 1s–kσ coupling. The consequence of this stronger coupling is that in the valence region there is only one strong shape resonance feature in the process. This resonant state then decays into the two possible asymptotic 3Π and 1Π channels. In the partial cross sections the resonance then appears at the same energy in each channel.
The core ionization of NO is a unique system. We observed that the MCSCI calculation for this system is very sensitive to the selection of the molecular basis set. Our calculations also showed that the cross sections are very sensitive to the N—O bond distance (Lin & Lucchese, 2000a). This geometry dependency of the may be the cause of the difference between theory and experiments.
In our current calculations, we excluded the many possible shake-up channels in the close-coupling expansion. However, we believe that the couplings between the core-hole states and shake-up channels cannot be neglected. The exclusion of the shake-up channels in the expansion may also be the cause of the disagreement.
3.2. N 1s ionization of N2
The calculation of N 1s of N2 is performed at its equilibrium geometry, in which RN–N = 1.09768 Å. An RHF calculation is first performed for the ground state of the N2 molecule. A CASSCF calculation including full valence orbitals is then performed, followed by a multi-reference (MRCI) calculation. A set of natural orbitals is generated after the MRCI calculation. The following natural orbitals are selected to form the set in the selected CI calculation,
Two restrictions are applied in the selected CI calculation, i.e. only two electrons are allowed in the last six orbitals, and only one electron is allowed in the last four orbitals. Total and partial cross sections are obtained by coupled-channel scattering calculations.
The 2 is 1σg21σu22σg22σu23σg21πu4. The N 1s ionization in N2 includes both 1σg-1 and 1σu-1 ionization channels. The calculated IP for N 1s of N2 is that of the 1σu-1 state, which equals 413.054 eV. The measured IP (Chen et al., 1989) is 409.938 eV. The calculated splitting between two core-hole states, 1σu-1 and 1σg-1, is 0.144 eV. The small splitting makes it difficult to distinguish their contributions to the N 1s However, the two channels did have very different dynamics in the core process due to their symmetry properties.
of the ground state of NThe total s ionization and the partial cross sections of the two channels are plotted in Fig. 2, together with the total measured by a high-resolution study (Kempgens et al., 1996). The experimentalists found a shape resonance at ∼419 eV. Many studies using a multiple-scattering model (MSM) (Dehmer & Dill, 1975), or the Stieltjes–Tchebycheff moment theory (STMT) (Rescigno & Langhoff, 1977), or the Schwinger method with frozen-core Hartree–Fock including core relaxation effects (RCHF) (Lynch & McKoy, 1984), were not able to reproduce both the magnitude and position of this resonance. Our calculation found the shape resonance located at about 420 eV with the maximum of about 2 Mb, which agrees very well with the experiments. Further investigation found that this is a σu*-type resonance, which is found in the components with a total Σu symmetry after The main contribution comes from the channel leading to the 1σg-1 ion state. By looking only at the components with final Σu symmetry, the single-channel calculations show that the contribution from the channel of the final 1σu-1 ion state is about 10% of the channel of the final 1σg-1 ion state at the peak of the shape resonance, and, in the coupled-channel calculation, the ratio only slightly increases. The result apparently suggests that there is little interchannel coupling effects between these two channels, which contradicts the recent random phase approximation (RPA) results that suggest strong coupling between these two channels (Cherepkov et al., 2000). We also compare the calculations using the RCHF approximation and using the CI approximation. The coupled-channel calculations using the RCHF approximation also suggest much stronger interchannel coupling effects (Lin & Lucchese, 2000b). In the study of C 1s of C2H2, we have already found that neglecting the correlation of the target ion states leads to exaggerated interchannel coupling effects in the coupled-channel scattering calculation (Lin & Lucchese, 2000a). Our results in the N 1s study of N2 suggest that only the simple treatment of the core relaxation effect is not sufficient in such a highly correlated system.
of the N 1New experimental developments also enable the measurement of the angular distribution of the photoelectron in the fixed-in-space molecular frame, which provides extra data and information about the dynamics of the core et al., 2000). Theoretical studies are in progress to calculate such angular distributions, and also to include the shake-up channels in the close-coupling expansion. Thus, we will be able to better understand the effects in the process, and to study the interchannel couplings between core-hole states and shake-up channels.
process (Cherepkov4. Conclusions
In this paper we studied the N 1s of NO and N2. There are several significant differences between these two processes. First of all, N2 is a close-shell molecule while NO is an open-shell molecule. The N 1s ionization of NO leads to target ion states with different multiplicities, and NO will have more channels leading to the target states due to different spin couplings. Also, N2 is a homonuclear molecule; both N atoms contribute to the N 1s Thus, there are more channels coming from the initial state of N2. The two molecules have different symmetry properties, which lead to different coupling and transition rules. Their spectra also show distinct and interesting characters of the dynamic processes.
Besides these differences, we also observed some similarities in the N 1s Both the N2 and NO molecules have shape resonances located near the threshold, which have the character of the antibonding orbital. We also find that there are only small couplings between channels of different final core-hole states when properly correlated ion states are used in the calculations.
Acknowledgements
We acknowledge support of this work by the Welch Foundation (Houston) under grant number A-1020. We also acknowledge the support of the Texas A&M University Supercomputing Facility.
References
Bandarage, G. & Lucchese, R. R. (1993). Phys. Rev. A, 47, 1989–2003. CrossRef CAS PubMed Web of Science
Botting, S. K. & Lucchese, R. R. (1997). Phys. Rev. A, 56, 3666–3674. CrossRef CAS Web of Science
Chen, C. T., Ma, Y. & Sette, F. (1989). Phys. Rev. A, 40(11), 6737–6740. CrossRef Web of Science
Cherepkov, N. A., Semenov, S. K., Hikosaka, Y., Ito, K., Motoki, S. & Yagishita, A. (2000). Phys. Rev. Lett. 84(2), 250–253. Web of Science CrossRef
Dehmer, J. L. & Dill, D. (1975). Phys. Rev. Lett. 35, 213–215. CrossRef CAS Web of Science
Domcke, W. & Cederbaum, L. S. (1977). Chem. Phys. 25, 189–196. CrossRef CAS Web of Science
Dunning, T. H. Jr (1989). J. Chem. Phys. 90(2), 1007–1023. CrossRef Web of Science
Hinze, J. (1981). Editor. The Unitary Group for the Evaluation of Electronic Energy Matrix Elements. Lecture Notes in Chemistry. Berlin: Springer.
Kempgens, B., Kivimäki, A., Neeb, M., Köppe, H. M., Bradshaw, A. M. & Feldhaus, J. (1996). J. Phys. B, 29, 5389–5402. CrossRef CAS Web of Science
Lin, P. & Lucchese, R. R. (2000a). J. Chem. Phys. 113(5), 1843–1851. Web of Science CrossRef
Lin, P. & Lucchese, R. R. (2000b). In preparation.
Lucchese, R. R. & McKoy, V. (1983). Phys. Rev. A, 28, 1382–1394. CrossRef CAS Web of Science
Lucchese, R. R., Raseev, G. & McKoy, V. (1982). Phys. Rev. A, 25(5), 2572–2587. CrossRef Web of Science
Lucchese, R. R., Takatsuka, K. & McKoy, V. (1986). Phys. Rep. 131(3), 147–221. CrossRef Web of Science
Lynch, D. L. & McKoy, V. (1984). Phys. Rev. A, 30(3), 1561–1564. CrossRef Web of Science
Natalense, A. P. P. & Lucchese, R. R. (1999). J. Chem. Phys. 111(12), 5344–5348. Web of Science CrossRef
Rescigno, T. N. & Langhoff, P. W. (1977). Chem. Phys. Lett. 51(1), 65–70. CrossRef Web of Science
Ruedel, A. & Hergenhahn, U. (2000). Private communication.
Siegbahn, K., Nordling, C., Johansson, G., Hedman, J., Hedin, P. F., Hamrin, K., Gelius, U., Bergmark, T., Werme, L. O., Manne, R. & Baer, Y. (1969). ESCA Applied to Free Molecules. Amsterdam: North-Holland.
Stratmann, R. E. & Lucchese, R. R. (1995). J. Chem. Phys. 102(21), 8493–8505. CrossRef Web of Science
Stratmann, R. E., Zurales, R. W. & Lucchese, R. R. (1996). J. Chem. Phys. 104(22), 8989–9000. CrossRef Web of Science
Thomas, T. D., Berrah, N., Bozek, J., Carroll, T. X., Hahne, J., Karlsen, T., Kukk, E. & Saethre, L. J. (1999). Phys. Rev. Lett. 82(6), 1120–1123. Web of Science CrossRef
Wallace, S., Dill, D. & Dehmer, J. L. (1982). J. Chem. Phys. 76(3), 1217–1222. CrossRef Web of Science
Wells, M. & Lucchese, R. R. (1999). J. Chem. Phys. 110(13), 6365–6380. Web of Science CrossRef
Werner, H.-J. & Knowles, P. J. (1996). MOLPRO. University of Birmingham, UK.
© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.