The Chemical Bond and s–d Hybridization in Coinage Metal(I) Cyanides

Metal-ligand bond in group-11 complexes and nanoclusters

Density functional theory is used to study geometric, energetic, and electronic properties of metal-ligand bonds in a series of group-11 metal complexes and ligand-protected metal clusters. We study complexes as the forms of M-L (L = SCH3, SC8H9, PPh3, NHCMe, NHCEt, NHCiPr, NHCBn, CCMe, CCPh) and L1-M-L2 (L1 = NHCBn, PPh3, and L2 = CCPh). Furthermore, we study clusters denoted as [M13L6Br6]- (L = PPh3, NHCMe, NHCEt, NHCiPr, NHCBn). The systems were studied at the standard GGA level using the PBE functional and including vdW corrections via BEEF-vdW. Generally, Au has the highest binding energies, followed by Cu and Ag. PBE and BEEF-vdW functionals show the order Ag-L > Au-L > Cu-L for bond lengths in both M-L complexes and metal clusters. In clusters, the smallest side group (CH3) in NHCs leads to the largest binding energy whereas no significant variations are seen concerning different side groups of NHC in M-L complexes. By analyzing the projected density of states and molecular orbitals in complexes and clusters, the M-thiolate bonds were shown to have σ and π bond characteristics whereas phosphines and carbenes were creating σ bonds to the transition metals. Interestingly, this analysis revealed divergent behavior for M-alkynyl complexes: while the CCMe group displayed both σ and π bonding features, the CCPh ligand was found to possess only σ bond properties in direct head-to-head binding configuration. Moreover, synergetic effects increase the average binding strength to the metal atom significantly in complexes of two different ligands and underline the potential of adding Cu to synthesize structurally richer cluster systems. This study helps in understanding the effects of different ligands on the stability of M-L complexes and clusters and suggests that PPh3 and NHCs-protected Cu clusters are most stable after Au clusters.

Impact of ligand fields on Kubas interaction of open copper sites in MOFs with hydrogen molecules: an electronic structural insight

We investigate hydrogen sorption on open copper sites in various ligand coordinations of metal-organic frameworks (MOFs), including the triangular T(CuL3) in MFU-4l, the linear L(CuL2) in NU2100, and the paddlewheel P(CuL4)2 in HKUST-1 from an electronic structure perspective using DFT calculations. The ligand-field-induced splitting of d states and spd hybridizations in copper are thoroughly examined. The hybridization between Cu s, p, and d orbitals occurs in various forms to optimize the Coulomb repulsion of different ligand fields. Despite the Cu+ oxidation state, which is typically conducive to strong Kubas interactions with hydrogen molecules, the vacant spd z 2 hybrid orbitals of the open copper site in the L(CuL2) coordination are unsuitable for facilitating electron forward donation, thereby preventing effective hydrogen adsorption. In contrast, the vacant spd z 2 hybrid orbitals in the T(CuL3) and P(CuL4)2 coordinations can engage in electron forward donations, forming bonding states between the Cu spd z 2 and H2 σ bonding orbitals. The forward donation in the T(CuL3) configuration is significantly stronger than in the P(CuL4)2 configuration due to both the lower energy of the vacant orbitals and the larger contributions of p and d z 2 characters to the hybrid orbital. Additionally, the occupied Cu pd xz/yz and pd x 2-y 2 hybrid orbitals in the T(CuL3) configuration promote electron back donation to the H2 σ* antibonding orbital, leading to the formation of π bonding states. In the P(CuL4)2 coordination, the repulsion from the electron density distributed over the surrounding ligands prevents the H2 molecule from approaching the copper center closely enough for the back donation to occur. The complete Kubas interaction, involving both forward and back electron donations, results in a large dihydrogen-copper binding energy of 37.6 kJ mol-1 in the T(CuL3) coordination. In contrast, the binding energy of 10.6 kJ mol-1 in the P(CuL4)2 coordination is primarily driven by electrostatic interactions with a minor contribution of the Kubas-like forward donation interaction. This analysis highlights the pivotal role of coordination environments in determining the hydrogen sorption properties of MOFs.

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.

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.

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 antitumor coinage metal-based N-heterocyclic carbene complexes with cysteine and selenocysteine protein sites

The reaction of the antitumor M(I)-bis-N-heterocyclic carbene (M(I)-NHC) complexes, M = Cu, Ag, and Au, with their potential protein binding sites, i.e. cysteine and selenocysteine, was investigated by means of density functional theory approaches. Capped cysteine and selenocysteine were employed to better model the corresponding residues environment within peptide structures. By assuming the neutral or deprotonated form of the side chains of these amino acids and by considering the possible assistance of an external proton donor such as an adjacent acidic residue or the acidic component of the surrounding buffer environment, we devised five possible routes leading to the binding of the investigated M(I)-NHC scaffolds to these protein sites, reflecting their different location in the protein structure and exposure to the bulk. The targeting of either cysteine or selenocysteine in their neutral forms is a kinetically unfavored process, expected to be quite slow if observable at all at physiological temperature. On the other hand, the reaction with the deprotonated forms is much more favored, even though an external proton source is required to assist the protonation of the leaving carbene. Our calculations also show that all coinage metals are characterized by a similar reactivity toward the binding of cysteine and selenocysteine sites, although the Au(I) complex has significantly higher reaction and activation free energies compared to Cu(I) and Ag(I).

Oxidation of Silver Cyanide Ag(CN)2- by the OH Radical: From Ab Initio Calculation to Molecular Simulation and to Experiment

We investigate the oxidation of silver cyanide AgI(CN)2- in water by the OH radical in order to compare this complex with the free cation Ag⁺ and to measure the influence of the ligands. High-level ab initio calculations of the model species AgII(CN)2· enable the calibration of molecular simulations and the prediction of the oxidized species: AgII(CN)2(H2O)2· and its absorption spectrum, with an intense band at 292 nm and a weaker one at 390 nm. Pulse radiolysis measurements of the oxidation of AgI(CN)2- by the OH radical in water yields a transient species with a broad, intense band at 290 nm and a weaker band at 410 nm at short times after the pulse and a blue shift of the spectrum at longer times. The prediction of the simulations, that the oxidized complex AgII(CN)2(H2O)2· is formed, is confirmed by thermochemistry. Our calculations also suggest that the formation of the OH-adduct is possible only in very basic solution and that the blue shift observed at long times after the pulse is due to disproportionation of the oxidized complex. We also perform molecular simulations of the oxidation of free Ag⁺ cations by the OH radical. The results are compared to that of the literature and to the results obtained with the AgI(CN)2- complex.

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.

Excimer and Exciplex Formation in Gold(I) Complexes Preconditioned by Aurophilic Interactions

Excimers and exciplexes are defined as assemblies of atoms or molecules A/A' where interatomic/intermolecular bonding appears only in excited states such as [A2 ]* (for excimers) and [AA']* (for exciplexes). Their formation has become widely known because of their role in gas-phase laser technologies, but their significance in general chemistry terms has been given little attention. Recent investigations in gold chemistry have opened up a new field of excimer and exciplex chemistry that relies largely on the preorganization of gold(I) compounds (electronic configuration AuI (5d10 )) through aurophilic contacts. In the corresponding excimers, a new type of Au⋅⋅⋅Au bonding arises, with bond energies and lengths approaching those of ground-state Au-Au bonds between metal atoms in the Au0 (5d10 6s1 ) and AuII (5d9 ) configurations. Excimer formation gives rise to a broad range of photophysical effects, for which some of the relaxation dynamics have recently been clarified. Excimers have also been shown to play an important role in photoredox binuclear gold catalysis.

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.