Charge-displacement analysis via natural orbitals for chemical valence: Charge transfer effects in coordination chemistry

The Chemical Bond at the Bottom of the Periodic Table: The Case of the Heavy Astatine and the Super-Heavy Tennessine

In this work, we study the chemical bond in molecules containing heavy and super-heavy elements according to the current state-of-the-art bonding models. An Energy Decomposition Analysis in combination with Natural Orbital for Chemical Valence (EDA-NOCV) within the relativistic four-component Dirac-Kohn-Sham (DKS) framework is employed, which allows to successfully include the spin-orbit coupling (SOC) effects on the chemical bond description. Simple halogen-bonded adducts ClX⋯L (X=At, Ts; L=NH3, Br-, H2O, CO) of astatine and tennessine have been selected to assess a trend on descending along a group, while modulating the ClX⋯L bond features through the different electronic nature of the ligand L. Interesting effects caused by SOC have been revealed: i) a huge increase of the ClTs dipole moment (which is almost twice as that of ClAt), ii) a lowering of the ClX⋯L bonding energy arising from different contributions to the ClX…L interaction energy strongly depending on the nature of L, iii) a quenching of one of the π back-donation components to the bond. In the ClTs(CO) adduct, the back-donation from ClTs to CO becomes the most important component. The analysis of the electronic structure of the ClX dimers allows for a clear interpretation of the SOC effects in these systems.

Conformation-Associated C···dz2-PtII Tetrel Bonding: The Case of Cyclometallated Platinum(II) Complex with 4-Cyanopyridyl Urea Ligand

The nucleophilic addition of 3-(4-cyanopyridin-2-yl)-1,1-dimethylurea (1) to cis-[Pt(CNXyl)2Cl2] (2) gave a new cyclometallated compound 3. It was characterized by NMR spectroscopy (1H, 13C, 195Pt) and high-resolution mass spectrometry, as well as crystallized to obtain two crystalline forms (3 and 3·2MeCN), whose structures were determined by X-ray diffraction. In the crystalline structure of 3, two conformers (3A and 3B) were identified, while the structure 3·2MeCN had only one conformer 3A. The conformers differed by orientation of the N,N-dimethylcarbamoyl moiety relative to the metallacycle plane. In both crystals 3 and 3·2MeCN, the molecules of the Pt(II) complex are associated into supramolecular dimers, either {3A}2 or {3B}2, via stacking interactions between the planes of two metal centers, which are additionally supported by hydrogen bonding. The theoretical consideration, utilizing a number of computational approaches, demonstrates that the C···dz2(Pt) interaction makes a significant contribution in the total stacking forces in the geometrically optimized dimer [3A]2 and reveals the dz2(Pt)→π*(PyCN) charge transfer (CT). The presence of such CT process allowed for marking the C···Pt contact as a new example of a rare studied phenomenon, namely, tetrel bonding, in which the metal site acts as a Lewis base (an acceptor of noncovalent interaction).

Monomeric gold hydrides for carbon dioxide reduction: ligand effect on the reactivity

We analyzed the ligand electronic effect in the reaction between a [LAu(I)H]0/- hydride species and CO2, leading to a coordinated formate [LAu(HCOO)]0/-. We explored 20 different ligands, such as carbenes, phosphines and others, carefully selected to cover a wide range of electron-donor and -acceptor properties. We included in the study the only ligand, an NHC-coordinated diphosphene, that, thus far, experimentally demonstrated facile and reversible reaction between the monomeric gold(I) hydride and carbon dioxide. We elucidated the previously unknown reaction mechanism, which resulted to be concerted and common to all the ligands: the gold-hydrogen bond attacks the carbon atom of CO2 with one oxygen atom coordinating to the gold center. A correlation between the ligand σ donor ability, which affects the electron density at the reactive site, and the kinetic activation barriers of the reaction has been found. This systematic study offers useful guidelines for the rational design of new ligands for this reaction, while suggesting a few promising and experimentally accessible potential candidates for the stoichiometric or catalytic CO2 activation.

Cooperative small molecule activation by apolar and weakly polar bonds through the lens of a suitable computational protocol

Small molecule activation processes are central in chemical research and cooperativity is a valuable tool for the fine-tuning of the efficiency of these reactions. In this contribution, we discuss recent and remarkable examples in which activation processes are mediated by bimetallic compounds featuring apolar or weakly polar metal-metal bonds. Relevant experimental breakthroughs are thoroughly analyzed from a computational perspective. We highlight how the rational and non-trivial application of selected computational approaches not only allows rationalization of the observed reactivities but also inferring of general principles applicable to activation processes, such as the breakdown of the structure-reactivity relationship in carbon dioxide activation in a cooperative framework. We finally provide a simple yet unbiased computational protocol to study these reactions, which can support experimental advances aimed at expanding the range of applications of apolar and weakly polar bonds as catalysts for small molecule activation.

Terminal end-on coordination of dinitrogen versus isoelectronic CO: A comparison using the charge displacement analysis

The metal dinitrogen bonding in a wide series of terminal end-on dinitrogen complexes is investigated with the charge displacement analysis based on natural orbitals of chemical valence (CD-NOCV). The effect of the σ donation and π backdonation on the NN bond are discussed and compared with the observations for a series of carbonyl complexes, published in 2016 by Tarantelli et al. The σ donation is relative invariant over the series of dinitrogen complexes and has no significant effect on the NN bond strength, whereas the π backdonation causes a considerable elongation of the NN bond. Some uncommon examples of weakly bound dinitrogen with blue-shifted stretching frequency compared to free N2 (ν = 2330 cm-1 ) are known. The dinitrogen bonding in these complexes is simulated with a point charge. Apparently, electrostatics account for the shortened N─N bond in these systems.

Radical-like reactivity for dihydrogen activation by coinage metal-aluminyl complexes: computational evidence inspired by experimental main group chemistry

The computational study of an unprecedented reactivity of coinage metal-aluminyl complexes with dihydrogen is reported. In close resemblance to group 14 dimetallenes and dimetallynes, the complexes are predicted to activate H₂ under mild conditions. Two different reaction pathways are found disclosing a common driving force, i.e., the nucleophilic behavior of the electron-sharing M-Al (M = Cu, Ag, Au) bond, which enables a cooperative and diradical-like mechanism. This mode of chemical reactivity emerges as a new paradigm for dihydrogen activation and calls for experimental feedback.

Widening the Landscape of Small Molecule Activation with Gold-Aluminyl Complexes: A Systematic Study of E-H (E=O, N) Bonds, SO2 and N2 O Activation

The electronic features of gold-aluminyl complexes have been thoroughly explored. Their similarity with Group 14 dimetallenes and other metal-aluminyl complexes suggests that their reactivity with small molecules beyond carbon dioxide could be accessed. In this work, the reactivity of the [^t Bu₃ PAuAl(NON)] (NON=4,5-bis(2,6 diisopropylanilido)-2,7-ditert-butyl-9,9-dimethylxanthene) complex towards water, ammonia, sulfur dioxide and nitrous oxide is computationally explored. The reaction mechanisms computed for each substrate strongly suggest that all activation processes are in principle experimentally feasible. Electronic structure analysis highlights that, in all cases, the reactivity is driven by the presence of the poorly polarized electron-sharing gold-aluminyl bond, which induces a radical-like reactivity of the complex towards all the substrates. A flat topology of the potential energy surface (PES) has been found for the reaction with N₂ O, where two almost isoenergetic transition states can be located along the same reaction coordinate with different geometries, suggesting that the N₂ O binding mode may not be a good indicator of the nature of N₂ O activation in a cooperative bimetallic reactivity. In addition, the catalytic potentialities of these complexes have been explored in the framework of nitrous oxide reduction. The study reveals that the [^t Bu₃ PAuAl(NON)] complex might be an efficient catalyst towards oxidation of phosphines (and boranes) via N₂ O reduction. These findings underline recurring trends in the novel chemistry of gold-aluminyl complexes and call for experimental feedbacks.

Mechanistic Study of Alkyne Insertion into Cu-Al and Au-Al Bonds: A Paradigm Shift for Coinage Metal Chemistry

In this work, the mechanism of the insertion reaction of 3-hexyne into Cu-Al and Au-Al bonds in M-aluminyl (M = Cu, Au) complexes is computationally elucidated. The mechanism is found to be radical-like, with the Cu-Al and Au-Al bonds acting as nucleophiles toward the alkyne, and predicts a less efficient reactivity for the gold-aluminyl complex. The proposed mechanism well rationalizes the kinetic (or thermodynamic) control on the formation of the syn (or anti) insertion product into the Cu-Al bond (i.e., dimetallated alkene) which has been recently reported. A comparative analysis of the electronic structure reveals that the reduced reactivity at the gold site─usually showing higher efficiency than copper as a "standard" electrophile in alkyne activation─arises from a common feature, i.e., the highly stable 6s Au orbital. The relativistic lowering of the 6s orbital, making it more suitable for accepting electron density and thus enhancing the electrophilicity of gold complexes, in the gold-aluminyl system is responsible for a less nucleophilic Au-Al bond and, consequently, a less efficient alkyne insertion. These findings demonstrate that the unconventional electronic structure and the electron-sharing nature of the M-Al bond induce a paradigm shift in the properties of coinage metal complexes. In particular, the peculiar radical-like reactivity, previously shown also with carbon dioxide, suggests that these complexes might efficiently insert/activate other small molecules, opening new and unexplored paths for their reactivity.

Local Charge-Displacement Analysis: Targeting Local Charge-Flows in Complex Intermolecular Interactions

Charge-displacement (CD) analysis has recently proven to be a simple and powerful scheme for quantitatively analyzing the profile of the charge redistribution occurring upon intermolecular interactions along a given interaction axis. However, when two molecular fragments bind through complex interactions involving multiple concurrent charge flows, ordinary CD analysis is capable of providing only an averaged picture of the related charge-flow profiles and no detailed information on each of them. In this article, we combine CD analysis with a Hirshfeld partitioning of the molecular charge redistribution for a local analysis on focused portions of the molecule, allowing for a detailed characterization of one charge flow at a time. The resulting scheme - the local charge-displacement (LCD) analysis - is tested on the intriguing case of the dimethyl sulfide (DMS)-sulfur dioxide (SO2) complex, characterized by concurrent charge flows relating to a sulfur-sulfur homo-chalcogen interaction and a pair of hydrogen bonds. The LCD scheme is then applied to the analysis of multiple hydrogen bonding in the acetic acid dimer, of base-pairing interactions in DNA, and of ambifunctional hydrogen bonding in the ammonia-pyridine complex.

Gold-Aluminyl and Gold-Diarylboryl Complexes: Bonding and Reactivity with Carbon Dioxide

The unconventional carbon dioxide insertion reaction of a gold-aluminyl [^tBu₃PAuAl(NON)] complex has been recently shown to be related to the electron-sharing character of the Au–Al bond that acts as a nucleophile and stabilizes the insertion product through a radical-like behavior. Since a gold-diarylboryl [IPrAuB(o-tol)₂] complex with similar reactivity features has been recently reported, in this work we computationally investigate the reaction of carbon dioxide with [LAuX] (L = phosphine, N-heterocyclic carbene (NHC); X = Al(NON), B(o-tol)₂) complexes to get insights into the Al/B anionic and gold ancillary ligand effects on the Au–Al/B bond nature, electronic structure, and reactivity of these compounds. We demonstrate that the Au–Al and Au–B bonds possess a similar electron-sharing nature, with diarylboryl complexes displaying a slightly more polarized bond as Au(δ⁺)–B(δ^–). This feature reduces the radical-like reactivity toward CO₂, and the Al/B anionic ligand effect is found to favor aluminyls over boryls, despite the greater oxophilicity of B. Remarkably, the ancillary ligand of gold has a negligible electronic trans effect on the Au–X bond and only a minor impact on the formation of the insertion product, which is slightly more stable with carbene ligands. Surprisingly, we find that the modification of the steric hindrance at the carbene site may exert a sizable control over the reaction, with more sterically hindered ligands thermodynamically disfavoring the formation of the CO₂ insertion product.

The Au−Al and Au−B bonds have both an electron-sharing nature, with the diarylboryl gold complexes displaying a more polarized Au^δ+−B^δ− bond. The gold ligand (phosphine or N-heterocyclic carbene) has a negligible electronic effect on the Au−X bond, consistently with a radical-like reactivity of the complexes with carbon dioxide, which favors the gold-aluminyl over the gold-diarylboryl complexes.

Unraveling differences in aluminyl and carbene coordination chemistry: bonding in gold complexes and reactivity with carbon dioxide

The electronic properties of aluminyl anions have been reported to be strictly related to those of carbenes, which are well-known to be easily tunable via selected structural modifications imposed on their backbone. Since peculiar reactivity of gold-aluminyl complexes towards carbon dioxide has been reported, leading to insertion of CO2 into the Au-Al bond, in this work the electronic structure and reactivity of Au-Al complexes with different aluminyl scaffolds have been systematically studied and compared to carbene analogues. The analyses reveal that, instead, aluminyls and carbenes display a very different behavior when bound to gold, with the aluminyls forming an electron-sharing and weakly polarized Au-Al bond, which turns out to be poorly modulated by structural modifications of the ligand. The reactivity of gold-aluminyl complexes towards CO2 shows, both qualitatively and quantitatively, similar reaction mechanisms, reflecting the scarce tunability of their electronic structure and bond nature. This work provides further insights and perspectives on the properties of the aluminyl anions and their behavior as coordination ligands.

What Singles out Aluminyl Anions? A Comparative Computational Study of the Carbon Dioxide Insertion Reaction in Gold-Aluminyl, -Gallyl, and -Indyl Complexes

Anionic aluminum(I) anions ("aluminyls") are the most recent discovery along Group 13 anions, and the understanding of the unconventional reactivity they are able to induce at a coordinated metal site is at an early stage. A striking example is the efficient insertion of carbon dioxide into the Au-Al bond of a gold-aluminyl complex. The reaction occurs via a cooperative mechanism, with the gold-aluminum bond being the actual nucleophile and the Al site also behaving as an electrophile. In the complex, the Au-Al bond has been shown to be mainly of an electron-sharing nature, with the two metal fragments displaying a diradical-like reactivity with CO2. In this work, the analogous reactivity with isostructural Au-X complexes (X = Al, Ga, and In) is computationally explored. We demonstrate that a kinetically and thermodynamically favorable reactivity with CO2 may only be expected for the gold-aluminyl complex. The Au-Al bond nature, which features the most (nonpolar) electron-sharing character among the Group 13 anions analyzed here, is responsible for its highest efficiency. The radical-like reactivity appears to be a key ingredient to stabilize the CO2 insertion product. This investigation elucidates the special role of Al in these hetero-binuclear compounds, providing new insights into the peculiar electronic structure of aluminyls, which may help for the rational control of their unprecedented reactivity toward carbon dioxide.

A Comprehensive Analysis of the Metal-Nitrile Bonding in an Organo-Diiron System

Nitriles (N≡CR) are ubiquitous in coordination chemistry, yet literature studies on metal-nitrile bonding based on a multi-technique approach are rare. We selected an easily-available di-organoiron framework, containing both π-acceptor (CO, aminocarbyne) and donor (Cp = η⁵-C₅H₅) ligands, as a suitable system to provide a comprehensive description of the iron-nitrile bond. Thus, the new nitrile (2-12)CF₃SO₃ and the related imine/amine complexes (8-9)CF₃SO₃ were synthesized in 58-83% yields from the respective tris-carbonyl precursors (1a-d)CF₃SO₃, using the TMNO strategy (TMNO = trimethylamine-N-oxide). The products were fully characterized by elemental analysis, IR (solution and solid state) and multinuclear NMR spectroscopy. In addition, the structures of (2)CF₃SO₃, (3)CF₃SO₃, (5)CF₃SO₃ and (11)CF₃SO₃ were ascertained by single crystal X-ray diffraction. Salient spectroscopic data of the nitrile complexes are coherent with the scale of electron-donor power of the R substituents; otherwise, this scale does not match the degree of Fe → N π-back-donation and the Fe-N bond energies, which were elucidated in (2-7)CF₃SO₃ by DFT calculations.

Charge-Flow Profiles along Curvilinear Paths: A Flexible Scheme for the Analysis of Charge Displacement upon Intermolecular Interactions

The Charge-Displacement (CD) analysis has proven to be a powerful tool for a quantitative characterization of the electron-density flow occurring upon chemical bonding along a suitably chosen interaction axis. In several classes of interesting intermolecular interactions, however, an interaction axis cannot be straightforwardly defined, and the CD analysis loses consistency and usefulness. In this article, we propose a general, flexible reformulation of the CD analysis capable of providing a quantitative view of the charge displacement along custom curvilinear paths. The new scheme naturally reduces to ordinary CD analysis if the path is chosen to be a straight line. An implementation based on a discrete sampling of the electron densities and a Voronoi space partitioning is described and shown in action on two test cases of a metal-carbonyl and a pyridine-ammonia complex.

Reactivity of a Gold-Aluminyl Complex with Carbon Dioxide: A Nucleophilic Gold?

A gold-aluminyl complex has been recently reported to feature an unconventional gold nucleophilic center, which was revealed through reactivity with carbon dioxide leading to the Au-CO₂ coordination mode. In this work, we computationally investigate the reaction mechanism, which is found to be cooperative, with the gold-aluminum bond being the actual nucleophile and Al also behaving as electrophile. The Au-Al bond is shown to be mainly of an electron-sharing nature, with the two metal fragments displaying a diradical-like reactivity with CO₂.

Assessing the Orbital Contribution in the "Spodium Bond" by Natural Orbital for Chemical Valence-Charge Displacement Analysis

The term "spodium bond" (SpB) has been recently proposed to describe the noncoordinative interaction that can be established between a polarized group 12 metal and a mild Lewis base (LB). Most of the systems showing short metal-donor distances compatible with SpB are characterized by the coexistence of multiple weak interactions, including hydrogen and halogen bonding, making the assessment of real importance of SpB difficult. Here, we show that the relative importance of each contribution can be probed by dissecting the orbital component of the interaction through the extended transition state-natural orbital for chemical valence-charge displacement analysis (ETS-NOCV-CD). The latter gives useful information about relative energies and electrons involved, for model systems ([(thiourea)₂MX₂]···LB; M = Zn, Cd, and Hg; X = Cl and I; and LB = CH₂S, CH₂O, CH₃CN, and CO) and a variety of structures extracted from experimentally characterized adducts, allowing us to demonstrate the lack of a direct correlation between a favorable metal-base distance and the presence of an orbital contribution for the SpB.

Spectroscopic/Bond Property Relationship in Group 11 Dihydrides via Relativistic Four-Component Methods

Group 11 dihydrides MH2- (M = Cu, Ag, Au, Rg) have been much less studied than the corresponding MH compounds, despite having potentially several interesting applications in chemical research. In this work, their main spectroscopic constants (bond lengths, dissociation energies, and force constants) have been evaluated by means of highly accurate relativistic four-component coupled cluster (4c-CCSD(T)) calculations in combination with large basis sets. Periodic trends have been quantitatively explained by the charge-displacement/natural orbitals for chemical valence (CD-NOCV) analysis based on the four-component relativistic Dirac-Kohn-Sham method, which allows a consistent picture of the nature of the M-H bond to be obtained on going down the periodic table in terms of Dewar-Chatt-Duncanson bonding components. A strong ligand-to-metal donation drives the M-H bond and it is responsible for the heterolytic (HM···H-) dissociation energies to increase monotonically from Cu to Rg, with RgH2- showing the strongest and most covalent M-H bond. The "V"-shaped trend observed for the bond lengths, dissociation energies, and stretching frequencies can be explained in terms of relativistic effects and, in particular, of the relativistically enhanced sd hybridization occurring at the metal, which affects the metal-ligand distances in heavy transition-metal complexes. The sd hybridization is very small for Cu and Ag, whereas it becomes increasingly important for Au and Rg, being responsible for the increasing covalent character of the bond, the sizable contraction of the Au-H and Rg-H bonds, and the observed trend. This work rationalizes the spectroscopic/bond property relationship in group 11 dihydrides within highly accurate relativistic quantum chemistry methods, paving the way for their applications in chemical bond investigations involving heavy and superheavy elements.

The role of the halogen bond in iodothyronine deiodinase: Dependence on chalcogen substitution in naphthyl-based mimetics

The effects on the activity of thyroxine (T4) due to the chalcogen replacement in a series of peri-substituted naphthalenes mimicking the catalytic function of deiodinase enzymes are computationally examined using density functional theory. In particular, T4 inner-ring deiodination pathways assisted by naphthyl-based models bearing two tellurols and a tellurol-thiol pair in peri-position are explored and compared with the analogous energy profiles for the naphthalene mimic having two selenols. The presence of a halogen bond (XB) in the intermediate formed in the first step and involved in the rate-determining step of the reaction is assumed to facilitate the process increasing the rate of the reaction. The rate-determining step calculated energy barrier heights allow rationalizing the experimentally observed superior catalytic activity of tellurium containing mimics. Charge displacement analysis is used to ascertain the presence and the role of the electron density charge transfer occurring in the rate-determining step of the reaction, suggesting the incipient formation or presence of a XB interaction. © 2019 Wiley Periodicals, Inc.

Disentanglement of orthogonal hydrogen and halogen bonds via natural orbital for chemical valence: A charge displacement analysis

As known, the electron density of covalently bound halogen atoms is anisotropically distributed, making them potentially able to establish many weak interactions, acting at the same time as halogen bond donors and hydrogen bond acceptors. Indeed, there are many examples in which the halogen and hydrogen bond coexist in the same structure and, if a correct bond analysis is required, their separation is mandatory. Here, the advantages and limitations of coupling the charge displacement analysis with natural orbital for chemical valence method (NOCV-CD) to separately analyze orthogonal weak interactions are shown, for both symmetric and asymmetric adducts. The methodology gives optimal results with intermolecular adducts but, in the presence of an organometallic complex, also intramolecular interactions can be correctly analyzed. Beyond the methodological aspects, it is shown that correctly separate and quantify the interactions can give interesting chemical insights about the systems.

Orbital Decomposition of the Carbon Chemical Shielding Tensor in Gold(I) N-Heterocyclic Carbene Complexes

The good performance of N‐heterocyclic carbenes (NHCs), in terms of versatility and selectivity, has called the attention of experimentalists and theoreticians attempting to understand their electronic properties. Analyses of the Au(I)–C bond in [(NHC)AuL]^+/0 (L stands for a neutral or negatively charged ligand), through the Dewar–Chatt–Duncanson model and the charge displacement function, have revealed that NHC is not purely a σ‐donor but may have a significant π‐acceptor character. It turns out, however, that only the σ‐donation bonding component strongly correlates with one specific component of the chemical shielding tensor. Here, in extension to earlier works, a current density analysis, based on the continuous transformation of the current density diamagnetic zero approach, along a series of [(NHC)AuL]^+/0 complexes is presented. The shielding tensor is decomposed into orbital contributions using symmetry considerations together with a spectral analysis in terms of occupied to virtual orbital transitions. Analysis of the orbital transitions shows that the induced current density is largely influenced by rotational transitions. The orbital decomposition of the shielding tensor leads to a deeper understanding of the ligand effect on the magnetic response properties and the electronic structure of (NHC)‐Au fragments. Such an orbital decomposition scheme may be extended to other magnetic properties and/or substrate‐metal complexes.

Chemical promenades: Exploring potential-energy surfaces with immersive virtual reality

The virtual-reality framework AVATAR (Advanced Virtual Approach to Topological Analysis of Reactivity) for the immersive exploration of potential-energy landscapes is presented. AVATAR is based on modern consumer-grade virtual-reality technology and builds on two key concepts: (a) the reduction of the dimensionality of the potential-energy surface to two process-tailored, physically meaningful generalized coordinates, and (b) the analogy between the evolution of a chemical process and a pathway through valleys (potential wells) and mountain passes (saddle points) of the associated potential energy landscape. Examples including the discovery of competitive reaction paths in simple A + BC collisional systems and the interconversion between conformers in ring-puckering motions of flexible rings highlight the innovation potential that augmented and virtual reality convey for teaching, training, and supporting research in chemistry.

Spin-orbit coupling is the key to unraveling intriguing features of the halogen bond involving astatine

The effect of spin-orbit coupling (SOC) on the halogen bond involving astatine has been investigated using state-of-the-art two- and four-component relativistic calculations. Adducts between Cl-X (X = Cl, Br, I and At) and ammonia have been selected to establish a trend on going down the periodic table. The SOC influence has been explored not only on the geometric and energetic features that can be used to characterize the halogen bond strength but also on the three main contributions to it that are the charge transfer, the "σ-hole" (i.e. the localized region with a net positive electrostatic potential at the halogen site) and the "polar flattening" (which is related to the effective shape of the halogen site). A surprisingly large increase of the Cl-At dipole moment, due to the inclusion of SOC, has been worked out using four-component CCSD(T) reference calculations, indicating that this bond is significantly more ionic than one may predict. Due to the SOC effect, which induces a peculiar charge accumulation on the At side in the Cl-At dimer, a weakening of the astatine-mediated halogen bond occurs arising from the (i) reduced amount of charge transfer, (ii) decrease of the polar flattening and (iii) lowering of the short-range Coulomb potential. The analysis of the electronic structure of the Cl-At moiety allows for a rationalization of the SOC effects on all the considered features of the halogen bond, including an unprecedented unsymmetrical charge back-donation from Cl-At to ammonia.

Theoretical study of glycine amino acid adsorption on graphene oxide

The non-dissociative and dissociative adsorptions of zwitterionic Gly on graphene oxide (GO) was studied in the framework of DFT using a cluster model approach. In this work, the interaction with an epoxy group of GO basal plane was mainly considered. As a comparison, the non-dissociative and dissociative adsorptions of neutral Gly were also taken into account. The non-dissociative adsorption modes for zwitterionic and neutral Gly conformers show binding energies of 12.2 and 14.4 kcal mol^-1, respectively. These molecules are thought to remain over the GO surface due to attractive noncovalent interactions. Two dissociative adsorption modes, for Z-Gly and N-Gly, show smaller binding energies of 7.2 and 8.4 kcal mol^-1, where the deprotonated species links strongly through a C-O or C-N covalent bond to the GO surface. The results obtained in the present theoretical approach to the glycine/graphene oxide system support the fact that glycine can be attached to epoxy groups of graphene oxide basal planes in addition to the anchoring on edge oxidation groups. In summary, we conclude that glycine can be used as a reducing agent as well as a functionalizer of GO sheets.

Charge Displacement Analysis-A Tool to Theoretically Characterize the Charge Transfer Contribution of Halogen Bonds

Theoretical bonding analysis is of prime importance for the deep understanding of the various chemical interactions, covalent or not. Among the various methods that have been developed in the last decades, the analysis of the Charge Displacement function (CD) demonstrated to be useful to reveal the charge transfer effects in many contexts, from weak hydrogen bonds, to the characterization of σ hole interactions, as halogen, chalcogen and pnictogen bonding or even in the decomposition of the metal-ligand bond. Quite often, the CD analysis has also been coupled with experimental techniques, in order to give a complete description of the system under study. In this review, we focus on the use of CD analysis on halogen bonded systems, describing the most relevant literature examples about gas phase and condensed phase systems. Chemical insights will be drawn about the nature of halogen bond, its cooperativity and its influence on metal-ligand bond components.

Chemical bonding in cuprous complexes with simple nitriles: octet rule and resonance concepts versus quantitative charge-redistribution analysis

Resonance structures for six cuprous complexes with simple nitriles are interpreted by means of a quantitative analysis of charge redistribution upon copper-nitrile bonding.

Rhenium(i) complexation–dissociation strategy for synthesising fluorine-18 labelled pyridine bidentate radiotracers

A novel fluorine-18 radiolabelling method employing rhenium(i) mediation is described herein. In less than 1 minute, fluorine-18 labelled complexes and ligands were synthesised in greater than 80% and 60% radiochemical yields (RCY), respectively.

Influence of halogen bonding on gold(i)-ligand bond components and DFT characterization of a gold-iodine halogen bond

A gold(i) complex bearing nitrogen acyclic carbene (NAC) and selenourea (SeU) has been used to verify whether the second-sphere SeI halogen bond (XB) is able to modify the Dewar-Chatt-Duncanson components of the Au-C and Au-Se bonds. The chosen system was found to be thermically unstable but it allowed an in-depth theoretical study by means of Energy Decomposition Analysis, Natural Bond Orbital and Natural Orbitals for Chemical Valence methods, coupled with Charge Displacement analysis. Indeed, in the presence of iodoperfluoroalkanes as XB donors, iodine interacts with the lone pair of the coordinated selenium, enhancing the Au ← C σ donation and depressing the Au → C π back-donation, as demonstrated also by the increase of the rotational barrier of the C-N bond of the NAC (see G. Ciancaleoni and others, Chem. - Eur. J., 2015, 21, 2467). On the other hand, in the presence of N-iodosuccinimide (NIS), the gold directly establishes a XB with the iodine by using its d lone pairs. This AuI XB is favored by the low steric hindrance of the ligands coordinated to the gold and the presence of the amino protons of SeU, which establish additional hydrogen bonds with the NIS. Also in this case, the effect is to increase the σ acidity and decrease the π basicity of the metal.

Cationic Gold(I) Diarylallenylidene Complexes: Bonding Features and Ligand Effects

Using computational approaches, we qualitatively and quantitatively assess the bonding components of a series of experimentally characterized Au(I) diarylallenylidene complexes (N.Kim, R.A.Widenhoefer, Angew. Chem. Int. Ed. 2018, 57, 4722-4726). Our results clearly demonstrate that Au(I) engages only weakly in π-backbonding, which is, however, a tunable bonding component. Computationally identified trends in bonding are clearly correlated with the substitution patterns of the aryl substituents in the Au(I) diarylallenylidene complexes and good agreement is found with the previously reported experimental data, such as IR spectra, 13 C NMR chemical shifts and rates of decomposition together with their corresponding barrier heights, further substantiating the computational findings. The description of the bonding patterns in these complexes allow predictions of their spectroscopic features, their reactivity and stability.

Alkyne Activation with Gold(III) Complexes: A Quantitative Assessment of the Ligand Effect by Charge-Displacement Analysis

A quantitative assessment of the Dewar-Chatt-Duncanson components of the Au(III)-alkyne bond in a series of cationic and dicationic bis- and monocyclometalated gold(III) complexes with 2-butyne via charge-displacement (CD) analysis is reported. Bonding between Au(III) and 2-butyne invariably shows a dominant σ donation component, a smaller, but significant, π back-donation, and a remarkable polarization of the alkyne CC triple bond toward the metal fragment. A very large net electron charge transfer from CC triple bond to the metal fragment results, which turns out to be unexpectedly insensitive to the charge of the complex and more strictly related to the nature of the ancillary ligand. The combination of σ donation, π back-donation, and polarization effects is in fact modulated by the different ligand frameworks, with ligands bearing atoms different from carbon in trans position with respect to the alkyne emerging as especially interesting for both imparting Au(III)-alkyne bond stability and inducing a more effective alkyne activation. A first attempt to figure out a rationale on the bonding/reactivity relationship for Au(III)-alkyne is made by performing a comparative study in a model nucleophilic attack of water to the alkyne triple bond. Smaller π back-donation facilitates alkyne slippage in the transition states, which is energetically less demanding for Au(III) than for Au(I), and suggests a greater propensity of Au(III) to facilitate the nucleophilic attack.

Interplay between Gold(I)-Ligand Bond Components and Hydrogen Bonding: A Combined Experimental/Computational Study

The influence of weak interactions on the donation/back-donation bond components in the complex [(NHC)Au(SeU)]⁺ (NHC = N-heterocyclic carbene; SeU = selenourea) has been studied by coupling experimental and theoretical techniques. In particular, NMR ¹H and pulsed-field gradient spin-echo titrations allowed us to characterize the hydrogen bond (HB) between the -NH₂ moieties of SeU and the anions PF₆ ^- and ClO₄ ^-, whereas ⁷⁷Se NMR spectroscopy allowed us to characterize the Au-Se bond. Theoretically, the Au-Se and Au-C orbital interactions have been decomposed using the natural orbital for the chemical valence framework and the bond components quantified through the charge displacement analysis. This methodology provides the quantification of the Dewar-Chatt-Duncanson (DCD) components for the Au-C and Au-Se bonds in the absence and presence of the second-sphere HB. The results presented here show that the anion has a dual mode action: it modifies the conformation of the cation by ion pairing (and this already influences the DCD components) and it induces new polarization effects that depend on the relative anion/cation relative orientation. The perchlorate polarizes SeU, enhancing the Se → Au σ donation and the Au → C back-donation and depressing the C → Au σ donation. On the contrary, the hexafluorophosphate depresses both the Se → Au and C → Au σ donations.

Ferrocenes with simple chiral substituents: an in-depth theoretical and experimental VCD and ECD study

Ferrocenes bearing chiral pendants are investigated through VCD and ECD. The VCD spectra are best interpreted by GVPT2-anharmonic DFT calculations. Diagnostic bands related to the absolute configuration of the title compounds are found in both kinds of spectra.

Effect of Rhenium(I) Complexation on Aza-Michael Additions to 5-Amino-1,10-Phenanthroline with [18F]Ethenesulfonyl Fluoride towards PET Optical Tracer Development

Conjugations with the recently developed [18F]ethenesulfonyl fluoride ([18F]ESF) were performed on 5-amino-1,10-phenanthroline, in its free form and coordinated to a rhenium(i) tricarbonyl complex, as a means of radiosynthesizing dual-modal optical and positron emission tomography (PET) tracers. The Michael-donating ability of the aromatic amine was noticeably perturbed on coordination with the rhenium(i) centre, resulting in decreased radiochemical yields from 34%, in the case of the free ligand, to 1%. We attribute the decreased nucleophilicity of the amine to metal deactivation from the electron-withdrawing feature of the rhenium(i) tricarbonyl centre, based on spectroscopic and computational evidence, thus highlighting this effect as a crucial parameter in designing late-stage metal coordination methods employing related aza-Michael additions. Photophysical analyses were also performed on the ESF-conjugated rhenium(i) complex, exhibiting a longer decay lifetime from the triplet metal-to-ligand charge transfer excited state when compared with the non-conjugated analogue.

Effect of Electron Correlation on Intermolecular Interactions: A Pair Natural Orbitals Coupled Cluster Based Local Energy Decomposition Study

The development of post-Hartree-Fock (post-HF) energy decomposition schemes that are able to decompose the HF and correlation components of the interaction energy into chemically meaningful contributions is a very active field of research. One of the challenges is to provide a clear-cut quantification to the elusive London dispersion component of the intermolecular interaction. London dispersion is well-known to be a pure correlation effect, and as such it is not properly described by mean field theories. In this context, we have recently developed the local energy decomposition (LED) analysis, which provides a chemically meaningful decomposition of the interaction energy between two or more fragments computed at the domain-based local pair natural orbitals coupled cluster (DLPNO-CCSD(T)) level of theory. In this work, this scheme is used in conjunction with other interpretation tools to study a series of molecular adducts held together by intermolecular interactions of different natures. The HF and correlation components of the interaction energy are thus decomposed into a series of chemically meaningful contributions. Emphasis is placed on discussing the physical effects associated with the inclusion of electron correlation. It is found that four distinct physical effects can contribute to the magnitude of the correlation part of intermolecular binding energies (Δ E_int^C): (i) London dispersion, (ii) the correlation correction to the reference induction energy, (iii) the correlation correction to the electron sharing process, and (iv) the correlation correction to the permanent electrostatics. As expected, the largest contribution to the correlation binding energy of neutral, apolar molecules is London dispersion, as in the argon dimer case. In contrast, the correction for the HF induction energy dominates Δ E_int^C in systems in which an apolar molecule interacts with charged or strongly polar species, as in Ar-Li⁺. This effect has its origin in the systematic underestimation of polarizabilities at the HF level of theory. For similar reasons, electron sharing largely contributes to the correlation binding energy of covalently bound molecules, as in the beryllium dimer case. Finally, the correction for HF permanent electrostatics significantly contributes to Δ E_int^C in molecules with strong dipoles, such as water and hydrogen fluoride dimers. This effect originates from the characteristic overestimation of dipole moments at the HF level of theory, leading in some cases to positive Δ E_int^C values. Our results are apparently in contrast to the widely accepted view that Δ E_int^C is typically dominated by London dispersion, at least, in the strongly interacting region. Clearly, post-HF energy decomposition schemes are very powerful tools to analyze, categorize, and understand the various contributions to the intermolecular interaction energy. Hopefully, this will eventually lead to insights that are helpful in designing systems with tailored properties. All analysis tools presented in this work will be available free of charge in the next release of the ORCA program package.

Ubiquity of cis-Halide → Isocyanide Direct Interligand Interaction in Organometallic Complexes

We recently reported a density functional theory (DFT) analysis of the Nb(V)-C bond in various NbCl₅(L) complexes, discovering that the carbon ligand L receives electronic density from the metal (classical back-donation) and from the chlorides in the cis position (direct interligand interaction). Here we report the synthesis and the structural characterization of two new coordination compounds of niobium pentahalides, i.e., NbX₅(CNXyl) (X = Cl, Br; Xyl = 2,6-C₆H₃Me₂), and the corresponding DFT analyses of the Nb(V)-C bond using the Natural Orbitals for Chemical Valence-Charge Displacement (NOCV-CD) approach, confirming the presence of a cis-halide → isocyanide direct interligand interaction. To verify whether the latter is limited to Nb complexes or not, we performed a NOCV-CD analysis on a series of several organometallic complexes based on Ti(IV), Nb(V), Ta(V), Rh(III), Pd(II), and Au(III), all of which bear one halide ligand and m-xylyl-isocyanide in a mutual cis position, revealing that the cis-halide → isocyanide interaction is always present.

Lewis Base Activation of Lewis Acid: A Detailed Bond Analysis

The effect of a Lewis base (LB) on the nucleophilic attack on chalcogeniranium (chalcogen = sulfur and selenium) cations, the so-called LB activation of a Lewis acid, has been studied coupling natural orbital for chemical valence decomposition of the orbital interaction energy with charge displacement analysis. This methodology provides a detailed and accurate description of all the interactions (LB···chalcogen, chalcogen···olefin and olefin···ammonia) present in the system and leads to a deeper understanding of how they influence each other at all stages of the reaction: reactant complex, transition state, and product complex. In particular, the bond between the chalcogen and the olefin has been decomposed in terms of σ donation/π back-donation and the bond components quantified. This allowed determination of a linear relationship between the activation barrier of the nucleophilic attack and the net amount of charge donated by the olefin to the chalcogen.

Unveiling the Sulfur-Sulfur Bridge: Accurate Structural and Energetic Characterization of a Homochalcogen Intermolecular Bond

By combining rotational spectroscopy in supersonic expansion with the capability of state-of-the-art quantum-chemical computations in accurately determining structural and energetic properties, the genuine nature of a sulfur-sulfur chalcogen bond between dimethyl sulfide and sulfur dioxide has been unveiled in a gas-jet environment free from collision, solvent and matrix perturbations. A SAPT analysis pointed out that electrostatic S⋅⋅⋅S interactions play the dominant role in determining the stability of the complex, largely overcoming dispersion and C-H⋅⋅⋅O hydrogen-bond contributions. Indeed, in agreement with the analysis of the quadrupole-coupling constants and of the methyl internal rotation barrier, the NBO and NOCV/CD approaches show a marked charge transfer between the sulfur atoms. Based on the assignment of the rotational spectra for 7 isotopologues, an accurate semi-experimental equilibrium structure for the heavy-atom backbone of the molecular complex has been determined, which is characterized by a S⋅⋅⋅S distance (2.947(3) Å) well below the sum of van der Waals radii.

Perspective: Found in translation: Quantum chemical tools for grasping non-covalent interactions

Today's quantum chemistry methods are extremely powerful but rely upon complex quantities such as the massively multidimensional wavefunction or even the simpler electron density. Consequently, chemical insight and a chemist's intuition are often lost in this complexity leaving the results obtained difficult to rationalize. To handle this overabundance of information, computational chemists have developed tools and methodologies that assist in composing a more intuitive picture that permits better understanding of the intricacies of chemical behavior. In particular, the fundamental comprehension of phenomena governed by non-covalent interactions is not easily achieved in terms of either the total wavefunction or the total electron density, but can be accomplished using more informative quantities. This perspective provides an overview of these tools and methods that have been specifically developed or used to analyze, identify, quantify, and visualize non-covalent interactions. These include the quantitative energy decomposition analysis schemes and the more qualitative class of approaches such as the Non-covalent Interaction index, the Density Overlap Region Indicator, or quantum theory of atoms in molecules. Aside from the enhanced knowledge gained from these schemes, their strengths, limitations, as well as a roadmap for expanding their capabilities are emphasized.

Ligand Effect on Bonding in Gold(III) Carbonyl Complexes

We quantitatively assess the Dewar-Chatt-Duncanson (DCD) components of the Au(III)-CO bond and the charge density polarization at the CO, in a series of neutral, cationic, and dicationic bis- and monocyclometalated gold(III) complexes via charge-displacement (CD) analysis. A striking feature concerns the very small net electron charge flux from CO to the metal fragment which is unexpectedly stable toward both the charge of the complex and the oxidation state of gold (I, III). All systems exhibit a similar trend for the σ charge rearrangement in the region of the carbonyl bond, where, by contrast, the π back-donation trend variation is large, which is strictly correlated to the change in CO bond distance and the shift in CO stretching frequencies, in close analogy with the gold(I) carbonyl complexes. In the whole series of gold(III) compounds, a large Au(III) ← CO σ donation is measured (from 0.19 to 0.31 electrons), as well as a significant Au(III) → CO π back-donation (from -0.09 to -0.22 electrons), which however is not generally able to completely balance the polarization of the CO π electrons in the direction from oxygen to carbon (C ← O) induced by the presence of the metal fragment [LAu(III)]^0/+1/+2. Surprisingly, all the gold(III) complexes in the series are characterized by a very small anisotropy in the Au(III) → CO in-plane and out-of-plane π back-donation components, in sharp contrast with the marked anisotropy found before for the experimentally characterized [(C^N^C)Au(III)CO]⁺ complex. A first attempt to figure out a rationale on the bonding/reactivity relationship for Au(III)-CO is made by performing a comparative study with an isostructural [(N^N^C)Pt(II)CO]⁺ complex in a model water-gas shift (WGS) reaction.

Exploiting coordination geometry to selectively predict the σ-donor and π-acceptor abilities of ligands: a back-and-forth journey between electronic properties and spectroscopy

Coordination geometry switches the carbonyl stretching frequency into a selective probe of the σ-donor and π-acceptor abilities of ligands.