A Core-Valence Separated Similarity Transformed EOM-CCSD Method for Core-Excitation Spectra

Core-Excited States for Open-Shell Systems in Similarity-Transformed Equation-of-Motion Theory

X-ray absorption spectroscopy (XAS) is a powerful method for exploring molecular electronic structure by exciting core electrons into higher unoccupied molecular orbitals. In this study, we present the first integration of the spin-unrestricted similarity-transformed equation-of-motion coupled cluster method (CVS-USTEOM-CCSD) for core-excited and core-ionized states into the ORCA quantum chemistry package. Using the core-valence separation (CVS) approach, we evaluate the accuracy of CVS-USTEOM-CCSD across 13 open-shell organic systems, covering over 20 core excitations with diverse spin multiplicities (doublet, triplet, and quartet). The implementation leverages automated active space selection, incorporating CIS natural orbitals to efficiently capture electronic transitions. We benchmark the predicted K- and L-edge spectra against experimental data, underscoring the accuracy of the CVS-USTEOM-CCSD method for high-precision core excitation studies.

Elucidating Anomalous Intensity Ratios in Chlorine L-Edge X-ray Absorption Spectroscopy: Multiplet Effects and Core Rydberg Transitions

A relativistic core-valence-separated equation-of-motion coupled cluster (CVS-EOM-CC) study of chlorine L2,3-edge X-ray absorption near-edge structure (XANES) spectra using CH3Cl and CH2ICl as representative molecules is reported. The nearly identical intensity for the main features in the L2- and L3-edge XANES spectra is attributed to multiplet effects and the overlap between core-valence and core Rydberg transitions. The multiplet effects originating from the interaction between the core hole and the C-Cl σ* orbitals account for around half of the deviation of the L3 and L2 intensity ratio from the 2:1 ratio of the numbers of 2p3/2 and 2p1/2 electrons. The 2p3/2 → 4s core Rydberg transitions are shown to overlap with the 2p1/2 → σ* transitions and contribute to the other half of the intensity anomaly. We demonstrate that triple excitations in CVS-EOM-CC calculations play important roles in accurate simulation of the overlap between the 2p1/2 → σ* and 2p3/2 → 4s transitions.

Long-range Corrected Density Functional Theory Including a Two-Gaussian Hartree-Fock Operator for High Accuracy Core-excitation Energy Calculations of Both the Second- and Third-Row Atoms (LC2gau-core-BOP)

In the previous work, LCgau-core-BOP, which includes the short-range interelectronic Gaussian attenuating Hartree-Fock (HF) exchange to the long-range HF exchange, showed high accuracy core-excitation energies from 1s orbitals of the second-row atoms (1s → π*, 1s → σ*, 1s → n*, and 1s → Rydberg), but underestimates the core-excitation energies from 1s orbitals of the third-row atoms. To improve this, we added one more Gaussian attenuating HF exchange to LCgau-core-BOP. We named it LC2gau-core-BOP, which achieves a mean absolute error (MAE) of 0.6 and 0.3 eV for core excitation energies of the second- and third-row atoms of the tested small molecules, respectively. We found that the inclusion of the short-range interelectronic HF exchange at a distance ranging from 0.2 to 0.6 a.u. contributes to the increase of performances on 1s orbital energy calculations of the second-row atoms, while the inclusion of more short-range interelectronic HF exchange at a distance ranging from 0 to 0.2 a.u. does to the increase of performance on 1s orbital energy calculations of the third-row atoms. It is notable that all of these improvements were accomplished using flexible Gaussian attenuating HF exchange inclusion. LC2gau-core-BOP shows deviations of less than 0.8 eV from experimental values for all of the core-excitation energies of the tested medium-size molecules consisting of thymine, oxazole, glycine, and dibenzothiophene sulfone. Moreover, by optimizing one parameter of the OP correlation functional, LC2gau-core-BOP provides atomization energies over the G3 test set with an accuracy comparable to that of B3LYP.

Efficient State-Specific Natural Orbital Based Equation of Motion Coupled Cluster Method for Core-Ionization Energies: Theory, Implementation, and Benchmark

We have implemented a reduced-cost partial triples correction scheme to the equation of motion coupled cluster method for core-ionization energy based on state-specific natural orbitals. The second-order Algebraic Diagrammatic Construction (ADC) method is used to generate the state-specific natural orbital, which provides quicker convergence of the core-IP value with respect to the size of the virtual space than that observed in standard MP2-based natural orbitals. The error due to truncation of the virtual orbital can be reduced by using a perturbative correction. The accuracy of the method can be controlled by a single threshold, and there is a black box to use. The inclusion of the partial triples correction in the natural orbital based EOM-CCSD method greatly improves the agreement of the results with the experiment. The efficiency of the present implementation is demonstrated by calculating the core-ionization energy of a molecule containing 60 atoms and more than 2000 basis functions.

Ab initio investigation of the geometrical behavior in solution and electronic structure of the anion complexes [bis(1,3-dithiole-2-thione-4,5-dithiolate)M], for M = Bi(III), Sb(III), and Zn(II)

CONTEXT

1,3-Dithiole-2-thione-4,5-dithiolate (dmit) ligands are known for their conductive and optical properties. Dmit compounds have been assessed for use in sensor devices, information storage, spintronics, and optical material applications. Associations with various metallic centers endow dmit complexes with magnetic, optical, conductive, and antioxidant properties. Optical doping can facilitate the fabrication of magnetic conductor materials from ground-state nonmagnetic cations. While most studied complexes involve transition-metal centers due to their diverse chemistry, compounds with representative elements are less explored in the literature. This study investigated the structural and electronic properties of bisdmit complexes with representative Bi(III), Sb(III), and Zn(II) cations. AIMD calculations revealed two new geometries for Bi(III) and Zn(II) complexes, diverging from the isolated geometry typically used in quantum chemical calculations. The coordination of acetonitrile molecules to the cationic centers of the complexes resulted in unstable structures, while the dimerization of the complexes was stable. SA-CASSCF/NEVPT2 calculations were applied to the structures of the isolated complexes and stable dimers, confirming the multireference character of the electronic structure of the three systems and the multiconfigurational character of the Bi(III) complex. The electronic spectra simulated by the STEOM-DLPNO-CCSD calculations accurately reproduced the experimental UV‒Vis spectra indicating the participation of the isolated Bi(III) dmit complex and its dimeric form in solution.

METHODOLOGY

AIMD calculations of the dmit salts were conducted using the GFN2-xTB method with 60 explicit acetonitrile molecules as the solvent at 300 K for a total simulation time of 50.0 ps, with printing intervals of 0.5 fs. The final geometries were optimized employing the PBEh-3c compound method, incorporating implicit conductor-like polarizable continuum model (CPCM) solvation for acetonitrile. Local energy decomposition (LED) analysis at the DLPNO-CCSD(T)/Def2-TZVP level of theory was utilized to investigate the stability of the complex geometries identified by AIMD. The electronic structures of the complexes were assessed using the SA-CASSCF/NEVPT2/Def2-TZVP method to confirm the multiconfigurational and multireference nature of their electronic structures. Electronic spectra were analyzed using the STEOM-DLPNO-CCSD/Def2-TZVP method, with CPCM used to simulate an acetonitrile medium.

Frequency-Dependent Quadratic Response Properties and Two-Photon Absorption from Relativistic Equation-of-Motion Coupled Cluster Theory

We present the implementation of quadratic response theory based upon the relativistic equation-of-motion coupled cluster method. We showcase our implementation, whose generality allows us to consider both time-dependent and time-independent electric and magnetic perturbations, by considering the static and frequency-dependent hyperpolarizability of hydrogen halides (HX, X = F-At), providing comprehensive insights into their electronic response characteristics. Additionally, we evaluated the Verdet constant for noble gases Xe and Rn and discussed the relative importance of relativistic and electron correlation effects for these magneto-optical properties. Finally, we calculate the two-photon absorption cross sections of transition [ns1S0 → (n + 1)s1S0] of Ga+ and In+, which are suggested as candidates for new ion clocks. As our implementation allows for the use of nonrelativistic Hamiltonians as well, we have compared our EOM-QRCC results to the QR-CC implementation in the DALTON code and show that the differences between CC and EOMCC response are in general smaller than 5% for the properties considered. Collectively, the results underscore the versatility of our implementation and its potential as a benchmark tool for other approximated models, such as density functional theory for higher-order properties.

Fast and Accurate Calculation of the UV-Vis Spectrum with the Modified Local Excitation Approximation

The local excitation approximation (LEA), a method for the calculation of electronic excitations localized in a specific region of a molecule, has been modified with new approaches to enhance the accuracy of the original method. The primary concept behind LEA involves isolating the region of interest as a submolecule from the full molecule using a localization method, followed by calculating electronic excitations solely within this submolecule. In this study, we examined approaches that improve the accuracy in describing the region of interest, particularly its molecular orbital energies. Additionally, the localization method was extended with a new projection technique to accelerate calculations. These approaches were studied in time-dependent density functional theory (TDDFT) calculations applied to four testing systems with a chromophore as the region of interest: two basic linear molecules, acrolein surrounded by 24 water molecules, and a model of a green fluorescent protein. For all studied systems, the results of TDDFT calculations combined with LEA exhibited near-zero error when groups of atoms adjacent to the chromophore were explicitly included in the submolecule. This was achieved with at least a quadratic speedup of the calculation time as a function of the submolecule size.

Excited-State Dynamics during Primary C-I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to nπC═O* states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed nπC═O* and nσC-I* character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin-orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C-I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N4,5 edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C-I homolysis. By examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C-I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental-theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin-orbit coupled system.

Theoretical Investigation of the X-ray Stark Effect in Small Molecules

We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.

Double-Hybrid Density Functional Theory for Core Excitations: Theory and Benchmark Calculations

The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.

Pair Natural Orbital Equation-of-Motion Coupled-Cluster Method for Core Binding Energies: Theory, Implementation, and Benchmark

We present the theory and implementation of a lower scaling core-valence separated equation-of-motion coupled-cluster approach based on domain-based local pair natural orbitals for core binding energies. The accuracy of the new method has been compared with that of the standard equation-of-motion coupled-cluster method and experimentally measured results. The use of pair natural orbitals significantly reduces the computation cost and can be applied to large molecules.

Benchmark relativistic delta-coupled-cluster calculations of K-edge core-ionization energies of third-row elements

A benchmark computational study of K-edge core-ionization energies of third-row elements using relativistic delta-coupled-cluster (ΔCC) methods and a revised core-valence separation (CVS) scheme is reported. High-level relativistic (HLR) corrections beyond the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), including the contributions from two-electron picture-change effects, spin-orbit coupling, the Breit term, and quantum electrodynamics effects, have been taken into account and demonstrated to play an important role. Relativistic ΔCC calculations are shown to provide accurate results for core-ionization energies of third-row elements. The SFX2C-1e-CVS-ΔCC results augmented with HLR corrections show a maximum deviation of less than 0.5 eV with respect to experimental values.

Accurate Core-Excited States via Inclusion of Core Triple Excitations in Similarity-Transformed Equation-of-Motion Theory

The phenomenon of orbital relaxation upon excitation of core electrons is a major problem in the linear-response treatment of core-hole spectroscopies. Rather than addressing relaxation through direct dynamical correlation of the excited state via equation-of-motion coupled cluster theory (EOMEE-CC), we extend the alternative similarity-transformed equation-of-motion coupled cluster theory (STEOMEE-CC) by including the core-valence separation (CVS) and correlation of triple excitations only within the calculation of core ionization energies. This new method, CVS-STEOMEE-CCSD+cT, significantly improves on CVS-EOMEE-CCSD and unmodified CVS-STEOMEE-CCSD when compared to full CVS-EOM-CCSDT for K-edge core-excitation energies of a set of small molecules. The improvement in both absolute and relative (shifted) peak positions is nearly as good as that for transition-potential coupled cluster (TP-CC), which includes an explicit treatment of orbital relaxation, and CVS-EOMEE-CCSD*, which includes a perturbative treatment of triple excitations.