Gauge origin independence in finite basis sets and perturbation theory

We show that origin independence in finite basis sets for the oscillator strengths is possibly in any gauge contrary to what is stated in literature. This is proved from a discussion of the consequences in perturbation theory when the exact eigenfunctions and eigenvalues to the zeroth order Hamiltonian H-0 cannot be found. We demonstrate that the erroneous conclusion for the lack of gauge origin independence in the length gauge stems from not transforming the magnetic terms in the multipole expansion leading to the use of a mixed gauge. Numerical examples of exact origin dependence are shown

Simulations of iron K pre-edge X-ray absorption spectra using the restricted active space method

The intensities and relative energies of metal K pre-edge features are sensitive to both geometric and electronic structures. With the possibility to collect high-resolution spectral data it is important to find theoretical methods that include all important spectral effects: ligand-field splitting, multiplet structures, 3d-4p orbital hybridization, and charge-transfer excitations. Here the restricted active space (RAS) method is used for the first time to calculate metal K pre-edge spectra of open-shell systems, and its performance is tested against on six iron complexes: [FeCl6](n-), [FeCl4](n-), and [Fe(CN)(6)](n-) in ferrous and ferric oxidation states. The method gives good descriptions of the spectral shapes for all six systems. The mean absolute deviation for the relative energies of different peaks is only 0.1 eV. For the two systems that lack centrosymmetry [FeCl4](2-/1-), the ratios between dipole and quadrupole intensity contributions are reproduced with an error of 10%, which leads to good descriptions of the integrated pre-edge intensities. To gain further chemical insight, the origins of the pre-edge features have been analyzed with a chemically intuitive molecular orbital picture that serves as a bridge between the spectra and the electronic structures. The pre-edges contain information about both ligand-field strengths and orbital covalencies, which can be understood by analyzing the RAS wavefunction. The RAS method can thus be used to predict and rationalize the effects of changes in both the oxidation state and ligand environment in a number of hard X-ray studies of small and medium-sized molecular systems.QC 20210215</p

Efficient calculations of a large number of highly excited states for multiconfigurational wavefunctions

Electronically excited states play important roles in many chemical reactions and spectroscopic techniques. In quantum chemistry, a common technique to solve excited states is the multiroot Davidson algorithm, but it is not designed for processes like X-ray spectroscopy that involves hundreds of highly excited states. We show how the use of a restricted active space wavefunction together with a projection operator to remove low-lying electronic states offers an efficient way to reach single and double-core-hole states. Additionally, several improvements to the stability and efficiency of the configuration interaction (CI) algorithm for a large number of states are suggested. When applied to a series of transition metal complexes the new CI algorithm does not only resolve divergence issues but also leads to typical reduction in computational time by 70%, with the largest savings for small molecules and large active spaces. Together, the projection operator and the improved CI algorithm now make it possible to simulate a wide range of single- and two-photon spectroscopies

Fingerprinting Electronic Structure of Heme Iron by Ab Initio Modeling of Metal L-Edge X-ray Absorption Spectra

The capability of the multiconfigurational restricted active space approach to identify electronic structure from spectral fingerprints is explored by applying it to iron L-edge X-ray absorption spectroscopy (XAS) of three heme systems that represent the limiting descriptions of iron in the Fe-O-2 bond, ferrous and ferric [Fe(P)(ImH)(2)](0/1+) (P = porphine, ImH = imidazole), and Fe-II(P). The level of agreement between experimental and simulated spectral shapes is calculated using the cosine similarity, which gives a quantitative and unbiased assignment. Further dimensions in fingerprinting are obtained from the L-edge branching ratio, the integrated absorption intensity, and the edge position. The results show how accurate ab initio simulations of metal L-edge XAS can complement calculations of relative energies to identify unknown species in chemical reactions