Search

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 1sC* 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, 1sC→π* and 1sC*→σ*(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 1sC→π* 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 1sC core orbital to the fully vacant π* orbitals. The 1sC*→σ* 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+.

We report high-level calculations of the excited states of [2,2]-paracyclophane (PCP), which was recently investigated experimentally by ultrafast pump–probe experiments on oriented single crystals [Haggag et al., ChemPhotoChem 6 e202200181 (2022)]. PCP, in which the orientation of the two benzene rings and their range of motion are constrained, serves as a model for studying benzene excimer formation. The character of the excimer state and the state responsible for the brightest transition are similar to those of the benzene dimer. The constrained structure of PCP allows one to focus on the most important degree of freedom, the inter-ring distance. The calculations explain the main features of the transient absorption spectral evolution. This brightest transition of the excimer is polarized along the inter-fragment axis. The absorption of the light polarized in the plane of the rings reveals the presence of other absorbing states of Rydberg character, with much weaker intensities. We also report new transient absorption data obtained by a broadband 8 fs pump, which time-resolve strong modulations of the excimer absorption. The combination of theory and experiment provides a detailed picture of the evolution of the electronic structure of the PCP excimer in the course of a single molecular vibration. [ABSTRACT FROM AUTHOR]

A roadmap is developed that integrates simulation methodology and data science methods to target new theories that traverse the multiple length- and time-scale features of many-body phenomena.

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 1B2u (ππ*) (S2) and 1B3u (nπ*) (S1) states, the participation of the optically dark 1Au (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 1Au (nπ*) and 1B3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The 1Au (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.

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 1sC → π 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 1sC → π* 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 1sC → π* 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.

Ultrafast table-top X-ray spectroscopy at the carbon K-edge is used to measure the X-ray spectral features of benzene radical cations (Bz+). The ground state of the cation is prepared selectively by two-photon ionization of neutral benzene, and the X-ray spectra are probed at early times after the ionization by transient absorption using X-rays produced by high harmonic generation (HHG). Bz+ is well-known to undergo Jahn-Teller distortion, leading to a lower symmetry and splitting of the π orbitals. Comparison of the X-ray absorption spectra of the neutral and the cation reveals a splitting of the two degenerate π* orbitals as well as an appearance of a new peak due to excitation to the partially occupied π-subshell. The π* orbital splitting of the cation, elucidated on the basis of high-level calculations in a companion theoretical paper [Vidal et al. J. Phys. Chem. A. http://dx.doi.org/10.1021/acs.jpca.0c08732], is discovered to be due to both the symmetry distortion and even more dominant spin coupling of the unpaired electron in the partially vacant π orbital (from ionization) with the unpaired electrons resulting from the transition from the 1sC core orbital to the fully vacant π* orbitals.

Application of photoionization mass spectroscopy, a technique capable of assessing protonation states in complex molecules in the gas phase, is challenging for arginine due to its fragility. We report photoionization efficiencies in the valence region of aqueous aerosol particles produced from arginine solutions under various pH and vaporization conditions. By using ab initio calculations, we investigate the stability of different conformers. Our results show that neutral arginine fragments upon ionization in the gas phase but solvation stabilizes the molecular ion, resulting in different photoionization dynamics. We also report the valence-band photoelectron spectra of the aerosol solutions obtained at different pH values.

Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention. [ABSTRACT FROM AUTHOR]

Coupled-cluster methods provide highly accurate models of molecular structure through explicit numerical calculation of tensors representing the correlation between electrons. These calculations are dominated by a sequence of tensor contractions, motivating the development of numerical libraries for such operations. While based on matrix–matrix multiplication, these libraries are specialized to exploit symmetries in the molecular structure and in electronic interactions, and thus reduce the size of the tensor representation and the complexity of contractions. The resulting algorithms are irregular and their parallelization has been previously achieved via the use of dynamic scheduling or specialized data decompositions. We introduce our efforts to extend the Libtensor framework to work in the distributed memory environment in a scalable and energy-efficient manner. We achieve up to 240× speedup compared with the optimized shared memory implementation of Libtensor. We attain scalability to hundreds of thousands of compute cores on three distributed-memory architectures (Cray XC30 and XC40, and IBM Blue Gene/Q), and on a heterogeneous GPU-CPU system (Cray XK7). As the bottlenecks shift from being compute-bound DGEMM's to communication-bound collectives as the size of the molecular system scales, we adopt two radically different parallelization approaches for handling load-imbalance, tasking and bulk synchronous models. Nevertheless, we preserve a unified interface to both programming models to maintain the productivity of computational quantum chemists.

Linear-scaling algorithms must be developed in order to extend the domain of applicability of electronic structure theory to molecules of any desired size. However, the increasing complexity of modern linear-scaling methods makes code development and maintenance a significant challenge. A major contributor to this difficulty is the lack of robust software abstractions for handling block-sparse tensor operations. We therefore report the development of a highly efficient symbolic block-sparse tensor library in order to provide access to high-level software constructs to treat such problems. Our implementation supports arbitrary multi-dimensional sparsity in all input and output tensors. We avoid cumbersome machine-generated code by implementing all functionality as a high-level symbolic C++ language library and demonstrate that our implementation attains very high performance for linear-scaling sparse tensor contractions.

This commentary documents how federal funding agencies are changing the criteria by which they distribute taxpayer money intended for scientific research. Increasingly, STEMM (Science, Technology, Engineering, Mathematics, and Medicine) funding agencies are requiring applicants for funding to include a plan to advance DEI (“Diversity, Equity, and Inclusion”) in their proposals and to dedicate a part of the research budget to its implementation. These mandates undermine the academic freedom of researchers and the unbiased generation of knowledge needed for a well-functioning democracy. Maintaining excellence in science is fundamental to the continuation of the U.S. as a global economic leader. Science provides a basis for solving important global challenges such as security, energy, climate, and health. Diverting funding from science into activities unrelated to the production of knowledge undermines science’s ability to serve humankind. When funding agencies politicize science by using their power to further a particular ideological agenda, they contribute to public mistrust in science. Hijacking science funding to promote DEI is thus a threat to our society. [ABSTRACT FROM AUTHOR]

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 1e1g and 3e2g 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. [ABSTRACT FROM AUTHOR]

The equation-of-motion coupled-cluster singles and doubles method with double electron attachment (EOM-DEA-CCSD) is capable of computing reliable energies, wave functions, and first-order properties of excited states in diradicals and polyenes that have a significant doubly excited character with respect to the ground state, without the need for including the computationally expensive triple excitations. Here, we extend the capabilities of the EOM-DEA-CCSD method to the calculations of a multiphoton property, two-photon absorption (2PA) cross sections. Closed-form expressions for the 2PA cross sections are derived within the expectation-value approach using response wave functions. We analyze the performance of this new implementation by comparing the EOM-DEA-CCSD energies and 2PA cross sections with those computed using the CC3 quadratic response theory approach. As benchmark systems, we consider transitions to the states with doubly excited character in twisted ethene and in polyenes, for which EOM-EE-CCSD (EOM-CCSD for excitation energies) performs poorly. The EOM-DEA-CCSD 2PA cross sections are comparable with the CC3 results for twisted ethene; however, the discrepancies between the two methods are large for hexatriene. The observed trends are explained by configurational analysis of the 2PA channels. [ABSTRACT FROM AUTHOR]

We present a new implementation for computing spin–orbit couplings (SOCs) within a time-dependent density-functional theory (TD-DFT) framework in the standard spin-conserving formulation as well in the spin–flip variant (SF-TD-DFT). This approach employs the Breit–Pauli Hamiltonian and Wigner–Eckart's theorem applied to the reduced one-particle transition density matrices, together with the spin–orbit mean-field treatment of the two-electron contributions. We use a state-interaction procedure and compute the SOC matrix elements using zero-order non-relativistic states. Benchmark calculations using several closed-shell organic molecules, diradicals, and a single-molecule magnet illustrate the efficiency of the SOC protocol. The results for organic molecules (described by standard TD-DFT) show that SOCs are insensitive to the choice of the functional or basis sets, as long as the states of the same characters are compared. In contrast, the SF-TD-DFT results for small diradicals (CH2, NH 2 + , SiH2, and PH 2 + ) show strong functional dependence. The spin-reversal energy barrier in a Fe(III) single-molecule magnet computed using non-collinear SF-TD-DFT (PBE0, ωPBEh/cc-pVDZ) agrees well with the experimental estimate. [ABSTRACT FROM AUTHOR]

Libtensor is a framework designed to implement the tensor contractions arising form the coupled cluster and equations of motion computational quantum chemistry equations. It has been optimized for symmetry and sparsity to be memory efficient. This allows it to run efficiently on the ubiquitous and cost-effective SMP architectures. Unfortunately, movement of memory controllers on chip has endowed these SMP systems with strong NUMA properties. Moreover, the manycore trend in processor architecture demands that the implementation be extremely thread-scalable on node. To date, Libtensor has been generally agnostic of these effects. To that end, in this paper, we explore a number of optimization techniques including a thread-friendly and NUMA-aware memory allocator and garbage collector, tuning the tensor tiling factor, and tuning the scheduling quanta. In the end, our optimizations can improve the performance of contractions implemented in Libtensor by up to 2× on representative Ivy Bridge, Nehalem, and Opteron SMPs.

We report a combined experimental and theoretical investigation of electron–molecule interactions using pyrrole as a model system. Experimental two-dimensional electron energy loss spectra (EELS) encode information about the vibrational states of the molecule as well as the position and structure of electronic resonances. The calculations using complex-valued extensions of equation-of-motion coupled-cluster theory (based on non-Hermitian quantum mechanics) facilitate the assignment of all major EELS features. We confirm the two previously described π resonances at about 2.5 and 3.5 eV (the calculations place these two states at 2.92 and 3.53 eV vertically and 2.63 and 3.27 eV adiabatically). The calculations also predict a low-lying resonance at 0.46 eV, which has a mixed character—of a dipole-bound state and σ* type. This resonance becomes stabilized at one quanta of the NH excitation, giving rise to the sharp feature at 0.9 eV in the corresponding EELS. Calculations of Franck–Condon factors explain the observed variations in the vibrational excitation patterns. The ability of theory to describe EELS provides a concrete illustration of the utility of non-Hermitian quantum chemistry, which extends such important concepts as potential energy surfaces and molecular orbitals to states embedded in the continuum. [ABSTRACT FROM AUTHOR]

A production-level implementation of equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) for electron attachment and excitation energies augmented by a complex absorbing potential (CAP) is presented. The new method enables the treatment of metastable states within the EOM-CC formalism in a similar manner as bound states. The numeric performance of the method and the sensitivity of resonance positions and lifetimes to the CAP parameters and the choice of one-electron basis set are investigated. A protocol for studying molecular shape resonances based on the use of standard basis sets and a universal criterion for choosing the CAP parameters are presented. Our results for a variety of π(*) shape resonances of small to medium-size molecules demonstrate that CAP-augmented EOM-CCSD is competitive relative to other theoretical approaches for the treatment of resonances and is often able to reproduce experimental results.

We present the two‐photon absorption (2PA) spectrum of aqueous thiocyanate calculated using high‐level quantum‐chemistry methods. The 2PA spectrum is compared to the one‐photon absorption (1PA) spectrum computed using the same computational protocol. Although the two spectra probe the same set of electronic states, the intensity patterns are different, leading to an apparent red‐shift of the 2PA spectrum relative to the 1PA spectrum. The presented analysis explains the intensity patterns and attributes the differences between the 1PA and 2PA spectra to the native symmetry of isolated SCN −, which influences the spectra in the low‐symmetry solvated environment. The native symmetry also manifests itself in variations of the polarization ratio (e.g., parallel vs. perpendicular cross sections) across the spectrum. The presented results highlight the potential of 2PA spectroscopy and high‐level quantum‐chemistry methods in studies of condensed‐phase phenomena. [ABSTRACT FROM AUTHOR]

The politicization of science – the infusion of ideology into the scientific enterprise – threatens the ability of science to serve humanity. Today, the greatest such threat comes from a set of ideological viewpoints collectively referred to as Critical Social Justice (CSJ). This contribution describes how CSJ has detrimentally affected scientific publishing by means of social engineering, censorship, and the suppression of scholarship. [ABSTRACT FROM AUTHOR]

We present an approach for obtaining a molecular orbital picture of the first dipole hyperpolarizability (β) from correlated many-body electronic structure methods. Ab initio calculations of β rely on quadratic response theory, which recasts the sum-over-all-states expression of β into a closed-form expression by calculating a handful of first- and second-order response states; for resonantly enhanced β, damped response theory is used. These response states are then used to construct second-order response reduced one-particle density matrices (1PDMs), which, upon visualization in terms of natural orbitals (NOs), facilitate a rigorous and black-box mapping of the underlying electronic structure with β. We explain the interpretation of different components of the response 1PDMs and the corresponding NOs within both the undamped and damped response theory framework. We illustrate the utility of this new tool by deconstructing β for cis-difluoroethene, para-nitroaniline, and hemibonded OH· + H2O complex, computed within the framework of coupled-cluster singles and doubles response theory, in terms of the underlying response 1PDMs and NOs for a range of frequencies. [ABSTRACT FROM AUTHOR]

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2p and 3p1h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections. [ABSTRACT FROM AUTHOR]

X-ray photon absorption leads to the creation of highly excited species, which often decay through the Auger process. The theoretical treatment of Auger decay is challenging because of the resonance nature of the initial core-excited or core-ionized states and the continuous nature of the ejected electron. In Paper I [W. Skomorowski and A. I. Krylov, J. Chem. Phys. 154, 084124 (2021)], we have introduced a theoretical framework for computing Auger rates based on the Feshbach–Fano approach and the equation-of-motion coupled-cluster ansätze augmented with core–valence separation. The outgoing Auger electron is described with a continuum orbital. We considered two approximate descriptions—a plane wave and a Coulomb wave with an effective charge. Here, we use the developed methodology to calculate Auger transition rates in core-ionized and core-excited benchmark systems (Ne, H2O, CH4, and CO2). Comparison with the available experimental spectra shows that the proposed computational scheme provides reliable ab initio predictions of the Auger spectra. The reliability, cost efficiency, and robust computational setup of this methodology offer advantages in applications to a large variety of systems. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]