A Unified Strategy for the Chemically Intuitive Interpretation of Molecular Optical Response Properties

Optical Activity of Nonactin and Its Cation Complexes

Nonactin is a non-enantiomorphous (S4 symmetric), optically active natural product with a specific rotation of zero in solutions at all frequencies and temperatures. All optically active, non-enantiomorphous natural products have specific rotations of zero as a consequence of the spatial average of bisignate chiroptical (magnetoelectric or gyration) tensors with equal and opposite eigenvalues. Zeros that arise in the spatial average are distinct in principle, though not necessarily in practice, from zeros that arise in optical inactivity-chiroptical tensors with zero values for all elements as in centric molecules. Nonactin would be measurably optically active when oriented. The anisotropy of the optical activity of nonactin and its cation complexes, likewise S4 symmetric, are studied here by computation to emphasize the infelicitous linkage between optical activity and chirality. Computations show that changes in the conformation of the nonactin macrocycle upon complexation principally are responsible for diminishing the computed optical activity; the metals are incidental.

The eXact integral simplified time-dependent density functional theory (XsTD-DFT)

In the framework of simplified quantum chemistry methods, we introduce the eXact integral simplified time-dependent density functional theory (XsTD-DFT). This method is based on the simplified time-dependent density functional theory (sTD-DFT), where all semi-empirical two-electron integrals are replaced by exact one- and two-center two-electron integrals, while other approximations from sTD-DFT are kept. The performance of this new parameter-free XsTD-DFT method was benchmarked on excited state and (non)linear response properties, including ultra-violet/visible absorption, first hyperpolarizability, and two-photon absorption (2PA). For a set of 77 molecules, the results from the XsTDA approach were compared to the TDA data. XsTDA/B3LYP excitation energies only deviate on average by 0.14 eV from TDA while drastically cutting computational costs by a factor of 20 or more depending on the energy threshold chosen. The absolute deviations of excitation energies with respect to the full scheme are decreasing with increasing system size, showing the suitability of XsTDA/XsTD-DFT to treat large systems. Comparing XsTDA and its predecessor sTDA, the new scheme generally improves excitation energies and oscillator strengths, in particular, for charge transfer states. TD-DFT first hyperpolarizability frequency dispersions for a set of push-pull π-conjugated molecules are faithfully reproduced by XsTD-DFT, while the previous sTD-DFT method provides redshifted resonance energy positions. Excellent performance with respect to the experiment is observed for the 2PA spectrum of the enhanced green fluorescent protein. The obtained robust accuracy similar to TD-DFT at a fraction of the computational cost opens the way for a plethora of applications for large systems and in high throughput screening studies.

Molecular-Orbital Framework of Two-Electron Processes: Application to Auger and Intermolecular Coulomb Decay

States with core- or inner-shell vacancies, which are commonly created by absorption of high-energy photons, can decay by a two-electron process in which one electron fills the core hole and the second one is ejected. These processes accompany many X-ray spectroscopies. Depending on the nature of the initial core- or inner-shell-hole state and the decay valence-hole state, these processes are called Auger decay, intermolecular Coulomb decay, or electron-transfer-mediated decay. To connect many-body wave functions of the initial and final states with the molecular orbital picture of the decay, we introduce the concept of natural Auger orbitals (NAOs). NAOs are obtained by a two-step singular value decomposition of the two-body Dyson orbitals, reduced quantities that enter the expression of the decay rate in the Feshbach-Fano treatment. NAOs afford chemical insight and interpretation of the high-level ab initio calculations of Auger decay and related two-electron relaxation processes.

Natural orbitals and two-particle correlators as tools for the analysis of effective exchange couplings in solids

Using generalizations of spin-averaged natural orbitals and two-particle charge correlators for solids, we investigate the electronic structure of antiferromagnetic transition-metal oxides with a fully self-consistent, imaginary-time GW method. Our findings disagree with the Goodenough-Kanamori (GK) rules that are commonly used for the qualitative interpretation of such solids. First, we found a strong dependence of the natural orbital occupancies on momenta, contradicting GK assumptions. Second, along the momentum path, the character of natural orbitals changes. In particular, the contributions of oxygen 2s orbitals are important, which has not been considered in the GK rules. To analyze the influence of the electronic correlation on the values of effective exchange coupling constants, we use both natural orbitals and two-particle correlators and show that electronic screening modulates the degree of superexchange by stabilizing the charge-transfer contributions, which greatly affects these coupling constants. Finally, we give a set of predictions and recommendations regarding the use of density functional, Green's function, and wave-function methods for evaluating effective magnetic couplings in molecules and solids.

Multi-State Second-Order Nonlinear Optical Switches Incorporating One to Three Benzazolo-Oxazolidine Units: A Quantum Chemistry Investigation

This contribution employs quantum chemistry methods to describe the variations of the second nonlinear optical responses of molecular switches based on benzazolo-oxazolidine (BOX) units, connected by π-linkers, along their successive opening/closing. Under the fully closed forms, all of them display negligible first hyperpolarizability (β) values. When one BOX is opened, which is sketched as C→O, a push-pull π-conjugated segment is formed, having the potential to enhance β and to set the depolarization ratio (DR) to its one-dimensional-like value (DR = 5). This is observed when only one BOX is open, either for the monoBOX species (C→O) or for the diBOX (CC→CO) and triBOX (CCC→CCO) compounds, i.e., when the remaining BOXs stay closed. The next BOX openings have much different effects. For the diBOXs, the second opening (CO→OO) is associated with a decrease of β, and this decrease is tuned by controlling the conformation of the π-linker, i.e., the centrosymmetry of the whole compound because β vanishes in centrosymmetric compounds. For the triBOXs, the second opening gives rise to a Λ-shape compound, with a negligible change of β, but a decrease of the DR whereas, along the third opening, β remains similar and the DR decreases to the typical value of octupolar systems (DR = 1.5).

Plasmonic Circular Dichroism in Chiral Gold Nanowire Dimers

We report a computational study at the time-dependent density functional theory (TDDFT) level of the chiro-optical spectra of chiral gold nanowires coupled in dimers. Our goal is to explore whether it is possible to overcome destructive interference in single nanowires that damp chiral response in these systems and to achieve intense plasmonic circular dichroism (CD) through a coupling between the nanostructures. We predict a huge enhancement of circular dichroism at the plasmon resonance when two chiral nanowires are intimately coupled in an achiral relative arrangement. Such an effect is even more pronounced when two chiral nanowires are coupled in a chiral relative arrangement. Individual component maps of rotator strength, partial contributions according to the magnetic dipole component, and induced densities allow us to fully rationalize these findings, thus opening the way to the field of plasmonic CD and its rational design.

Spin-Orbit Natural Transition Orbitals and Spin-Forbidden Transitions

Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. The present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be nonzero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T₂-S₁ state mixing in pyrazine, spin-forbidden versus spin-allowed 4f-5d transitions in the Tb³⁺ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy)₃].

All-Atom Quantum Mechanical Calculation of the Second-Harmonic Generation of Fluorescent Proteins

Fluorescent proteins (FPs) are biotags of choice for second-harmonic imaging microscopy (SHIM). Because of their large size, computing their second-harmonic generation (SHG) response represents a great challenge for quantum chemistry. In this contribution, we propose a new all-atom quantum mechanics methodology to compute SHG of large systems. This is now possible because of two recent implementations: the tight-binding GFN2-xTB method to optimize geometries and a related version of the simplified time-dependent density functional theory (sTD-DFT-xTB) to evaluate quadratic response functions. In addition, a new dual-threshold configuration selection scheme is introduced to reduce the computational costs while retaining overall similar accuracy. This methodology was tested to evaluate the SHG of the proteins iLOV and bacteriorhodopsin (bR). In the case of bR, quantitative agreement with respect to experiment was reached for the out-of-resonance low-energy part of the β_HRS frequency dispersion. This work paves the way toward an accurate prediction of the SHG of large structures-a requirement for the design of new and improved SHIM biotags.

Perspective on Simplified Quantum Chemistry Methods for Excited States and Response Properties

We review recent developments in the framework of simplified quantum chemistry for excited state and optical response properties (sTD-DFT) and present future challenges for new method developments to improve accuracy and extend the range of application. In recent years, the scope of sTD-DFT was extended to molecular response calculations of the polarizability, optical rotation, first hyperpolarizability, two-photon absorption (2PA), and excited-state absorption for large systems with hundreds to thousands of atoms. The recently introduced spin-flip simplified time-dependent density functional theory (SF-sTD-DFT) variant enables an ultrafast treatment for diradicals and related strongly correlated systems. A few drawbacks were also identified, specifically for the computation of 2PA cross sections. We propose solutions to this problem and how to generally improve the accuracy of simplified schemes. New possible simplified schemes are also introduced for strongly correlated systems, e.g., with a second-order perturbative correlation correction. Interpretation tools that can extract chemical structure-property relationships from excited state or response calculations are also discussed. In particular, the recently introduced method-agnostic RespA approach based on natural response orbitals (NROs) as the key concept is employed.