The one-particle Green’s function method in the Dirac–Hartree–Fock framework. II. Third-order valence ionization energies of the noble gases, CO and ICN

Interatomic and Intermolecular Coulombic Decay

Interatomic or intermolecular Coulombic decay (ICD) is a nonlocal electronic decay mechanism occurring in weakly bound matter. In an ICD process, energy released by electronic relaxation of an excited atom or molecule leads to ionization of a neighboring one via Coulombic electron interactions. ICD has been predicted theoretically in the mid nineties of the last century, and its existence has been confirmed experimentally approximately ten years later. Since then, a number of fundamental and applied aspects have been studied in this quickly growing field of research. This review provides an introduction to ICD and draws the connection to related energy transfer and ionization processes. The theoretical approaches for the description of ICD as well as the experimental techniques developed and employed for its investigation are described. The existing body of literature on experimental and theoretical studies of ICD processes in different atomic and molecular systems is reviewed.

Effect of spin-orbit coupling on decay widths of electronic decay processes

Auger-Meitner processes are electronic decay processes of energetically low-lying vacancies. In these processes, the vacancy is filled by an electron of an energetically higher lying orbital, while another electron is simultaneously emitted to the continuum. In low-lying orbitals, relativistic effects can not, even for light elements, be neglected. At the same time, lifetime calculations are computationally expensive. In this context, we investigate which effect spin-orbit coupling has on Auger-Meitner decay widths and aim for a rule of thumb for the relative decay widths of initial states split by spin-orbit coupling. We base this rule of thumb on Auger-Meitner decay widths of Sr4p^-1 and Ra6p^-1 obtained by relativistic FanoADC-Stieltjes calculations and validate it against Auger-Meitner decay widths from the literature.

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

Investigation of the valence ionization spectrum of chromium carbonyl using an ab initio quantum dynamical approach

The nuclear dynamics of the chromium carbonyl cation following an ionization process corresponding to the ²T_2g ← ¹A_1g transition is studied theoretically, for the first time, using a fully quantal approach as well as high levels of the ab initio electronic structure and semiempirical density functional theory (DFT) methods. The photoelectron spectrum is calculated by the construction of a Hamiltonian model, in which the two totally symmetric modes ν₁₉ (the Cr-C stretching mode) and ν₃₉ (the C-O stretching mode) together with the spin-orbit (SO) coupling up to the zeroth-order SO splitting are treated. The potential energy curves along these two vibrational modes are computed by using the DFT. The simulated photoelectron spectrum is found to be in good agreement with the corresponding experimental one, leading to the conclusion that the potential energy surfaces and the diabatic population analysis are accurately determined. Our calculation confirms that the vibrational modes ν₁₉ and ν₃₉ are the vibrational progression of the valence ionization spectrum of the chromium carbonyl cation.

Fragmentation of mercury compounds under ultraviolet light irradiation

Ultraviolet light induced photofragmentation of mercury compounds is studied experimentally with electron energy resolved photoelectron-photoion coincidence techniques and theoretically with computational quantum chemical methods. A high resolution photoelectron spectrum using synchrotron radiation is presented. Fragmentation of the molecule is studied subsequent to ionization to the atomic-mercury-like d orbitals. State dependent fragmentation behaviour is presented and specific reactions for dissociation pathways are given. The fragmentation is found to differ distinctly in similar orbitals of different mercury compounds.

Communication: electron transfer mediated decay enabled by spin-orbit interaction in small krypton/xenon clusters

In this work we study the influence of relativistic effects, in particular spin-orbit coupling, on electronic decay processes in KrXe2 clusters of various geometries. For the first time it is shown that inclusion of spin-orbit coupling has decisive influence on the accessibility of a specific decay pathway in these clusters. The radiationless relaxation process is initiated by a Kr 4s ionization followed by an electron transfer from xenon to krypton and a final second ionization of the system. We demonstrate the existence of competing electronic decay pathways depending in a subtle way on the geometry and level of theory. For our calculations a fully relativistic framework was employed where omission of spin-orbit coupling leads to closing of two decay pathways. These findings stress the relevance of an adequate relativistic description for clusters with heavy elements and their fragmentation dynamics.

The band 12 issue of norbornane: a study of higher shake-up states

In line with a recent study of the electronic structure of the cage compound norbornane (J. Chem. Phys. 121 (2004), 10525; J. Phys. Chem. A 109 (2005), 4267), symmetry adapted cluster expansion configuration interaction (SAC-CI) general R calculations have been performed and compared with results obtained by the third order algebraic diagrammatic construction scheme [ADC(3)]. Comparison has been made with previously performed electron momentum spectroscopy (EMS) and ultraviolet photo-electron measurements. The region around ~25 eV (band 12), characterized by an elaborated band in the EMS spectrum which is missing in previous Green's function and ADC calculations, is investigated. This study is completed with outer-valence Green's function (OVGF) and SAC-CI/SD-R calculations, and results are obtained by employing (single and double) ionization extended second order ADC [ADC(2)-x]. Since ADC(3) only includes 2h-1p shake-up states, while SAC-CI general-R also includes higher order states, the agreement between both methods assures that the higher order shake-up states do not play an important role in the ionization spectrum of norbornane. While the band-12 issue of norbornane is therefore still open for further discussion, a tentative description in terms of ultrafast nuclear dynamical effects and autoionization processes has become more plausible.

Toward improved density functionals for the correlation energy

Eleven density functionals, including some of the most widely used ones, are tested on their ability to predict nonrelativistic, electron correlation energies for the 17 atoms from He to Ar, the 17 cations from Li(+) to K(+), and 11 (1)S state atoms from Ca to Rn. They all lead to relatively poor predictions for the heavier atoms. Reparametrization of these functionals improves their performance for light atoms but does not alleviate their problems with the heavier, closed-shell atoms. Several novel, few-parameter, density functionals for the correlation energy are developed heuristically. Four new functionals lead to qualitatively improved predictions for the heavier atoms without unreasonably compromising accuracy for the lighter atoms. Further progress would be facilitated by reliable estimates of electron correlation energies for more atoms, particularly heavy ones.

Implementation and use of a direct, partially integral-driven non-Dyson propagator method for molecular ionization

The Green's function ADC(3) scheme has been for many years a successful method to predict theoretically the ionization (and electron affinity) spectrum of molecules. However, a dramatic enhancement of the method's power has come only recently, with the development of an approximation method to the one-particle Green's function which does not make direct use of the Dyson equation. In the present work, we present an efficient computer implementation of this novel approach, with first comparative tests demonstrating its enormous computational advantage over the conventional approach.

Ionization spectra and electronic decay in small iodide clusters: fully relativistic results

Singly ionized systems in high-lying energetic final states can stabilize themselves via various electronic decay mechanisms. With increasing system size interatomic and intermolecular processes dominate over intra-atomic (Auger) decay channels. For the small (HI)(2) and (LiI)(2) clusters fully relativistic ionization spectra are calculated and the subsequent electronic decay of the cations is investigated. Due to the presence of the iodine atom a fully relativistic description is mandatory and was performed by the algebraic diagrammatic construction technique in its four-component form. The lifetimes of the singly ionized final states are estimated by the application of Weisskopf-Wigner [Z. Phys. 63, 54 (1930)] theory.

PtF62− dianion and its detachment spectrum: A fully relativistic study

In this work we calculate the photoelectron spectrum of the PtF(6)2- dianion by application of the third-order Dirac-Hartree-Fock one-particle propagator technique. Relativistic effects and electron correlation are hereby treated on a consistent theoretical basis which is mandatory for systems containing heavy elements. A PtF6(2-) gas phase photoelectron spectrum is not yet available and our calculations therefore have predictive character. As it is characteristic for dianionic systems a strong dependence on basis set size and molecular geometry is observed. In contrast to the already calculated PtCl(6)2- photoelectron spectrum no valence orbital inversion due to strong interplay of spin-orbit coupling and electron correlation is observed. Furthermore an unusually strong spin-orbit splitting was found for the sigma-type subvalence 1t1u molecular spinor despite its very small platinum p population. The double ionization threshold is strongly lowered by relativistic effects now enabling an interatomic Coulombic decay process after ionization from the sigma-bonding orbitals. The results stress the importance of spin-orbit coupling for the understanding of the spectral structure which cannot be reproduced by a scalar-relativistic treatment only.

Molecular ionization energies and ground- and ionic-state properties using a non-Dyson electron propagator approach

An earlier proposed propagator method for the treatment of molecular ionization is tested in first applications. The method referred to as the non-Dyson third-order algebraic-diagrammatic construction [nD-ADC(3)] approximation for the electron propagator represents a computationally promising alternative to the existing Dyson ADC(3) method. The advantage of the nD-ADC(3) scheme is that the (N+/-1)-electronic parts of the one-particle Green's function are decoupled from each other and the corresponding equations can be solved separately. For a test of the method the nD-ADC(3) results for the vertical ionization transitions in C(2)H(4), CO, CS, F(2), H(2)CO, H(2)O, HF, N(2), and Ne are compared with available experimental and theoretical data including results of full configuration interaction (FCI) and coupled cluster computations. The mean error of the nD-ADC(3) ionization energies relative to the experimental and FCI results is about 0.2 eV. The nD-ADC(3) method, scaling as n(5) with the number of orbitals, requires the solution of a relatively simple Hermitian eigenvalue problem. The method renders access to ground-state properties such as dipole moments. Moreover, also one-electron properties of (N+/-1) electron states can now be studied as a consequence of a specific intermediate-state representation (ISR) formulation of the nD-ADC approach. Corresponding second-order ISR equations are presented.

Effect of relativity on the ionization spectra of the xenon fluorides XeFn (n=2, 4, 6)

Noble gas compounds exhibit special chemical bonding situations and have been investigated by various spectroscopic and theoretical techniques. In this work we calculate the ionization spectra of the xenon fluorides (XeF2,XeF4, and XeF6) in the valence and subvalence (down to Xe 4d) areas by application of the recently developed Dirac-Hartree-Fock one-particle propagator technique. In this technique, the relativistic (four-component) and electron correlation effects are computed simultaneously. The xenon compounds show considerable spin-orbit splitting strongly influencing the photoelectron spectrum not reproducible in prior calculations. Comparison to one-component methods is made and the occurring satellite structures are interpreted. The satellite structures can be attributed either to the breakdown of the one-particle picture or to a reflection of intra-atomic and interatomic Auger decay processes within the molecule.

Remarkable interplay of electron correlation and relativity in the photodetachment spectrum of PtCl62−

In this work we calculate the photoelectron spectrum of the PtCl(6) (2-) dianion by application of the recently developed third-order Dirac-Hartree-Fock implementation of the one-particle propagator technique allowing for a consistent treatment of spin-orbit and scalar relativistic effects together with electron correlation. For PtCl(6) (2-) a gas phase photoelectron spectrum is available showing clearly discernible structures not reproducible by a nonrelativistic or purely scalar-relativistic computation. A population analysis of the valence orbitals allows for an assignment of the photoelectron peaks and reveals the strong influence of relativity in combination with electron correlation.