Structure and Stability of M−CO, M = First-Transition-Row Metal: An Application of Density Functional Theory and Topological Approaches

The Scandium(II) Carbonyl Complex (C5H2tBu3)2Sc(CO) and Its Isocyanide Analog (C5H2tBu3)2Sc(CNC6H3Me2-2,6)

Treatment of the scandium(II) metallocene Cpttt2Sc (Cpttt = C5H2tBu3) with CO or the isocyanide CNXyl (Xyl = C6H3Me2-2,6) yields the carbonyl complex Cpttt2Sc(CO), 1, or the isocyanide complex Cpttt2Sc(CNXyl), 2, which were identified by X-ray crystallography. Isotopic labeling with 13CO shows the CO stretch of 1 at 1875 cm-1 shifts to 1838 cm-1 in 1-13CO. The CN stretch in 2 is shifted to 1939 cm-1 compared to 2118 cm-1 for the free isocyanide. The 80.1 MHz (28.7 G) 45Sc hyperfine coupling in 1 and 74.7 MHz (26.8 G) in 2 are similar to the 82.6 MHz (29.6 G) coupling constant in Cpttt2Sc and indicate that 1 and 2 are Sc(II) complexes. A comprehensive analysis of the electronic structures of 1 and 2 using DFT calculations is reported.

Microsolvation of cobalt, nickel, and copper atoms with ammonia: a theoretical study of the solvated electron precursors

CONTEXT

The s-block metals dissolved in ammonia form metal-ammonia complexes with diffuse electrons which could be used for redox catalysis. In this theoretical paper, we investigated the possibility of the d-bloc transition metals (Mn, Fe, Co, Ni, and Cu) solvated by ammonia. It has been demonstrated that both Mn and Fe atoms undergo into an oxidative reaction with NH3 forming an inserted species, HMNH2. On the contrary, the Co, Ni, and Cu atoms can accommodate four NH3, via the coordination bond, to form the first solvation sphere within C2v, D2d, and Td point groups, respectively. Addition of a fifth NH3 constitute the second solvation shell by forming hydrogen bond with the other NH3s. Interestingly, M(NH3)4 (M = Co, Ni, and Cu) is a so-called solvated electron precursor and should be considered as a monocation M(NH3)4+ kernel in tight contact with one electron distributed over its periphery. This nearly free electron could be used to capture a CO2 molecule and engages in a reduction reaction.

METHODS

Geometry optimization of the stationary points on the potential energy surface was performed using density functional theory - CAM-B3LYP functional including the GD3BJ dispersion contribution - in combination with the 6-311 +  + G(2d, 2p) basis set for all the atoms. All first-principles calculations were performed using the Gaussian 09 quantum chemical packages. The natural electron configuration of transition atom engaged in the compounds has been found using the natural bond orbital (NBO) method. We used the EDR (electron delocalization range) approach to analyze the structure of solvated electrons in real space. We also used the electron localization function (ELF) to measure the degree of electronic localization within a chemical compound. The EDR and ELF analyses are done using the TopMod and Multiwfn packages, respectively.

Spectroscopic Characterization of Highly Excited Neutral Chromium Tricarbonyl

Spectroscopic characterization of highly excited neutral transition-metal complexes is important for understanding the multifaceted reaction mechanisms between metals and ligands. In this work, the reactions of neutral chromium atoms with carbon monoxide were probed by size-specific infrared spectroscopy. Interestingly, Cr(CO)3 was found to have an unprecedented 7A2″ septet excited state rather than the singlet ground state. A combination of experiment and theory shows that the gas-phase formation of this highly excited Cr(CO)3 is facile both thermodynamically and kinetically. Electronic structure and bonding analyses indicate that the valence electrons of Cr atoms in the septet Cr(CO)3 are in a relatively stable configuration, which facilitate the highly excited structure and the planar geometric shape (D3h symmetry). The observed septet Cr(CO)3 affords a paradigm for exploring the structure, properties, and formation mechanism of a large variety of excited neutral compounds.

Exploring the Reactivity of Donor-Acceptor Systems through a Combined Conceptual and Constrained DFT Approach

In the context of the conceptual density functional theory (cDFT) and based on the computational efficiency of the constrained DFT (CDFT), we demonstrate that chemical reactivity can be governed by the difference between the local interacting chemical potentials of the reactants (referred as Edual), in agreement with Sanderson's equalization principle. In a proof-of-concept study, we investigated illustrative examples involving typical non-covalent donor-acceptor systems and reactive systems are provided. For the selected systems, our approach reveals significant mimicking between Edual and the DFT-computed intermolecular interaction energy profiles. We further evaluate the influence of the Coulomb and exchange-correlation contributions in Edual. These latter results suggest that numerous potential energy surfaces of clusters can be explored using a Sanderson-like model only based on classical interactions between molecular orbitals domains. To conclude, this study achieved a deeper understanding of the principles of cDFT and assessed, in a wider context, its efficiency in predicting the chemical reactivity.

Spectroscopic Signature of the Carbon-Carbon Coupling Reaction between Carbon Monoxide and Nickel Carbide

Spectroscopic characterization of ketenylidene complexes is of essential importance for understanding the structure-reactivity relationships of the catalytic sites. Here, we report a size-specific photoelectron velocity map imaging spectroscopic study of the reactions of carbon monoxide with nickel carbide. Quantum chemical calculations have been conducted to search for the energetically favorable isomers and to recognize the experimental spectra. The target products with the chemical formula of NiC(CO)n- (n = 3-5) are characterized to have an intriguing ketenylidene CCO unit. The evolution from NiC(CO)3- to NiC(CO)4- involves the breaking and formation of the Ni-C bond and the coordination conversion between the terminal and bridging carbonyls. Experimental and theoretical analyses reveal an efficient C-C bond formation process within the reactions of carbon monoxide and laser-vaporized nickel carbide. This work highlights the pivotal roles played by metal carbides in the C-C bond formation and also proposes new ideas for the design and chemical control of a broad class of complexes with unique physical and chemical properties.

Metal-CO Bonding in Mononuclear Transition Metal Carbonyl Complexes

DFT calculations have been carried out for coordinatively saturated neutral and charged carbonyl complexes [M(CO) _n ] ^q where M is a metal atom of groups 2-10. The model compounds M(CO)₂ (M = Ca, Sr, Ba) and the experimentally observed [Ba(CO)]⁺ were also studied. The bonding situation has been analyzed with a variety of charge and energy partitioning approaches. It is shown that the Dewar-Chatt-Duncanson model in terms of M ← CO σ-donation and M → CO π-backdonation is a valid approach to explain the M-CO bonds and the trend of the CO stretching frequencies. The carbonyl ligands of the neutral complexes carry a negative charge, and the polarity of the M-CO bonds increases for the less electronegative metals, which is particularly strong for the group 4 and group 2 atoms. The NBO method delivers an unrealistic charge distribution in the carbonyl complexes, while the AIM approach gives physically reasonable partial charges that are consistent with the EDA-NOCV calculations and with the trend of the C-O stretching frequencies. The AdNDP method provides delocalized MOs which are very useful models for the carbonyl complexes. Deep insight into the nature of the metal-CO bonds and quantitative information about the strength of the [M] ← (CO)₈ σ-donation and [M(d)] → (CO)₈ π-backdonation visualized by the deformation densities are provided by the EDA-NOCV method. The large polarity of the M-CO π orbitals toward the CO end in the alkaline earth octacarbonyls M(CO)₈ (M = Ca, Sr, Ba) leads to small values for the delocalization indices δ(M-C) and δ(M···O) and significant overlap between adjacent CO groups, but the origin of the charge migration and the associated red-shift of the C-O stretching frequencies is the [M(d)] → (CO)₈ π-backdonation. The heavier alkaline earth metals calcium, strontium and barium use their s/d valence orbitals for covalent bonding. They are therefore to be assigned to the transition metals.

Thermodynamics and Kinetics of Gas-Phase CO Oxidation on the Scandium Monoxide Carbonyl Complexes

The CO chemisorption onto the ScO⁺ cation was investigated using infrared photodissociation spectroscopy combined with density functional theory calculations. The spectra were recorded in the CO stretching vibrational region for the OSc(CO)_n⁺ (n = 4-6) complex series. Comparisons of the experimental spectra with the simulated ones have established the geometries and present strong evidence that all of the CO ligands are chemisorbed, which could not be readily oxidized by scandium monoxide core into CO₂. Complementary calculations demonstrate that, regardless of the thermodynamic feasibility, the CO oxidation on the scandium monoxide carbonyl complexes is kinetically unfavorable due to the significant barriers involved in the CO oxidation process. Nevertheless, the consecutive CO adsorption has a positive influence on the Sc-O bond activation.

Electronic structure calculations permit identification of the driving forces behind frequency shifts in transition metal monocarbonyls

Our direct DFT decomposition of CO frequency shifts updates the paradigm for metal carbonyl binding.

Carbon monoxide adsorption at forsterite surfaces as models of interstellar dust grains: An unexpected bathochromic (red) shift of the CO stretching frequency

Carbon monoxide (CO) is one of the most abundant species in the interstellar medium (ISM). In the colder regions of the ISM, it can directly adsorb onto exposed Mg cations of forsterite (Fo, Mg₂SiO₄), one of the main constituents of the dust grains. Its energetic of adsorption can strongly influence the chemico-physical evolution of cold interstellar clouds; thus, a detailed description of this process is desirable. We recently simulated the CO adsorption on crystalline Fo surfaces by computer ab initio methods and, surprisingly, reported cases where the CO stretching frequency underwent a bathochromic (red) shift (i.e., it is lowered with respect to the CO gas phase frequency), usually not experimentally observed for CO adsorbed onto oxides with non-d cations, like the present case. Here, we elucidate in deep when and under which conditions this case may happen and concluded that this red shift may be related to peculiar surface sites occurring at the morphologically complex Fo surfaces. The reasons for the red shift are linked to both the quadrupolar nature of the CO molecule and the role of dispersion interactions with surfaces of complex morphology. The present work, albeit speculative, suggests that, at variance with CO adsorption on simple oxides like MgO, the CO spectrum may exhibit features at lower frequencies than the reference gas frequency when CO is adsorbed on complex oxides, even in the absence of transition metal ions.

Revisiting π backbonding: the influence of d orbitals on metal–CO bonds and ligand red shifts

Stronger C–O bonds are found if metal d-functions are present in carbonyl complexes, contrary to the common understanding of π-backbonding.

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.

Understanding the CO Dissociation in [Fe(CN)2(CO)2(dithiolate)]2- Complexes with Quantum Chemical Topology Tools

The active site of the [NiFe]-hydrogenase contains a pentacoordinated iron atom; therefore, a vacant coordination site is available for substrate binding. Nonetheless, most organometallic models of the [NiFe]-hydrogenase failed to reproduce this key feature of the active site. In order to rationalize such behavior, quantum chemical calculations were carried out on a series of [Fe(CN)₂(CO)_n(dithiolate)]^2- n = 1,2 complexes, where dithiolate denotes the ligands (CF₃)₂C₂S₂^2-, (CO₂Me)₂C₂S₂^2-, Ph₂C₂S₂^2-, C₆Cl₂H₂S₂^2-, C₆H₄S₂^2-, C₂H₄S₂^2-, and C₃H₆S₂^2-. Structural and energetic features are discussed, and a topological analysis based on two scalar fields, the one-electron density and the electron localization function (ELF), has been attempted to describe the nature of the metal-ligand bonds. The present approach contributes to better understand the ability of noninnocent dithiolene to strongly labilize one CO whereas innocent dithiolate cannot. The methodology developed throughout the paper could be useful in the field of the CO-releasing molecules.

How π back-donation quantitatively controls the CO stretching response in classical and non-classical metal carbonyl complexes

The CO stretching response upon coordination to a metal M to form [(L) n M(CO)] m complexes (L is an auxiliary ligand) is investigated in relation to the σ donation and π back-donation components of the M-CO bond and to the electrostatic effect exerted by the ligand-metal fragment. Our analysis encompasses over 30 carbonyls, in which the relative importance of donation, back-donation and electrostatics are varied either through the ligand in a series of [(L)Au(CO)]0/+ gold(i) complexes, or through the metal in a series of anionic, neutral and cationic homoleptic carbonyls. Charge-displacement analysis is used to obtain well-defined, consistent measures of σ donation and π back-donation charges, as well as to quantify the σ and π components of CO polarization. It is found that all complexes feature a comparable charge flow of σ symmetry (both in the M-CO bonding region and in the CO fragment itself), which is therefore largely uncorrelated to CO response. By contrast, π back-donation is exceptionally variable and is found to correlate tightly with the change in CO bond distance, with the shift in CO stretching frequency, and with the extent and direction (C → O or C ← O) of the CO π polarization. As a result, we conclusively show that π back-donation can be an important bond component also in non-classical carbonyls and we provide the framework in which the spectroscopic data on coordinated CO can be used to extract quantitative information on the π donor properties of metal-ligand moieties.

Benchmark calculations of metal carbonyl cations: relativistic vs. electron correlation effects

In this paper we present benchmark results for isoelectronic metal carbonyl complexes of the groups 11 and 12 of the periodic table. The focus is on the geometry, vibrational frequencies, bond dissociation energy and chemical bonding. The description of these complexes requires a good balance between electron correlation and relativistic effects. Our results demonstrate that the combination of the effective core potential and the MP2 method gives quantitative results for the first- and the second-row transition metal complexes and only qualitative agreement for the third-row complexes. In order to obtain quantitative results for the whole series the use of four-component or X2C methods is mandatory. The fourth-row transition metal carbonyl complexes from groups 11 and 12 have been studied for the first time. The metal-carbon bond strength pattern along the group is shown to be highly dependent on the correct description of the relativistic effects. Finally, the relativistic effects on the bonding are studied by means of electron density difference maps, the analysis of the bond critical points of the electron density and the mechanism for σ-donation and π-backdonation. Our analysis indicates that the fourth-row complexes exhibit a strong covalent character induced by relativistic effects.

Comparative density functional theory and post-Hartree-Fock (CCSD, CASSCF) studies on the electronic structure of halogen nitrites ClONO and BrONO using quantum chemical topology

In this paper, the electronic structures of cis- and trans-ClONO and BrONO are studied at the CCSD∕aug-cc-pVTZ, CASSCF(14,12)/aug-cc-pVTZ, and B3LYP/aug-cc-pVTZ computational levels. For the Cl-O bond, topological analysis of the electron density field, ρ(r), shows the prevalence of the shared-electron type bond (∇(2)ρ((3,-1)) < 0). The Br-O bond, however, represents the closed-shell interaction (∇(2)ρ((3,-1)) > 0). Topological analysis of the electron localization function, η(r), and electron localizability indicator (ELI-D), (D) (σ)(r), shows that the electronic structure of the central N-O bond is very sensitive to both electron correlation improvements (coupled-cluster single double (CCSD), CASSCF, density functional theory (DFT)) and bond length alteration. Depending on the method used, the N-O bond can be characterized as a "normal" N-O bond with a disynaptic V(N,O) basin (DFT); a protocovalent N-O bond with two monosynaptic, V(N) and V(O), basins (CCSD, CASSCF); or a new type, first discovered for FONO, characterized by a single monosynaptic, V(N) basin (CCSD, DFT). The total basin population oscillates between 0.46-0.96 e (CCSD) and 0.86-1.02 e (CASSCF). The X-O bond is described by the single disynaptic basin, V(X,O), with a basin population between 0.76 and 0.81 e (CCSD) or 0.77 and 0.85 e (CASSCF). Analysis of the localized electron detector distribution for the cis-Cl-O1-N=O2 shows a manifold in the Cl···O2 region, associated with decreased electron density.

Nonlinear time-dependent density functional theory studies of the ionization of CO2 by ultrashort intense laser pulses

Time-dependent density functional theory (TDDFT) studies of the ionization of CO2 by intense laser pulses (3.50 × 1014, 1.40 × 1015, 2.99 × 1015, and 1.25 × 1016 W/cm2) at 800 nm (ω = 0.0584 au) are presented in the nonlinear nonpertubative regime. Special emphasis is placed on elucidating molecular orbital orientation and various peak-intensities effects on the ionization processes. The results reveal that molecular orbital ionizations are strongly sensitive to their symmetry and the laser intensities. Most notably, we found that with a proper choice of the laser intensity (3.5 × 1014 W/cm2), the sensitivity is strong enough such that the nature and symmetry of the highest occupied molecular orbital (HOMO) can be directly probed and visualized from the angular dependence of laser-induced ionization. At higher intensities, ionization is found to occur also from inner orbitals, thus complicating the imaging of simple orbitals. A time-dependent electron-localization function (TDELF) is used to get a visual insight on the time evolution process of the electron density.

Advancing beyond charge analysis using the electronic localization function: Chemically intuitive distribution of electrostatic moments

We propose here an evaluation of chemically intuitive distributed electrostatic moments using the topological analysis of the electron localization function (ELF). As this partition of the total charge density provides an accurate representation of the molecular dipole, the distributed electrostatic moments based on the ELF partition (DEMEP) allows computing of local moments located at non atomic centers such as lone pairs, sigma bonds and pi systems. As the local dipole contribution can be decomposed in polarization and charge transfer components, our results indicate that local dipolar polarization of the lone pairs and chemical reactivity are closely related whereas the charge transfer contribution is the key factor driving the local bond dipole. Results on relevant molecules show that local dipole contributions can be used to rationalize inductive polarization effects in alcohols derivatives and typical hydrogen bond interactions. Moreover, bond quadrupole polarization moments being related to a pi character enable to discuss bond multiplicities, and to sort families of molecules according to their bond order. That way, the nature of the C-O bond has been revisited for several typical systems by means of the DEMEP analysis which appears also helpful to discuss aromaticity. Special attention has been given to the carbon monoxide molecule, to the CuCO complex and to a weak intramolecular N|-CO interaction involved in several biological systems. In this latter case, it is confirmed that the bond formation is mainly linked to the CO bond polarization. Transferability tests show that the approach is suitable for the design of advanced force fields.

The performance of the new 6-31G## basis set: Molecular structures and vibrational frequencies of transition metal carbonyls

The performance of the newly proposed 6-31G(##) basis set for calculating the equilibrium structure and vibrational frequencies of transition metal carbonyl complexes has been studied at the HF and DFT levels of theory. The 6-31G(##) basis set has been constructed by augmentation of the 6-31G basis set by diffuse and polarization functions, which are generated from the corresponding 6-31G basis AOs response functions obtained in the frame of propagator approach. The predicted values of bond distances and vibrational frequencies for the title compounds are in good agreement with the experimental data. The relative energies and HOMO-LUMO gaps were also estimated for the series of MCO complexes.

Activation of Methane by the Iron Dimer Cation. A Theoretical Study

A detailed investigation of the reaction mechanisms underlying the observed reactivity of the iron dimer cation with respect to methane has been performed by density functional hybrid (B3LYP) and nonhybrid (BPW91) calculations. Minima and transition states have been fully optimized and characterized along the potential energy surfaces leading to three different exit channels for both the ground and the first excited states of the dimer. A comparison with our previous work covering the reactivity of the Fe(+) monomer was made to underline similarities and differences of the reaction mechanisms. Results show that geometric arrangements corresponding to bridged positions of the ligands with respect to iron atoms are always favored and stabilize intermediates, transition states and products, facilitating their formation. Binding energies of reaction products have been computed and compared with experimental measurements, and ELF analysis of the bond has been performed to rationalize trends as a function of the structure.

Intrinsic Dinitrogen Activation at Bare Metal Atoms

There is ongoing interest in metal complexes which bind dinitrogen and facilitate either its reduction or oxidation under mild conditions. In nature, the enzyme nitrogenase catalyzes this process, and dinitrogen fixation occurs under mild and ambient conditions at a metal-sulfur cluster in the center of the MoFe protein, but the mechanism of this process remains largely unknown. In the last few years, new important discoveries have been made in this field. In this review are discussed recent findings on the interaction of N(2) with metal atoms and metal-atom dimers from all groups of the periodic table as provided by gas-phase as well as matrix-isolation experiments. Intrinsic dinitrogen activation at such bare metal atoms is then related to corresponding processes at complexes, clusters, and surfaces.

Comments on the nature of the bonding in oxygenated dinuclear copper enzyme models

The nature of the bonding in model complexes of di-copper metalloenzymes has been analyzed by means of the electronic localization function (ELF) and by the quantum theory of atoms in molecules (QTAIM). The constrained space orbital variations (CSOV) approach has also been used. Density functional theory (DFT) and CASSCF calculations have been carried out on several models of tyrosinase such as the sole Cu2O22+ central core, the Cu2O2(NH3)62+ complex and the Cu2O2(Imidazol)62+ complex. The influence on the central Cu(2)O(2) moiety of both levels of calculation and ligand environment have been discussed. The distinct bonding modes have been characterized for the two major known structures: [Cu(2)(mu-eta(2): eta(2)-O(2))](2+) and [Cu(2)(mu-O(2))](2+). Particular attention has been given to the analysis of the O-O and Cu-O bonds and the nature of the bonding modes has also been analyzed in terms of mesomeric structures. The ELF topological approach shows a significant conservation of the topology between the DFT and CASSCF approaches. Particularly, three-center Cu-O-Cu bonds are observed when the ligands are attached to the central core. At the DFT level, the importance of self interaction effects are emphasized. Although, the DFT approach does not appear to be suitable for the computation of the electronic structure of the isolated Cu(2)O(2) central core, competitive self interaction mechanisms lead to an imperfect but acceptable model when using imidazol ligands. Our results confirm to a certain extent the observations of [M.F. Rode, H.J. Werner, Theoretical Chemistry Accounts 4-5 (2005) 247.] who found a qualitative agreement between B3LYP and localized MRCI calculations when dealing with the Cu(2)O(2) central core with six ammonia ligands.

Density Functional Study of the Interaction of Carbon Monoxide with Small Neutral and Charged Silver Clusters

CO adsorption on small neutral, anionic, and cationic silver clusters Ag(n) (n = 1-7) has been studied with use of the PW91PW91 density functional theory (DFT) method. The adsorption of CO on-top site, among various possible sites, is energetically preferred irrespective of the charge state of the silver cluster. The cationic silver clusters generally have a greater tendency to adsorb CO than the anionic and neutral silver ones, except for n = 3 and 4, and the binding energies reach a local minimum at n = 5. The binding energies on the neutral clusters, instead, reach a local maximum at n = 3, which is about 0.87 eV, probably large enough to be captured in the experiments. Binding of CO to the silver clusters is generally weaker than that to the copper and gold counterparts at the same size and charge state. This is due to the weaker orbital interaction between silver and CO, which is caused by the larger atomic radius of the silver atom. In contrast, Au atoms with a larger nuclear charge but a similar atomic radius to silver owing to the lanthanide contraction are able to have a stronger interaction with CO.

The Electronic Structure of Alkali Aurides. A Four-Component Dirac−Kohn−Sham Study

Spectroscopic constants, including dissociation energies, harmonic and anharmonic vibrational frequencies, and dipole moments, are calculated for the complete alkali auride series (LiAu, NaAu, KAu, RbAu, CsAu). The four-component formulation of relativistic density functional theory has been employed in this study, using the G-spinor basis sets implemented recently in the program BERTHA. The performance of four standard nonrelativistic density functionals employed is investigated by comparing the results with the best available theoretical and experimental data. The present work provides the first theoretical predictions on the molecular properties of RbAu. The intermetallic bond that occurs in the alkali auride series is highly polar and is characterized by a large charge transfer from the alkali metals to gold. The extent of this electron transfer has been investigated using several different charge analysis methods, enabling us to reach some general conclusions on their relative performance. We further report a detailed analysis of the topological properties of relativistic electron density in the bonding region, discussing the features of this approach which characterize the nature of the chemical bond. We have also computed the fully relativistic density for the alkali halides MBr and MI (M = Li, Na, K, Rb, and Cs). The comparative study shows that, on the basis of several topological properties and the variation in bond lengths, the gold atom behaves similarly to a halogen intermediate between Br and I.