Unraveling charge transfer processes with the quantum theory of atoms-in-molecules

Seeking for Optimal Excited States in Photoinduced Electron-Transfer Processes─The Case Study of Brooker's Merocyanine

Material design enters an era in which control of electrons in atoms, molecules, and materials is an essential property to be predicted and thoroughly understood in view of discovering new compounds with properties optimized toward specific optical/optoelectronic applications. π-electronic delocalization and charge separation/recombination enter notably into the set of features that are highly desirable to tailor. Diverse domains are particularly relying on photoinduced electron-transfer (PET), including fields of paramount importance such as energy production through light-harvesting, efficient chemoreceptive sensors, or organic field-effect transistors. In view of completing the arsenal of strategies in this area, we selected Brooker's merocyanine─a typical [D-π-A] compound─as the case study and examined from time-dependent density functional theory the opportunity offered by selected excited states to reach a suited manipulation of the charge transfer (CT) extent. In addition to the consideration of diagnostic tools able to spot the charge amount (i.e., magnitude of electron fraction) transferred upon excitation (q_CT), the spatial extent associated with such an electronic transition or CT length (D_CT), as well as the corresponding variation in dipole moment between the ground and the excited states (μ_CT), further analysis of the excitation process was undertaken. The advantage of going beyond the above-mentioned molecular indicators─which can be considered as PET global indices─was explored on the basis of a partitioning of the electron density. Relevant insight was gained on the relation these global indices have with the evolution of (local) features characterizing either chemical bond or electron delocalization upon vertical excitations.

Halogen Bonding Involving Gold Nucleophiles in Different Oxidation States

A single-crystal X-ray diffraction (XRD) study of diaryliodonium tetrachloroaurates (or, in the recent terminology, tetrachloridoaurates), [(p-XC₆H₄)₂I][AuCl₄] (X = Cl, 1; Br, 2), was performed for 1 (the structure is denoted as 1a to show similarity with the isomorphic structure 2a) and two polymorphs─2a (obtained from MeOH) and 2b (from 1,2-C₂H₄Cl₂). Examination of the XRD data for these three structures revealed 2-center C-X···Au^III (X = Cl and Br) and 3-center bifurcated C-Br···(Cl-Au) halogen bonding (abbreviated as XB) between the p-Cl or p-Br atoms of the diaryliodonium cations and the gold(III) atom of [AuCl₄]^-. The noncovalent nature of Au^III-involving interactions, the nucleophilicity of the gold(III) atoms, and the electrophilic role of p-X atoms of the diaryliodonium cations in the XBs were studied by a set of complementary computational methods. Combined experimental and theoretical studies allowed the recognition of the d-nucleophilicity of the [d⁸Au^III] atom which, regardless of its rather substantial formal 3+ charge, can function as a d-nucleophilic partner of XB. This conclusion was also supported by theoretical calculations performed for the structures' refcodes BINXOM and ICSD 62511; the obtained data verified the nucleophilicity of Au^III toward a K⁺ ions or a σ-(Cl)-hole, respectively. All our results, together with consideration of relevant literature, indicate that gold atoms in the three oxidation states (0, I, and even III) exhibit nucleophilicity in XBs.

Partition of the electronic energy of the PM7 method via the interacting quantum atoms approach

Partitions of the electronic energy such as that provided by the Interacting Quantum Atoms (IQA) approach have given valuable insights for numerous chemical systems and processes. Unfortunately, this kind of analysis may involve the integration of scalar fields over very irregular volumes, a condition which leads to a large and often prohibitive computational effort. These circumstances have limited the use of these energy partitions to systems comprising a few tens of atoms at most. On the other hand, semiempirical methods have proved useful in the study of systems of several thousands of atoms. Therefore, the goal of this work is to carry out partitions of the semiempirical method PM7 in compliance with the IQA approach. For this purpose, we computed one- and two-atomic energetic contributions whose sum equals the PM7 electronic energy. We illustrate how one might exploit the partition of electronic energies computed via the PM7 method by considering small organic and inorganic molecules and the energetics of individual hydrogen bond interactions within several water clusters which include (H₂O)₃₀, (H₂O)₅₀ and (H₂O)₁₀₀. We also considered the solvation of the amphiphilic caprylate anion to exemplify how to exploit the energy partition proposed in this paper. Overall, this investigation shows how the approach put forward herein might give further insights of the interactions occurring within complex systems in physical and biological chemistry.

Inorganic–Organic {dz2-MIIS4}···π-Hole Stacking in Reverse Sandwich Structures. The Case of Cocrystals of Group 10 Metal Dithiocarbamates with Electron-deficient Arenes

Cocrystallization of the dithiocarbamate complexes [M(S2CNEt2)2] (M = Ni 1, Pd 2, Pt 3) and X-substituted perfluoroarenes (X = I, Br; 1,2-dibromoperfluorobenzene FBrB and 1,2-diiodoperfluorobenzene FIB) gives isomorphous cocrystals of...

Spin-orbit coupling as a probe to decipher halogen bonding

The nature of halogen-bond interactions is scrutinized from the perspective of astatine, the heaviest halogen element. Potentially the strongest halogen-bond donor, its ability is shown to be deeply affected by relativistic effects and especially by the spin-orbit coupling. Complexes between a series of XY dihalogens (X, Y = At, I, Br, Cl and F) and ammonia are studied with two-component relativistic quantum calculations, revealing that the spin-orbit interaction leads to a weaker halogen-bond donating ability of the diastatine species with respect to diiodine. In addition, the donating ability of the lighter halogen elements, iodine and bromine, in the AtI and AtBr species is more decreased by the spin-orbit coupling than that of astatine. This can only be rationalized from the evolution of a charge-transfer descriptor, the local electrophilicity ω+S,max, determined for the pre-reactive XY species. Finally, the investigation of the spin-orbit coupling effects by means of quantum chemical topology methods allows us to unveil the connection between the astatine propensity to form charge-shift bonds and the astatine ability to engage in halogen bonds.

Detection of simple inorganic and organic molecules over Cu-decorated circumcoronene: a combined DFT and QTAIM study

Nowadays graphene materials have attracted a considerable attention because of their potential utilization as gas sensors, biosensors, or adsorbents. Doping or decorating the graphene surface with transition metals can significantly tune its electronic properties and chemical reactivity. Circumcoronene, being a polyaromatic hydrocarbon composed of 19 benzene rings, can be used as a model system of a tiny graphene quantum dot. The adsorption of a set of small molecules (water, hydrogen peroxide, methanol, ethanol, and oxygen) over the copper-decorated circumcoronene was theoretically investigated using density functional theory (DFT) and Bader's quantum theory of atoms in molecules (QTAIM). Following the obtained B3LYP energies, the adsorption of O2 and the chemisorption of H2O2 were found to be energetically the most favorable, with energetic outcomes of -3.6 eV and -3.7 eV, respectively. Moreover, an H2O2 molecule was decomposed during the chemisorption on the Cu atom to form a neutral Cu(OH)2 molecule. Changes in the electronic structure of the studied systems, in particular the oxidation of copper, after the adsorption were investigated within the framework of QTAIM (e.g., charges, bond critical points, and delocalization indices) and partial density of states (PDOS) analysis. The results of this study suggest the suitability of the Cu-decorated graphene materials as adsorbents and/or gas sensors in practical applications.

Hydrogen detachment driven by a repulsive 1πσ* state – an electron localization function study of 3-amino-1,2,4-triazole

Electron localization function analysis reveals the details of a charge induced hydrogen detachment mechanism of 3-amino-1,2,4-triazole, identified recently as responsible for phototautomerization of the molecule.

A quantitative analysis of light-driven charge transfer processes using voronoi partitioning of time dependent DFT-derived electron densities

An analytical method is presented that provides quantitative insight into light-driven electron density rearrangement using the output of standard time-dependent density functional theory (TD-DFT) computations on molecular compounds. Using final and initial electron densities for photochemical processes, the subtraction of summed electron density in each atom-centered Voronoi polyhedron yields the electronic charge difference, QVECD . This subtractive method can also be used with Bader, Mulliken and Hirshfeld charges. A validation study shows QVECD to have the most consistent performance across basis sets and good conservation of charge between electronic states. Besides vertical transitions, relaxation processes can be investigated as well. Significant electron transfer is computed for isomerization on the excited state energy surface of azobenzene. A number of linear anilinepyridinium donor-bridge-acceptor chromophores was examined using QVECD to unravel the influence of its pi-conjugated bridge on charge separation. Finally, the usefulness of the presented method as a tool in optimizing charge transfer is shown for a homologous series of organometallic pigments. The presented work allows facile calculation of a novel, relevant quantity describing charge transfer processes at the atomic level. © 2017 Wiley Periodicals, Inc.

Driving Forces in the Sharpless Epoxidation Reaction: A Coupled AIMD/QTAIM Study

In order to better understand the epoxide-formation step of the Sharpless epoxidation process, a set of 263 oxygen-transfer reactions reflecting the complexity of the Sharpless epoxidation process were studied using density functional theory (DFT) and Bader's quantum theory of atoms in molecules (QTAIM). The diversity within these reactions reflects the different ligands in the coordination sphere of vanadium and also different substrates (alkene and an allylic alcohol both free and in the form of an alcoxo ligand). The transition states for 76 of these reactions were also characterized using DFT and QTAIM, allowing for an estimation of the impact of the different ligands and substrates on the activation barriers. A smaller subset of the latter was further subjected to an ab initio molecular dynamics (AIMD) simulation coupled to QTAIM analysis. The results show that the type of active catalyst plays an important role in the thermodinamic outcome of these reactions, with vanadium(V) tert-butylhydroperoxide adducts being responsible for the most exoenergetic reactions. On the other hand, the different ligands tested play only a limited role in modulating the thermodynamics and kinetics of these reactions. Moreover, no evidence was found to support a thermodynamic or kinetic preference for the epoxidation of an allylic alcohol over that of an unfunctionalized alkene. However, the results suggest that the reaction path is strongly influenced by the orientation of the substrate upon approximation to the active catalyst, confirming the well-known regioselectivity of the Sharpless epoxidation process.

Molecular Relaxations in Supercooled Liquid and Glassy States of Amorphous Quinidine: Dielectric Spectroscopy and Density Functional Theory Approaches

In this article, we conduct a comprehensive molecular relaxation study of amorphous Quinidine above and below the glass-transition temperature (Tg) through broadband dielectric relaxation spectroscopy (BDS) experiments and theoretical density functional theory (DFT) calculations, as one major issue with the amorphous state of pharmaceuticals is life expectancy. These techniques enabled us to determine what kind of molecular motions are responsible, or not, for the devitrification of Quinidine. Parameters describing the complex molecular dynamics of amorphous Quinidine, such as Tg, the width of the α relaxation (βKWW), the temperature dependence of α-relaxation times (τα), the fragility index (m), and the apparent activation energy of secondary γ relaxation (Ea-γ), were characterized. Above Tg (> 60 °C), a medium degree of nonexponentiality (βKWW = 0.5) was evidenced. An intermediate value of the fragility index (m = 86) enabled us to consider Quinidine as a glass former of medium fragility. Below Tg (< 60 °C), one well-defined secondary γ relaxation, with an apparent activation energy of Ea-γ = 53.8 kJ/mol, was reported. From theoretical DFT calculations, we identified the most reactive part of Quinidine moieties through exploration of the potential energy surface. We evidenced that the clearly visible γ process has an intramolecular origin coming from the rotation of the CH(OH)C9H14N end group. An excess wing observed in amorphous Quinidine was found to be an unresolved Johari-Goldstein relaxation. These studies were supplemented by sub-Tg experimental evaluations of the life expectancy of amorphous Quinidine by X-ray powder diffraction and differential scanning calorimetry. We show that the difference between Tg and the onset temperature for crystallization, Tc, which is 30 K, is sufficiently large to avoid recrystallization of amorphous Quinidine during 16 months of storage under ambient conditions.