How Electronic Excitation Can be Used to Inhibit Some Mechanisms Associated to Substituent Effects

Explainable chemical artificial intelligence from accurate machine learning of real-space chemical descriptors

Machine-learned computational chemistry has led to a paradoxical situation in which molecular properties can be accurately predicted, but they are difficult to interpret. Explainable AI (XAI) tools can be used to analyze complex models, but they are highly dependent on the AI technique and the origin of the reference data. Alternatively, interpretable real-space tools can be employed directly, but they are often expensive to compute. To address this dilemma between explainability and accuracy, we developed SchNet4AIM, a SchNet-based architecture capable of dealing with local one-body (atomic) and two-body (interatomic) descriptors. The performance of SchNet4AIM is tested by predicting a wide collection of real-space quantities ranging from atomic charges and delocalization indices to pairwise interaction energies. The accuracy and speed of SchNet4AIM breaks the bottleneck that has prevented the use of real-space chemical descriptors in complex systems. We show that the group delocalization indices, arising from our physically rigorous atomistic predictions, provide reliable indicators of supramolecular binding events, thus contributing to the development of Explainable Chemical Artificial Intelligence (XCAI) models.

Calculation of the ELF in the excited state with single-determinant methods

Since its first definition, back in 1990, the electron localization function (ELF) has settled as one of the most commonly employed techniques to characterize the nature of the chemical bond in real space. Although most of the work using the ELF has focused on the study of ground-state chemical reactivity, a growing interest has blossomed to apply these techniques to the nearly unexplored realm of excited states and photochemistry. Since accurate excited electronic states usually require to account appropriately for electron correlation, the standard single-determinant ELF formulation cannot be blindly applied to them, and it is necessary to turn to correlated ELF descriptions based on the two-particle density matrix (2-PDM). The latter requires costly wavefunction approaches, unaffordable for most of the systems of current photochemical interest. Here, we compare the exact, 2-PDM-based ELF results with those of approximate 2-PDM reconstructions taken from reduced density matrix functional theory. Our approach is put to the test in a wide variety of representative scenarios, such as those provided by the lowest-lying excited electronic states of simple diatomic and polyatomic molecules. Altogether, our results suggest that even approximate 2-PDMs are able to accurately reproduce, on a general basis, the topological and statistical features of the ELF scalar field, paving the way toward the application of cost-effective methodologies, such as time-dependent-Hartree-Fock or time-dependent density functional theory, in the accurate description of the chemical bonding in excited states of photochemical relevance.

Interacting Quantum Atoms-A Review

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

Does the Intra-Atomic Deformation Energy of Interacting Quantum Atoms Represent Steric Energy?

We show that the mutual, through-space compression of atomic volume experienced by approaching topological atoms causes an exponential increase in the intra-atomic energy of those atoms, regardless of approach orientation. This insight was obtained using the modern energy partitioning method called interacting quantum atoms (IQA). This behaviour is consistent for all atoms except hydrogen, which can behave differently depending on its environment. Whilst all atoms experience charge transfer when they interact, the intra-atomic energy of the hydrogen atom is more vulnerable to these changes than larger atoms. The difference in behaviour is found to be due to hydrogen's lack of a core of electrons, which, in heavier atoms, consistently provide repulsion when compressed. As such, hydrogen atoms do not always provide steric hindrance. In accounting for hydrogen's unusual behaviour and demonstrating the exponential character of the intra-atomic energy in all other atoms, we provide evidence for IQA's intra-atomic energy as a quantitative description of steric energy.

Unveiling the role of intra and interatomic interactions in the energetics of reaction schemes: a quantum chemical topology analysis

In this work we present a detailed analysis of selected reaction schemes in terms of the atomic components of the electronic energy defined by the quantum theory of atoms in molecules and the interacting quantum atoms method. The aim is to provide an interpretation tool for the energy change involved in a chemical reaction by means of the atomic and interaction contributions to the energies of the molecules involved. Ring strain in cyclic alkanes, the resonance energy of aromatic and antiaromatic molecules, local aromaticity in polycyclic aromatic hydrocarbons, intermolecular bonding in hydrogen fluoride clusters, and hydration of d-block metal dications were selected for the study. It was found that in addition to the changes in the strong C-C interactions in the carbon skeleton of the organic molecular rings, other contributions not usually considered to be important such as those between C and H atoms (either bonded or not) need to be considered in order to account for the net energy changes. The analysis unveils the role of the ionic and covalent contributions to the hydrogen bonding in HF clusters and the energetic origin and extent of cooperative effects involved. Moreover, the "double-hump" behavior observed for the hydration energy trend of [M(H2O)6]2+ complexes is explained in terms of the deformation energy of the metal cation and the increasingly covalent metal-water interactions. In addition, proper comparisons with the description provided by other methodologies are briefly discussed. The topological approach proposed in this contribution proves to be useful for the description of energy changes of apposite reaction schemes in chemically meaningful terms.

Energetic Analysis of Conjugated Hydrocarbons Using the Interacting Quantum Atoms Method

A number of aromatic, antiaromatic, and nonaromatic organic molecules was analyzed in terms of the contributions to the electronic energy defined in the quantum theory of atoms in molecules and the interacting quantum atoms method. Regularities were found in the exchange and electrostatic interatomic energies showing trends that are closely related to those of the delocalization indices defined in the theory. In particular, the CC interaction energies between bonded atoms allow to rationalize the energetic stabilization associated with the bond length alternation in conjugated polyenes. This approach also provides support to Clar's sextet rules devised for aromatic systems. In addition, the H⋯H bonding found in some of the aromatic molecules studied was of an attractive nature, according to the stabilizing exchange interaction between the bonded H atoms. © 2017 Wiley Periodicals, Inc.

The ANANKE relative energy gradient (REG) method to automate IQA analysis over configurational change

Much chemical insight ultimately comes down to finding out which fragment of a total system behaves like the total system, in terms of an energy profile. A simple example is that of the water dimer, where this system is regarded as held together by a hydrogen bond. The hydrogen bond consists of two atoms (H···O), which energetically behave similarly to the total system (H₂O)₂. However, from a quantum mechanical point of view, each atom in the total system interacts with any other atom. Thus, the view that the hydrogen bond by itself governs the energetic stability of the water dimer needs rigorous justification. In this work, we propose a method that provides such a justification, in general, but only illustrated on the water dimer here. This method is based on the topological energy partitioning method called interacting quantum atoms (IQA). The method is implemented in the program ANANKE, which calculates correlations between the energy profile of the total system and those of subsystems (or fragments). ANANKE acts on the IQA energy contributions obtained for a sequence of full-system geometries controlled by a coordinate of interest (e.g. the O···H distance in the water dimer). Although applied only for the water dimer in this work, the method is general and able to explain the gauche effect, the torsional barrier in biphenyl, the arrow-pushing scheme of an enzymatic reaction (peptide hydrolysis in the HIV-1 Protease active site), and halogen-alkane nucleophilic substitution (S_N2) reactions. Those applications will appear elsewhere as separate and elaborated case studies; here we focus on the details of the ANANKE method and its justification, using the water dimer as a concrete case.