Time dependent density functional study of the photoionization dynamics of SF6

Photoionization and core resonances from range-separated density-functional theory: General formalism and example of the beryllium atom

We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely, the time-dependent range-separated hybrid (TDRSH) scheme, which uses a global range-separation parameter, and the time-dependent locally range-separated hybrid (TDLRSH), which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a non-orthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing, in particular, on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum.

Photoionization Observables from Multi-Reference Dyson Orbitals Coupled to B-Spline DFT and TD-DFT Continuum

We present a theoretical model to compute the accurate photoionization dynamical parameters (cross-sections, asymmetry parameters and orbital, or cross-section, ratios) from Dyson orbitals obtained with the multi-state complete active space perturbation theory to the second order (MS-CASPT2) method. Our new implementation of Dyson orbitals in OpenMolcas takes advantage of the full Abelian symmetry point group and has the corrected normalization. The Dyson orbitals are coupled to an accurate description of the electronic continuum obtained with a multicentric B-spline basis at the DFT and TD-DFT levels. Two prototype diatomic molecules, i.e., CS and SiS, have been chosen due to their smallness, which hides important correlation effects. These effects manifest themselves in the appearance of well-characterized isolated satellite bands in the middle of the valence region. The rich satellite structures make CS and SiS the perfect candidates for a computational study based on our highly accurate MS-CASPT2/B-spline TD-DFT protocol.

Cation and Anion Formation in F 1s-Excited SF6-Clusters

Cation and anion formation from molecular SF6 and SF6 clusters is reported in the F 1s-excitation regime (680–720 eV) using tunable synchrotron radiation. Stable products are identified by mass spectrometry, while the tunablility of the soft X-rays source is used to determine mechanistic details on the formation of the detected charged products. Characteristic differences are observed between product formation from the isolated molecule as well as from clusters for both, cations and anions. Specifically, dication formation, such as SF2 2+, is only observed for molecular precursors, which are formed with high efficiency at the F 1s → t2g-transition. SF6 − is efficiently formed from clusters upon F 1s → t1u-transition. Mechanistic details on the formation of stable products are discussed.

Photoionization Dynamics of the Tetraoxo Complexes OsO4 and RuO4

The photoionization dynamics of OsO₄ and RuO₄, chosen as model systems of small-size mononuclear heavy-metal complexes, has been theoretically studied by the time-dependent density functional theory (TDDFT). Accurate experimental measurements of photoionization dynamics as a benchmarking test for the theory are reported for the photoelectron asymmetry parameters of outer valence ionizations of OsO₄, measured in the 17–90 eV photon energy range. The theoretical results are in good agreement with the available experimental data. The observed dynamical behavior of partial cross sections and asymmetry parameters has been related to both the coupling to the continuum of discrete excited states, giving strong modulations in the photon energy dependency, and the atomic composition of the initial ionized states, which determines the rate of decay of ionization probability for increasing excitation energies. Overall, an extensive analysis of the photoionization dynamics for valence and core orbitals is presented, showing good agreement with all the available experimental data. This provides confidence for the validity of the TDDFT approach in describing photoionization of heavy transition element compounds, with the perspective of being used for larger systems. Further experimental work is suggested for RuO₄ to gather evidence of the sensitivity of the theoretical method to the nature of the metal atom.

In this work, the time-dependent density functional theory is used to calculate the photoionization dynamics of the valence and core ionizations of OsO₄ and RuO₄ complexes, which are chosen as model systems of small-size mononuclear heavy-metal complexes. Accurate experimental measurements of the branching ratios and photoelectron asymmetry parameters of outer valence ionizations of OsO₄ provide a sound benchmark of the computational methodology.

Ultrafast charge dynamics in glycine induced by attosecond pulses

The combination of attosecond pump-probe techniques with mass spectrometry methods has recently led to the first experimental demonstration of ultrafast charge dynamics in a biomolecule, the amino acid phenylalanine [Calegari et al., Science, 2014, 346, 336]. Using an extension of the static-exchange density functional theory (DFT) method, the observed dynamics was explained as resulting from the coherent superposition of ionic states produced by the broadband attosecond pulse. Here, we have used the static-exchange DFT method to investigate charge migration induced by attosecond pulses in the glycine molecule. We show that the observed dynamics follows patterns similar to those previously found in phenylalanine, namely that charge fluctuations occur all over the molecule and that they can be explained in terms of a few typical frequencies of the system. We have checked the validity of our approach by explicitly comparing with the photoelectron spectra obtained in synchrotron radiation experiments and with the charge dynamics that follows the removal of an electron from a given molecular orbital, for which fully correlated ab initio results are available in the literature. From this comparison, we conclude that our method provides an accurate description of both the coherent superposition of cationic states generated by the attosecond pulse and its subsequent time evolution. Hence, we expect that the static-exchange DFT method should perform equally well for other medium-size and large molecules, for which the use of fully correlated ab initio methods is not possible.

Interchannel coupling effects in the valence photoionization of SF6

The complex Kohn and polyatomic Schwinger variational techniques have been employed to illustrate the interchannel coupling correlation effects in the valence photoionization dynamics of SF6. Partial photoionization cross sections and asymmetry parameters of six valence subshells (1t1g, 5t1u, 1t2u, 3eg, 1t2g, 4t1u) are discussed in the framework of several theoretical and experimental studies. The complex Kohn results are in rather good agreement with experimental results, indicative of the fact that the interchannel coupling effects alter the photoionization dynamics significantly. We find that the dominant effect of interchannel coupling is to reduce the magnitude of shape resonant cross sections near the threshold and to induce resonant features in other channels to which resonances are coupled. The long-standing issue concerning ordering of the valence orbitals is addressed and confirmed 4t1u (6)1t2g (6)3eg (4)(5t1u (6)+1t2u (6)) 1t1g (6) as the most likely ordering.

Multi-channel electronic and vibrational dynamics in polyatomic resonant high-order harmonic generation

High-order harmonic generation in polyatomic molecules generally involves multiple channels of ionization. Their relative contribution can be strongly influenced by the presence of resonances, whose assignment remains a major challenge for high-harmonic spectroscopy. Here we present a multi-modal approach for the investigation of unaligned polyatomic molecules, using SF6 as an example. We combine methods from extreme-ultraviolet spectroscopy, above-threshold ionization and attosecond metrology. Fragment-resolved above-threshold ionization measurements reveal that strong-field ionization opens at least three channels. A shape resonance in one of them is found to dominate the signal in the 20-26 eV range. This resonance induces a phase jump in the harmonic emission, a switch in the polarization state and different dynamical responses to molecular vibrations. This study demonstrates a method for extending high-harmonic spectroscopy to polyatomic molecules, where complex attosecond dynamics are expected.

Determination of chiral asymmetries in the valence photoionization of camphor enantiomers by photoelectron imaging using tunable circularly polarized light

An electron imaging technique has been used to study the full angular distribution of valence photoelectrons produced from enantiomerically pure molecular beams of camphor when these are photoionized with circularly polarized light. In addition to the familiar beta parameter, this provides a new chiral term, taking the form of an additional cosine function in the angular distribution which consequently displays a forward-backward electron ejection asymmetry. Several ionization channels have been studied using synchrotron radiation in the 8.85-26 eV photon energy range. With alternating left and right circularly polarized radiations the photoelectron circular dichroism (PECD) in the angular distribution can be measured and shows some strong dynamical variations with the photon energy, depending in sign and intensity on the ionized orbital. For all orbitals the measured PECD has a quite perfect antisymmetry when switching between R and S enantiomers, as expected from theory. In the HOMO(-1) channel the PECD chiral asymmetry curves show a double maxima reaching nearly 10% close to threshold, and peaking again at approximately 20% some 11 eV above threshold. This is attributed to a resonance that is also visible in the beta parameter curve. Newly optimized CMS-Xalpha photoionization dynamics calculations are also presented. They are in reasonably good agreement with the experimental data, including in the very challenging threshold regions. These calculations show that PECD in such randomly oriented samples can be understood in the electric dipole approximation and that, unlike the case pertaining in core-shell ionization-where a highly localized achiral initial orbital means that the dichroism arises purely as a final state scattering effect-in valence shell ionization there is a significant additional influence contributed by the initial orbital density.