Relativistic effects in gold chemistry. V. Group 11 dipole polarizabilities and weak bonding in monocarbonyl compounds

Diffuse Basis Functions for Relativistic s and d Block Gaussian Basis Sets

Diffuse s, p, and d functions have been optimized for use with previously reported relativistic basis sets for the s and d blocks of the periodic table. The functions were optimized on the 4:1 weighted average of the s² and p² configurations of the anion, with the d shell in the dⁿ⁺¹ configuration for the d blocks. Exponents were extrapolated for groups 2 and 12, which have unstable or weakly bound anions. The diffuse basis sets have been tested by application to calculations of electron affinities of the group 11 elements (Cu, Ag, and Au), double electron affinities of the group 11 monocations, and potential energy curves of Mg₂ and Ca₂ van der Waals dimers, as well as some response properties of the group 1 anions (Rb^-, Cs^-, and Fr^-), the group 2 elements (Sr, Ba, and Ra), and RbLi, CsLi, and FrLi molecules.

Anion Photoelectron Spectroscopy and Theoretical Study of HAuCN and [HAuCN]-: Spin-Orbit Coupling and Low-Lying Excited States

We report a study of the electronic structures and chemical bonding of gaseous [HAuCN]^- and the corresponding neutral molecule using photoelectron spectroscopy and relativistic quantum chemistry calculations. The electron affinity of the neutral HAuCN is reported to be 4.75 eV for the first time. The low-lying excited states of neutral molecule are observed and assigned according to the calculations utilizing a sophisticated electron correlation method incorporating both the scalar and spin-orbit relativistic effects. Our theoretical calculations suggest the geometry will be distorted from linear structure to the bent during the process of detaching one electron from the anion. Various chemical bonding analyses based on theoretical calculations have been performed for the titled complexes, and the apparent covalent natures of interactions between gold and the studied ligands have been verified.

Insight into induced charges at metal surfaces and biointerfaces using a polarizable Lennard-Jones potential

Metallic nanostructures have become popular for applications in therapeutics, catalysts, imaging, and gene delivery. Molecular dynamics simulations are gaining influence to predict nanostructure assembly and performance; however, instantaneous polarization effects due to induced charges in the free electron gas are not routinely included. Here we present a simple, compatible, and accurate polarizable potential for gold that consists of a Lennard–Jones potential and a harmonically coupled core-shell charge pair for every metal atom. The model reproduces the classical image potential of adsorbed ions as well as surface, bulk, and aqueous interfacial properties in excellent agreement with experiment. Induced charges affect the adsorption of ions onto gold surfaces in the gas phase at a strength similar to chemical bonds while ions and charged peptides in solution are influenced at a strength similar to intermolecular bonds. The proposed model can be applied to complex gold interfaces, electrode processes, and extended to other metals.

Molecular dynamics models for predicting the behavior of metallic nanostructures typically do not take into account polarization effects in metals. Here, the authors introduce a polarizable Lennard–Jones potential that provides quantitative insight into the role of induced charges at metal surfaces and related complex material interfaces.

Aun (n = 1,11) Clusters Interacting With Lone-Pair Ligands

We analyze the pattern of binding energies (BEs) of small Aun clusters (n = 1-7, 11) with lone-pair ligands (L = H2O, SH2, NH3, PH3, PF3, PCl3, and PMe3) employing the density functional theory. We use PBE0 functional with the dispersion correction and scalar relativistic effective core potential. This approach provides correct BEs when compared with benchmark CCSD(T) calculations for Au-L and Au2-L complexes. The pattern of BEs of Aun-L complexes is irregular with BE for Au3 ≈ Au4 > Au2 > Au7 > Au5 > Au11 > Au6 > Au1. Electron affinities (EAs) of Aun clusters exhibit oscillatory pattern with the cluster size. Binding energies of Aun-L complexes are oscillatory as well following EAs of Aun clusters. BEs of odd and even Aun-L complexes were analyzed separately. The bonding mechanism in odd Aun-L complexes is dominated by the lone pair → metal electron donation to the singly occupied valence Aun orbital accompanied by the back-donation. Even Aun clusters create covalent Aun-L bonds with BEs higher than those in odd Aun-L complexes. The BEs pattern and optimized geometries of Aun-L complexes correspond to the picture of creating the gold-ligand bond through the lone pair of a ligand interacting with the singly occupied molecular orbital in odd clusters or lowest unoccupied molecular orbital in even clusters of Aun. Ligands in both odd and even Aun-L complexes form three groups with binding energies that correlate with their ionization energies. The lowest BE is calculated for H2O as a ligand, followed by SH2 and NH3. PX3 ligands exhibit highest BEs.

On the gold-ligand covalency in linear [AuX2](-) complexes

Gold compounds, clusters, and nanoparticles are widely used as catalysts and therapeutic medicines; the interactions between gold and its ligands in these systems play important roles in their chemical properties and functionalities. In order to elucidate the nature of the chemical interactions between Au(I) and its ligands, herein we use several theoretical methods to study the chemical bonding in a variety of linear [AuX2](-) complexes, where X = halogen atoms (F, Cl, Br, I, At and Uus), H, OH, SH, OCH3, SCH3, CN and SCN. It is shown that the most important bonding orbitals in these systems have significant contributions from the Au sd hybridized atomic orbitals. The ubiquitous linear or quasi-linear structures of [AuX2](-) are attributed to the well-balanced optimal overlap in both σ and π bonding orbitals and minimal repulsion between the two negatively charged ligands. The stability of these complexes is related to the covalency of the Au-X bond and a periodic trend is found in the evolution of covalency along the halogen group ligands. The special stability of [Au(CN)2](-) is a result of strong covalent and ionic interactions. For the superheavy element Uus, the covalency of Au-Uus is enhanced through the spin-orbit interactions.

Adsorption of Pd, Pt, Cu, Ag, and Au monomers on NiAl(110) surface: a comparative study from DFT calculations

First principles calculations based on periodic density functional theory (DFT) have been used to investigate the structural, energetic and electronic properties of different transition metal atoms (Pd, Pt, Cu, Ag, and Au) on the NiAl(110) surface at low coverages (0.08 and 0.25 monolayer). All adatoms prefer to adsorb on 4-fold coordinated sites interacting with two Al and two Ni atoms and forming polar and covalent bonds, respectively. The calculated negative work function changes are explained by the effect of positive surface image created after adsorption, which induces the polarization of the negatively charged adsorbates. Consequently, for metals with similar electronegativity as Ni (Ag and Cu), this polarization effect becomes more significant and leads to larger negative work function changes, but the charge transferred is small.

The pseudopotential approximation in electronic structure theory

A short review is presented on one of the most successful theories for electronic structure calculations, the pseudopotential approximation, originally introduced by Hans G. A. Hellmann in 1934. Recent developments in relativistic quantum theory allow for the accurate adjustment of pseudopotential parameters to valence spectra, producing results for properties of atoms, molecules, and the solid-state in excellent agreement with more accurate all-electron results if a small-core definition is used. Thus the relativistic pseudopotential approximation is now the most widely applied method for systems containing heavy elements.

Detecting positron-atom bound states through resonant annihilation

A method is proposed for detecting positron-atom bound states by observing enhanced positron annihilation due to electronic Feshbach resonances at electron-volt energies. The method is applicable to a range of open-shell transition-metal atoms which are likely to bind the positron: Fe, Co, Ni, Tc, Ru, Rh, Sn, Sb, Ta, W, Os, Ir, and Pt. Estimates of their binding energies are provided.

Interaction of coinage metal clusters with chalcogen dihydrides

The interaction of chalcogen dihydrides (H 2E; E = O, S, and Se) with small coinage metal clusters (M n ; M = Cu, Ag, and Au, n = 3 and 4) is studied based on density functional theory, with a focus on the nature of chalcogen-metal bonds. A newly developed pseudopotential-based correlation-consistent basis set is used for metal clusters together with the 6-311++G** basis set for the remaining atoms. Geometrical data identified that no significant deviation has been observed for molecules before and after complexation. For these three metals, binding energy calculations indicate that gold has the highest and silver has the lowest affinities for interaction with H 2E. In comparison with gold and copper, complexation between silver and chalcogen dihydrides is significantly weaker. It is found that interaction of H 2E molecules with the coinage metals have the order of H 2Se > H 2S > H 2O. Therefore, in agreement with experimental works, our calculations confirm that the gold-selenium bond is the most stable. The nature of M-E bonds is also interpreted by means of the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses. According to the QTAIM results, the bonds are found to be partially ionic and partially covalent. Natural resonance theory (NRT) is used to calculate natural bond order and bond polarity. The NRT result indicates that the percentage of polarity of M-E bonds is affected by coinage metals.

Theoretical study of binding interactions and vibrational Raman spectra of water in hydrogen-bonded anionic complexes: (H2O)n- (n = 2 and 3), H2O...X- (X = F, Cl, Br, and I), and H2O...M- (M = Cu, Ag, and Au)

Binding interactions and Raman spectra of water in hydrogen-bonded anionic complexes have been studied by using the hybrid density functional theory method (B3LYP) and ab initio (MP2) method. In order to explore the influence of hydrogen bond interactions and the anionic effect on the Raman intensities of water, model complexes, such as the negatively charged water clusters ((H2O)n-, n = 2 and 3), the water...halide anions (H2O...X-, X = F, Cl, Br, and I), and the water-metal atom anionic complexes (H2O...M-, M = Cu, Ag, and Au), have been employed in the present calculations. These model complexes contained different types of hydrogen bonds, such as O-H...X-, O-H...M-, O-H...O, and O-H...e-. In particular, the last one is a dipole-bound electron involved in the anionic water clusters. Our results showed that there exists a large enhancement in the off-resonance Raman intensities of both the H-O-H bending mode and the hydrogen-bonded O-H stretching mode, and the enhancement factor is more significant for the former than for the latter. The reasons for these spectral properties can be attributed to the strong polarization effect of the proton acceptors (X-, M-, O, and e-) in these hydrogen-bonded complexes. We proposed that the strong Raman signal of the H-O-H bending mode may be used as a fingerprint to address the local microstructures of water molecules in the chemical and biological systems.

The effects of Group 11 transition metals, including gold, on ethylene binding to the ETR1 receptor and growth of Arabidopsis thaliana

It has been previously shown that Cu(I) and the ethylene response antagonist, Ag(I), support ethylene binding to exogenously expressed ETR1 ethylene receptors. Both are Group 11 transition metals that also include gold. We compared the effects of gold ions with those of Cu(I) and Ag(I) on ethylene binding in exogenously expressed ETR1 receptors and on ethylene growth responses in etiolated Arabidopsis seedlings. We find that gold ions also support ethylene binding but, unlike Ag(I), do not block ethylene action on plants. Instead, like Cu(I), gold ions affect seedlings independently of ethylene signaling.

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.

Vibrational and electron paramagnetic resonance properties of free and MgO supported AuCO complexes

The bonding, spin density related properties, and vibrational frequency of CO bound to single Au atom in the gas-phase or supported on MgO surfaces have been investigated with a variety of computational methods and models: periodic plane waves calculations have been compared with molecular approaches based on atomic orbital basis sets; pseudopotential methods with all electron fully relativistic calculations; various density functional theory (DFT) exchange-correlation functionals with the unrestricted coupled-cluster singles and doubles with perturbative connected triples [CCSD(T)]. AuCO is a bent molecule but the potential for bending is very soft, and small changes in the bond angle result in large changes in the CO gas-phase vibrational frequency. At the equilibrium geometry the DFT calculated vibrational shift of CO with respect to the free molecule is about -150 cm(-1), whereas smaller values -60-70 cm(-1) are predicted by the more accurate CCSD(T) method. These relatively large differences are due to the weak and nonclassic bonding in this complex. Upon adsorption on MgO, the CO vibrational shift becomes much larger, about -290 cm(-1), due to charge transfer from the basic surface oxide anion to AuCO. This large redshift is predicted by all methods, and is fully consistent with that measured for MgOAuCO complexes. The strong influence of the support on the AuCO bonding is equally well described by all different approaches.

Oxidation of Formic Acid and Carbon Monoxide on Gold Electrodes Studied by Surface-Enhanced Raman Spectroscopy and DFT

The oxidation of formic acid and carbon monoxide was studied at a gold electrode by a combination of electrochemistry, in situ surface-enhanced Raman spectroscopy (SERS), differential electrochemical mass spectrometry, and first-principles DFT calculations. Comparison of the SERS results and the (field-dependent) DFT calculations strongly suggests that the relevant surface-bonded intermediate during oxidation of formic acid on gold is formate HCOO- ad*. Formate reacts to form carbon dioxide via two pathways: at low potentials, with a nearby water to produce carbon dioxide and a hydronium ion; at higher potentials, with surface-bonded hydroxyl (or oxide) to give carbon dioxide and water. In the former pathway, the rate-determining step is probably related to the reaction of surface-bonded formate with water, as measurements of the reaction order imply a surface almost completely saturated with adsorbate. The potential dependence of the rate of the low-potential pathway is presumably governed by the potential dependence of formate coverage. There is no evidence for CO formation on gold during oxidation of formic acid. The oxidation of carbon monoxide must involve the carboxyhydroxyl intermediate, but SERS measurements do not reveal this intermediate during CO oxidation, most likely because of its low surface coverage, as it is formed after the rate-determining step. Based on inconclusive spectroscopic evidence for the formation of surface-bonded OH at potentials substantially below the surface oxidation region, the question whether surface-bonded carbon monoxide reacts with surface hydroxyl or with water to form carboxyhydroxyl and carbon dioxide remains open. The SERS measurements show the existence of both atop and bridge-bonded CO on gold from two distinguishable low-frequency modes that agree very well with DFT calculations.

Relativistic effects on the electronic structure and volatility of group-8 tetroxides MO4, where M=Ru, Os, and element 108, Hs

The influence of relativistic effects on properties and volatility of the group-8 tetroxides MO4, where M=Ru, Os, and element 108, Hs, was studied on the basis of results of the fully relativistic (four component) and nonrelativistic density functional theory calculations. Relativistic effects were shown to increase bond strengths and decrease bond lengths in these molecules. They are responsible for a decrease in molecular polarizabilities and an increase in ionization potentials. The effects are much stronger in HsO4 than in the lighter congeners. Relativistic effects were also shown to slightly decrease dispersion interaction energies of RuO4, OsO4, and HsO4 with an inert (quartz or silicon nitride) surface, i.e., they increase volatility of these compounds as studied in the "one-atom-at-a-time" gas-phase chromatography experiments. They do, however, not influence the trend in group 8: both relativistically and nonrelativistically, volatility should change as RuO4<OsO4<HsO4. The reason for that is identical trends in the relativistic and nonrelativistic space distributions of the valence d electrons.

A Comparison of Structure and Stability between the Group 11 Halide Tetramers M4X4 (M = Cu, Ag, or Au; X = F, Cl, Br, or I) and the Group 11 Chloride and Bromide Phosphanes (XMPH3)4

The tetramers of the group 11 (I) halides, M(4)X(4) (M = Cu, Ag, or Au; X = F, Cl, Br, or I), and corresponding group 11 (I) phosphanes, chloride and bromide (XMPH(3))(4) (X = Cl or Br), are investigated by the density functional theory. All coinage metal(I) halide tetramers adopt squarelike ring structures with an out-of-plane distorted (butterfly) D(2d) symmetry. These structures are much lower in energy than the more compact cubelike T(d) arrangements, which maximize dipole-dipole interactions and more closely resemble the solid-state structures of the copper and silver halides. Phosphine coordination completely changes the structures of these M(4)X(4) clusters. The copper(I) and silver(I) phosphane chloride and bromide tetramers adopt a heterocubane structure, slightly preferred over a step (ladder-type)-cluster structure well-known in the coordination chemistry of such compounds. In stark contrast, gold(I) phosphane chloride and bromide tetramers prefer assemblies of linear XAuPH(3) units with direct gold-gold contacts, resulting in a square planar, centered trigonal planar, or tetrahedral gold core.

Localized molecular constraint on electron delocalization in a metallic chain

An artificial quantum structure consisting of a single CO molecule adsorbed on a Au chain was assembled by manipulating single Au atoms on NiAl(110) at 12 K with a scanning tunneling microscope (STM). The CO disrupts the delocalization of electron density waves in the chain, as it suppresses the coupling between neighboring chain atoms. The possibility to specify the CO position on the chain allows controlled modification of the electronic properties in a quantum system. Inelastic electron tunneling spectroscopy with the STM provides vibrational characterization of the adsorbed CO.

The stability of gold iodides in the gas phase and the solid state

The stability of gold iodides in the oxidation state +I and +III is investigated at the ab initio and density functional level using relativistic and nonrelativistic energy-adjusted pseudopotentials for gold and iodine. The calculations reveal that relativistic effects stabilize the higher oxidation state of gold as expected, that is Au2I6 is thermodynamically stable at the relativistic level, whilst at the nonrelativistic level the complex of two iodine molecules weakly bound to both gold atoms in Au2I2 is energetically preferred. The rather low stability of AuI3 with respect to dissociation into AuI and I2 will make it difficult to isolate this species in the solid state as (possibly) Au2I6 or detect it by matrix-isolation techniques. The monomer AuI3 is Jahn-Teller distorted from the ideal trigonal planar (D3h) form, but adopts a Y-shaped structure (in contrast to AuF3 and AuCl3), and in the nonrelativistic case can be described as I2 weakly bound to AuI. Relativistic effects turn AuI3 from a static Jahn-Teller system to a dynamic one. For the yet undetected gas-phase species AuI accurate coupled-cluster calculations for the potential energy curve are used to predict vibrational-rotational constants. Solid-state density functional calculations are performed for AuI and Au2I6 in order to predict cohesive energies.

Spectroscopic and theoretical investigations of vibrational frequencies in binary unsaturated transition-metal carbonyl cations, neutrals, and anions

Figure 18 presents the C-O stretching vibrational frequencies of the first-row transition-metal monocarbonyl cations, neutrals, and anions in solid neon; similar diagrams have been reported for neutral MCO species in solid argon, but three of the early assignments have been changed by recent work and one new assignment added. The laser-ablation method produces mostly neutral atoms with a few percent cations and electrons for capture to make anions; in contrast, thermal evaporation gives only neutral species. Hence, the very recent neon matrix investigations in our laboratory provide carbonyl cations and anions for comparison to neutrals on a level playing field. Several trends are very interesting. First, for all metals, the C-O stretching frequencies follow the order cations > neutrals > anions with large diagnostic 100-200 cm-1 separations, which is consistent with the magnitude of the metal d to CO pi * donation. Second, for a given charge, there is a general increase in C-O stretching vibrational frequencies with increasing metal atomic number, which demonstrates the expected decrease in the metal to CO pi * donation with increasing metal ionization potential. Some of the structure in this plot arises from the extra stability of the filled and half-filled d shell and from the electron pairing that occurs at the middle of the TM row; the plot resembles the "double-humped" graph found for the variation in properties across a row of transition metals. For the anions, the variation with metal atom is the smallest since all of the metals can easily donate charge to the CO ligand. Third, for the early transition-metal Ti, V, and Cr families, the C-O stretching frequencies decrease when going down the family, but the reverse relationship is observed for the late transition-metal Fe, Co, and Ni families. In most of the present discussion, we have referred to neon matrix frequencies; however, the argon matrix frequencies are complementary, and useful information can be obtained from comparison of the two matrix hosts. In most cases, the neon-to-argon red shift for neutral carbonyls is from 11 to 26 cm-1, but a few (CrCO) lie outside of this range. In the case of FeCO and Fe(CO)2, it appears that neon and argon trap different low-lying electronic states. In general, the carbonyl neutrals and anions have similar shifts but carbonyl cations have larger matrix shifts. For example, the FeCO+ fundamental is at 2123.0 cm-1 in neon and 2081.5 cm-1 in argon, a 42.5 cm-1 shift, which is larger than those found for FeCO- (11.7 cm-1) and FeCO (11.7 cm-1). It is unusual for different low-lying electronic states to be trapped in different matrices, but CUO provides another example. The linear singlet state (1047.3, 872.2 cm-1) is trapped in solid neon, and a calculated 1.2 kcal/mol higher triplet state is trapped in solid argon (852.5, 804.3 cm-1) and stabilized by a specific interaction with argon. The bonding trends are well described by theoretical calculations of vibrational frequencies. Table 5 compares the scale factors (observed neon matrix/calculated) for the C-O stretching modes of the monocarbonyl cations, neutrals, and anions of the first-row transition metals observed in a neon matrix using the B3LYP and BP86 density functionals. Most of the calculated carbonyl harmonic stretching frequencies are within 1% of the experimental fundamentals at the BP86 level of theory, while calculations using the B3LYP functional give frequencies that are 3-4% higher as expected for these density functionals and calculations on saturated TM-carbonyls. For second- and third-row carbonyls using the BP86 density functional and the LANL effective core potential in conjunction with the DZ basis set, the agreement between theory and experiment is just as good. For example, the 16 M(CO)1-4 neutral and anion and 2 MCO+ cation (M = Ru, Os) carbonyl frequencies are fit within 1.5%. The 16 species (M = Rh, Ir) are fit within 1%, but the Rh(CO)1-4+ calculations are 2-3% too low and Ir(CO)1-4+ computations are 1-2% too low. In addition to predicting the vibrational frequencies, DFT can be used to calculate different isotopic frequencies, and isotopic frequency ratios can be computed as a measure of the normal vibrational mode in the molecule for an additional diagnostic. For diatomic CO, the 12CO/13CO ratio 1.0225 and C16O/C18O ratio 1.0244 characterize a pure C-O stretching mode. In a series of molecules such as RhCO+, RhCO, and RhCO-, where the metal-CO bonding varies, the Rh-C, C-O vibrational interaction is different and the unique isotopic ratios for the carbonyl vibration are characteristic of that particular molecule. Table 6 summarizes the isotopic ratios observed and calculated for the RhCO+,0,- species. Note that RhCO+ exhibits slightly more carbon-13 and less oxygen-18 involvement in the C-O vibration than CO itself and that this trend increases to RhCO and to RhCO- as the Rh-C bond becomes shorter and stronger. Note also how closely the calculated and observed ratios both follow this trend. In a molecule with two C-O stretching modes, for example, bent Ni(CO)2 exhibits a strong b2 mode at 1978.9 cm-1 and a weak a1 mode at 2089.7 cm-1 in solid neon, and these two modes involve different C and O participations. The symmetric mode shows substantially more C (1.0242) and less O (1.0217) participation than does the antisymmetric mode with C (1.0228) and O (1.0238) involvement, based on the given isotopic frequency ratios, which are nicely matched by DFT calculations (a1 1.0244, 1.0224 and b2 1.0232, 1.0241, respectively). These investigations of vibrational frequencies in unsaturated transition-metal carbonyl cations, neutrals, and anions clearly demonstrate the value of a close working relationship between experiment and theory to identify and characterize new molecular species.