A methodology is developed to compute photoionization cross sections beyond the electric dipole approximation from response theory, using Gaussian type orbitals and plane waves for the initial and final states, respectively. The methodology is applied to compute photoionization cross sections of atoms and ions from the first four rows of the periodic table. Analyzing the error due to the plane wave description of the photoelectron, we find kinetic energy and concomitant photon energy thresholds above which the plane wave approximation becomes applicable. The correction introduced by going beyond the electric dipole approximation increases with photon energy and depends on the spatial extension of the initial state. In general, the corrections are below 10% for most elements, at a photon energy reaching up to 12 keV.

Photoionization is the physical process in which atoms, molecules, or solids emit electrons upon irradiation. It constitutes the basis for photoelectron spectroscopy (PES), also known as photoemission, an experimental technique routinely used for the characterization of materials.1,2 By recording the kinetic energy of emitted electrons, PES provides a description of the occupied electronic structure of the sample. Besides information on the binding energies of the electronic levels, PE spectra measured using different photon energies can be used to distinguish between states with different atomic compositions.3,4

Commonly, the comparison between experimental PE spectra and calculated electronic structures or total densities of states (DOS) is used to interpret the quality of a particular electronic structure method.5–13 A step beyond such a simple comparison is to include in the calculation the probability of ionization of a particular state. This can be performed by multiplying the projected density of states with the corresponding photoionization cross section calculated for atomic orbitals within the dipole approximation, as performed, for example, in Refs. 14–18. Another possibility is to compute the dipole transition matrix elements using the Fermi golden rule, as in Refs. 19–23.

Photoionization cross sections have been calculated within the dipole approximation for different ranges of photon energies and using different mathematical descriptions for the initial and final states. Most notably, Yeh and Lindau computed the total and partial ionization cross sections for all atoms in the periodic table for photon energies between 10 and 1500 eV.24,25 They used an algorithm introduced by Cooper and Manson,26,27 where both initial and final states were described on a radial grid and the final state was relaxed in the presence of the positive ion.28 Other attempts have also been made using plane waves,21,22,29–31 orthogonalized plane waves,22,29–31 Coulomb waves,21,32 or B-splines33,34 to describe the final states. In the solid state community, time-reversed (TR) low energy electron diffraction (LEED) states20,35,36 for the final state are used in a Green’s function method approach to obtain angle-resolved photoemission spectra (ARPES) of 3d- and 4f-materials.20,37–39

The electric dipole (ED) approximation is based on the assumption that the electromagnetic field can be considered uniform across the spatial extent of the initial state wave function. This allows for the truncation of the series expansion of the spatial part of the field at zeroth order, i.e., eik·r ≈ 1. Here, k represents the wavevector of the electromagnetic radiation with wavelength λ, and we have k = 2π/λ. Thus, the dipole approximation is not valid when the wavelength λ is comparable to or smaller than the spatial extent of the initial-state wave function (or orbital), i.e., when the term k · r becomes comparable to or larger than 1. Going beyond, the electric dipole approximation becomes especially important when the ionizing photon energy is large and the ionized electronic states are highly delocalized, such as the valence orbitals of conjugated organic molecules or the valence bands of solids.

Non-dipole (ND) effects have been theoretically described and experimentally observed especially in terms of the angular distributions of the emitted photoelectrons.23,40–47 In general, non-dipole photoionization cross sections and photoelectron angular distributions are calculated by including higher order terms, usually up to the first or second order, in the series expansion of eik·r as performed, for example, in Refs. 45 and 47–53. For a review on the topic, we refer the reader to Ref. 54. Not many studies include the full field operator in the photoionization cross section, but we note the work of Demekhin41 for atoms and of Seabra and co-workers23 for molecules. Demekhin included the full field operator to calculate non-dipole cross sections for H, He, Be, Ne, Ar, and Kr using atomic orbitals relaxed in a scalar relativistic framework, while Seabra et al.23 used Dyson orbitals and a plane wave (PW) final state to compute the error introduced by the dipole approximation for HF, H2O, NH3, and CH4.

In the present study, we analyze the correction introduced by using the full field operator in the expression for the atomic photoionization cross section, following the theoretical derivation based on response theory from Ref. 55. By using Gaussian type orbitals (GTOs) and plane waves (PWs) for the initial and final states, respectively, we compute ND cross sections (denoted σND throughout this article) for the atoms in the first four rows of the periodic table, alongside their most common ions, and estimate the magnitude of the ND correction as a function of photon energy for each system. Furthermore, we analyze the error introduced by the PW description for the outgoing photoelectron and find a kinetic energy threshold above which this approximation becomes applicable. It is often assumed that a PW description for the final state is correct at large photoelectron kinetic energies, but, to our knowledge, there is no systematic study of how the PW approximation performs for different kinetic energies, atoms, and atomic shells.

The article is organized as follows: the derivation of the expression for σND is presented in Sec. II, the computational details are described in Sec. III, while results for selected atoms are presented in Sec. IV. The complete set of data for all atoms and ions included in this study is provided at https://www.theochem.kth.se/BED_atomic_cross_sections.

In this section, we present the derivation of the equation used to calculate atomic and molecular photoionization cross sections beyond the electric dipole approximation. We start from the expression of the absorption cross section derived from complex response theory by List and co-workers55 and apply it to the photoionization process.

The photoabsorption cross section obtained as the linear response of the system of electrons to a time-dependent electromagnetic field is55 

(1)

where , e, ε0, c, and me are the reduced Planck constant, elementary charge, vacuum permittivity, speed of light in vacuum, and electron mass, respectively. ω is the angular frequency of the electromagnetic field, k is the associated wavevector, and rj is the position vector corresponding to electron j. |0 and |n represent the many-body ground state and excited state, respectively. The summation is performed over all excited states and ωn0 designates the angular frequency associated with the energy difference between the excited state and the ground state.

Equation (1) has been derived from response theory by List and co-workers,55 but for the sake of completeness, we include the intermediate steps required to arrive at this expression in the supplementary material.

Equation (1) is a general expression applicable to single photon absorption where the initial and final states are described by the many-body wave functions |0⟩ and |n⟩, respectively. Considering that the electronic structure of the atoms we are interested in may be approximated by a single Slater determinant, the matrix element of the one-electron operator eikrjϵj can be re-written in terms of spin-orbitals,56 

(2)

where χi is an occupied atomic orbital and χf is the wave function corresponding to the photoelectron. We note that the same expression is also valid in the case of molecular orbitals.

Considering that ND corrections become more important at high photon energies, where the wavelength of the ionizing electromagnetic radiation becomes comparable to the spatial extent of (valence) atomic orbitals, it is reasonable to assume that the photoelectron, with a high kinetic energy, can be described as a periodic plane wave normalized over an arbitrary volume Ω,

(3)

with density of states

(4)

Using PWs and the corresponding free electron density of states, in the limit of a large volume Ω, the sum over final states in Eq. (2) becomes an integral over a continuum,

(5)

where the kinetic energy of the photoelectron (Ef = ℏωf) satisfies energy conservation ωf = ωi + ω and the photoelectron wavevector ke is correspondingly chosen. The advantage of using a PW for the final state is the fact that PWs are eigenfunctions of the momentum operator, and the expression for the photoionization cross section may be, therefore, further simplified,

(6)

Finally, by expressing the initial state wave function χi as a linear combination of Gaussian type orbitals (GTOs), the photoionization cross section becomes

(7)

where the summation is over all primitives used to describe the initial state and gl is a GTO with l = lx + ly + lz,

(8)

with Rx, Ry, and Rz the coordinates of the atom where the GTO is centred.

Equation (7) can be straightforwardly generalized to molecular orbitals, by replacing χi with a linear combination of atomic orbitals.

Finally, the matrix element to be calculated is

(9)

where r=xe^x+ye^y+ze^z and K=Kxe^x+Kye^y+Kze^z determined from the input parameters ke and k, K = kek.

The advantage of working with PWs and GTOs is the fact that the integrals in Eq. (9) can be computed analytically.57 We provide analytic formulas up to l = 5 in the supplementary material. Furthermore, the dipole approximation is obtained simply by setting the photon wavevector k = 0, and hence, the dipole and ND photoionization cross sections are treated with the same formalism, and no further variables are introduced in the comparison.

The electronic structure of the atoms from the first four rows of the periodic table (H–Kr) was calculated with the Gaussian 1658 quantum chemistry software, at the B3LYP level of theory59 and using the universal Gaussian basis set (UGBS) by de Castro and co-workers.60–68 The electronic structure calculations were performed using Cartesian Gaussian functions because the expression for the matrix elements involved in the photoionization cross section, Eq. (7), has been derived for this GTO variant. For the N atom, we have additionally computed cross sections using five other basis sets, namely, STO-3G,69,70 4-31G,71–74 6-31G,71–80 cc-pVTZ,81 and cc-pVQZ.82 The basis set results are shown in Fig. S1 of the supplementary material.

In addition to the neutral atoms, we have also computed the most common ions corresponding to each atom. The atoms and ions included are listed in Table I, alongside their electronic configurations. The B3LYP-optimized electronic structures were used to compute the dipole and beyond dipole photoionization cross sections for a range of photon energies between 20 eV and 12 keV, with an energy step of 10 eV. A photoionization cross section was calculated for each atomic orbital in the relaxed electronic configuration using Eq. (7) and the computational setup shown in Fig. 1. The wavevector of the photon, k, was set in the y-direction, while the polarization vector was chosen in the z-direction. The norm of k is calculated from the photon energy ℏω. The norm of the photoelectron wavevector ke is calculated from the kinetic energy KE, in turn, computed from energy conservation,

(10)

where εAO represents the computed negative eigenvalue of a particular atomic orbital.

TABLE I.

Electronic configuration of the atoms and ions included in this study.

ZElementConfigurationElementConfiguration
1s1   
He 1s2   
Li 1s22s1 Li1+ 1s2 
Be 1s22s2 Be2+ 1s2 
1s22s22p1 B3+ 1s2 
1s22s22p2   
1s22s22p3   
1s22s22p4 O2− 1s22s22p6 
1s22s22p5 F1− 1s22s22p6 
10 Ne 1s22s22p6   
11 Na [Ne] 3s1 Na1+ [Ne] 
12 Mg [Ne] 3s2 Mg2+ [Ne] 
13 Al [Ne] 3s23p1 Al3+ [Ne] 
14 Si [Ne] 3s23p2   
15 [Ne] 3s23p3   
16 [Ne] 3s23p4 S2− [Ne] 3s23p6 
17 Cl [Ne] 3s23p5 Cl1− [Ne] 3s23p6 
18 Ar [Ne] 3s23p6   
19 [Ar] 4s1 K1+ [Ar] 
20 Ca [Ar] 4s2 Ca2+ [Ar] 
21 Sc [Ar] 4s23d1 Sc3+ [Ar] 
22 Ti [Ar] 4s23d2 Ti4+ [Ar] 
23 [Ar] 4s23d3 V5+ [Ar] 
24 Cr [Ar] 4s13d5 Cr6+ [Ar] 
25 Mn [Ar] 4s23d5 Mn2+ [Ar] 3d5 
26 Fe [Ar] 4s23d6 Fe2+ [Ar] 3d6 
27 Co [Ar] 4s23d7 Co2+ [Ar] 3d7 
28 Ni [Ar] 4s23d8 Ni2+ [Ar] 3d8 
29 Cu [Ar] 4s13d10 Cu2+ [Ar] 3d9 
30 Zn [Ar] 4s23d10 Zn2+ [Ar] 3d10 
31 Ga [Ar] 4s23d104p1 Ga3+ [Ar] 3d10 
32 Ge [Ar] 4s23d104p2   
33 As [Ar] 4s23d104p3   
34 Se [Ar] 4s23d104p4 Se2− [Ar] 4s23d104p6 
35 Br [Ar] 4s23d104p5 Br1− [Ar] 4s23d104p6 
36 Kr [Ar] 4s23d104p6   
ZElementConfigurationElementConfiguration
1s1   
He 1s2   
Li 1s22s1 Li1+ 1s2 
Be 1s22s2 Be2+ 1s2 
1s22s22p1 B3+ 1s2 
1s22s22p2   
1s22s22p3   
1s22s22p4 O2− 1s22s22p6 
1s22s22p5 F1− 1s22s22p6 
10 Ne 1s22s22p6   
11 Na [Ne] 3s1 Na1+ [Ne] 
12 Mg [Ne] 3s2 Mg2+ [Ne] 
13 Al [Ne] 3s23p1 Al3+ [Ne] 
14 Si [Ne] 3s23p2   
15 [Ne] 3s23p3   
16 [Ne] 3s23p4 S2− [Ne] 3s23p6 
17 Cl [Ne] 3s23p5 Cl1− [Ne] 3s23p6 
18 Ar [Ne] 3s23p6   
19 [Ar] 4s1 K1+ [Ar] 
20 Ca [Ar] 4s2 Ca2+ [Ar] 
21 Sc [Ar] 4s23d1 Sc3+ [Ar] 
22 Ti [Ar] 4s23d2 Ti4+ [Ar] 
23 [Ar] 4s23d3 V5+ [Ar] 
24 Cr [Ar] 4s13d5 Cr6+ [Ar] 
25 Mn [Ar] 4s23d5 Mn2+ [Ar] 3d5 
26 Fe [Ar] 4s23d6 Fe2+ [Ar] 3d6 
27 Co [Ar] 4s23d7 Co2+ [Ar] 3d7 
28 Ni [Ar] 4s23d8 Ni2+ [Ar] 3d8 
29 Cu [Ar] 4s13d10 Cu2+ [Ar] 3d9 
30 Zn [Ar] 4s23d10 Zn2+ [Ar] 3d10 
31 Ga [Ar] 4s23d104p1 Ga3+ [Ar] 3d10 
32 Ge [Ar] 4s23d104p2   
33 As [Ar] 4s23d104p3   
34 Se [Ar] 4s23d104p4 Se2− [Ar] 4s23d104p6 
35 Br [Ar] 4s23d104p5 Br1− [Ar] 4s23d104p6 
36 Kr [Ar] 4s23d104p6   
FIG. 1.

Schematic representation of the computational setup. The atom is placed at the origin of the coordinate system, the photon wavevector is directed along the y-axis, the polarization direction is along the z-axis, and the direction of the photoelectron wavevector is generated using the Lebedev quadrature. The points generated by the quadrature at the 6th order are marked with yellow dots.

FIG. 1.

Schematic representation of the computational setup. The atom is placed at the origin of the coordinate system, the photon wavevector is directed along the y-axis, the polarization direction is along the z-axis, and the direction of the photoelectron wavevector is generated using the Lebedev quadrature. The points generated by the quadrature at the 6th order are marked with yellow dots.

Close modal

The ND photoionization cross section is finally computed by averaging over all possible photoelectron directions generated on a sphere using the Lebedev quadrature.83 We tested the Lebedev quadrature at different orders. Due to the spherical symmetry of the atomic system, the lowest order of the quadrature (i.e., 6) sufficed in the case of the dipole photoionization cross section. For ND, the order 50 was instead required to fully converge the cross sections. We have therefore computed both ND and dipole cross sections using the order 50 for the Lebedev quadrature.

In the case of open-shell atoms, we additionally averaged over the possible initial state configurations of the partially filled shell. For example, in the case of the C atom, the 2p cross section is computed as follows:

(11)

where the index pxα means that the spin-up 2px orbital is occupied.

The total atomic photoionization cross section was calculated by summing over the cross sections of the occupied atomic shells. Dipole photoionization cross sections were calculated by the same algorithm with the only difference that the wavevector k was set to zero.

In the following (Sec. IV A), we analyze the performance of the PW approximation for the final state by comparing our calculated photoionization cross sections to calculated results from the literature and experimentally measured data (Sec. IV B). We then compare the dipole and ND cross sections to each other and calculate the magnitude of the relative correction introduced by going beyond the electric dipole approximation as a function of photon energy.

Figure 2 shows the total dipole photoionization cross section calculated using a PW final state (σPW) in comparison with the cross section calculated by Yeh and Lindau24,25,84 using a final state relaxed in the presence of the positive ion (σYL). The comparison is shown as a function of photon energy, from 20 eV up to 1500 eV, for Li [Fig. 2(a)], C [Fig. 2(b)], Mn [Fig. 2(c)], and O [Fig. 2(d)]. The figure also depicts the ratio σPW/σYL as a function of photon energy.

FIG. 2.

Comparison between total dipole photoionization cross sections calculated using a PW final state (σPW, in dark green) and total dipole photoionization cross sections calculated by Yeh and Lindau24,25,84 (σYL, in red) for (a) lithium, (b) carbon, (c) manganese, and (d) oxygen. The ratio between σPW and σYL is also depicted (dotted line). The labels 1s, 2s, and 2p mark the ionization thresholds for these orbitals. The PW cross sections have been multiplied by a factor of 1/3 to average over the directions of the incoming photon, as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

FIG. 2.

Comparison between total dipole photoionization cross sections calculated using a PW final state (σPW, in dark green) and total dipole photoionization cross sections calculated by Yeh and Lindau24,25,84 (σYL, in red) for (a) lithium, (b) carbon, (c) manganese, and (d) oxygen. The ratio between σPW and σYL is also depicted (dotted line). The labels 1s, 2s, and 2p mark the ionization thresholds for these orbitals. The PW cross sections have been multiplied by a factor of 1/3 to average over the directions of the incoming photon, as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

Close modal

What becomes clear from the figure is that σPW converges toward σYL as the photon energy is increased, i.e., as the kinetic energy of the photoelectron increases. In the vicinity of the ionization threshold, the PW approximation performs poorly, but it improves as the photon energy is increased. The photon energy above which the PW approximation becomes applicable depends on the particular element, becoming larger as the atomic number increases. For example, going approximately 100 eV above the 1s threshold is required in the case of the Li atom, but more than 700 eV above the 2s threshold is instead required for the Mn atom. The poor performance in the vicinity of the ionization threshold is expected since the outgoing electron, in this case, has low kinetic energy and is not free electron like.

In addition, we have also analyzed selected partial dipole photoionization cross sections and how they compare to the results obtained by Yeh and Lindau.24,25,84 These comparisons are shown in Fig. 3 for C [Fig. 3(a)] and O [Fig. 3(b)] in the photon energy window between 20 eV and 1500 eV. As it can be seen from Fig. 3, the PW approximation becomes better further away from the ionization threshold, but different atomic orbitals behave differently. The agreement to the Yeh and Lindau results24,25,84 is better in the case of s orbitals than in the case of p orbitals.

FIG. 3.

Comparison between partial dipole photoionization cross sections calculated using a PW final state (blue solid lines) and partial dipole photoionization cross sections calculated by Yeh and Lindau24,25,84 (red symbols) for (a) carbon and (b) oxygen. The ratio between σPW and σYL is also depicted (blue dotted lines). The PW cross sections have been multiplied by a factor of 1/3 to average over the directions of the incoming photon as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

FIG. 3.

Comparison between partial dipole photoionization cross sections calculated using a PW final state (blue solid lines) and partial dipole photoionization cross sections calculated by Yeh and Lindau24,25,84 (red symbols) for (a) carbon and (b) oxygen. The ratio between σPW and σYL is also depicted (blue dotted lines). The PW cross sections have been multiplied by a factor of 1/3 to average over the directions of the incoming photon as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

Close modal

It is interesting to note that the cross sections calculated for the 2s orbitals of both the C and O atoms present a pronounced dip (Cooper minimum), at ℏω ≈ 90 eV and ℏω ≈ 160 eV, respectively (marked by arrows in Fig. 3). The reason for this minimum is that the transition matrix element is zero at that particular photon energy. Because the 2s atomic orbital has one node, it is expressed as a linear combination of GTOs with both positive and negative coefficients. Since the integrals involving s-type GTOs and a PW are always real and positive (see the analytical formulas for the integrals in the supplementary material), there are both negative and positive contributions to the transition matrix element. For a particular value of the photoelectron wavevector, the positive and negative contributions cancel out, giving rise to the minimum in the photoionization cross section. All other orbitals which have nodes present similar minima. General equations for different types of atomic orbitals described as linear combinations of GTOs may be worked out. We include one example for a very simple 2s orbital in the supplementary material.

Minima arising from the cancellation of the positive and negative contributions to the integrand have been first described by Cooper,26 using final states relaxed in the presence of the positive ion. Subsequently, Cooper minima have been observed experimentally, most notably in the case of Ar, due to the ionization of the 3s subshell.85 

We should note that a perfect cancellation between the positive and negative parts of the transition matrix element does not arise for Slater type orbitals (STOs) in combination with a PW final state, as exemplified in the supplementary material for a 2s STO.

To further elucidate if the PW approximation for the final state is applicable, we have compared the non-dipole cross section (σND) to experimental data and calculations performed over a much larger photon energy range. Figure 4 shows this comparison for photon energies in the range of 20 eV to 12 keV for Li [Fig. 4(a)], N [Fig. 4(b)], Ar [Fig. 4(c)], and O [Fig. 4(d)]. The experimental data were compiled for nine atoms from selected reference sources by Berkowitz.86 The measurements are typically performed in gas phase and at room temperature, as described in Refs. 89–91. The calculated data were obtained within the dipole approximation using the Hartree-Fock-Slater method (the same algorithm used by Yeh and Lindau) and are stored on the National Institute of Standards and Technologies (NIST) database.88 

FIG. 4.

Comparison of the total ND photoionization cross sections calculated using a PW final state (σND, black line) with measured photoionization cross sections compiled by Berkowitz86 (stored on the NIFS database,87σexp, blue squares) and with total dipole cross sections calculated using the Hartree-Fock-Slater method (σNIST, red circles) from the NIST database.88 The comparison is shown for (a) Li, (b) N, (c) Ar, and (d) O. σND has been multiplied by a factor of 1/3 to average over the directions of the incoming photon, as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

FIG. 4.

Comparison of the total ND photoionization cross sections calculated using a PW final state (σND, black line) with measured photoionization cross sections compiled by Berkowitz86 (stored on the NIFS database,87σexp, blue squares) and with total dipole cross sections calculated using the Hartree-Fock-Slater method (σNIST, red circles) from the NIST database.88 The comparison is shown for (a) Li, (b) N, (c) Ar, and (d) O. σND has been multiplied by a factor of 1/3 to average over the directions of the incoming photon, as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

Close modal

The figure shows that far from the ionization threshold the cross sections calculated using a PW final state (beyond the electric dipole approximation) compare very well with the experimental data. Additionally, the values of σND are very close to the calculated dipole cross sections σNIST. As it will be discussed in Sec. IV B, the reason for this similarity is that the relative correction introduced by ND to the total photoionization cross section does not amount to more than ∼5%.

Figure 5 shows the comparison between partial cross sections calculated beyond the electric dipole using a PW final state and the first order corrected photoionization cross sections computed by Cooper.47 The two results are in good agreement in the case of the s orbitals, once far enough from the ionization threshold. In the case of the 4s orbital of Kr, the comparison is hindered by the presence of one of the Cooper minima in σND at approximately 5 keV photon energy. The cross sections calculated using a PW final state and the full field operator are larger than the first order corrected cross sections for the p and d orbitals. This difference may partly be assigned to the larger error introduced by the PW approximation for these types of orbitals (see Fig. 3) and partly to the use of the full field operator in σNDvs. the first order correction in Ref. 47.53,54

FIG. 5.

Comparison of partial ND photoionization cross sections calculated using a PW final state (solid lines) with first order corrected non-dipole photoionization cross sections computed by Cooper47 (circles). The comparison is shown for (a) He, (b) Ne, (c) Ar, and (d) Kr. σND has been multiplied by a factor of 1/3 to average over the directions of the incoming photon, as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

FIG. 5.

Comparison of partial ND photoionization cross sections calculated using a PW final state (solid lines) with first order corrected non-dipole photoionization cross sections computed by Cooper47 (circles). The comparison is shown for (a) He, (b) Ne, (c) Ar, and (d) Kr. σND has been multiplied by a factor of 1/3 to average over the directions of the incoming photon, as performed in Refs. 24 and 25. Note that the vertical axis is logarithmic.

Close modal

We have also computed the photoelectron angular distribution at a photon energy ℏω = 5206 eV, for the Ar 1s orbital, shown in Fig. 6. The angular distribution of photoelectrons is typically measured using linearly polarized light by placing an electron analyser at a fixed polar angle θ, as depicted in Fig. 6(b). The number of electrons emitted at different azimuthal angles is then recorded by rotating the electron analyser around the z-axis and performing measurements at different values of φ. We note that the polar angle θ is defined as the angle between the polarization direction of the photon (ϵ) and the photoelectron direction (ke), while the azimuthal angle φ [shaded in gray in Fig. 6(b)] is the angle between the projection of ke in the xy-plane and the photon direction (k).47,92,93 If the photoionization process would take place according to the dipole approximation, the number of emitted photoelectrons would not depend on φ and the graph in Fig. 6(a) would be a straight horizontal line. The deviation of the curve from a straight line is therefore a measure of nondipolar asymmetry. The differential cross section for photoionization (dσ/dΩ) using linearly polarized light, taking into account quadrupole terms, may be expressed as47,92

(12)

where Ω is the solid angle, β is the anisotropy parameter, δ and γ are the asymmetry parameters which quantify the nondipolar effects, and P2(cos θ) is the second Legendre polynomial.

FIG. 6.

(a) Comparison between the angular distributions of photoelectrons calculated using the ND cross sections (black solid line) and the experimental photoelectron yield (red circles) as a function of the azimuthal angle φ. The angular distribution is determined for the 1s orbital of Ar and a photon energy ℏω = 5206 eV. The experimental data are from Ref. 92. The values on the y-axis have been normalized to the cross section calculated at φ = 0°. (b) depicts the experimental setup, where the analyser is placed at a fixed polar angle θ = 54.7° and rotated on the yellow circle to vary φ.

FIG. 6.

(a) Comparison between the angular distributions of photoelectrons calculated using the ND cross sections (black solid line) and the experimental photoelectron yield (red circles) as a function of the azimuthal angle φ. The angular distribution is determined for the 1s orbital of Ar and a photon energy ℏω = 5206 eV. The experimental data are from Ref. 92. The values on the y-axis have been normalized to the cross section calculated at φ = 0°. (b) depicts the experimental setup, where the analyser is placed at a fixed polar angle θ = 54.7° and rotated on the yellow circle to vary φ.

Close modal

At the magic angles θ = 54.7° and 180° − θ, P2(cos θ) is zero and the differential cross section then depends only on the azimuthal angle φ,93 

(13)

We have calculated σND for the 1s atomic orbital of Ar at different φ angles by rotating the ke vector around the z-axis, as depicted in Fig. 6(b). The photon polarization was kept in the z-direction and the wavevector of the photon was kept in the y-direction. We normalized the value of the calculated cross section to the value calculated for φ = 0°, as performed for the experimental data.92 The experimental and calculated curves match rather well, and our calculated angle-dependent σND cross section is able to reproduce the magnitude of nondipolar effects observed experimentally.

To summarize this section, by comparing our atomic photoionization cross sections calculated using a PW final state to previously calculated results from the literature and experimental data, we can conclude that the PW approximation for the final state is applicable at high kinetic energies, i.e., far from the ionization threshold. However, the kinetic energy above which the approximation is applicable depends on the atom and type of orbitals since the comparison becomes worse for larger atoms and for orbitals with larger l quantum numbers. Using the PW final state, the cross sections computed for orbitals with one or more nodes present minima at photon energies (generally close to the ionization threshold) where the positive and negative contributions to the transition matrix element cancel out. Finally, our calculated σND is able to reproduce the nondipolar angular distribution of photoelectrons at high photon energies.

Having analyzed the advantages and limitations of a PW final state, we now move to examine the magnitude of the correction introduced by going beyond the electric dipole approximation. Figure 7 shows the correction obtained for total and partial photoionization cross sections of selected atoms as a function of photon energy, between 20 eV and 12 keV. The correction is represented as

(14)

The partial cross sections are grouped together based on the orbital type (s, p, d) because the values computed from different shells (1s, 2s, etc.), but the same l quantum number, are very similar to each other.

FIG. 7.

The magnitude of the correction introduced by going beyond the dipole approximation as a function of photon energy for (a) Li, (b) N, (c) Mn, and (d) O. The corrections to partial cross sections are represented in pink (s orbitals), blue (p orbitals), and brown (d orbitals).

FIG. 7.

The magnitude of the correction introduced by going beyond the dipole approximation as a function of photon energy for (a) Li, (b) N, (c) Mn, and (d) O. The corrections to partial cross sections are represented in pink (s orbitals), blue (p orbitals), and brown (d orbitals).

Close modal

As expected, the ND correction increases as the photon energy is increased. In the case of the s partial cross sections, the values reach only up to 5% for a photon energy ℏω = 12 keV. At this photon energy the wavelength of the photon (1.95 a.u.) begins to become comparable to the size of the orbitals (∼2 a.u., for 2s orbitals). The total cross sections (not shown) follow closely the behavior of the s-type orbitals. This is expected since the last shell to ionize will be the one which will dominate, as illustrated in Fig. 3, where the 1s cross section is 2-3 orders of magnitude larger than the cross sections of the 2s and 2p orbitals. For all the elements of the first four rows of the periodic table, at 12 keV, the last shell to ionize is either a 1s or a 2s shell, resulting in the high similarity between the correction for s partial photoionization cross section and the total one.

In the case of p and d atomic orbitals, the correction introduced by ND is slightly larger, between 5% and 10% for the three atoms shown in Fig. 7. One last aspect to note is the appearance of sharp peaks, generally at low photon energies, notably visible in the case of the Mn atom [marked with arrows in Fig. 7(c)]. These are related to the Cooper minima previously discussed, resulting from the cancellation between the positive and negative contributions to the transition matrix element which involves an orbital with at least one node.

We further compare our calculated ND corrections to the results of Demekhin.41Figure 8 shows the comparison between the ratio σED/σND calculated using a PW final state and the same ratio obtained in Ref. 41, using atomic orbitals within a scalar relativistic framework. We obtain the same ratio as Demekhin41 for the s orbitals of H, He, and Be, and for the 2p orbitals of Ne. Our calculated ratio is larger than the one in Ref. 41 in the case of Kr 2p, Ar 1s > Ar 2s > Ar 2p, and Ne 1s > N 2s (where the atoms and subshells are ordered from the largest difference to the lowest). Considering this ordering, the difference between our results and those of Demekhin41 may be attributed partly to the fact that our electronic structure calculation does not include relativistic effects.

FIG. 8.

Comparison of the ratio σED/σND calculated using a PW final state (solid lines) to the ratio calculated by Demekhin41 (symbols) as a function of photon energy. The comparison is shown for (a) the 1s shells of H, He, Be (circles), Ne (diamonds), and Ar (triangles) and (b) the 2s and 2p shells of Ne (diamonds), Ar (triangles), and Kr (squares, 2p only).

FIG. 8.

Comparison of the ratio σED/σND calculated using a PW final state (solid lines) to the ratio calculated by Demekhin41 (symbols) as a function of photon energy. The comparison is shown for (a) the 1s shells of H, He, Be (circles), Ne (diamonds), and Ar (triangles) and (b) the 2s and 2p shells of Ne (diamonds), Ar (triangles), and Kr (squares, 2p only).

Close modal

Finally, we collect the corrections introduced by going beyond the electric dipole approximation for all atoms at the largest photon energy included in the study, i.e., ℏω = 12 keV. The results are represented as a function of atomic number in Fig. 9. The corrections to the total cross section are depicted in gray [Fig. 9(b)], while the l-shells with the largest correction to the cross sections are depicted in color [Fig. 9(a)]. The atomic number at which a new shell starts to fill is marked by a dotted line. Note that the maximum correction may be obtained for a different orbital than the last to be filled. For example, the largest correction obtained for the Na atom corresponds to a photoionization of the 2p orbitals, rather than the 3s which is the last to be filled.

FIG. 9.

The magnitude of the correction introduced by going beyond the dipole approximation as a function of the atomic number, at a photon energy ℏω = 12keV, for (a) the maximum calculated corrections to partial cross sections and (b) the corrections to total cross sections. For the points marked with dotted lines, the data shown were computed for the ions instead of the neutral atoms. (For the Ni and Co atoms, the B3LYP electronic configuration relaxed to a situation where the 4s orbital remained unoccupied irrespective of the initial guess. We have, therefore, used the calculated data for the ions instead.) The corrections to partial cross sections in (a) are coloured according to the orbital shell they were calculated for (red for s orbitals, blue for p orbitals, and brown for d orbitals). (c) The computational setup used to calculate the partial cross sections and (d) the radial part of the 2pz orbital as calculated using the B3LYP functional and the UGBS basis set for atoms B to Ne.

FIG. 9.

The magnitude of the correction introduced by going beyond the dipole approximation as a function of the atomic number, at a photon energy ℏω = 12keV, for (a) the maximum calculated corrections to partial cross sections and (b) the corrections to total cross sections. For the points marked with dotted lines, the data shown were computed for the ions instead of the neutral atoms. (For the Ni and Co atoms, the B3LYP electronic configuration relaxed to a situation where the 4s orbital remained unoccupied irrespective of the initial guess. We have, therefore, used the calculated data for the ions instead.) The corrections to partial cross sections in (a) are coloured according to the orbital shell they were calculated for (red for s orbitals, blue for p orbitals, and brown for d orbitals). (c) The computational setup used to calculate the partial cross sections and (d) the radial part of the 2pz orbital as calculated using the B3LYP functional and the UGBS basis set for atoms B to Ne.

Close modal

The first feature to notice in Fig. 9(a) is that there are three “edges”: the first at Z = 1, the second at Z = 5, and the third at Z = 21. Interestingly, these edges arise when making a transition from l to l + 1, i.e., going from s to p, and, respectively, from p to d, where the corrections to the p and d cross sections are larger.

The larger sensitivity of the p and d orbitals to the ND correction may be understood by visualizing the effect of the ND correction in the computational setup, as represented in Fig. 9(c). The dipole approximation is obtained by setting the photon wavevector to zero, meaning the transition matrix elements that have to be calculated involve the initial state (s, p, or d orbital) and a PW of wavevector ke directed along the z-axis. Going beyond the dipole approximation (k ≠ 0) translates into having to calculate transition matrix elements between the initial state and a new PW now directed along the z′-axis, with a new K = kek. This new PW has both a different magnitude and a different direction, compared to ke. In the case of the s orbital, which is spherically symmetric, only the difference in magnitude affects the transition matrix element, while in the case of the p and d orbitals, both magnitude and direction affect the transition matrix element, making this type of orbitals more sensitive to the dipole approximation.

Another feature of the graph in Fig. 9(a) is that each edge is followed by a linear decrease in the magnitude of the ND correction. Considering that all these are valence orbitals of relatively similar binding energies (between 5 and 20 eV for the majority of elements) and that relativistic effects may be neglected, this decrease may be explained by the decrease in the spatial extent of the particular type of orbitals when going through the series. This is illustrated for the 2p elements in Fig. 9(d).

In summary, we have derived an equation to compute photoionization cross sections beyond the electric dipole approximation from response theory, using Gaussian type orbitals for the initial state and a plane wave with appropriate kinetic energy for the final state. We compared our calculated cross sections to the experimental data and to the calculated dipole photoionization cross sections from the literature. This comparison showed that the PW approximation for the final state is applicable far from the ionization threshold. However, how far is sufficient depends on the atom and the type of atomic orbital, the comparison becoming worse for the larger atoms and for orbitals with larger l quantum numbers.

Given that the dipole and beyond dipole photoionization cross sections can be computed analytically and on the same footing, we could estimate the relative correction introduced by going beyond the electric dipole approximation. As expected, this correction increases with photon energy, but it is, in general, rather small, ∼5%–10% at the largest photon energy we included (12 keV). The fact that we obtain such a small correction is not surprising, considering the small spatial extent of the atomic orbitals. Even though the ND corrections for the atoms are quite small, we expect larger corrections for more delocalized states as, for example, molecular orbitals of conjugated molecules or valence states of solids.

Finally, the description for the atomic photoionization cross sections may be improved by replacing the PW final state with a linear combination of PWs relaxed in the presence of the positive ion. This is, however, rather complicated because the relaxation should be performed for a small kinetic energy window centred around the kinetic energy obtained from energy conservation. Considering all states up to this kinetic energy is simply not feasible.

Supplementary material contains (a) the complete derivation of the linear response function corresponding to a time-dependent perturbation V^ and a generic operator Ô; (b) derivation of the equation for the photoabsorption cross section from complex response theory; (c) analytical formulas for the integrals in Eq. (9); (d) photoionization cross sections calculated using different basis sets for N; (e) the analytical form of a transition matrix element between a simplified 2s orbital (described as a linear combination of two GTOs) and a plane wave.

Financial support from the Knut and Alice Wallenberg Foundation (Grant No. KAW-2013.0020) and the Swedish Research Council (Grant No. 621-2014-4646) is acknowledged. The computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at NSC and HPC2N. I.E.B. is very grateful to Faris Gel’mukhanov, Nanna Holmgaard List, Barbara Brena, and Igor Di Marco for all the insightful discussions. O.E. also acknowledges support from eSSENCE.

1.
H.
Stefan
,
Photoelectron Spectroscopy: Principles and Applications
(
Springer
,
2003
), Vol. 3.
2.
A.
Damascelli
,
Z.
Hussain
, and
Z.-X.
Shen
,
Rev. Mod. Phys.
75
,
473
(
2003
).
3.
J. C.
Green
,
Acc. Chem. Res.
27
,
131
(
1994
).
4.
J. C.
Green
and
P.
Decleva
,
Coord. Chem. Rev.
249
,
209
(
2005
).
5.
S.
Kümmel
and
L.
Kronik
,
Rev. Mod. Phys.
80
,
3
(
2008
).
6.
L.
Kronik
and
Y.
Morikawa
, “
Understanding the metal-molecule interface from first principles
,” in
The Molecule-Metal Interface
, edited by
N.
Koch
,
N.
Ueno
, and
A. T. S.
Wee
(
John Wiley & Sons, Ltd.
,
2013
), Chap. 3, pp.
51
89
.
7.
L.
Kronik
and
S.
Kümmel
, “
Gas-phase valence-electron photoemission spectroscopy using density functional theory
,” in
First Principles Approaches to Spectroscopic Properties of Complex Materials
, edited by
C.
De Valentin
,
S.
Botti
, and
M.
Cococcioni
(
Springer Berlin Heidelberg
,
2014
), pp.
137
191
.
8.
L.
Kronik
and
J. B.
Neaton
,
Annu. Rev. Phys. Chem.
67
,
587
(
2016
).
9.
D. A.
Egger
,
S.
Weissman
,
S.
Refaely-Abramson
,
S.
Sharifzadeh
,
M.
Dauth
,
R.
Baer
,
S.
Kümmel
,
J. B.
Neaton
,
E.
Zojer
, and
L.
Kronik
,
J. Chem. Theory Comput.
10
,
1934
(
2014
).
10.
N.
Marom
,
F.
Caruso
,
X.
Ren
,
O. T.
Hofmann
,
T.
Körzdörfer
,
J. R.
Chelikowsky
,
A.
Rubio
,
M.
Scheffler
, and
P.
Rinke
,
Phys. Rev. B
86
,
245127
(
2012
).
11.
G.
Kotliar
,
S.
Savrasov
,
K.
Haule
,
V.
Oudovenko
,
O.
Parcollet
, and
C.
Marianetti
,
Rev. Mod. Phys.
78
,
865
(
2006
).
13.
F.
Manghi
,
V.
Bellini
,
J.
Osterwalder
,
T.
Kreutz
,
P.
Aebi
, and
C.
Arcangeli
,
Phys. Rev. B
59
,
R10409
(
1999
).
14.
S. K.
Panda
,
B.
Pal
,
S.
Mandal
,
M.
Gorgoi
,
S.
Das
,
I.
Sarkar
,
W.
Drube
,
W.
Sun
,
I.
Di Marco
,
A.
Lindblad
,
P.
Thunström
,
A.
Delin
,
O.
Karis
,
Y. O.
Kvashnin
,
M.
van Schilfgaarde
,
O.
Eriksson
, and
D. D.
Sarma
,
Phys. Rev. B
93
,
235138
(
2016
).
15.
B.
Brena
,
C.
Puglia
,
M.
de Simone
,
M.
Coreno
,
K.
Tarafder
,
V.
Feyer
,
R.
Banerjee
,
E.
Gthelid
,
B.
Sanyal
,
P. M.
Oppeneer
, and
O.
Eriksson
,
J. Chem. Phys.
134
,
074312
(
2011
).
16.
I. E.
Brumboiu
,
R.
Totani
,
M.
de Simone
,
M.
Coreno
,
C.
Grazioli
,
L.
Lozzi
,
H. C.
Herper
,
B.
Sanyal
,
O.
Eriksson
,
C.
Puglia
, and
B.
Brena
,
J. Phys. Chem. A
118
,
927
(
2014
).
17.
I.
Bidermane
,
I.
Brumboiu
,
R.
Totani
,
C.
Grazioli
,
M.
Shariati-Nilsson
,
H.
Herper
,
O.
Eriksson
,
B.
Sanyal
,
B.
Ressel
,
M.
de Simone
,
L.
Lozzi
,
B.
Brena
, and
C.
Puglia
,
J. Electron Spectrosc.
205
,
92
(
2015
).
18.
T.
Zhang
,
I. E.
Brumboiu
,
V.
Lanzilotto
,
J.
Lüder
,
C.
Grazioli
,
E.
Giangrisostomi
,
R.
Ovsyannikov
,
Y.
Sassa
,
I.
Bidermane
,
M.
Stupar
,
M.
de Simone
,
M.
Coreno
,
B.
Ressel
,
M.
Pedio
,
P.
Rudolf
,
B.
Brena
, and
C.
Puglia
,
J. Phys. Chem. C
121
,
26372
(
2017
).
19.
N. L.
Nguyen
,
G.
Borghi
,
A.
Ferretti
,
I.
Dabo
, and
N.
Marzari
,
Phys. Rev. Lett.
114
,
166405
(
2015
).
20.
J.
Minár
, “
Theoretical description of arpes: The one-step model
,” in
DMFT at 25: Infinite Dimensions
, Volume 4 of Modeling and Simulation, edited by
E.
Pavarini
,
E.
Koch
,
D.
Vollhardt
, and
A.
Lichtenstein
(
Verlag des Forschungszentrum Jülich
,
2014
), Chap. 13.
21.
S.
Gozem
,
A. O.
Gunina
,
T.
Ichino
,
D. L.
Osborn
,
J. F.
Stanton
, and
A. I.
Krylov
,
J. Phys. Chem. Lett.
6
,
4532
(
2015
).
22.
G.
Seabra
,
I.
Kaplan
,
V.
Zakrzewski
, and
J.
Ortiz
,
J. Chem. Phys.
121
,
4143
(
2004
).
23.
G. M.
Seabra
,
I. G.
Kaplan
, and
J. V.
Ortiz
,
J. Chem. Phys.
123
,
114105
(
2005
).
24.
J.
Yeh
and
I.
Lindau
,
At. Data Nucl. Data Tables
32
,
1
(
1985
).
25.
J.
Yeh
,
Atomic Calculation of Photoionization Cross-Sections and Asymmetry Parameters
(
Gordon & Breach Science Publishers
,
1993
).
26.
J. W.
Cooper
,
Phys. Rev.
128
,
681
(
1962
).
27.
S. T.
Manson
and
J. W.
Cooper
,
Phys. Rev.
165
,
126
(
1968
).
28.
U.
Fano
and
J. W.
Cooper
,
Rev. Mod. Phys.
40
,
441
(
1968
).
29.
F. O.
Ellison
,
J. Chem. Phys.
61
,
507
(
1974
).
30.
F. O.
Ellison
,
J. Chem. Phys.
62
,
4587
(
1975
).
31.
M.
Deleuze
,
B. T.
Pickup
, and
J.
Delhalle
,
Mol. Phys.
83
,
655
(
1994
).
32.
B.
Ritchie
,
J. Chem. Phys.
60
,
898
(
1974
).
33.
H.
Bachau
,
E.
Cormier
,
P.
Decleva
,
J.
Hansen
, and
F.
Martín
,
Rep. Prog. Phys.
64
,
1815
(
2001
).
34.
M.
Ruberti
,
V.
Averbukh
, and
P.
Decleva
,
J. Chem. Phys.
141
,
164126
(
2014
).
36.
J.
Braun
,
J.
Minár
, and
H.
Ebert
, “
One step model description of HARPES: Inclusion of disorder and temperature effects
,” in
Hard X-Ray Photoelectron Spectroscopy (HAXPES)
, edited by
J.
Woicik
(
Springer International Publishing
,
2016
), pp.
159
174
.
37.
A.
Keqi
,
M.
Gehlmann
,
G.
Conti
,
S.
Nemšák
,
A.
Rattanachata
,
J.
Minár
,
L.
Plucinski
,
J.
Rault
,
J.
Rueff
,
M.
Scarpulla
 et al.,
Phys. Rev. B
97
,
155149
(
2018
).
38.
M.
Veis
,
J.
Minár
,
G.
Steciuk
,
L.
Palatinus
,
C.
Rinaldi
,
M.
Cantoni
,
D.
Kriegner
,
K.
Tikuišis
,
J.
Hamrle
,
M.
Zahradník
 et al.,
Phys. Rev. B
97
,
125109
(
2018
).
39.
A.
Gray
,
C.
Papp
,
S.
Ueda
,
B.
Balke
,
Y.
Yamashita
,
L.
Plucinski
,
J.
Minár
,
J.
Braun
,
E.
Ylvisaker
,
C.
Schneider
 et al.,
Nat. Mater.
10
,
759
(
2011
).
40.
R.
Guillemin
,
O.
Hemmers
,
D. W.
Lindle
, and
S. T.
Manson
,
Radiat. Phys. Chem.
75
,
2258
(
2006
).
41.
P. V.
Demekhin
,
J. Phys. B: At., Mol. Opt. Phys.
47
,
025602
(
2014
).
42.
D. W.
Lindle
and
O.
Hemmers
,
J. Electron Spectrosc.
100
,
297
(
1999
).
43.
P.-S.
Shaw
,
U.
Arp
, and
S.
Southworth
,
Phys. Rev. A
54
,
1463
(
1996
).
44.
A. N.
Grum-Grzhimailo
,
J. Phys. B: At., Mol. Opt. Phys.
36
,
2385
(
2003
).
45.
A.
Derevianko
,
W. R.
Johnson
, and
K. T.
Cheng
,
At. Data Nucl. Data Tables
73
,
153
(
1999
).
46.
O.
Hemmers
,
G.
Fisher
,
P.
Glans
,
D.
Hansen
,
H.
Wang
,
S.
Whitfield
,
R.
Wehlitz
,
J.
Levin
,
I.
Sellin
,
R. C.
Perera
 et al.,
J. Phys. B: At., Mol. Opt. Phys.
30
,
L727
(
1997
).
48.
A.
Bechler
and
R. H.
Pratt
,
Phys. Rev. A
39
,
1774
(
1989
).
49.
A.
Bechler
and
R. H.
Pratt
,
Phys. Rev. A
42
,
6400
(
1990
).
50.
V. I.
Nefedov
,
V. G.
Yarzhemsky
,
I. S.
Nefedova
,
M. B.
Trzhaskovskaya
, and
I. M.
Band
,
J. Electron Spectrosc.
107
,
123
(
2000
).
51.
M. Y.
Amusia
,
A. S.
Baltenkov
,
L. V.
Chernysheva
,
Z.
Felfli
, and
A. Z.
Msezane
,
Phys. Rev. A
63
,
052506
(
2001
).
52.
V. K.
Dolmatov
and
S. T.
Manson
,
Phys. Rev. A
63
,
022704
(
2001
).
53.
A.
Derevianko
,
O.
Hemmers
,
S.
Oblad
,
P.
Glans
,
H.
Wang
,
S. B.
Whitfield
,
R.
Wehlitz
,
I. A.
Sellin
,
W. R.
Johnson
, and
D. W.
Lindle
,
Phys. Rev. Lett.
84
,
2116
(
2000
).
54.
M.
Simon
,
M. N.
Piancastelli
, and
D. W.
Lindle
, “
Hard-x-ray photoelectron spectroscopy of atoms and molecules
,” in
Hard X-Ray Photoelectron Spectroscopy (HAXPES)
, edited by
J.
Woicik
(
Springer International Publishing
,
2016
), pp.
65
110
.
55.
N. H.
List
,
J.
Kauczor
,
T.
Saue
,
H. J. A.
Jensen
, and
P.
Norman
,
J. Chem. Phys.
142
,
244111
(
2015
).
56.
A.
Szabo
and
N. S.
Ostlund
,
Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory
(
Courier Corporation
,
2012
).
57.
W. O.
Straub
, “
An elementary primer on Gaussian integrals
,” viXra preprint, viXra:1404.0026,
2014
.
58.
M. J.
Frisch
,
G. W.
Trucks
,
H. B.
Schlegel
,
G. E.
Scuseria
,
M. A.
Robb
,
J. R.
Cheeseman
,
G.
Scalmani
,
V.
Barone
,
B.
Mennucci
,
G. A.
Petersson
,
H.
Nakatsuji
,
M.
Caricato
,
X.
Li
,
H. P.
Hratchian
,
A. F.
Izmaylov
,
J.
Bloino
,
G.
Zheng
,
J. L.
Sonnenberg
,
M.
Hada
,
M.
Ehara
,
K.
Toyota
,
R.
Fukuda
,
J.
Hasegawa
,
M.
Ishida
,
T.
Nakajima
,
Y.
Honda
,
O.
Kitao
,
H.
Nakai
,
T.
Vreven
,
J. J. A.
Montgomery
,
J. E.
Peralta
,
F.
Ogliaro
,
M.
Bearpark
,
J. J.
Heyd
,
E.
Brothers
,
K. N.
Kudin
,
V. N.
Staroverov
,
R.
Kobayashi
,
J.
Normand
,
K.
Raghavachari
,
A.
Rendell
,
J. C.
Burant
,
S. S.
Iyengar
,
J.
Tomasi
,
M.
Cossi
,
N.
Rega
,
J. M.
Millam
,
M.
Klene
,
J. E.
Knox
,
J. B.
Cross
,
V.
Bakken
,
C.
Adamo
,
J.
Jaramillo
,
R.
Gomperts
,
R. E.
Stratmann
,
O.
Yazyev
,
A. J.
Austin
,
R.
Cammi
,
C.
Pomelli
,
J. W.
Ochterski
,
R. L.
Martin
,
K.
Morokuma
,
V. G.
Zakrzewski
,
G. A.
Voth
,
P.
Salvador
,
J. J.
Dannenberg
,
S.
Dapprich
,
A. D.
Daniels
,
Ö.
Farkas
,
J. B.
Foresman
,
J. V.
Ortiz
,
J.
Cioslowski
, and
D. J.
Fox
, gaussian 16, Revision A.03,
Gaussian, Inc.
,
Wallingford, CT
,
2016
.
59.
A. D.
Becke
,
J. Chem. Phys.
98
,
1372
(
1993
).
60.
D. M.
Silver
,
S.
Wilson
, and
W.
Nieuwpoort
,
Int. J. Quantum Chem.
14
,
635
(
1978
).
61.
D. M.
Silver
and
W.
Nieuwpoort
,
Chem. Phys. Lett.
57
,
421
(
1978
).
62.
J.
Mohallem
,
R.
Dreizler
, and
M.
Trsic
,
Int. J. Quantum Chem.
30
,
45
(
1986
).
63.
J.
Mohallem
and
M.
Trsic
,
J. Chem. Phys.
86
,
5043
(
1987
).
64.
H. F. M.
da Costa
,
M.
Trsic
, and
J.
Mohallem
,
Mol. Phys.
62
,
91
(
1987
).
65.
A. B. F.
da Silva
,
H. F. M.
Da Costa
, and
M.
Trsic
,
Mol. Phys.
68
,
433
(
1989
).
66.
F. E.
Jorge
,
E. V. R.
De Castro
, and
A. B. F.
Da Silva
,
J. Comput. Chem.
18
,
1565
(
1997
).
67.
F. E.
Jorge
,
E. V. R.
De Castro
, and
A. B. F.
Da Silva
,
Chem. Phys.
216
,
317
(
1997
).
68.
E. V. R.
de Castro
and
F. E.
Jorge
,
J. Chem. Phys.
108
,
5225
(
1998
).
69.
W. J.
Hehre
,
R. F.
Stewart
, and
J. A.
Pople
,
J. Chem. Phys.
51
,
2657
(
1969
).
70.
J. B.
Collins
,
P. v. R.
Schleyer
,
J. S.
Binkley
, and
J. A.
Pople
,
J. Chem. Phys.
64
,
5142
(
1976
).
71.
R.
Ditchfield
,
W. J.
Hehre
, and
J. A.
Pople
,
J. Chem. Phys.
54
,
724
(
1971
).
72.
W. J.
Hehre
,
R.
Ditchfield
, and
J. A.
Pople
,
J. Chem. Phys.
56
,
2257
(
1972
).
73.
P.
Hariharan
and
J. A.
Pople
,
Mol. Phys.
27
,
209
(
1974
).
74.
M. S.
Gordon
,
Chem. Phys. Lett.
76
,
163
(
1980
).
75.
P. C.
Hariharan
and
J. A.
Pople
,
Theor. Chim. Acta
28
,
213
(
1973
).
76.
M. M.
Francl
,
W. J.
Pietro
,
W. J.
Hehre
,
J. S.
Binkley
,
M. S.
Gordon
,
D. J.
DeFrees
, and
J. A.
Pople
,
J. Chem. Phys.
77
,
3654
(
1982
).
77.
R.
Binning
, Jr.
and
L.
Curtiss
,
J. Comput. Chem.
11
,
1206
(
1990
).
78.
J.-P.
Blaudeau
,
M. P.
McGrath
,
L. A.
Curtiss
, and
L.
Radom
,
J. Chem. Phys.
107
,
5016
(
1997
).
79.
V.
Rassolov
,
J.
Pople
,
M.
Ratner
, and
T.
Windus
,
J. Chem. Phys.
109
,
1223
(
1998
).
80.
V. A.
Rassolov
,
M. A.
Ratner
,
J. A.
Pople
,
P. C.
Redfern
, and
L. A.
Curtiss
,
J. Comput. Chem.
22
,
976
(
2001
).
81.
R. A.
Kendall
,
T. H.
Dunning
, Jr.
, and
R. J.
Harrison
,
J. Chem. Phys.
96
,
6796
(
1992
).
82.
D. E.
Woon
and
T. H.
Dunning
, Jr.
,
J. Chem. Phys.
98
,
1358
(
1993
).
83.
V.
Lebedev
,
USSR Comput. Math. Math. Phys.
16
,
10
(
1976
).
84.
See https://vuo.elettra.eu/services/elements/WebElements.html for atomic calculation of photoionization cross-sections and asymmetry parameters; accessed 13 September 2018 (archived by WebCite® at http://www.webcitation.org/72OwdTOR5).
85.
J. A. R.
Samson
and
W. C.
Stolte
,
J. Electron Spectrosc.
123
,
265
(
2002
).
86.
J.
Berkowitz
,
Atomic and Molecular Photoabsorption: Absolute Total Cross Sections
(
Academic Press
,
2002
).
87.
N.
Sakamoto
,
H.
Tsuchida
,
T.
Kato
,
S.
Kamakura
,
H.
Ogawa
,
K.
Ishii
,
I.
Murakami
,
D.
Kato
,
H.
Sakaue
, and
M.
Kato
, “
Atomic and molecular numerical databases, photoabsorption database
,” http://dpc.nifs.ac.jp/photoab; accessed 13 September 2018 (archived by WebCite® at http://www.webcitation.org/72Outg34w).
88.
See https://physics.nist.gov/PhysRefData/Xcom/Text/intro.html for Nist xcom, photon cross sections database; accessed 13 September 2018 (archived by WebCite® at http://www.webcitation.org/72OwB0BWM).
89.
V.
Schmidt
,
Electron Spectrometry of Atoms Using Synchrotron Radiation
(
Cambridge University Press
,
2005
), Vol. 6.
90.
J. A. R.
Samson
,
Z.
He
,
L.
Yin
, and
G.
Haddad
,
J. Phys. B: At., Mol. Opt. Phys.
27
,
887
(
1994
).
91.
J.
Berkowitz
,
J. Chem. Phys.
70
,
2819
(
1979
).
92.
B.
Krässig
,
M.
Jung
,
D.
Gemmell
,
E.
Kanter
,
T.
LeBrun
,
S.
Southworth
, and
L.
Young
,
Phys. Rev. Lett.
75
,
4736
(
1995
).
93.
M.
Jung
,
B.
Krässig
,
D.
Gemmell
,
E.
Kanter
,
T.
LeBrun
,
S.
Southworth
, and
L.
Young
,
Phys. Rev. A
54
,
2127
(
1996
).

Supplementary Material