# Probing Chirality across the Electromagnetic Spectrum with the Full Semi-classical Light-Matter Interaction

Martin van Horn, Trond Saue, Nanna List

► To cite this version:

Martin van Horn, Trond Saue, Nanna List. Probing Chirality across the Electromagnetic Spectrum with the Full Semi-classical Light-Matter Interaction. The Journal of Chemical Physics, 2022, 156 (5), pp.054113. 10.1063/5.0077502. hal-03414476

HAL Id: hal-03414476

https://hal.science/hal-03414476

Submitted on 20 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
The Journal of
Chemical Physics

RESEARCH ARTICLE | FEBRUARY 03 2022

# Probing chirality across the electromagnetic spectrum with the full semi-classical light-matter interaction

Martin van Horn ; Trond Saue ; Nanna Holmgaard List

<figure><img src="image_2.png" /></figure>

Check for updates

J. Chem. Phys. 156, 054113 (2022)

https://doi.org/10.1063/5.0077502

<figure><img src="image_3.png" /></figure>

View
Online

<figure><img src="image_4.png" /></figure>

Export
Citation

CrossMark

<figure><img src="image_5.png" /></figure>

AIP
Publishing

<figure><img src="image_6.png" /></figure>

The Journal of Chemical Physics

Special Topic: Adhesion and Friction

Submit Today!

<figure><img src="image_7.png" /></figure>

AIP
Publishing
# Probing chirality across the electromagnetic spectrum with the full semi-classical light-matter interaction

Cite as: J. Chem. Phys. 156, 054113 (2022); doi: 10.1063/5.0077502

Submitted: 3 November 2021 • Accepted: 10 January 2022 •

Published Online: 3 February 2022

<figure><img src="image_8.png" /></figure>

View Online

<figure><img src="image_9.png" /></figure>

Export Citation

<figure><img src="image_10.png" /></figure>

CrossMark

Martin van Horn,<fnref n="1" />,ᵃ⁾ Trond Saue,<fnref n="1" />,ᵇ⁾ and Nanna Holmgaard List<fnref n="2" />,ᶜ⁾

## AFFILIATIONS

## ABSTRACT

We present a formulation and implementation of anisotropic and isotropic electronic circular dichroism (ECD) using the *full* semi-classical light-matter interaction operator within a four-component relativistic framework. Our treatment uniquely accounts for *both* beyond-first-order light-matter interactions and relativistic effects, enabling us to investigate the ECD response across the electromagnetic spectrum from optical to x-ray wavelengths where relativistic selection rules and spatial field variations gain increasing importance. We consider the isotropic and oriented ECD across the valence transition and sulfur *L*- and *K*-edge transitions in the simplest disulfides, $H_2S_2$ and $(CH_3S)_2$, and evaluate the influence of the full interaction by comparing to a traditional truncated formulation in the Coulomb gauge (velocity representation). Additionally, we demonstrate that in the relativistic formalism, it is possible to work in the velocity representation, hence keeping order-by-order gauge-origin invariance, contrary to the multipolar gauge, yet being able to distinguish electric and magnetic multipole contributions. Going beyond a first-order treatment in the wave vector is mandatory in the higher-energy end of the soft x-ray region and beyond where the consequent intensity redistribution becomes significant. While the sulfur *K*-edge absorption spectrum is essentially unaffected by this redistribution, the signed differential counterpart is not: At least third-order contributions are required to describe the differential absorption profile that is otherwise overestimated by a factor of about two. The first-order description deteriorates at higher transition energies (beyond ~1000 eV) where it may even fail to predict the sign of individual differential oscillator strengths.

© 2022 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/5.0077502

## I. INTRODUCTION

Chirality is fundamental to life on Earth, and the origin of homochirality of natural amino acids and sugars remains an unresolved question.<fnref n="1" /> Enantiomers of a chiral molecule interact differently with chiral objects, such as chiral receptors or left and right circularly polarized light. The former is central in chiral molecular recognition, e.g., drug-receptor binding, while the latter enables spectroscopic detection of chirality. Electronic circular dichroism

(ECD) is the lowest-order natural chiroptical response: It measures the differential linear absorption between left and right circularly polarized light in the spectral region of electronic transitions.<fnref n="2" /> Many biomolecules are optically active either intrinsically due to the presence of a chiral center or induced through structural and electronic perturbations by a chiral environment such as through electronic coupling within a chiral arrangement of achiral chromophores.³,⁴ ECD in the valence region is therefore extensively used to provide stereochemical, conformational, and binding information of

<footnote n="1">¹Laboratoire de Chimie et Physique Quantiques, UMR 5626 CNRS—Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France</footnote>

<footnote n="2">²School of Engineering Sciences in Chemistry, Biotechnology and Health (CBH), Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Stockholm, Sweden</footnote>

<footnote>ᵃ⁾Electronic mail: mvanhorn@irsamc.ups-tlse.fr</footnote>

<footnote>ᵇ⁾Electronic mail: trond.saue@irsamc.ups-tlse.fr</footnote>

<footnote>ᶜ⁾Author to whom correspondence should be addressed: nalist@kth.se</footnote>
molecular systems.⁷ Anisotropic circular dichroism allows for additional information, which would otherwise be lost in the isotropic signal. As an example, anisotropic ECD was recently proposed as a way to distinguish different transitions in light-harvesting complex II oriented by its embedding in a membrane.⁶,⁷

While the discovery of the natural optical activity in the visible region dates back more than two hundred years to the work of Arago⁸ and Biot,⁹ its existence in the x-ray range was only established in 1998: The detection of x-ray natural circular dichroism (XNCD) in non-centrosymmetric crystals¹⁰⁻¹³ was enabled by the advent of intense third-generation synchrotron x-ray sources with full polarization control. Like other resonant core-level spectroscopies, XNCD exploits the site- and element-specificity of x-ray probes to provide information about the local symmetry around the absorbing atom(s). In addition to chiral crystals, XNCD has also been measured for small organic molecules including several amino acids (see, e.g., Refs. 14–17).

From a multipole expansion of the light-matter interactions in orders of the wave vector, the leading-order contributions to ECD (i.e., first order) are governed by the electric-dipole/magnetic-dipole (E1-M1) and electric-dipole/electric-quadrupole (E1-E2) interference terms. The latter is traceless and hence vanishes for isotropic samples. While the electric-dipole/magnetic-dipole interference dominates in the optical region (in the center-of-mass frame), it is difficult to observe at $K(L_1)$-edges¹⁸,¹⁹ because the magnetic transition-dipole moment involving the atomic-like $1(2)s$-orbitals becomes forbidden within a zeroth-order approximation in the non-relativistic domain.¹³,²⁰ The electric-dipole/electric-quadrupole interference term, on the other hand, gains intensity in the x-ray region due to its dependency on the spatial variations of the electromagnetic field.²⁰,²¹ For these reasons, substantial XNCD signals are expected only for oriented samples, although weak electric-dipole/magnetic-dipole contributions have been reported in the absence of orientational order (e.g., in solution or powdered samples¹⁹,²²).

Theoretical simulations play a key role in the assignment of experimental ECD spectra. Traditionally, these rely on the aforementioned multipole expansion of the ECD signal through first order in the wave vector. However, caution has to be exercised to ensure gauge-origin independent results, e.g., by performing calculations in the velocity representation²³ or by employing gauge-including atomic orbitals.²⁴,²⁵ An additional aspect arising in the x-ray regime, where the wavelength of the electromagnetic field approaches molecular dimensions, is the importance of higher-order terms in the expansion. To address these issues in the context of linear absorption of linearly polarized light, we recently proposed using the full semi-classical light-matter interaction operator in both the non-relativistic²⁶,²⁷ and relativistic regimes.²⁸ In addition to providing compatible transformation properties for both light and matter, the relativistic formulation is particularly simple because the light-matter interaction is linear in both the scalar and vector potentials.

In this work, we present a four-component relativistic formulation and implementation of isotropic and anisotropic ECD using the full semi-classical light-matter interaction. This allows us to investigate the ECD response across the electromagnetic spectrum, from optical to x-ray regimes. So far, previous theoretical studies of the isotropic XNCD of molecules have only considered

the first-order truncated interaction in a non-relativistic framework.¹⁸,²²,²⁹⁻³⁵ In these studies, the only source of non-vanishing response at the $L_1$- and $K$-edges is the polarization away from the atomic symmetry of the core orbitals. Our implementation accounts for two additional possible contributions: (i) effects of beyond first-order light-matter interactions and (ii) inclusion of relativistic effects, notably spin-orbit coupling that modifies selection rules (in particular, the magnetic transition-dipole selection rule³⁶,³⁷). Consequently, this allows, for the first time, to realistically examine the ECD response of molecules across the valence and x-ray regimes.

As test systems, we consider the simplest disulfide chromophore models, dihydrogen disulfide $H_2S_2$ and dimethyl disulfide $(CH_3S)_2$. Because of the low disulfide torsional barriers (~6–11 kcal/mol³⁸⁻⁴¹), the two enantiomeric forms (P- and M-helix) cannot be resolved experimentally. However, the disulfide bridge is an important structural element in proteins,⁴² where it preferentially occurs in non-planar, chiral conformations ($C_2$ symmetry) and hence displays structurally induced axial chirality. An interesting perspective for complex systems (e.g., proteins) is the potential use of XNCD as a local probe of chirality.³⁰,³⁵ This could potentially complement the delocalized conformational information encoded in valence ECD. Because of its computational tractability, $H_2S_2$ has been widely used to benchmark electronic structure methods for the calculation of chiroptical properties.⁴³⁻⁴⁹ For the same reason, Goulon et al. also used it to estimate the relative magnitudes of XNCD responses within the first-order truncated interaction and non-relativistic framework, reporting values below the experimental detection limits.¹⁸ Here, we revisit the ECD of the disulfide chromophore across the valence and $L$- and $K$-edges, going beyond these approximations.

## II. THEORY

### A. Full light-matter interaction

We begin with a description of our ECD formulation using the full semi-classical light-matter interaction, which has been implemented in a development version of the DIRAC electronic structure program.⁵⁰,⁵¹ We consider a circularly polarized monochromatic (angular frequency $\omega$) electromagnetic field with electric and magnetic components given by

$$ 
\begin{aligned}
\mathbf{E}^{(M)}(\mathbf{r}, t) &= \frac{1}{\sqrt{2}} E_{\omega} [\boldsymbol{\epsilon}_1 \sin(\mathbf{k} \cdot \mathbf{r} - \omega t) + (-1)^M \boldsymbol{\epsilon}_2 \cos(\mathbf{k} \cdot \mathbf{r} - \omega t)], \\
\mathbf{B}^{(M)}(\mathbf{r}, t) &= \frac{1}{\sqrt{2}} \frac{E_{\omega}}{c} [\boldsymbol{\epsilon}_2 \sin(\mathbf{k} \cdot \mathbf{r} - \omega t) - (-1)^M \boldsymbol{\epsilon}_1 \cos(\mathbf{k} \cdot \mathbf{r} - \omega t)],
\end{aligned}
\quad (1) $$

where appears the polarization vectors $\boldsymbol{\epsilon}_1$ and $\boldsymbol{\epsilon}_2$ that, together with the wave vector $\mathbf{k} = \boldsymbol{k}_\mathbf{k}$ ($k = \frac{\omega}{c}$), form a right-handed coordinate system; the field amplitude $E_\omega$; and the integer $M$ that specifies the handedness of the field, i.e., left (even) and right (odd) circularly polarized light, respectively, following the International Union of Pure and Applied Chemistry (IUPAC) recommendation.⁵² Such an electromagnetic wave is conventionally represented in the Coulomb gauge by the scalar and vector potentials,
$$ \phi^{(M)}(\mathbf{r}, t) = 0, $$

$$ \mathbf{A}^{(M)}(\mathbf{r}, t) = -\frac{1}{\sqrt{2}} \frac{E_{\omega}}{\omega} [\boldsymbol{\epsilon}_1 \cos(\mathbf{k} \cdot \mathbf{r} - \omega t) - (-1)^M \boldsymbol{\epsilon}_2 \sin(\mathbf{k} \cdot \mathbf{r} - \omega t)]. \quad (2) $$

In the relativistic domain, this leads to an interaction operator for left/right-handed (L/R) circularly polarized light of the form

$$ \hat{V}_{L/R}(t) = e c (\boldsymbol{\alpha} \cdot \mathbf{A}_{L/R}) = -\frac{1}{2} E_{\omega} \sum_{\pm \omega} \hat{T}_{L/R}(\omega) e^{-i\omega t}, \quad (3) $$

$$ \hat{T}_{L/R}(\omega) = \frac{e}{\omega} (c\boldsymbol{\alpha} \cdot \boldsymbol{\epsilon}_{L/R}) e^{i\mathbf{k} \cdot \mathbf{r}}, $$

where $\hat{T}_{L/R}(\omega)$ is the effective interaction operator (cf. Ref. 28), $\boldsymbol{\alpha}$ are Dirac matrices, and

$$ \boldsymbol{\epsilon}_{L/R} = \frac{1}{\sqrt{2}} (\boldsymbol{\epsilon}_1 \pm i \boldsymbol{\epsilon}_2) = \boldsymbol{\epsilon}_{R/L}^*. \quad (4) $$

For notational simplicity, we tacitly assume the summation over electrons in the effective interaction operator in Eq. (3). Considerable simplifications arise from the fact that the relativistic interaction operator is linear in the vector potential and there is no separation into spin and spatial parts.

The linear absorption cross section between initial and final states, $|i\rangle$ and $|f\rangle$, respectively, is then given by

$$ \sigma_{L/R}(\omega) = \frac{\pi \omega}{\epsilon_0 \hbar c} |\boldsymbol{\epsilon}_{L/R} \cdot \mathbf{T}|^2 f(\omega, \omega_{fi}, \gamma_{fi}), \quad \mathbf{T} = \langle f | \frac{e}{\omega} c \boldsymbol{\alpha} e^{i \mathbf{k} \cdot \mathbf{r}} | i \rangle. \quad (5) $$

The use of Fermi's golden rule for a transition between discrete states is justified by the inclusion of a (Lorentzian) line shape function $f(\omega, \omega_{fi}, \gamma_{fi})$ expressing the finite lifetime of excited states (see the discussion in Ref. 28). Moreover, the line shape function, having the dimension of time, has to be included in order for the absorption cross section to have the correct dimension of area. On the other hand, the dimensionless differential oscillator strength is obtained from the absorption cross section by the pre-factor substitution $\frac{\pi \omega}{\epsilon_0 \hbar c} \rightarrow \frac{2m_e \omega}{\hbar c^2}$ and removal of the line shape function. In the interest of compact expressions, we shall in the following carry out our derivations using oscillator strengths rather than absorption cross sections.

Evaluating the differential oscillator strengths between left and right circularly polarized light starting from Eq. (5) provides the expression for the anisotropic differential oscillator strength using the full semi-classical interaction operator,

$$ \begin{aligned} \Delta f = f_L - f_R &= \frac{2m_e \omega}{\hbar e^2} \{ (\boldsymbol{\epsilon}_L \cdot \mathbf{T})(\boldsymbol{\epsilon}_R \cdot \mathbf{T}^*) - (\boldsymbol{\epsilon}_L \cdot \mathbf{T}^*)(\boldsymbol{\epsilon}_R \cdot \mathbf{T}) \} \\ &= \frac{2m_e \omega}{\hbar e^2} (\boldsymbol{\epsilon}_L \times \boldsymbol{\epsilon}_R) \cdot (\mathbf{T} \times \mathbf{T}^*) \\ &= -i \frac{2m_e \omega}{\hbar e^2} \boldsymbol{\epsilon}_k \cdot (\mathbf{T} \times \mathbf{T}^*), \end{aligned} \quad (6) $$

where we have used Eq. (4) and the vector relation

$$ (\mathbf{a} \times \mathbf{b}) \cdot (\mathbf{c} \times \mathbf{d}) = (\mathbf{a} \cdot \mathbf{c})(\mathbf{b} \cdot \mathbf{d}) - (\mathbf{a} \cdot \mathbf{d})(\mathbf{b} \cdot \mathbf{c}). \quad (7) $$

This final, compact expression represents the relativistic extension of the expression previously reported by Hansen and Avery.⁵³

It is instructive to compare the differential oscillator strength in Eq. (6) (ECD) with the oscillator strength associated with linearly polarized light. The latter is obtained by replacing the complex polarization vector $\boldsymbol{\epsilon}_{L/R}$ in Eq. (5) by the polarization vector $\boldsymbol{\epsilon}$ giving the direction of the electric field vector in the plane perpendicular to the wave vector. One immediately notes that the differential oscillator strength contains no reference to polarization vectors at all, which can be understood as a manifestation of axial symmetry. This simplifies the rotational averaging appropriate for samples such as liquids and gases, with freely tumbling molecules. Specifically, the rotational average involves three angles $(\theta, \phi, \chi)$: The first two are associated with an average over propagation direction, whereas the polarization angle $\chi$ now becomes redundant.

## B. Truncated interaction

### 1. Multipoles in the velocity representation

In our previous work on the absorption of linearly polarized light, we considered two schemes involving truncated light-matter interaction.²⁸ In the first scheme, referred to as the generalized velocity representation, the truncated interaction is obtained via a Taylor expansion of the oscillator strength in orders of the wave vector. This is a relativistic extension of the theory developed by Bernadotte et al.⁵⁴ Applying this to Eq. (6), we obtain the following expressions for the differential oscillator strength to even and odd order in the wave vector:

$$ \Delta f^{[2n]} = \frac{2m_e \omega}{\hbar e^2} \boldsymbol{\epsilon}_k \cdot \sum_{m=0}^{n} (-1)^m (2 - \delta_{m0}) \\ \times \mathrm{Im} \left\{ (\mathbf{T}^{[n+m]} \times \mathbf{T}^{[n-m]*}) \right\} = 0, \quad (8) $$

$$ \Delta f^{[2n+1]} = \frac{2m_e \omega}{\hbar e^2} \boldsymbol{\epsilon}_k \cdot \sum_{m=0}^{n} (-1)^m 2 \mathrm{Re} \{ \mathbf{T}^{[n+m+1]} \times \mathbf{T}^{[n-m]*} \}, \quad (9) $$

where

$$ \mathbf{T}^{[n]} \cdot \boldsymbol{\epsilon}_{L/R} = \langle f | \hat{T}_{L/R}^{[n]} | i \rangle, \quad \hat{T}_{L/R}^{[n]} = \frac{e}{\omega} \frac{1}{n!} (i \boldsymbol{\alpha} \cdot \boldsymbol{\epsilon}_{L/R}) (\mathbf{k} \cdot \mathbf{r})^n. \quad (10) $$

Note that only odd-order terms are sensitive to the handedness and hence contribute to the expansion. This is a consequence of the time-reversal symmetry of the truncated interaction operators, making the corresponding transition moments real.²⁸ It should be stressed that the (differential) oscillator strength is gauge-origin invariant in each order, whereas the transition moments are intrinsically gauge-origin dependent.²⁸,⁵⁴

In the second scheme, referred to as generalized length representation, we transform to multipolar gauge (mg).²⁸,⁵⁵⁻⁵⁹ Using the Einstein summation convention, the full interaction operator is then expressed as

$$ \hat{V}_{\text{mg}}(\omega) = \sum_{n=0}^{\infty} \hat{V}_{\text{mg}}^{[n]}(\omega), \quad \left\{ \begin{array}{ll} \hat{V}_{\text{mg}}^{[0]}(\omega) & = -\hat{Q}_{p}^{[1]} E_{p}^{[0]}, \\ \hat{V}_{\text{mg}}^{[n \neq 0]}(\omega) & = -\frac{1}{(n+1)!} \hat{Q}_{j_1 \dots j_n p}^{[n+1]} E_{p j_1 \dots j_n}^{[n]} \\ & \qquad -\frac{1}{n!} \hat{m}_{j_1 \dots j_{n-1}; p} B_{r j_1 \dots j_{n-1}}^{[n-1]} \end{array} \right. \quad (11) $$
where appears electric and magnetic multipole operators,

$$ \hat{Q}_{j_1 \dots j_n p}^{[n+1]} = -e r_{j_1} r_{j_2} \dots r_{j_n} r_p, \quad (12) $$

$$ \hat{m}_{j_1 \dots j_{n-1}; r}^{[n]} = \frac{n}{n+1} r_{j_1} r_{j_2} \dots r_{j_{n-1}} (\mathbf{r} \times \hat{\mathbf{j}})_r, \quad \hat{\mathbf{j}} = -e\mathbf{c}\alpha, \quad (13) $$

which couple to the electric and magnetic fields and their derivatives,

$$ F_{p; j_1 \dots j_n}^{[n]} = \left. \frac{\partial^n F_p}{\partial r_{j_1} \partial r_{j_2} \dots \partial r_{j_n}} \right|_{\mathbf{r}=0}, \quad \mathbf{F} = \mathbf{E}, \mathbf{B}, \quad (14) $$

at a selected expansion point, here set to the origin. The corresponding nth-order effective interaction operator is given by

$$ \hat{T}_{\text{mg;L/R}}^{[n]}(\omega) = -e \left[ \frac{1}{(n+1)!} (\mathbf{r} \cdot \boldsymbol{\epsilon}_{\text{L/R}})(\mathbf{i}\mathbf{k} \cdot \mathbf{r})^n - \frac{i}{\omega} \frac{n}{(n+1)!} \times (\mathbf{i}\mathbf{k} \times \boldsymbol{\epsilon}_{\text{L/R}}) \cdot (\mathbf{r} \times \boldsymbol{c}\boldsymbol{\alpha}) (\mathbf{i}\mathbf{k} \cdot \mathbf{r})^{n-1} \right]. \quad (15) $$

In the component form, the effective interaction operator is conveniently expressed as

$$ \hat{T}_{\text{mg;L/R}}^{[n]}(\omega) = i^n \epsilon_{\text{L/R}} k_{j_1} k_{j_2} \dots k_{j_n} \hat{X}_{j_1 \dots j_{n-1}; r}^{[n]}(\omega), \quad (16) $$

in terms of a multipole operator,²⁸

$$ \hat{X}_{j_1 \dots j_n; p}^{[n]}(\omega) = \frac{1}{(n+1)!} \hat{Q}_{j_1 \dots j_n; p}^{[n+1]} - \frac{i}{\omega} \frac{1}{n!} \hat{m}_{j_1 \dots j_{n-1}; r}^{[n]} \epsilon_{r_{j_n} p}(\omega), \quad (17) $$

($\epsilon_{r_{j_n} p}$ is the Levi-Civita symbol). However, in contrast to the generalized velocity representation, we now lose the origin independence of oscillator strengths.²⁸

It is possible, though, to combine the best of both worlds. In our previous work,²⁸ we have shown that it is possible within the generalized velocity representation to distinguish electric and magnetic multipole components, while maintaining origin independence of the oscillator strengths, although its individual components still exhibit this dependency. The operator manipulations required for this separation do not involve commutators with the Hamiltonian and are therefore also valid when using finite basis sets. We here review and expand this result. Using the vector relation in Eq. (7), the effective interaction operator of Eq. (10) can be rewritten as

$$ \hat{T}_{\text{L/R}}^{[n]} = \frac{e}{\omega} \left[ \frac{1}{(n+1)!} \left\{ (ic\boldsymbol{\alpha} \cdot \boldsymbol{\epsilon}_{\text{L/R}})(\mathbf{k} \cdot \mathbf{r})^n + n(\mathbf{k} \cdot ic\boldsymbol{\alpha}) \times (\boldsymbol{\epsilon}_{\text{L/R}} \cdot \mathbf{r})(\mathbf{k} \cdot \mathbf{r})^{n-1} \right\} + \frac{n}{(n+1)!} (\mathbf{k} \times \boldsymbol{\epsilon}_{\text{L/R}}) \cdot (\mathbf{r} \times ic\boldsymbol{\alpha})(\mathbf{k} \cdot \mathbf{r})^{n-1} \right]. \quad (18) $$

As evident upon comparison to Eq. (15), we have effectively recovered the magnetic multipole part of the operator, and, in fact,

writing the corresponding nth-order interaction operator $\hat{V}_{\text{L/R}}^{[n]}$ out in the component form, we recover the expression of Eq. (11), with the important modification that the electric multipole operator in the length representation Eq. (12) is replaced by its counterpart in the velocity representation

$$ \hat{Q}_{j_1 \dots j_n p}^{[n+1]} = \frac{ie}{\omega} r_{j_1} \dots r_{j_{n-1}} ((c\alpha_p)r_{j_n} + n(c\alpha_{j_n})r_p). \quad (19) $$

This somewhat curious form is perhaps best understood by noting that the operator can equally well be expressed as

$$ Q_{j_1 \dots j_n p}^{[n+1]} = \frac{ie}{\omega} \underbrace{(c \alpha_{j_1} r_{j_2} \dots r_{j_n} r_p + \dots + r_{j_1} r_{j_2} \dots c \alpha_{j_n} r_p + r_{j_1} r_{j_2} \dots r_{j_n} c \alpha_p)}_{(n+1) \text{ terms}} \quad (20) $$

being perfectly symmetric in all indices {$j_1 \dots j_n p$}, just as Eq. (12). The more compact expression in Eq. (19) is possible because the relativistic velocity operator $c\alpha$ commutes with the coordinates, contrary to the corresponding non-relativistic operator $\mathbf{p}/m$. Finally, using the commutator relation with the Hamiltonian

$$ -\frac{i}{\hbar} \left[ (\mathbf{r} \cdot \boldsymbol{\epsilon}_{\text{L/R}})(\mathbf{k} \cdot \mathbf{r})^n, \hat{\mathbf{h}} \right] = (c\alpha \cdot \boldsymbol{\epsilon}_{\text{L/R}})(\mathbf{k} \cdot \mathbf{r})^n + n(\mathbf{r} \cdot \boldsymbol{\epsilon}_{\text{L/R}}) \times (\mathbf{k} \cdot c\alpha)(\mathbf{k} \cdot \mathbf{r})^{n-1}, \quad (21) $$

one can show that transition moments of the electric multipole operator in the length and velocity representation will be identical.²⁸ However, very importantly, this generally only holds in a complete basis.

The odd-order oscillator strength of Eq. (9) may now be expressed as

$$ \Delta f^{[2n+1]} = \frac{2m_e \omega}{\hbar e^2} k^{2n+1} e_{k; j_1} \dots e_{k; j_{2n+1}} \mathbf{e}_k \cdot \sum_{m=0}^{n} (-1)^m 2 \left( \mathbf{X}_{j_1 \dots j_{n+m+1}}^{[n+m+1]} \times \mathbf{X}_{j_{n+m+2} \dots j_{2n+1}}^{[n-m]} \right), \quad (22) $$

where the components of the real vector $\mathbf{X}_{j_1 \dots j_n}^{[n]}$ are given as the transition moments of the multipole operator in Eq. (17) in the velocity representation, that is, using Eq. (19) for the electric multipole operator.

## 2. Angular dependence of truncated ECD

In the anisotropic case, both the full and truncated differential oscillator strengths [Eqs. (6) and (22), respectively] are functions of the orientation of the incident radiation. In the case of truncated interaction, we make this angular dependence explicit by writing the unit wave vector as

$$ \mathbf{e}_k = (\sin \theta \cos \phi, \sin \theta \sin \phi, \cos \theta). \quad (23) $$

The first-order differential oscillator strength can then be expanded as
$$
\begin{align}
\Delta f^{[1]}(\theta, \phi) &= e_{k;\ell} R_{ij}^{[1]} e_{k;\ell} \nonumber \\
&= R_{\text{iso}}^{[1]} + \frac{1}{2} (R_{zz}^{[1]} - R_{\text{iso}}^{[1]}) (3\cos^2\theta - 1) \nonumber \\
&\quad + \frac{1}{2} (R_{xz}^{[1]} + R_{zx}^{[1]}) \sin 2\theta \cos \phi \nonumber \\
&\quad + \frac{1}{2} (R_{yz}^{[1]} + R_{zy}^{[1]}) \sin 2\theta \sin \phi \nonumber \\
&\quad + \frac{1}{2} (R_{xx}^{[1]} - R_{yy}^{[1]}) \sin^2 \theta \cos 2\phi \nonumber \\
&\quad + \frac{1}{2} (R_{xy}^{[1]} + R_{yx}^{[1]}) \sin^2 \theta \sin 2\phi, \tag{24}
\end{align}
$$

where we have introduced the rotational strength tensor $R_{ij}^{[1]}$ (cf. Ref. 60) and its isotropic part $R_{\text{iso}}^{[1]}$, known as the rotational strength,

$$R_{ij}^{[1]} = \frac{4m_e \omega^2}{\hbar c e^2} (\mathbf{X}_j^{[1]} \times \mathbf{X}^{[0]})_i, \quad R_{\text{iso}}^{[1]} = \frac{1}{3} (R_{xx}^{[1]} + R_{yy}^{[1]} + R_{zz}^{[1]}). \quad (25)$$

From Eq. (24), we see that the anisotropic part of $\Delta f^{[1]}(\theta, \phi)$ is a linear combination of d-orbitals weighted by the relevant elements of the rotatory strength tensor. This can be made more explicit by rewriting the first-order differential oscillator strength as

$$
\begin{align}
\Delta f^{[1]}(\theta, \phi) ={}& \sqrt{15} R_{\text{iso}}^{[1]} s + \sqrt{3} (R_{zz}^{[1]} - R_{\text{iso}}^{[1]}) d_{z2} \nonumber \\
& + (R_{xx}^{[1]} - R_{yy}^{[1]}) d_{x^2-y^2} + (R_{xz}^{[1]} + R_{zx}^{[1]}) d_{xz} \nonumber \\
& + (R_{yz}^{[1]} + R_{zy}^{[1]}) d_{yz} + (R_{xy}^{[1]} + R_{yx}^{[1]}) d_{xy}, \tag{26}
\end{align}
$$

where we have used a common normalization for all solid harmonics.

From Eq. (22), we see that rotational strength tensors can be generalized to arbitrary odd orders. It will also be seen that the differential oscillator strength $\Delta f^{[2n+1]}$ contains products of $2(n+1)$ components of the unit wave vector and is therefore spanned by spherical harmonics of even $\ell = 0, 2, \dots, 2(n+1)$.

In the *supplementary material*, the (first-order) rotational strength tensor is analyzed in terms of symmetry for the $C_2$ point group of the disulfide configurations considered, assuming the rotation axis to coincide with the z-axis. We also show that the rotational strength tensor in the generalized velocity representation is invariant under a shift of origin, also in a finite basis, which generally does not hold true in the generalized length representation. As pointed out above, the separation into electric and magnetic multipole contributions generally depends on the choice of gauge origin and, thus, is not unique. However, the presence of symmetry may introduce additional gauge-origin invariances.⁶¹ The E1–M1 and E1–E2 contributions to the $R_{z*}^{[1]}$ elements are separately invariant under origin shifts along the rotation axes (Sec. S3); in A symmetry, this holds for all elements of the *full* rotational strength tensor.

In the case of the full interaction, the factorization of radial and angular parts is possible using the plane wave expansion

$$e^{\mathbf{k} \cdot \mathbf{r}} = 4\pi \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{\ell} i^{\ell} f_{\ell}(kr) Y_{\ell}^{m}(\mathbf{e}_{k}) Y_{\ell}^{m*}(\mathbf{e}_{r}) \quad (27)$$

in terms of spherical harmonics $Y_{\ell}^{m}$ and spherical Bessel functions $j_{\ell}$. However, it is more cumbersome since it involves infinite sums over angular momentum $\ell$ and will not be pursued further here.

### C. Rotational average

As pointed out above, rotational averaging is simplified for the differential oscillator strength because of the absence of reference to the polarization vectors. It is therefore limited to the solid angle indicating the propagation direction and is defined as

$$\langle g(\mathbf{r}) \rangle_{\theta, \phi} = \frac{1}{4\pi} \int_{0}^{2\pi} \int_{0}^{\pi} g(\mathbf{r}) \sin \theta d\theta d\phi. \quad (28)$$

In the case of the full light–matter interaction, the rotational average is handled numerically using a Lebedev grid.²⁷,²⁸ In the case of the truncated interaction, the rotational average is expressed as

$$
\begin{equation}
\begin{split}
\langle \Delta f^{[2n+1]} \rangle_{\theta, \phi} ={}& \Delta f_{\text{iso}}^{[2n+1]} = \frac{2m_e \omega}{\hbar c e^2} k^{2n+1} \langle e_{k; j_1} \dots e_{k; j_{2n+1}} \mathbf{e}_k \rangle_{\theta, \phi} \\
& \cdot \sum_{m=0}^{n} (-1)^m \binom{2n}{2} \left\{ \mathbf{X}_{j_1 \dots j_{n+m+1}}^{[n+m+1]} \times \mathbf{X}_{j_{n+m+2} \dots j_{2n}}^{[n-m]} \right\}.
\end{split}
\tag{29}
\end{equation}
$$

In the lowest order, the isotropic differential oscillator strength becomes

$$
\begin{align}
\langle \Delta f^{[1]} \rangle_{\theta, \phi} ={}& \frac{4m_e \omega^2}{3\hbar c e^2} \sum_j (\mathbf{x}_j^{[1]} \times \mathbf{x}^{[0]})_j \nonumber \\
& - \frac{8m_e \omega^2}{3\hbar c e^2} \sum_j \langle f | \hat{\mathcal{O}}_j^{[1]} | i \rangle \langle f | -\frac{i}{\omega} \hat{m}_j^{[1]} | i \rangle = R_{\text{iso}}^{[1]}, \tag{30}
\end{align}
$$

where we have used that $\langle e_{k;\ell} \mathbf{e}_k \rangle_{\theta,\phi} = \frac{1}{c} \mathbf{e}_j$ and that the E1–E2 contribution vanishes since its elements form a traceless matrix.

### III. IMPLEMENTATION DETAILS

In practice, the interaction operator is separated into Hermitian and anti-Hermitian parts for compatibility with the quaternion symmetry scheme of DIRAC (see Refs. 28 and 62). The implementation of Eq. (6) requires the same integrals as the linear oscillator strength, and thus, our previous integral implementation could be readily extended to ECD.²⁶,²⁷ The integral evaluation over Lebedev points has been parallelized using MPI (Message Passing Interface).

### IV. COMPUTATIONAL DETAILS

The geometries of $H_2S_2$ and $(CH_3S)_2$ were obtained using the B3LYP⁶³⁻⁶⁶ exchange-correlation functional and the cc-pVTZ⁶⁷,⁶⁸ basis set. Geometry optimizations were performed in Gaussian 16.⁶⁹ To mimic the $\chi_3$ disulfide angle typical for protein structures,⁷⁰ we performed constrained geometry optimizations for $\chi_3 = -87^\circ$, corresponding to M-helical chirality.⁷¹ The restricted excitation window approach⁷²,⁷³ was used to selectively target the sulfur L- and K-edges. This also eliminates the issue of artificial transitions to quasi-continuum orbitals caused by finite basis set effects that otherwise often interfere with simulations at the L-edge.⁷⁴,⁷⁵ A Gaussian model was employed for the nuclear charge distribution,⁷⁶ and an 86-point Lebedev grid ($L_{\text{max}} = 12$) was used for the isotropic
averaging of the differential linear absorption based on the full inter-
action operator. The gauge origin was placed in the center-of-mass
(COM), and spatial symmetry was invoked in all cases except for the
gauge-origin dependence calculations.

Excitation energies, linear and differential absorption cross
sections for the full interaction operator, and the multipole expan-
sions within the generalized velocity gauge were computed using
the PBE0<sup>77,78</sup> exchange-correlation functional and the uncontracted
aug-pcX-3<sup>79</sup> and aug-pc-3<sup>80-82</sup> basis sets for sulfur and hydrogen,
respectively. The pcX-n basis set series was developed for describ-
ing core-excitation processes using the ΔSCF (Self-Consistent Field)
approach at both the nonrelativistic and relativistic levels. The
small component basis sets were generated within the condition
of restricted kinetic balance. The relativistic calculations were per-
formed using a Dirac-Coulomb Hamiltonian in which the (SS|SS)
integrals are replaced by an interatomic SS energy correction.⁸³
The gauge-origin invariance of our implementation of the full
semi-classical formulation of the isotropic and anisotropic rota-
tory strengths [Eq. (6)] and its first-order truncated counterpart
was confirmed numerically by shifting the gauge-origin (from 0
to 100 a₀) along the C₂ axis. This leads to a redistribution of the

E1–E2 and E1–M1 contributions to $R_{xx}^{[1]}$ and $R_{yy}^{[1]}$ for transitions of *B* symmetry. As expected, the results remained unchanged for both the full and truncated formulations (data not shown). Simulated spectra were obtained by convolving the stick spectrum with Gaussian line shape functions with a full width at half maximum (FWHM) of 0.4 eV, and those for $(CH_3S)_2$ were shifted by different offsets for each absorption edge to match their experimental counterparts.

V. RESULTS AND DISCUSSION

Before considering the ECD response of the disulfides, we
assign the linear absorption features across the valence and x-ray
regions. We initially focus on (CH3S)2 for which experimental gas-
phase absorption spectra are available.⁸⁴⁻⁸⁸ Expectedly, and as shown
in Fig. S1, the spectral profiles for H2S2 are similar, and because
of its greater computational tractability, we consider this minimal
disulfide in subsequent analyses.

Figure 1 displays the rotationally averaged linear and differ-
ential absorption spectra for valence transition and sulfur L- and
K-edge transitions of (CH3S)2, computed using the full interaction

<figure><img src="image_11.png" /><figcaption>FIG. 1. Isotropic linear (top) and differential (bottom) absorption spectra of (CH₃S)₂: (a) valence, (b) L₂,₃-edge, (c) L₁-edge, and (d) K-edge spectra using the full interaction operator in Eq. (3) (green shadings) or the lowest non-vanishing generalized velocity representation (orange lines). The left axes correspond to (differential) absorption cross sections, whereas (differential) oscillator strengths (sticks) are shown on the right axes. The black sticks indicate the location of all computed transitions, whereas the black lines are experimental spectra.⁸⁴⁻⁸⁸ The stick spectra were convolved with a Gaussian line shape with a FWHM of 0.4 eV. The theoretical absorption spectra have been uniformly shifted to align with the experimental counterparts (shift values indicated in the top panels). The same shifts were applied to the ECD spectra.</figcaption></figure>
operator (green shading) with corresponding oscillator strengths indicated as green sticks. Hereafter, we explicitly indicate the results of the full interaction with the superscript “full.” For comparison, we also provide the lowest non-vanishing terms in the truncated generalized velocity representation (orange lines) for linear and differential absorption, i.e., zeroth- and first-order in the magnitude of the wave vector, respectively. The black sticks at the top indicate the location of the underlying electronic transitions, and black lines (solid and dashed) indicate the experimental absorption spectra.⁸⁴⁻⁸⁸ Apart from uniform shifts necessary to align the lowest-energy band to the respective experimental spectrum, the theoretical spectra capture well both relative intensities and peak splittings.

The first valence band [Fig. 1(a)] is dominated by the two excitations from the symmetric and antisymmetric combinations of the non-bonding 3p orbital on each sulfur to the lowest unoccupied $\sigma_{SS}^*$ orbital (b symmetry). This assignment is consistent with the analysis by Lindberg and Michl on H₂S₂.⁸⁹ Since the valence orbitals are found to have well-defined spin, we adopt a non-relativistic state notation for the valence states, i.e., 1¹B and 2¹A, respectively. They are separated by ~ 0.15 eV. The second valence band originates from transitions into the $\sigma_{CS}^*$ orbital (3¹A and 2¹B).

Turning to the x-ray region, the first two bands at the sulfur L₂,₃-edge [Fig. 1(b)] are dominated by transitions from the symmetric and antisymmetric combinations of the S 2p₃/₂ and 2p₁/₂ core orbitals to the $\sigma_{SS}^*$ orbital. We note that the A/B pair of associated transitions are essentially degenerate (less than 0.1 meV splitting) because of the limited overlap between the core orbitals on each of the sulfur atoms. The calculated spin-orbit splitting of ~1.35 eV between the L₃- and L₂-branches is comparable to the experimental splitting reported for dihydrogen sulfide (~ 1.2 eV)⁹⁰,⁹¹. The L₃/L₂ branching ratio of ~1.4 : 1 (obtained by summing the underlying oscillator strengths) deviates significantly from its statistical value, which is obtained only in the limit of j-j coupling.⁹²,⁹³ The second peak in each branch is dominated by excitations to the $\sigma_{CS}^*$ orbital and is separated from the first peak by ~1.4–1.5 eV. Consequently, the second band in the spectrum contains contributions from both branches, whereas the third band is associated with the second peak in the L₂ branch. These assignments agree with previous studies.⁸⁶,⁹⁴ The energy range considered as well as the basis set used in our calculations does not cover the fourth band in the experimental spectrum that, according to previous work,⁸⁶ originates from excitations to higher-lying orbitals of mixed

<table><caption>**TABLE I.** Beyond first-order effects in the isotropically averaged and anisotropic rotatory strengths for the lowest intense valence transition and L- and K-edge transitions of (CH₃S)₂ for the full semi-classical light-matter interaction operator and in first-order within the Coulomb gauge, computed at the 4c-TD-PBE0 level of theory and the uncontracted aug-pcx-3/aug-pc3 basis set. For the oriented case, $\theta = \phi = \pi/2$ for the valence and L₂,₃-edge transitions and $\theta = \pi/2, \phi = \pi/4$ for the L₁- and K-edge transitions [cf. Eq. (23)]. For each spectral region, the accumulated (differential) oscillator strength (ac) across the two pairs of transitions (A/B symmetry) is given in the third line. The numbers in parentheses are exponents of 10.</caption><thead><tr><th>Transition</th><th>ΔE (eV)</th><th>f<sup>full</sup><sub>iso</sub></th><th>f<sup>[0]</sup><sub>iso</sub></th><th>Δf<sup>full</sup><sub>iso</sub></th><th>Δf<sup>[1]</sup><sub>iso</sub></th><th>Δf<sup>full</sup><sub>k</sub></th><th>Δf<sup>[1]</sup><sub>k</sub></th></tr></thead><tbody><tr><td colspan="8" style="text-align: center; font-weight: bold;">Valence</td></tr><tr><td>1<sup>1</sup>B (n<sub>S</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>4.695 14</td><td>7.7977(-03)</td><td>7.7974(-03)</td><td>-1.0856(-04)</td><td>-1.0856(-04)</td><td>-6.1154(-05)</td><td>-6.1147(-05)</td></tr><tr><td>2<sup>1</sup>A (n<sub>S</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>4.770 69</td><td>2.0797(-03)</td><td>2.0793(-03)</td><td>8.1247(-05)</td><td>8.1244(-05)</td><td>1.1759(-04)</td><td>1.1758(-04)</td></tr><tr><td></td><td></td><td>ac</td><td>9.8773(-03)</td><td>9.8767(-03)</td><td>-2.7316(-05)</td><td>-2.7317(-05)</td><td>5.6436(-05)</td><td>5.6435(-05)</td></tr><tr><td colspan="8" style="text-align: center; font-weight: bold;">L<sub>2,3</sub>-edge</td></tr><tr><td>6B (2p<sub>3/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>158.087 21</td><td>3.1232(-03)</td><td>3.1230(-03)</td><td>5.0231(-04)</td><td>5.0326(-04)</td><td>-1.1661(-03)</td><td>-1.169 2(-03)</td></tr><tr><td>6A (2p<sub>3/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>158.087 23</td><td>2.6587(-03)</td><td>2.6566(-03)</td><td>-4.9798(-04)</td><td>-4.9892(-04)</td><td>1.1726(-03)</td><td>1.1757(-03)</td></tr><tr><td></td><td></td><td>ac</td><td>5.7818(-03)</td><td>5.7796(-03)</td><td>4.3272(-06)</td><td>4.3363(-06)</td><td>6.4960(-06)</td><td>6.5232(-06)</td></tr><tr><td>10A (2p<sub>1/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>159.293 75</td><td>1.2239(-03)</td><td>1.2227(-03)</td><td>-2.5299(-04)</td><td>-2.5345(-04)</td><td>-5.7854(-04)</td><td>-5.7988(-04)</td></tr><tr><td>10B (2p<sub>1/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>159.293 87</td><td>1.3957(-03)</td><td>1.3953(-03)</td><td>2.5605(-04)</td><td>2.5651(-04)</td><td>5.8340(-04)</td><td>5.8476(-04)</td></tr><tr><td></td><td></td><td>ac</td><td>2.6196(-03)</td><td>2.6180(-03)</td><td>3.0568(-06)</td><td>3.0622(-06)</td><td>4.8616(-06)</td><td>4.8783(-06)</td></tr><tr><td colspan="8" style="text-align: center; font-weight: bold;">L<sub>1</sub>-edge</td></tr><tr><td>4B (2s<sub>1/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>216.393 99</td><td>3.6733(-02)</td><td>3.6872(-02)</td><td>-3.1538(-05)</td><td>-3.1716(-05)</td><td>-1.0027(-03)</td><td>-1.008 5(-03)</td></tr><tr><td>4A (2s<sub>1/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>216.395 99</td><td>2.9891(-04)</td><td>1.7881(-04)</td><td>3.1469(-05)</td><td>3.1646(-05)</td><td>1.0279(-03)</td><td>1.033 8(-03)</td></tr><tr><td></td><td></td><td>ac</td><td>3.7032(-02)</td><td>3.7051(-02)</td><td>-6.9310(-08)</td><td>-6.9637(-08)</td><td>2.5196(-05)</td><td>2.521 2(-05)</td></tr><tr><td colspan="8" style="text-align: center; font-weight: bold;">K-edge</td></tr><tr><td>4B (1s<sub>1/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>2427.833 34</td><td>7.7555(-03)</td><td>1.0438(-02)</td><td>-2.2670(-05)</td><td>-4.6773(-05)</td><td>-1.0771(-03)</td><td>-1.984 2(-03)</td></tr><tr><td>4A (1s<sub>1/2</sub> → σ<sup>*</sup><sub>SS</sub>)</td><td>2427.833 49</td><td>2.6746(-03)</td><td>2.4290(-05)</td><td>2.2665(-05)</td><td>4.6761(-05)</td><td>1.0754(-03)</td><td>1.982 3(-03)</td></tr><tr><td></td><td></td><td>ac</td><td>1.0430(-02)</td><td>1.0462(-02)</td><td>-5.6000(-09)</td><td>-1.1871(-08)</td><td>-1.6290(-06)</td><td>-1.821 6(-06)</td></tr></tbody></table>
$\sigma_{\text{CS}}^{*}$/Rydberg character. Not surprisingly, the L₁- and K-edge spectra bear strong resemblance [Figs. 1(c) and 1(d)]: They display two pre-edge features, separated by ~1.5 eV, which originate from pairs of near-degenerate excitations from the bonding and antibonding combinations of sulfur s-orbitals into the $\sigma_{\text{SS}}^{*}$ and $\sigma_{\text{CS}}^{*}$ orbitals, respectively.$^{87,95,96}$

A non-vanishing ECD response in these minimal disulfides results from axial chirality caused by trapping the disulfide bridge in a non-planar (i.e., C₂) conformation. As described above, the two lowest-energy transitions in each spectral domain are dominated by an excitation from the bonding or antibonding combinations of the relevant atomic orbitals on the sulfurs into the $\sigma_{\text{SS}}^{*}$ orbital. This pairing of transitions manifests as bisignate features in the low-energy region of the ECD spectra. On the basis of the simple Bergson model for the low-energy transitions in the disulfide chromophore,$^{97,98}$ Lindberg and Michl$^{89}$ formulated a quadrant rule for the optical activity of the two low-energy valence transitions (dominated by excitations from the symmetric and antisymmetric combinations of non-bonding 3p orbitals on sulfurs to the $\sigma_{\text{SS}}^{*}$ orbital) in organic disulfides. This rule relates the sign of the long-wavelength Cotton effect across the four dihedral quadrants and is a specific case of the C₂-rule for general chromophores of effective C₂ symmetry.$^{99}$ Woody extended the theoretical analysis to also include the absolute sign of the lowest ECD band within each quadrant,$^{100}$ providing

predictions in agreement with experimental results across different dihedral angles.$^{101-104}$

For the M-helical form considered here, we find a negative-first Cotton effect, consistent with the quadrant rule.$^{100}$ The intensity asymmetry of the lowest-energy valence couplet is attributed to different intrafragment (i.e., CH₃S-) contributions to the ECD signal of each transition. At higher energies, the electronic coupling between the (core) orbitals on different fragments decreases, reducing both the energetic splitting between symmetric and antisymmetric combinations as well as the intrafragment ECD contribution itself, since the environment is increasingly achiral on each fragment. As a consequence, the paired core transitions become near-degenerate (energy splitting of a few meV or less) with weak rotational strengths of almost equal magnitudes but opposite signs (see Table I). After additionally accounting for sources of broadening, including finite core-hole lifetimes (~0.1 and ~0.5 eV at the sulfur L- and K-edges,$^{105,106}$ respectively), the bisignate feature may coalesce into a single signal with the absolute sign given by the most intense of the pair of transitions (this will be seen in Fig. S2). In sum, the effective signals are reduced by orders of magnitude in the x-ray region (see the accumulated values in Table I).

The linear absorption profiles with the full interaction and the electric-dipole approximation essentially coincide across the four spectral regions. However, as shown by the underlying

<figure><img src="image_12.png" /><figcaption>**FIG. 2.** Comparison of full and truncated differential oscillator strength $\Delta f(\theta, \phi)$ across the spectral regions (valence, L₃-edge, L₁-edge, and K-edge). Transitions of [(a)–(d)] A and [(e)–(h)] B symmetry in H₂S₂. The black arrow points along the direction of the wave vector for the anisotropic ECD intensity given in Table I. The truncated ECD is represented by the smooth surface, whereas the full ECD is shown as individual points generated with a 5810-point Lebedev grid ($L_{\text{max}}$ = 131). Blue: negative and red: positive ECD signal. Note that different scaling factors (upper right corner) have been applied. The corresponding isotropic differential oscillator strengths are indicated in each inset.</figcaption></figure>
oscillator strengths in Table I, non-dipolar effects at the K-edge lead to intensity redistribution among the underlying near-degenerate transitions (i.e., unrelated to the arbitrary mixing allowed for degenerate states). Nonetheless, the overall spectral profiles within and beyond the electric-dipole approximation are essentially identical because of the nearly overlapping transitions. This is consistent with our previous findings for Cl K-edge absorption in $TiCl_4$.²⁸ In contrast, the beyond-first-order effects become evident in the differential K-edge absorption profile because of the signed nature of the underlying quantities. This leads to a factor-of-two overestimation of the ECD within the conventional first-order treatment. Introducing third-order contributions largely corrects this discrepancy at the sulfur K-edge, but going to higher orders in the expansion is not a general remedy, as will be discussed below.

To better understand the nature of the chiral response across spectral regions, we computed the underlying anisotropic differential oscillator strength distributions, considering now the smaller $H_2S_2$. Figure 2 shows the full ECD distributions (points), compared with the first-order truncated counterparts (surfaces). The solid angle represents the propagation direction, the distance from the origin (COM) indicates the magnitude of the associated signal, and the color indicates its sign. Note that different scaling factors have been applied across the transitions (see the upper right corner of each inset). The $C_2$-rotation axis coincides with the z-axis, whereas the disulfide bond is along the x-axis. The shapes of the anisotropic distributions can be understood by decomposing the first-order signals into isotropic and d-orbital contributions [Eq. (26)]. The resulting orbital weights of the excitations plotted in Fig. 2 are reported in Table II.

From symmetry considerations detailed in the supplementary material, we find that contributions from $d_{xz}$ and $d_{yz}$ vanish for excitations of both A and B symmetries. In A symmetry, $R_{z2}^{[1]}$ is also zero by symmetry, such that the s- and $d_{z^2}$-contributions come in a fixed ratio, giving a toroid in the xy-plane. For valence and $L_3$-edge excitations of A symmetry, this shape is modulated by the $d_{x^2-y^2}$-contribution, giving the shape of a biconcave disk elongated along the y- and x-axis, respectively, depending on its relative sign. For the $L_1$- and K-edge excitations, the $d_{xy}$-contribution completely dominates. This is also the case for the corresponding excitations in the

B symmetry. It may be noted that in the A symmetry, the E1–M1 contribution to the $d_{xy}$ term is zero by symmetry, whereas in the B symmetry, the E1–E2 and E1–M1 contributions are of similar magnitude and the same sign. Continuing to the $L_3$ excitation of the B symmetry, the angular plot resembles that of its counterpart in the A symmetry, albeit with opposite overall sign and rotated $\pi/2$ about the molecular axes. The latter can be understood from the relative weights of the $d_{z^2}$- and $d_{x^2-y^2}$-contributions, as seen in Table II. Finally, the angular plot of the valence excitation in the B symmetry is a biconcave disk elongated along the z-axis, arising from the positive interference between s- and $d_{z^2}$-contributions contrary to the negative interference observed for its counterpart of the A symmetry.

Next, we compare these conventional first-order truncated ECD distributions with their full counterparts. For the valence, $L_3$-edge, and $L_1$-edge transitions, the anisotropic distributions virtually coincide, thereby confirming the validity of the first-order truncated description also for the anisotropic signal in these energy regimes. On the other hand, the full and truncated ECD distribution for the K-edge transitions is seen to have the same overall shape, but markedly different size, such that the factor-of-two overestimation at first order of the isotropic response (Tables I and S1) arises largely from an overall scaling. Closer inspection of the angular distribution of the full ECD distribution reveals that the lobes are not strictly perpendicular as in a $d_{xy}$-orbital, hence indicating the contributions from solid harmonics of higher even angular momentum. To investigate this further, we provide the order-by-order contributions together with the full anisotropic ECD distribution in Fig. 3. Although $\Delta f^{[3]}$ ($\ell = 0, 2, 4$) resembles $\Delta f^{[1]}$, the inclusion of higher-order solid harmonics in the former leads to non-orthogonal lobes. Furthermore, the two distributions differ by an overall sign, such

<table><caption>**TABLE II.** Weights of contributions from solid harmonics ([24]) to the differential cross section $\Delta f(\theta, \phi)$ of selected transitions in $H_2S_2$. The weights have been scaled by the absolute value of the s-contribution. The numbers in parentheses are exponents of 10.</caption><thead><tr><th>Irrep</th><th>Excitation</th><th>Δf<sub>iso</sub><sup>[1]</sup></th><th>s</th><th>d<sub>z<sup>2</sup></sub></th><th>d<sub>x<sup>2</sup>-y<sup>2</sup></sub></th><th>d<sub>xy</sub></th></tr></thead><tbody><tr><td rowspan="5">A</td><td>Valence</td><td>5.176(-05)</td><td>1.000</td><td>-0.447</td><td>0.300</td><td>-0.056</td></tr><tr><td>L<sub>3</sub>-edge</td><td>-9.118(-05)</td><td>-1.000</td><td>0.447</td><td>0.414</td><td>-0.031</td></tr><tr><td>L<sub>2</sub>-edge</td><td>-2.824(-04)</td><td>-1.000</td><td>0.447</td><td>0.408</td><td>-0.024</td></tr><tr><td>L<sub>1</sub>-edge</td><td>7.316(-05)</td><td>1.000</td><td>-0.447</td><td>-0.862</td><td>8.817</td></tr><tr><td>K-edge</td><td>1.657(-04)</td><td>1.000</td><td>-0.447</td><td>-0.772</td><td>9.108</td></tr><tr><td rowspan="5">B</td><td>Valence</td><td>-7.624(-05)</td><td>-1.000</td><td>-0.328</td><td>-0.090</td><td>0.314</td></tr><tr><td>L<sub>3</sub>-edge</td><td>9.186(-05)</td><td>1.000</td><td>-0.132</td><td>-0.592</td><td>-0.005</td></tr><tr><td>L<sub>2</sub>-edge</td><td>2.883(-04)</td><td>1.000</td><td>-0.128</td><td>-0.591</td><td>-0.008</td></tr><tr><td>L<sub>1</sub>-edge</td><td>-7.312(-05)</td><td>-1.000</td><td>2.503</td><td>-0.347</td><td>-8.262</td></tr><tr><td>K-edge</td><td>-1.657(-04)</td><td>-1.000</td><td>0.377</td><td>0.813</td><td>-9.125</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>**FIG. 3.** Convergence of the truncated anisotropic ECD distributions, i.e., $\Delta f^{[2n+1]}(\theta, \phi)$, for the K-edge transitions of (a) A and (b) B symmetry in $H_2S_2$. The superscript notation $\Delta f^{[m]}$ indicates accumulated contributions up through mth order. The black arrows point along the direction of the wave vector for the anisotropic ECD intensity given in Table S1.</figcaption></figure>
that the inclusion of $\Delta f^{[3]}$ decreases the ECD signal. A possible issue is the rate of convergence of the truncated interaction toward the full one. Our implementation of truncated interaction for linear absorption allows us to go to arbitrary order, a unique functionality of the DIRAC code. In our previous work on the absorption of linearly polarized light,²⁸ we demonstrated that the truncated treatment converges to the full interaction upon inclusion of higher-order terms, but for higher-energy transitions (photon energies beyond ~1000 eV), the convergence behavior was too slow for practical applications. In the present case, we assessed the convergence of the isotropic differential oscillator strength expansion at the sulfur K-edge and found that the relative error is below the threshold of 1% at seventh order (Table III). Indeed, as shown in Fig. 3, the $\Delta f^{[5]}$ distribution ($\ell = 0, 2, 4, 6$) is minute. A similar convergence rate is found for the isotropic linear oscillator strength. From these data, it follows that the rate of convergence is acceptable for the application to the sulfur K-edge.

To assess whether the truncated formalism is applicable at higher transition energies, we have performed an additional series of calculations of the K- and L₁-edge isotropic differential oscillator strength for the heavier H₂X₂ analogs (X = Se and Te). As shown in Sec. S2, we find that only the selenide L₁-edge converges at a sufficient rate, whereas it is too slow to practically converge for all other edges. In both cases, the first-order treatment not only overestimates the rotatory strengths of individual transitions but also most

critically incorrectly predicts the signs of (the intense) pairs of near-degenerate transitions. It thereby appears that the breakdown of the first-order description at the sulfur K-edge, and not at the sulfur L-edges, is a consequence of its order-of-magnitude higher transition energy (~215 and 2427 eV at the L₁- and K-edges).

## VI. CONCLUSION

We have reported the first implementation and application of the anisotropic and isotropic ECD signal using the full semi-classical light-matter interaction operator within a four-component relativistic framework. This simultaneous account of beyond-first-order light-matter interactions and relativistic effects provides two additional sources of ECD, which become increasingly important at high photon energies. The linear form of the light-matter interaction operator in the relativistic domain further enabled straightforward extension to a multipole-based scheme in the velocity representation that allows for the traditional (albeit, in general, ambiguous) decomposition into electric and magnetic contributions while retaining order-by-order gauge-origin independence.

The presented approach was used to investigate the ECD response of two prototypical disulfides, H₂S₂ and (CH₃S)₂, across the electromagnetic spectrum, from valence to core transitions. To quantify the implications of higher-order effects, we compared the results of the full interaction to those obtained within the traditional lowest-order non-vanishing (i.e., first-order) truncated generalized velocity representation. Going beyond the electric-dipole approximation at the sulfur K-edge leads to non-negligible intensity redistribution among near-degenerate transitions but with no visible implications on the linear absorption profile. On the other hand, the differential absorption profile is affected by such redistribution because of its signed nature. This leads to an overall factor-of-two overestimation. By examining the shapes of the underlying anisotropic ECD distributions, we find this discrepancy to largely originate from an overall scaling that is corrected upon introducing third-order contributions. Critically, the first-order treatment deteriorates at higher transition energies (beyond ~1000 eV) where it may even fail to predict the sign of individual differential oscillator strengths. At such energies, going to higher orders is not a practical remedy because of the slow convergence of the truncated interaction—the full interaction is a must.

From a practical point of view, the full interaction operator provides a compact and inherently gauge-origin invariant treatment of light-matter interactions, which is both implementation- and computation-wise competitive with traditional multipole-based schemes for oriented and isotropic linear spectroscopies. For this reason, we believe that the use of the full light-matter interaction will become the new standard for theoretical x-ray spectroscopies. The question of how to efficiently handle the rotational averaging for nonlinear light-matter interactions will be explored in future work.

## SUPPLEMENTARY MATERIAL

See **supplementary material** for linear and differential absorption spectra and corresponding tabulated oscillator strengths for H₂S₂, The analysis of the L₁- and K-edges of the H₂X₂, X = S, Se, and Te series, and the analysis of the structure and origin dependence of the rotational strength tensor.

<table><caption>**TABLE III.** Contributions to the isotropic linear and differential oscillator strength ($\Delta f^{[2n+1]}$ and $f^{[2n]}$) at various orders for H₂S₂, $n = 0, 1, 2, 3$, compared to the result of the full interaction for the two $1s_{1/2} \rightarrow \sigma_{SS}^*$ transitions of the A/B symmetry. The superscript notation $\square^{(-m)}$ indicates accumulated contributions up through mth order. In particular, the “ac”-labeled row indicates the (differential) oscillator strengths accumulated through 8th (9th) order. The errors upon truncation are defined as $\%\delta\Delta f_{iso}^{(-2n+1)} = |(\Delta f_{iso}^{(-2n+1)} - \Delta f_{iso}^{\text{full}})/\Delta f_{iso}^{\text{full}}| \times 100\%$. The results were obtained using the 4c-TD-PBE0 level of theory and the uncontracted aug-pcx-3/aug-pc3 basis set. The numbers in parentheses are exponents of 10.</caption><thead><tr><th colspan="5">4A ($1s_{1/2} \rightarrow \sigma_{SS}^*$); $\omega_A = 2427.883$ 49 (eV)</th></tr><tr><th>n</th><th>$f_{iso}^{[2n]}$</th><th>$\Delta f_{iso}^{[2n+1]}$</th><th>$\%\delta f_{iso}^{(-2n)}$</th><th>$\%\delta\Delta f_{iso}^{(-2n+1)}$</th></tr></thead><tbody><tr><td>0</td><td>5.5576(-05)</td><td>1.6569(-04)</td><td>97.99</td><td>109.10</td></tr><tr><td>1</td><td>3.4276(-03)</td><td>-1.0818(-04)</td><td>26.07</td><td>27.42</td></tr><tr><td>2</td><td>-8.1305(-04)</td><td>2.4151(-05)</td><td>3.36</td><td>3.06</td></tr><tr><td>3</td><td>9.9900(-05)</td><td>-2.2416(-06)</td><td>2.57(-1)</td><td>2.32(-1)</td></tr><tr><td>4</td><td>-7.5661(-06)</td><td>1.2938(-07)</td><td>1.69(-2)</td><td>3.96(-1)</td></tr><tr><td>ac</td><td>2.7625(-03)</td><td>7.9240(-05)</td><td>...</td><td>...</td></tr><tr><td>full</td><td>2.7630(-03)</td><td>7.9240(-05)</td><td>...</td><td>...</td></tr></tbody><tbody><tr><th colspan="5">4B ($1s_{1/2} \rightarrow \sigma_{SS}^*$); $\omega_B = 2427.833$ 34 (eV)</th></tr><tr><td>0</td><td>1.0512(-02)</td><td>-1.6570(-04)</td><td>35.24</td><td>109.10</td></tr><tr><td>1</td><td>-3.4601(-03)</td><td>1.0818(-04)</td><td>9.27</td><td>27.42</td></tr><tr><td>2</td><td>8.1318(-04)</td><td>-2.4153(-05)</td><td>1.19</td><td>3.06</td></tr><tr><td>3</td><td>-9.9897(-05)</td><td>2.2413(-06)</td><td>9.60(-2)</td><td>2.24(-1)</td></tr><tr><td>4</td><td>7.5513(-06)</td><td>-1.2878(-07)</td><td>1.43(-3)</td><td>3.98(-1)</td></tr><tr><td>ac</td><td>7.7729(-03)</td><td>-7.9433(-05)</td><td>...</td><td>...</td></tr><tr><td>full</td><td>7.7728(-03)</td><td>-7.9246(-05)</td><td>...</td><td>...</td></tr></tbody></table>
ACKNOWLEDGMENTS

N.H.L. acknowledges financial support from the Carlsberg Foundation (Grant No. CF16-0290). M.v.H. acknowledges funding from the French Ministry of Higher Education and Research. Computing time from CALMIP (Calcul en Midi-Pyrenées) and SNIC (Swedish National Infrastructure for Computing) at the National Supercomputer Center (NSC) are gratefully acknowledged.

AUTHOR DECLARATIONS

Conflict of Interest

There are no conflicts to declare.

DATA AVAILABILITY

The data that support the findings of this study are available
from the corresponding author upon reasonable request.

REFERENCES

¹A. Dorta-Urra and P. Bargueño, *Symmetry* **11**, 661 (2019).

²L. D. Barron, *Molecular Light Scattering and Optical Activity* (Cambridge University Press, 2004).

³P. E. Schipper and A. Rodger, *J. Am. Chem. Soc.* **105**, 4541 (1983).

⁴N. H. List, J. Knoops, J. Rubio-Magnieto, J. Idé, D. Beljonne, P. Norman, M. Surin, and M. Linares, *J. Am. Chem. Soc.* **139**, 14947 (2017).

⁵N. Berova, K. Nakanishi, and R. W. Woody, *Circular Dichroism: Principles and Applications* (John Wiley & Sons, 2000).

⁶D. Lindorfer and T. Renger, *J. Phys. Chem. B* **122**, 2747 (2018).

⁷P. Akhtar, D. Lindorfer, M. Lingvay, K. Pawlak, O. Zsiros, G. Siligardi, T. Jávorfi, M. Dorogi, B. Ughy, G. Garab, T. Renger, and P. H. Lambrev, *J. Phys. Chem. B* **123**, 1090 (2019).

⁸M. Arago, Mem. Inst. Fr. **1**, 93 (1811); available at [https://babel.hathitrust.org/cgi/pt?id=ucm.5326746608 & view=image & seq=105](https://babel.hathitrust.org/cgi/pt?id=ucm.5326746608 & view=image & seq=105).

⁹J.-B. Biot, Mem. Inst. Fr. **1**, 1 (1812); available at [https://babel.hathitrust.org/cgi/pt?id=ucm.5326746626 & view=page & seq=11](https://babel.hathitrust.org/cgi/pt?id=ucm.5326746626 & view=page & seq=11).

¹⁰L. Alagna, T. Prosperi, S. Turchini, J. Goulon, A. Rogalev, C. Goulon-Ginet, C. R. Natoli, R. D. Peacock, and B. Stewart, *Phys. Rev. Lett.* **80**, 4799 (1998).

¹¹J. Goulon, C. Goulon-Ginet, A. Rogalev, V. Gotte, C. Malgrange, C. Brouder, and C. R. Natoli, *J. Chem. Phys.* **108**, 6394 (1998).

¹²J. Goulon, C. Goulon-Ginet, A. Rogalev, V. Gotte, C. Malgrange, and C. Brouder, *J. Synchrotron Radiat.* **6**, 673 (1999).

¹³B. Stewart, R. D. Peacock, L. Alagna, T. Prosperi, S. Turchini, J. Goulon, A. Rogalev, and C. Goulon-Ginet, *J. Am. Chem. Soc.* **121**, 10233 (1999).

¹⁴M. Tanaka, K. Nakagawa, A. Agui, K. Fujii, and A. Yokoya, *Phys. Scr.* **2005**, 873.

¹⁵K. Nakagawa, T. Matsui, Y. Izumi, A. Agui, M. Tanaka, and T. Muro, *Radiat. Phys. Chem.* **78**, 1198 (2009).

¹⁶Y. Izumi, A. Imazu, A. Mimoto, M. Tanaka, K. Nakagawa, M. Tanaka, A. Agui, and T. Muro, *J. Phys.: Conf. Ser.* **190**, 012209 (2009).

¹⁷Y. Izumi, M. Tanabe, A. Imazu, A. Mimoto, M. Tanaka, A. Agui, T. Muro, and K. Nakagawa, *J. Chem. Phys.* **138**, 074305 (2013).

¹⁸J. Goulon, F. Sette, C. Moise, A. Fontaine, D. Perey, P. Rudolf, and F. Baudelet, *Jpn. J. Appl. Phys., Part 1* **32**, 284 (1993).

¹⁹A. Rogalev, J. Goulon, F. Wilhelm, and A. Bosak, in *Magnetism and Synchrotron Radiation*, edited by E. Beaurepaire, H. Bulou, F. Scheurer, and K. Jean-Paul (Springer Berlin Heidelberg, Berlin, Heidelberg, 2010), pp. 169–190.

²⁰C. Brouder, *J. Phys.: Condens. Matter* **2**, 701 (1990).

²¹J. Goulon, in *Rayonnement Synchrotron Polarisé, Electrons Polarisés et Magnetisme*, edited by E. Beaurepaire, B. Carrière, and J.-P. Kappler (IPCMS, Mittelwihr, Haut-Rhin, 1990), pp. 333–386.

²²S. Turchini, N. Zema, S. Zennaray, L. Alagna, B. Stewart, R. D. Peacock, and T. Prosperi, *J. Am. Chem. Soc.* **126**, 4532 (2004).

²³A. E. Hansen and T. D. Bouman, "Natural chiroptical spectroscopy: Theory and computations," in *Advances in Chemical Physics* (John Wiley & Sons, 2007), pp. 545–644.

²⁴F. London, *J. Phys. Radium* **8**, 397 (1937).

²⁵K. L. Bak, A. E. Hansen, K. Ruud, T. Helgaker, J. Olsen, and P. Jørgensen, *Theor. Chim. Acta* **90**, 441 (1995).

²⁶N. H. List, J. Kauczor, T. Saue, H. J. A. Jensen, and P. Norman, *J. Chem. Phys.* **142**, 244111 (2015).

²⁷N. H. List, T. Saue, and P. Norman, *Mol. Phys.* **115**, 63 (2017).

²⁸N. H. List, T. R. L. Melin, M. van Horn, and T. Saue, *J. Chem. Phys.* **152**, 184110 (2020).

²⁹S. Villaume and P. Norman, *Chirality* **21**, E13 (2009).

³⁰Y. Zhang, J. R. Rouxel, J. Autschbach, N. Govind, and S. Mukamel, *Chem. Sci.* **8**, 5969 (2017).

³¹O. Takahashi, M. Kimoto, and L. G. M. Pettersson, *Chem. Phys.* **450-451**, 109 (2015).

³²V. Kimberg and N. Kosugi, *J. Chem. Phys.* **126**, 245101 (2007).

³³A. Jiemchooroj, U. Ekström, and P. Norman, *J. Chem. Phys.* **127**, 165104 (2007).

³⁴V. Carravetta, O. Plachkevytch, O. Vahtras, and H. Ågren, *Chem. Phys. Lett.* **275**, 70 (1997).

³⁵O. Plashkevych, V. Carravetta, O. Vahtras, and H. Ågren, *Chem. Phys.* **232**, 49 (1998).

³⁶J. Kaniauskas, Liet. Fiz. Rink **14**, 463 (1974) (in Russian; English translation available from the present authors); available at [http://www.itpa.lt/LFD/profile.html](http://www.itpa.lt/LFD/profile.html).

³⁷I. P. Grant, *Relativistic Quantum Theory of Atoms and Molecules: Theory and Computation* (Springer Science & Business Media, 2007), Vol. 40, pp. 673–713.

³⁸S. Urban, E. Herbst, P. Mittler, G. Winnewisser, K. M. T. Yamada, and M. Winnewisser, *J. Mol. Spectrosc.* **137**, 327 (1989).

³⁹E. Herbst and G. Winnewisser, *Chem. Phys. Lett.* **155**, 572 (1989).

⁴⁰G. S. Maciel, P. R. P. Barreto, F. Palazzetti, A. Lombardi, and V. Aquilanti, *J. Chem. Phys.* **129**, 164302 (2008).

⁴¹V. Aquilanti, M. Ragni, A. C. P. Bitencourt, G. S. Maciel, and F. V. Prudente, *J. Phys. Chem. A* **113**, 3804 (2009).

⁴²J. M. Thornton, *J. Mol. Biol.* **151**, 261 (1981).

⁴³A. Rauk, *J. Am. Chem. Soc.* **106**, 6517 (1984).

⁴⁴T.-K. Ha and W. Cencek, *Chem. Phys. Lett.* **182**, 519 (1991).

⁴⁵M. Pericou-Cayere, M. Rerat, and A. Dargelos, *Chem. Phys.* **226**, 297 (1998).

⁴⁶C. Diedrich and S. Grimme, *J. Phys. Chem. A* **107**, 2524 (2003).

⁴⁷M. Pecul and K. Ruud, in *Response Theory and Molecular Properties (A Tribute to Jan Lindenberg and Poul Jørgensen)*, Advances in Quantum Chemistry, edited by H. J. Aa. Jensen (Academic Press, 2005), Vol. 50, pp. 185–212.

⁴⁸L. Bednárová, P. Bouř, and P. Maloň, *Chirality* **22**, 514 (2010).

⁴⁹M. Scott, D. R. Rehn, S. Coriani, P. Norman, and A. Dreuw, *J. Chem. Phys.* **154**, 064107 (2021).

⁵⁰DIRAC, a relativistic *ab initio* electronic structure program, Release DIRAC21 (2021), written by R. Bast, A. S. P. Gomes, T. Saue, L. Visscher, and H. J. Aa. Jensen, with contributions from I. A. Aucar, V. Bakken, K. G. Dyall, S. Dubilard, U. Ekström, E. Eliav, T. Enevoldsen, E. Fašhauer, T. Fleig, O. Fossgaard, L. Halbert, E. D. Hedegård, T. Helgaker, B. Helmich-Paris, J. Henriksson, M. Iliaš, C. R. Jacob, S. Knecht, S. Komorovsky, O. Kullie, J. K. Laerdahl, C. V. Larsen, Y. S. Lee, N. H. List, H. S. Nataraj, M. K. Nayak, P. Norman, G. Olejniczak, J. Olsen, J. M. H. Olsen, A. Papadopoulos, Y. C. Park, J. K. Pedersen, M. Pernpointner, J. V. Pototschnig, R. Di Remigio, M. Repisky, K. Ruud, P. Salek, B. Schimmelpfennig, B. Senjean, A. Shee, J. Sikkema, A. Sunaga, A. J. Thorvaldsen, J. Thyssen, J. van Stralen, M. L. Vidal, S. Villaume, O. Visser, T. Winther, and S. Yamamoto (available at [http://dx.doi.org/10.5281/zenodo.4836496](http://dx.doi.org/10.5281/zenodo.4836496), see also [http://www.diracprogram.org](http://www.diracprogram.org)).

⁵¹T. Saue, R. Bast, A. S. P. Gomes, H. J. Aa. Jensen, L. Visscher, I. A. Aucar, R. Di Remigio, K. G. Dyall, E. Eliav, E. Fasshauer, T. Fleig, L. Halbert, E. D. Hedegård, B. Helmich-Paris, M. Iliaš, C. R. Jacob, S. Knecht, J. K. Laerdahl, M. L. Vidal, M. K. Nayak, M. Olejniczak, J. M. H. Olsen, M. Pernpointner, B. Senjean, A. Shee, A. Sunaga, and J. N. P. van Stralen, *J. Chem. Phys.* **152**, 204104 (2020).

⁵²S. E. Braslavsky, *Pure Appl. Chem.* **79**, 293 (2007).
<footnote n="53">⁵³A. E. Hansen and J. Avery, *Chem. Phys. Lett.* **13**, 396 (1972).</footnote>

<footnote n="54">⁵⁴S. Bernadotte, A. J. Atkins, and C. R. Jacob, *J. Chem. Phys.* **137**, 204106 (2012).</footnote>

<footnote n="55">⁵⁵F. Bloch, in *W. Heisenberg und die Physik unserer Zeit*, edited by F. Bopp (Vieweg & Sohn, Braunschweig, 1961).</footnote>

<footnote n="56">⁵⁶D. H. Kobe, *Am. J. Phys.* **50**, 128 (1982).</footnote>

<footnote n="57">⁵⁷A. M. Stewart, *J. Phys. A: Math. Gen.* **32**, 6091 (1999).</footnote>

<footnote n="58">⁵⁸T. Saue, in *Relativistic Electronic Structure Theory. Part 1. Fundamentals*, edited by P. Schwerdtfeger (Elsevier, Amsterdam, 2002), p. 332.</footnote>

<footnote n="59">⁵⁹P. Norman, K. Ruud, and T. Saue, *Principles and Practices of Molecular Properties: Theory, Modeling, and Simulations* (Wiley, Hoboken, NJ, 2018).</footnote>

<footnote n="60">⁶⁰T. B. Pedersen and A. E. Hansen, *Chem. Phys. Lett.* **246**, 1 (1995).</footnote>

<footnote n="61">⁶¹H.-G. Kuball, G. Sieber, S. Neubrech, B. Schultheis, and A. Schönhofer, *Chem. Phys.* **169**, 335 (1993).</footnote>

<footnote n="62">⁶²T. Saue and H. J. Aa. Jensen, *J. Chem. Phys.* **111**, 6211 (1999).</footnote>

<footnote n="63">⁶³P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, *J. Phys. Chem.* **98**, 11623 (1994).</footnote>

<footnote n="64">⁶⁴A. D. Becke, *J. Chem. Phys.* **98**, 5648 (1993).</footnote>

<footnote n="65">⁶⁵S. H. Vosko, L. Wilk, and M. Nusair, *Can. J. Phys.* **58**, 1200 (1980).</footnote>

<footnote n="66">⁶⁶C. Lee, W. Yang, and R. G. Parr, *Phys. Rev. B* **37**, 785 (1988).</footnote>

<footnote n="67">⁶⁷T. H. Dunning, *J. Chem. Phys.* **90**, 1007 (1989).</footnote>

<footnote n="68">⁶⁸D. E. Woon and T. H. Dunning, *J. Chem. Phys.* **98**, 1358 (1993).</footnote>

<footnote n="69">⁶⁹M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throllssell, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman, and D. J. Fox, *GAUSSIAN 16*, Revision C.01, Gaussian, Inc., Wallingford, CT, 2016.</footnote>

<footnote n="70">⁷⁰D. B. Craig and A. A. Dombkowski, *BMC Bioinf.* **14**, 346 (2013).</footnote>

<footnote n="71">⁷¹G. Pescitelli, L. Di Bari, and N. Berova, *Chem. Soc. Rev.* **40**, 4603 (2011).</footnote>

<footnote n="72">⁷²M. Stener, G. Fronzoni, and M. de Simone, *Chem. Phys. Lett.* **373**, 115 (2003).</footnote>

<footnote n="73">⁷³C. South, A. Shee, D. Mukherjee, A. K. Wilson, and T. Saue, *Phys. Chem. Chem. Phys.* **18**, 21010 (2016).</footnote>

<footnote n="74">⁷⁴M. Kadek, L. Konecny, B. Gao, M. Repisky, and K. Ruud, *Phys. Chem. Chem. Phys.* **17**, 22566 (2015).</footnote>

<footnote n="75">⁷⁵T. Fransson, D. Burdakova, and P. Norman, *Phys. Chem. Chem. Phys.* **18**, 13591 (2016).</footnote>

<footnote n="76">⁷⁶L. Visscher and K. G. Dyall, *At. Data Nucl. Data Tables* **67**, 207 (1997).</footnote>

<footnote n="77">⁷⁷J. P. Perdew, K. Burke, and M. Ernzerhof, *Phys. Rev. Lett.* **77**, 3865 (1996).</footnote>

<footnote n="78">⁷⁸C. Adamo and V. Barone, *J. Chem. Phys.* **110**, 6158 (1999).</footnote>

<footnote n="79">⁷⁹M. A. Ambroise and F. Jensen, *J. Chem. Theory Comput.* **15**, 325 (2019).</footnote>

<footnote n="80">⁸⁰F. Jensen, *J. Chem. Phys.* **115**, 9113 (2001).</footnote>

<footnote n="81">⁸¹F. Jensen, *J. Chem. Phys.* **116**, 7372 (2002).</footnote>

<footnote n="82">⁸²F. Jensen, *J. Chem. Phys.* **117**, 9234 (2002).</footnote>

<footnote n="83">⁸³L. Visscher, *Theor. Chem. Acc.* **98**, 68 (1997).</footnote>

<footnote n="84">⁸⁴I. Tokue, A. Hiraya, and K. Shobatake, *Chem. Phys.* **130**, 401 (1989).</footnote>

<footnote n="85">⁸⁵C. H. Hearn, E. Turcu, and J. A. Joens, *Atmos. Environ.*, Part A **24**, 1939 (1990).</footnote>

<footnote n="86">⁸⁶K. Schnorr, A. Bhattachacherjee, K. J. Oosterbaan, M. G. Delcey, Z. Yang, T. Xue, A. R. Attar, A. S. Chatterley, M. Head-Gordon, S. R. Leone, and O. Gessner, *J. Chem. Phys. Lett.* **10**, 1382 (2019).</footnote>

<footnote n="87">⁸⁷M. Ochmann, A. Hussain, I. von Ahnen, A. A. Cordones, K. Hong, J. H. Lee, R. Ma, K. Adamczyk, T. K. Kim, R. W. Schoenlein, O. Vendrell, and N. Huse, *J. Am. Chem. Soc.* **140**, 6554 (2018).</footnote>

<footnote n="88">⁸⁸A. P. Hitchcock, S. Bodeur, and M. Tronc, *Physica B* **158**, 257 (1989).</footnote>

<footnote n="89">⁸⁹J. Linderberg and J. Michl, *J. Am. Chem. Soc.* **92**, 2619 (1970).</footnote>

<footnote n="90">⁹⁰S. Svensson, A. Naves de Brito, M. P. Keane, N. Correia, and L. Karlsson, *Phys. Rev. A* **43**, 6441 (1991).</footnote>

<footnote n="91">⁹¹S. Svensson, A. Ausmees, S. J. Osborne, G. Bray, F. Gel'mukhanov, H. Ågren, A. Naves de Brito, O.-P. Sairanen, A. Kivimäki, E. Nömmiste, H. Aksela, and S. Aksela, *Phys. Rev. Lett.* **72**, 3021 (1994).</footnote>

<footnote n="92">⁹²B. T. Thole and G. van der Laan, *Phys. Rev. A* **38**, 1943 (1988).</footnote>

<footnote n="93">⁹³P. S. Bagus, H. Freund, H. Kuhlenbeck, and E. S. Ilton, *Chem. Phys. Lett.* **455**, 331 (2008).</footnote>

<footnote n="94">⁹⁴R. B. Bernini, L. B. G. da Silva, F. N. Rodrigues, L. H. Coutinho, A. B. Rocha, and G. G. B. de Souza, *J. Chem. Phys.* **136**, 144307 (2012).</footnote>

<footnote n="95">⁹⁵Y. Baba, K. Yoshii, and T. A. Sasaki, *J. Chem. Phys.* **105**, 8858 (1996).</footnote>

<footnote n="96">⁹⁶I. J. Pickering, M. Barney, J. J. H. Cotelesage, L. Vogt, M. J. Pushie, A. Nissan, R. C. Prince, and G. N. George, *J. Phys. Chem. A* **120**, 7279 (2016).</footnote>

<footnote n="97">⁹⁷G. Bergson, Ark. Kemi **12**, 233 (1958).</footnote>

<footnote n="98">⁹⁸G. Bergson, Ark. Kemi **18**, 409 (1961).</footnote>

<footnote n="99">⁹⁹G. Wagnière and W. Hug, *Tetrahedron Lett.* **11**, 4765 (1970).</footnote>

<footnote n="100">¹⁰⁰R. W. Woody, *Tetrahedron* **29**, 1273 (1973).</footnote>

<footnote n="101">¹⁰¹M. Carmack and L. A. Neubert, *J. Am. Chem. Soc.* **89**, 7134 (1967).</footnote>

<footnote n="102">¹⁰²R. M. Dodson and V. C. Nelson, *J. Org. Chem.* **33**, 3966 (1968).</footnote>

<footnote n="103">¹⁰³G. Claeson, *Acta Chem. Scand.* **22**, 2429 (1968).</footnote>

<footnote n="104">¹⁰⁴U. Ludescher and R. Schwyzer, *Helv. Chim. Acta* **54**, 1637 (1971).</footnote>

<footnote n="105">¹⁰⁵J. L. Campbell and T. Papp, *At. Data Nucl. Data Tables* **77**, 1 (2001).</footnote>

<footnote n="106">¹⁰⁶C. Nicolas and C. Miron, *J. Electron Spectrosc. Relat. Phenom.* **185**, 267 (2012).</footnote>