Table-Top X-ray Spectroscopy of Benzene Radical Cation

The x-ray absorption spectrum of the tert-butyl radical: An experimental and computational investigation

We report the x-ray absorption spectrum (XAS) of the tert-butyl radical, C4H9. The radical was generated pyrolytically from azo-tert-butane, and the XAS of the pure radical was obtained by subtraction of spectra recorded at different temperatures. The bands in the XAS were assigned by ab initio calculations that are in very good agreement with the experimental data. The lowest energy signal in the XAS is assigned to the C1s electron transition from the central carbon atom to the singly occupied molecular orbital (SOMO), while higher transitions correspond to C1s excitations from terminal carbon atoms. Furthermore, we investigated the fragmentation of the radical following resonant C1s excitation by electron-ion-coincidence spectroscopy. Several fragmentation channels were identified. The C1s excitation of the terminal carbons is associated with a stronger fragmentation tendency compared to the lowest C1s excitation of the central carbon into the SOMO. For this core excited state, we still observe an intact parent ion, C4H9+, and a comparatively higher tendency to dissociate into CH3+ + C3H6+.

Ultrafast dynamics of fluorene initiated by highly intense laser fields

We present an investigation of the ultrafast dynamics of the polycyclic aromatic hydrocarbon fluorene initiated by an intense femtosecond near-infrared laser pulse (810 nm) and probed by a weak visible pulse (405 nm). Using a multichannel detection scheme (mass spectra, electron and ion velocity-map imaging), we provide a full disentanglement of the complex dynamics of the vibronically excited parent molecule, its excited ionic states, and fragments. We observed various channels resulting from the strong-field ionization regime. In particular, we observed the formation of the unstable tetracation of fluorene, above-threshold ionization features in the photoelectron spectra, and evidence of ubiquitous secondary fragmentation. We produced a global fit of all observed time-dependent photoelectron and photoion channels. This global fit includes four parent ions extracted from the mass spectra, 15 kinetic-energy-resolved ionic fragments extracted from ion velocity map imaging, and five photoelectron channels obtained from electron velocity map imaging. The fit allowed for the extraction of 60 lifetimes of various metastable photoinduced intermediates.

Reinvented: An Attosecond Chemist

Attosecond science requires a substantial rethinking of how to make measurements on very short timescales; how to acquire the necessary equipment, technology, and personnel; and how to build a set of laboratories for such experiments. This entails a rejuvenation of the author in many respects, in the laboratory itself, with regard to students and postdocs, and in generating funding for research. It also brings up questions of what it means to do attosecond science, and the discovery of the power of X-ray spectroscopy itself, which complements the short timescales addressed. The lessons learned, expressed in the meanderings of this autobiographical article, may be of benefit to others who try to reinvent themselves.

Isomer-specific photofragmentation of C3H3+ at the carbon K-edge

Individual fingerprints of different isomers of C3H3+ cations have been identified by studying photoionization, photoexcitation, and photofragmentation of C3H3+ near the carbon K-edge. The experiment was performed employing the photon-ion merged-beams technique at the photon-ion spectrometer at PETRA III (PIPE). This technique is a variant of near-edge X-ray absorption fine-structure spectroscopy, which is particularly sensitive to the 1s → π* excitation. The C3H3+ primary ions were generated by an electron cyclotron resonance ion source. C3Hn2+ product ions with n = 0, 1, 2, and 3 were observed for photon energies in the range of 279.0 eV to 295.2 eV. The experimental spectra are interpreted with the aid of theoretical calculations within the framework of time-dependent density functional theory. To this end, absorption spectra have been calculated for three different constitutional isomers of C3H3+. We find that our experimental approach offers a new possibility to study at the same time details of the electronic structure and of the geometry of molecular ions such as C3H3+.

Spin Coupling Effect on Geometry-Dependent X-Ray Absorption of Diradicals

We theoretically investigate the influence of diradical electron spin coupling on the time-resolved X-ray absorption spectra of the photochemical ring opening of furanone. We predict geometry-dependent carbon K-edge signals involving transitions from core orbitals to both singly and unoccupied molecular orbitals. The most obvious features of the ring opening come from the carbon atom directly involved in the bond breaking through its transition to both the newly formed singly occupied and the available lowest unoccupied molecular orbitals (SOMO and LUMO, respectively). In addition to this primary feature, the singlet spin coupling of four unpaired electrons that arises in the core-to-LUMO states creates additional geometry dependence in some spectral features with both oscillator strengths and relative excitation energies varying observably as a function of the ring opening. We attribute this behavior to a spin-occupancy-induced selection rule, which occurs when singlet spin coupling is enforced in the diradical state. Notably, one of these geometry-sensitive core-to-LUMO transitions excites core electrons from a backbone carbon not involved in the bond breaking, providing a novel nonlocal X-ray probe of chemical dynamics arising from electron spin coupling.

Describing nuclear quantum effects in vibrational properties using molecular dynamics with Wigner sampling

In this work we discuss the generally applicable Wigner sampling and introduce a new, simplified Wigner sampling method, for computationally effective modeling of molecular properties containing nuclear quantum effects and vibrational anharmonicity. For various molecular systems test calculations of (a) vibrationally averaged rotational constants, (b) vibrational IR spectra and (c) photoelectron spectra have been performed. The performance of Wigner sampling has been assessed by comparing with experimental data and with results of other theoretical models, including harmonic and VPT2 approximations. The developed simplified Wigner sampling method shows advantages in application to large and flexible molecules.

Signatures of diradicals in x-ray absorption spectroscopy

Theoretical simulations are critical to analyze and interpret the x-ray absorption spectrum of transient open-shell species. In this work, we propose a model of the many-body core-excited states of symmetric diradicals. We apply this model to analyze the carbon K-edge transitions of o-, m-, and p-benzyne, three organic diradicals with diverse and unusual electronic structures. The predictions of our model are compared with high-level multireference computations of the K-edge spectrum of the benzynes obtained with the driven similarity renormalization group truncated to third order. Our model shows the importance of a many-body treatment of the core-excited states of the benzynes and provides a theoretical framework to understand which properties of the ground state of these diradicals can be extracted from their x-ray spectrum.

A study of core-excited states of organic molecules computed with the generalized active space driven similarity renormalization group

This work examines the accuracy and precision of x-ray absorption spectra computed with a multireference approach that combines generalized active space (GAS) references with the driven similarity renormalization group (DSRG). We employ the x-ray absorption benchmark of organic molecule (XABOOM) set, consisting of 116 transitions from mostly organic molecules [Fransson et al., J. Chem. Theory Comput. 17, 1618 (2021)]. Several approximations to a full-valence active space are examined and benchmarked. Absolute excitation energies and intensities computed with the GAS-DSRG truncated to second-order in perturbation theory are found to systematically underestimate experimental and reference theoretical values. Third-order perturbative corrections significantly improve the accuracy of GAS-DSRG absolute excitation energies, bringing the mean absolute deviation from experimental values down to 0.32 eV. The ozone molecule and glyoxylic acid are particularly challenging for second-order perturbation theory and are examined in detail to assess the importance of active space truncation and intruder states.

Real-time observation of the Woodward-Hoffmann rule for 1,3-cyclohexadiene by femtosecond soft X-ray transient absorption

The stereochemistry of pericyclic reactions is explained by orbital symmetry conservation, referred to as the Woodward-Hoffmann (WH) rule. Although this rule has been verified using the structures of reactants and products, the temporal evolution of the orbital symmetry during the reaction has not been clarified. Herein, we used femtosecond soft X-ray transient absorption spectroscopy to elucidate the thermal pericyclic reaction of 1,3-cyclohexadiene (CHD) molecules, i.e., their isomerization to 1,3,5-hexatriene. In the present experimental scheme, the ring-opening reaction is driven by the thermal vibrational energy induced by photoexcitation to the Rydberg states at 6.2 eV and subsequent femtosecond relaxation to the ground state of CHD molecules. The direction of the ring opening, which can be conrotatory or disrotatory, was the primary focus, and the WH rule predicts the disrotatory pathway in the thermal process. We observed the shifts in K-edge absorption of the carbon atom from the 1s orbital to vacant molecular orbitals around 285 eV at a delay between 340 and 600 fs. Furthermore, a theoretical investigation predicts that the shifts depend on the molecular structures along the reaction pathways and the observed shifts in induced absorption are attributed to the structural change in the disrotatory pathway. This confirms that the orbital symmetry is dynamically conserved in the ring-opening reaction of CHD molecules as predicted using the WH rule.

Jahn-Teller effects in initial and final states: high-resolution X-ray absorption, photoelectron and Auger spectroscopy of allene

Carbon K-edge resonant Auger spectra of gas-phase allene following excitation of the pre-edge 1s → π* transitions are presented and analysed with the support of EOM-CCSD/cc-pVTZ calculations. X-Ray absorption (XAS), X-ray photoelectron (XPS), valence band and non-resonant Auger spectra are also reanalysed with a series of computational approaches. The results presented demonstrate the importance of including nuclear ensemble effects for simulating X-ray observables and as an effective strategy for capturing Jahn-Teller effects in spectra.

Signatures of the Bromine Atom and Open-Shell Spin Coupling in the X-ray Spectrum of the Bromobenzene Cation

Tabletop X-ray spectroscopy measurements at the carbon K-edge complemented by ab initio calculations are used to investigate the influence of the bromine atom on the carbon core-valence transitions in the bromobenzene cation (BrBz⁺). The electronic ground state of the cation is prepared by resonance-enhanced two-photon ionization of neutral bromobenzene (BrBz) and probed by X-rays produced by high-harmonic generation (HHG). Replacing one of the hydrogen atoms in benzene with a bromine atom shifts the transition from the 1s_C* orbital of the carbon atom (C*) bonded to bromine by ∼1 eV to higher energy in the X-ray spectrum compared to the other carbon atoms (C). Moreover, in BrBz⁺, the X-ray spectrum is dominated by two relatively intense transitions, 1s_C→π* and 1s_C*→σ*(C*-Br), where the second transition is enhanced relative to the neutral BrBz. In addition, a doublet peak shape for these two transitions is observed in the experiment. The 1s_C→π* doublet peak shape arises due to the spin coupling of the unpaired electron in the partially vacant π orbital (from ionization) with the two other unpaired electrons resulting from the transition from the 1s_C core orbital to the fully vacant π* orbitals. The 1s_C*→σ* doublet peak shape results from several transitions involving σ* and vibrational C*-Br mode activations following the UV ionization, which demonstrates the impact of the C*-Br bond length on the core-valence transition as well as on the relaxation geometry of BrBz⁺.

The Auger spectrum of benzene

We present an ab initio computational study of the Auger electron spectrum of benzene. Auger electron spectroscopy exploits the Auger-Meitner effect, and although it is established as an analytic technique, the theoretical modeling of molecular Auger spectra from first principles remains challenging. Here, we use coupled-cluster theory and equation-of-motion coupled-cluster theory combined with two approaches to describe the decaying nature of core-ionized states: (i) Feshbach-Fano resonance theory and (ii) the method of complex basis functions. The spectra computed with these two approaches are in excellent agreement with each other and also agree well with experimental Auger spectra of benzene. The Auger spectrum of benzene features two well-resolved peaks at Auger electron energies above 260 eV, which correspond to final states with two electrons removed from the 1e_1g and 3e_2g highest occupied molecular orbitals. At lower Auger electron energies, the spectrum is less well resolved, and the peaks comprise multiple final states of the benzene dication. In line with theoretical considerations, singlet decay channels contribute more to the total Auger intensity than the corresponding triplet decay channels.

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the ¹ππ* (S₂) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the ¹nπ* (S₁) and ³ππ* (T₁) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576-16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core-valence separated (FC-CVS) EOM-CCSD approach - a successful variant of EOM-CCSD tailored towards core excitations.

X-ray absorption of molecular cations-a new challenge for electronic structure theory

In this paper we put forward some historical notes on the development of computational chemistry toward applications of x-ray spectroscopies. We highlight some of the important contributions by Enrico Clementi as method and program developer and as a supporter of this branch of computational research. We bring up a modern example based on the very recent experimental development of x-ray absorption of cationic molecules. As we show this spectroscopy poses new challenges for electronic structure theory and the electron correlation problem.

Jahn-Teller Distortion and Dissociation of CCl4+ by Transient X-ray Spectroscopy Simultaneously at the Carbon K- and Chlorine L-Edge

X-ray Transient Absorption Spectroscopy (XTAS) and theoretical calculations are used to study CCl₄⁺ prepared by 800 nm strong-field ionization. XTAS simultaneously probes atoms at the carbon K-edge (280–300 eV) and chlorine L-edge (195–220 eV). Comparison of experiment to X-ray spectra computed by orbital-optimized density functional theory (OO-DFT) indicates that after ionization, CCl₄⁺ undergoes symmetry breaking driven by Jahn–Teller distortion away from the initial tetrahedral structure (T_d) in 6 ± 2 fs. The resultant symmetry-broken covalently bonded form subsequently separates to a noncovalently bound complex between CCl₃⁺ and Cl over 90 ± 10 fs, which is again predicted by theory. Finally, after more than 800 fs, L-edge signals for atomic Cl are observed, indicating dissociation to free CCl₃⁺ and Cl. The results for Jahn–Teller distortion to the symmetry-broken form of CCl₄⁺ and formation of the Cl–CCl⁺₃ complex characterize previously unobserved new species along the route to dissociation.

Dynamics of CCl₄⁺ prepared by 800 nm strong-field ionization, as studied with X-ray transient absorption spectroscopy (XTAS) and quantum chemical calculations.

Exploring vibronic coupling in the benzene radical cation and anion with different levels of the GW approximation

The linear vibronic coupling constants of the benzene radical cation and anion have been obtained with different levels of the GW approximation, including G₀W₀, eigenvalue self-consistent GW, and quasiparticle self-consistent GW, as well as DFT with the following exchange-correlation functionals: BLYP, B3LYP, CAM-B3LYP, tuned CAM-B3LYP, and an IP-tuned CAM-B3LYP functional. The vibronic coupling constants were calculated numerically using the gradients of the eigenvalues of the degenerate HOMOs and LUMOs of the neutral benzene molecule for DFT, while the numerical gradients of the quasiparticle energies were used in the case of GW. The results were evaluated against those of high level wave function methods in the literature, and the approximate self-consistent GW methods and G₀W₀ with long-range corrected functionals were found to yield the best results on the whole.

o-Semiquinone radical anion isolated as an amorphous porous solid

The use of metal cations is a commonly applied strategy to create S > 1/2 stable molecular systems containing semiquinone radicals. Persistent mono-semiquinonato complexes of diamagnetic metal ions (S = 1/2) have been hitherto less common and mostly limited to the complexes of heavy metal ions. In this work, a mono-semiquinonato complex of aluminum, derived from 1,2-dihydroxybenzene, is obtained using a surprisingly short and uncomplicated procedure. The isolated product is an amorphous and porous solid that exhibits very good stability under ambient conditions. To characterise its molecular and electronic structure, 9.7, 34 and 406 GHz EPR spectroscopy was used in concert with computational techniques (DFT and DLPNO-CCSD). It was revealed that the radical complex is composed of two chemically equivalent aluminum cations and two catechol-like ligands with the unpaired electron uniformly distributed between the two organic molecules. The good stability and porous structure make this complex applicable in heterogeneous aerobic reactions.

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation

Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized ¹B_2u (ππ*) (S₂) and ¹B_3u (nπ*) (S₁) states, the participation of the optically dark ¹A_u (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite ¹A_u (nπ*) and ¹B_3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The ¹A_u (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules

Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the M_S = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.

Substituent Effect in the Cation Radicals of Monosubstituted Benzenes

In 30 monosubstituted benzene cation radicals, studied at the ωB97XD/aug-cc-pVTZ level, the phenyl rings usually adopt a compressed form, but a differently compressed form-equivalent to an elongated one-may coexist. The computational and literature ionization potentials are well correlated. The geometrical and magnetic aromaticity, estimated using HOMA and NICS indices, show the systems to be structurally aromatic but magnetically antiaromatic or only weakly aromatic. The partial charge is split between the substituent and ring and varies the most at C(ipso). In the ring, the spin is 70%, concentrated equally at the C(ipso) and C(p) atoms. The sEDA(D) and pEDA(D) descriptors of the substituent effect in cation radicals, respectively, were determined. In cation radicals, the substituent effect on the σ-electron system is like that in the ground state. The effect on the π-electron systems is long-range, and its propagation in the radical quinone-like ring is unlike that in the neutral molecules. The pEDA(D) descriptor correlates well with the partial spin at C(ipso) and C(p) and weakly with the HOMA(D) index. The correlation of the spin at the ring π-electron system and the pEDA(D) descriptor shows that the electron charge supplied to the ring π-electron system and the spin flow oppositely.

Feshbach-Fano approach for calculation of Auger decay rates using equation-of-motion coupled-cluster wave functions. I. Theory and implementation

X-ray absorption creates electron vacancies in the core shell. These highly excited states often relax by Auger decay-an autoionization process in which one valence electron fills the core hole and another valence electron is ejected into the ionization continuum. Despite the important role of Auger processes in many experimental settings, their first-principles modeling is challenging, even for small systems. The difficulty stems from the need to describe many-electron continuum (unbound) states, which cannot be tackled with standard quantum-chemistry methods. We present a novel approach to calculate Auger decay rates by combining Feshbach-Fano resonance theory with the equation-of-motion coupled-cluster single double (EOM-CCSD) framework. We use the core-valence separation scheme to define projectors into the bound (square-integrable) and unbound (continuum) subspaces of the full function space. The continuum many-body decay states are represented by products of an appropriate EOM-CCSD state and a free-electron state, described by a continuum orbital. The Auger rates are expressed in terms of reduced quantities, two-body Dyson amplitudes (objects analogous to the two-particle transition density matrix), contracted with two-electron bound-continuum integrals. Here, we consider two approximate treatments of the free electron: a plane wave and a Coulomb wave with an effective charge, which allow us to evaluate all requisite integrals analytically; however, the theory can be extended to incorporate a more sophisticated description of the continuum orbital.

Interplay of Open-Shell Spin-Coupling and Jahn-Teller Distortion in Benzene Radical Cation Probed by X-ray Spectroscopy

We report a theoretical investigation and elucidation of the X-ray absorption spectra of neutral benzene and of the benzene cation. The generation of the cation by multiphoton ultraviolet (UV) ionization and the measurement of the carbon K-edge spectra of both species using a table-top high-harmonic generation source are described in the companion experimental paper [Epshtein, M.; et al. J. Phys. Chem. A http://dx.doi.org/10.1021/acs.jpca.0c08736]. We show that the 1s_C → π transition serves as a sensitive signature of the transient cation formation, as it occurs outside of the spectral window of the parent neutral species. Moreover, the presence of the unpaired (spectator) electron in the π-subshell of the cation and the high symmetry of the system result in significant differences relative to neutral benzene in the spectral features associated with the 1s_C → π* transitions. High-level calculations using equation-of-motion coupled-cluster theory provide the interpretation of the experimental spectra and insight into the electronic structure of benzene and its cation. The prominent split structure of the 1s_C → π* band of the cation is attributed to the interplay between the coupling of the core → π* excitation with the unpaired electron in the π-subshell and the Jahn-Teller distortion. The calculations attribute most of the splitting (∼1-1.2 eV) to the spin coupling, which is visible already at the Franck-Condon structure, and we estimate the additional splitting due to structural relaxation to be around ∼0.1-0.2 eV. These results suggest that X-ray absorption with increased resolution might be able to disentangle electronic and structural aspects of the Jahn-Teller effect in the benzene cation.