Periodic table of 3d-metal dimers and their ions

Antiferromagnetic Chromium-Doped Tin Clusters

The trend toward further miniaturization of micronano antiferromagnetic (AFM) spintronic devices has led to a strong demand for low-dimensional materials. The assembly of AFM clusters to produce such materials is a potential pathway that promotes studies on such clusters. In this work, we report on the discovery of the AFM Cr2Snx (x = 3-20) clusters with a stepwise growth at the density functional theory (DFT) level. In comparison, the two Cr atoms tend to stay together and be buried by Sn atoms, forming endohedral structures with one Cr atom encapsulated at size 9 and finally forming a full-encapsulated structure at size 17. Each successive cluster size is composed of its predecessor with an extra Sn atom adsorbed onto the face, giving evidence of stepwise growth. All these Cr2Snx (x = 3-20) clusters are antiferromagnets, except for the triplet-state ferrimagnetic Cr2Sn11, and all their singly negatively and positively charged ions are ferromagnets. The found stable Cr2Sn17 cluster can dimerize, yielding dimers and trimers without noticeably distorting the geometrical structure and magnetic properties of each of its constituent cluster monomers, making it possible as a building block for AFM materials.

Discovery of a series of silicon-based ferrimagnets in CrMnSin (n = 4-20) clusters

Herein, the structural evolution, electronic and magnetic properties of silicon clusters with two different dopants, CrMnSin (n = 4-20) clusters were investigated at density functional theory (DFT) level. Small-sized CrMnSin (n = 4-9) clusters tend to adopt bipyramid-based geometries, while clusters with sizes n = 10 and 11 prefer to opening cage-like structures. For sizes n = 12 to 14, the half-encapsulated structures gradually transform into closed-cage Cr@Sin structures, with the Mn atom exposed outside. Starting from size 15, both the Cr and Mn atoms are completely encapsulated by silicon atoms. Meanwhile, the Cr and Mn atoms in smaller-sized CrMnSin (n = 4-7) clusters tend to be separated, while they prefer to stay together for larger sizes. Cr atom always acts as electron donor, but not for Mn atom. From the average binding energies, one can conclude that it is easier to form larger size clusters. Smaller and larger sized CrMnSin (n = 4-9 and 19-20) clusters prefer to exhibit ferromagnetic Cr-Mn coupling, while sizes n = 10-18 always exhibit ferrimagnetic state. To our knowledge, the CrMnSin clusters is the first kind of neutral transition-metal doped semiconductor clusters that show ferrimagnetic state within a wide size range.

Bond dissociation energy of FeCr+ determined through threshold photodissociation in a cryogenic ion trap

The bond dissociation energy of FeCr+ is measured using resonance enhanced photodissociation spectroscopy in a cryogenic ion trap. The onset for FeCr+ → Fe + Cr+ photodissociation occurs well above the lowest Cr+(6S, 3d5) + Fe(5D, 3d64s2) dissociation limit. In contrast, the higher energy FeCr+ → Fe+ + Cr photodissociation process exhibits an abrupt onset at the energy of the Cr(7S, 3d54s1) + Fe+(6D, 3d64s1) limit, enabling accurate dissociation energies to be extracted: D(Fe-Cr+) = 1.655 ± 0.006 eV and D(Fe+-Cr) = 2.791 ± 0.006 eV. The measured D(Fe-Cr+) bond energy is 10%-20% larger than predictions from accompanying CAM (Coulomb Attenuated Method)-B3LYP and NEVPT2 and coupled cluster singles, doubles, and perturbative triples electronic structure calculations, which give D(Fe-Cr+) = 1.48, 1.40, and 1.35 eV, respectively. The study emphasizes that an abrupt increase in the photodissociation yield at threshold requires that the molecule possesses a dense manifold of optically accessible, coupled electronic states adjacent to the dissociation asymptote. This condition is not met for the lowest Cr+(6S, 3d5) + Fe(5D, 3d64s2) dissociation limit of FeCr+ but is satisfied for the higher energy Cr(7S, 3d54s1) + Fe+(6D, 3d64s1) dissociation limit.

A Study of NbMo and NbMo- by Anion Photoelectron Spectroscopy

Photoelectron spectra of the niobium-molybdenum diatomic anion, obtained at 488 and 514 nm, display vibrationally resolved transitions from the ground state and one excited electronic state of the anion to the ground state and one excited electronic state of the neutral molecule. The electron affinity of NbMo is measured to be 1.130 ± 0.005 eV. Its ²Δ_3/2 spin-orbit component is observed to lie 870 ± 20 cm^-1 above its previously identified ²Δ_5/2 ground state. For ⁹³Nb⁹⁸Mo, vibrational energies measured for levels up to v = 4 for the ²Δ_5/2 and ²Δ_3/2 states give harmonic frequency and anharmonicity constant values of ω_e = 492 ± 12 cm^-1 and ω_ex_e = 8.0 ± 3.2 cm^-1, the former value corresponding to a force constant of 6.80 ± 0.35 mdyn/Å. These two vibrational parameters suggest a bond dissociation energy that is too low by at least a factor of 3, indicating that the ground state potential energy curve of NbMo deviates markedly from a Morse potential at higher energies. An excited electronic state of NbMo, assigned as a ²Σ⁺ state, is observed at 2900 ± 25 cm^-1 (T₀). Vibrational energies up to v = 8 in this excited state give values of ω_e = 544 ± 8 cm^-1 and ω_ex_e = 1.9 ± 1.2 cm^-1 for ⁹³Nb⁹⁸Mo. The former value corresponds to a high vibrational force constant of 8.30 ± 0.25 mdyn/Å. Both doublet states of the neutral molecule are accessed from the anion ground state, which is assigned as ¹Σ⁺. For the ⁹³Nb⁹⁸Mo^- anion, the fundamental vibrational frequency (ΔG_1/2) is 484 ± 15 cm^-1. Electron affinity data indicate that the bond dissociation energy of NbMo^- is 0.213 ± 0.005 eV greater than that of neutral NbMo, whose previously reported value then gives D₀ = 4.85 ± 0.27 eV for the anion. An excited state of the anion lying 3050 ± 25 cm^-1 (T₀) above its ground state is assigned as ³Δ, and the energies of its spin-orbit components above the ³Δ₃ lowest energy level are measured to be 450 ± 20 cm^-1 (³Δ₂) and 1100 ± 20 cm^-1 (³Δ₁). Their uneven spacing suggests that the energy of the ³Δ₂ level is lowered by interaction with a higher energy Ω = 2 anion state. The vibrational frequency (ΔG_1/2) for the ³Δ₁ and ³Δ₂ states is measured to be 433 ± 20 cm^-1. Bond length differences among the observed states are estimated from Franck-Condon fits to vibrational band intensity profiles. When combined with the previously reported NbMo bond length, these provide bond length estimates for the ground state of the anion (1.940 ± 0.025 Å) and for the observed excited states. These species offer extreme examples of multiple metal-metal bonding, with formal bond orders of 5¹/₂ for the ²Δ ground and ²Σ⁺ excited doublet states of NbMo, 6 for the singlet ground state of the anion, and 5 for its low-lying triplet state. The relationships among their bonding properties and those of related multiply bonded transition metal dimers are discussed.

Asymmetric Solvation of the Zinc Dimer Cation Revealed by Infrared Multiple Photon Dissociation Spectroscopy of Zn2+(H2O)n (n = 1-20)

Investigating metal-ion solvation—in particular, the fundamental binding interactions—enhances the understanding of many processes, including hydrogen production via catalysis at metal centers and metal corrosion. Infrared spectra of the hydrated zinc dimer (Zn₂⁺(H₂O)_n; n = 1–20) were measured in the O–H stretching region, using infrared multiple photon dissociation (IRMPD) spectroscopy. These spectra were then compared with those calculated by using density functional theory. For all cluster sizes, calculated structures adopting asymmetric solvation to one Zn atom in the dimer were found to lie lower in energy than structures adopting symmetric solvation to both Zn atoms. Combining experiment and theory, the spectra show that water molecules preferentially bind to one Zn atom, adopting water binding motifs similar to the Zn⁺(H₂O)_n complexes studied previously. A lower coordination number of 2 was observed for Zn₂⁺(H₂O)₃, evident from the highly red-shifted band in the hydrogen bonding region. Photodissociation leading to loss of a neutral Zn atom was observed only for n = 3, attributed to a particularly low calculated Zn binding energy for this cluster size.

Revealing noncollinear magnetic ordering at the atomic scale via XMCD

Mass-selected V and Fe monomers, as well as the heterodimer \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{Fe}}_1{\text{V}}_1$$\end{document}Fe1V1, were deposited on a Cu(001) surface. Their electronic and magnetic properties were investigated via X-ray absorption (XAS) and X-ray magnetic circular dichroism (XMCD) spectroscopy. Anisotropies in the magnetic moments of the deposited species could be examined by means of angle resolving XMCD, i.e. changing the X-ray angle of incidence. A weak adatom-substrate-coupling was found for both elements and, using group theoretical arguments, the ground state symmetries of the adatoms were determined. For the dimer, a switching from antiparallel to parallel orientation of the respective magnetic moments was observed. We show that this is due to the existence of a noncollinear spin-flop phase in the deposited dimers, which could be observed for the first time in such a small system. Making use of the two magnetic sublattices model, we were able to find the relative orientations for the dimer magnetic moments for different incidence angles.

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pK_a values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

SCAN over localizes orbitals leading to spin symmetry broken ground states in graphene and benzene.

A comparative study of CunX (X = Sc, Y; n = 1–10) clusters based on the structures, and electronic and aromatic properties

The MO diagrams and orbital contributions of the HOMO and LUMO for the Cu₇Sc and Cu₇Y clusters.

Comparative Study of Single and Double Hybrid Density Functionals for the Prediction of 3d Transition Metal Thermochemistry

The performance of 13 density functionals, including hybrid-GGA, hybrid-meta-GGA, and double-hybrid functionals, in combination with the correlation consistent basis sets, has been evaluated for the prediction of gas phase enthalpies of formation for a large set of 3d transition-metal-containing molecules with versatile bonding features. Of the methods studied, the hybrid B97-1 functional and the double hybrid functional mPW2-PLYP exhibit the best overall performance with mean absolute deviations (MAD) from experimental data of 7.2 and 7.3 kcal mol(-1), respectively. For single reference molecules, where dynamic correlation predominates, the results of the hybrid functionals B97-1, B98, and ωB97X and the double hybrid functionals B2-PLYP, B2GP-PLYP, and mPW2-PLYP yield the smallest deviations from the experimental enthalpies of formation. For the prediction of thermodynamic properties of coordination complexes including metal carbonyls, B97-1 and mPW2-PLYP are the most promising functionals of those investigated. When the size of the molecule is considered, B97-1 and B98 outperform mPW2-PLYP for diatomics and triatomics, while mPW2-PLYP yields the lowest MAD for larger molecules.

Toward accurate theoretical thermochemistry of first row transition metal complexes

The recently developed correlation consistent Composite Approach for transition metals (ccCA-TM) was utilized to compute the thermochemical properties for a collection of 225 inorganic molecules containing first row (3d) transition metals, ranging from the monohydrides to larger organometallics such as Sc(C(5)H(5))(3) and clusters such as (CrO(3))(3). Ostentatiously large deviations of ccCA-TM predictions stem mainly from aging and unreliable experimental data. For a subset of 70 molecules with reported experimental uncertainties less than or equal to 2.0 kcal mol(-1), regardless of the presence of moderate multireference character in some molecules, ccCA-TM achieves transition metal chemical accuracy of ±3.0 kcal mol(-1) as defined in our earlier work [J. Phys. Chem. A2007, 111, 11269-11277] by giving a mean absolute deviation of 2.90 kcal mol(-1) and a root-mean-square deviation of 3.91 kcal mol(-1). As subsets are constructed with decreasing upper limits of reported experimental uncertainties (5.0, 4.0, 3.0, 2.0, and 1.0 kcal mol(-1)), the ccCA-TM mean absolute deviations were observed to monotonically drop off from 4.35 to 2.37 kcal mol(-1). In contrast, such a trend is missing for DFT methods as exemplified by B3LYP and M06 with mean absolute deviations in the range 12.9-14.1 and 10.5-11.0 kcal mol(-1), respectively. Salient multireference character, as demonstrated by the T(1)/D(1) diagnostics and the weights (C(0)(2)) of leading electron configuration in the complete active self-consistent field wave function, was found in a significant amount of molecules, which can still be accurately described by the single reference ccCA-TM. The ccCA-TM algorithm has been demonstrated as an accurate, robust, and widely applicable model chemistry for 3d transition metal-containing species with versatile bonding features.

Growth mechanism and chemical bonding in scandium-doped copper clusters: experimental and theoretical study in concert

Size matters! The electronic structure and size-dependent stability of neutral and cationic scandium-doped copper clusters have been investigated by mass spectrometric studies (for the cations) and also quantum chemical computations. The proposed reaction paths ultimately lead to the most stable Frank-Kasper-shaped Cu(16)Sc(+) cluster (shown here), which could be the germ of a new crystallization process.Electronic structure and size-dependent stability of scandium-doped copper cluster cations, Cu(n)Sc(+), were investigated by using a dual-target dual-laser vaporization production scheme followed by mass spectrometric studies and also quantum chemical computations in the density functional theory framework. The neutral species also were studied by using computational methods. Enhanced abundances and dissociation energies were measured in the case of Cu(n)Sc(+) for n=4, 6, 8, 10 and 16, the last of these identified as being extraordinary stable. Neutral clusters are stable with n=5, 7, 9 and 15, which are isoelectronic with respect to the number of the valence s electrons with the stable cationic clusters; hence a simple electron count determines cluster properties to a great extent. The Cu(17)Sc cluster was found to be a superatomic molecule, containing Cu(16)Sc(+) and Cu(-) units; however, the charge separation is not as pronounced as in the case of CuLi. Cu(15)Sc was found to be a stable cluster with a large dissociation energy and a closed electronic structure; hence this can be regarded as a superatom, analogous to the noble gases. The main factors determining the growth patterns of these clusters are the central position of the scandium atom and the successive filling of the shell orbitals. For smaller clusters, the reaction paths appear to diverge yielding various products; however all paths ultimately lead to the most stable Frank-Kasper shaped Cu(16)Sc cluster, which in turn can be the germ of the crystallization process.

Density functional localized orbital corrections for transition metals

This paper describes the development of the B3LYP localized orbital correction model which improves the accuracy of the B3LYP thermochemical predictions for compounds containing transition metals. The development of this model employs a large data set containing 36 experimental atomic energies and 71 bond dissociation energies. B3LYP calculations were carried out on these systems with different basis sets. Based on an electronic structure analysis and physical arguments, we built a set of 10 parameters to correct atomic data and a set of 21 parameters to correct bond dissociation energies. Using the results from our biggest basis set, the model was shown to reduce the mean absolute deviation from 7.7 to 0.4 kcalmol for the atomic data and from 5.3 to 1.7 kcalmol for the bond dissociation energies. The model was also tested using a second basis set and was shown to give relatively accurate results too. The model was also able to predict an outlier in the experimental data that was further investigated with high level coupled-cluster calculations.

Hartree-Fock complete basis set limit properties for transition metal diatomics

Numerical Hartree-Fock (HF) energies accurate to at least 1 microhartree are reported for 27 diatomic transition-metal-containing species. The convergence of HF energies toward this numerical limit upon increasing the basis set size has been investigated, where standard nonrelativistic all-electron correlation consistent basis sets and augmented basis sets, developed by Balabanov and Peterson [J. Chem. Phys. 123, 064107 (2005)], were employed. Several schemes which enable the complete basis set (CBS) limit to be determined have been investigated, and the resulting energies have been compared to the numerical Hartree-Fock energies. When comparing basis set extrapolation schemes, those in the form of exponential functions perform well for our test set, with mean absolute deviations from numerical HF energies of 234 and 153 microE(h), when the CBS limit has been determined using a two-point fit as proposed by Halkier et al. [Chem. Phys. Lett. 302, 437 (1999)] on calculations of triple- and quadruple-zeta basis set qualities and calculations of quadruple- and quintuple-zeta basis set qualities, respectively. Overall, extrapolation schemes in the form of a power series are not recommended for the extrapolation of transition metal HF energies. The impact of basis set superposition error has also been examined.

First principle study of AlX (X=3d, 4d, 5d elements and Lu) dimer

The ground state equilibrium bond length, harmonic vibrational frequency, and dissociation energy of AlX (X=3d,4d,5d elements and Lu) dimers are investigated by density functional method B3LYP. The present results are in good agreement with the available experimental and other theoretical values except the dissociation energy of AlCr. The present calculations show that the late transition metal can combine strongly with aluminum compared with the former transition metal. The present calculation also indicates that it is more reasonable to replace La with Lu in the Periodic Table and that the bonding strengths of zinc, cadmium, and mercury with aluminum are very weak.

1 Pi<--X1 Sigma+ band systems of jet-cooled ScCo and YCo

Rotationally resolved resonant two-photon ionization (R2PI) spectra of ScCo and YCo are reported. The measured spectra reveal that these molecules possess ground electronic states of (1)Sigma(+) symmetry, as previously found in the isoelectronic Cr(2) and CrMo molecules. The ground state rotational constants for ScCo and YCo are B(0)(")=0.201 31(22) cm(-1) and B(0) (")=0.120 96(10) cm(-1), corresponding to ground state bond lengths of r(0) (")=1.812 1(10) A and r(0) (")=1.983 0(8) A, respectively. A single electronic band system, assigned as a (1)Pi<--X (1)Sigma(+) transition, has been identified in both molecules. In ScCo, the (1)Pi state is characterized by T(0)=15,428.8, omega(e)(')=246.7, and omega(e)(')x(e)(')=0.73 cm(-1). In YCo, the (1)Pi state has T(0)=13 951.3, omega(e)(')=231.3, and omega(e)(')x(e) (')=2.27 cm(-1). For YCo, hot bands originating from levels up to v(")=3 are observed, allowing the ground state vibrational constants omega(e)(")=369.8, omega(e)(")x(e)(")=1.47, and Delta G(12)(")=365.7 cm(-1) to be deduced. The bond energy of ScCo has been measured as 2.45 eV from the onset of predissociation in a congested vibronic spectrum. A comparison of the chemical bonding in these molecules to related molecules is presented.

Structure and Properties of Mnn, Mnn-, and Mnn+ Clusters (n = 3−10)

Electronic and geometrical structures of Mn(3)-Mn(10) together with their singly negatively and positively charged ions are computed using density functional theory with generalized gradient approximation. The ground-state spin multiplicities in the neutral series are 16, 21, 4, 9, 6, 5, 2, and 5, for Mn(3)-Mn(10), respectively. Thus, there is a transition from a ferromagnetic ground state to a ferrimagnetic ground state at Mn(5). The energy difference between ferrimagnetic and ferromagnetic states in Mn(n) grows rapidly with increasing n and exceeds 2 eV in Mn(10). The corresponding change from ferro- to ferrimagnetic ground state occurs at Mn(6)(-) and Mn(3)(+) in the anionic and cationic series, respectively. Beginning with Mn(6), the ion spin multiplicities differ from that of the neutral by +/-1 (i.e., they obey the empirical "+/-1 rule"). We found that the energy required to remove an Mn atom is nearly independent of the charge state of an Mn(n) cluster and the number of atoms in the cluster, except for Mn(3). The results of our calculations are in reasonable agreement with experiment, except for the experimental data on the magnetic moments per atom, where, in general, we predict smaller values than the experiment.

Origin of Trans-Bent Geometries in Maximally Bonded Transition Metal and Main Group Molecules

Recent crystallographic data unambiguously demonstrate that neither Ar'GeGeAr' nor Ar'CrCrAr' molecules adopt the expected linear (VSEPR-like) geometries. Does the adoption of trans-bent geometries indicate that Ar'MMAr' molecules are not "maximally bonded" (i.e., bond order of three for M = Ge and five for M = Cr)? We employ theoretical hybrid density functional (B3LYP/6-311++G) computations and natural bond orbital-based analysis to quantify molecular bond orders and to elucidate the electronic origin of such unintuitive structures. Resonance structures based on quintuple M-M bonding dominate for the transition metal compounds, especially for molybdenum and tungsten. For the main group, M-M bonding consists of three shared electron pairs, except for M = Pb. For both d- and p-block compounds, the M-M bond orders are reflected in torsional barriers, bond-antibond splittings, and heats of hydrogenation in a qualitatively intuitive way. Trans-bent structures arise primarily from hybridization tendencies that yield the strongest sigma-bonds. For transition metals, the strong tendency toward sd-hybridization in making covalent bonds naturally results in bent ligand arrangements about the metal. In the p-block, hybridization tendencies favor high p-character, with increasing avidity as one moves down the Group 14 column, and nonlinear structures result. In both the p-block and the d-block, bonding schemes have easily identifiable Lewis-like character but adopt somewhat unconventional orbital interactions. For more common metal-metal multiply bonded compounds such as [Re2Cl8]2-, the core Lewis-like fragment [Re2Cl4]2+ is modified by four hypervalent three-center/four-electron additions.

Characterization of synthetic oxomanganese complexes and the inorganic core of the O2-evolving complex in photosystem II: Evaluation of the DFT/B3LYP level of theory

The capabilities and limitations of the Becke-3-Lee-Yang-Parr (B3LYP) hybrid density functional are investigated as applied to studies of mixed-valent multinuclear oxomanganese complexes. Benchmark calculations involve the analysis of structural, electronic and magnetic properties of di-, tri- and tetra-nuclear Mn complexes, previously characterized both chemically and spectroscopically, including the di-mu-oxo bridged dimers [Mn(III)Mn(IV)(mu-O)(2)(H(2)O)(2)(terpy)(2)](3+) (terpy=2,2':6,2''-terpyridine) and [Mn(III)Mn(IV)(mu-O)(2)(phen)(4)](3+) (phen=1,10-phenanthroline), the Mn trimer [Mn(3)O(4)(bpy)(4)(H(2)O)(2)](4+) (bpy=2,2'-bipyridine), and the tetramer [Mn(4)O(4)L(6)](+) with L=Ph(2)PO(2)(-). Furthermore, the density functional theory (DFT) B3LYP level is applied to analyze the hydrated Mn(3)O(4)CaMn cluster completely ligated by water, OH(-), Cl(-), carboxylate and imidazole ligands, analogous to the '3+1 Mn tetramer' of the oxygen-evolving complex of photosystem II. It is found that DFT/B3LYP predicts structural and electronic properties of oxomanganese complexes in pre-selected spin-electronic states in very good agreement with X-ray and magnetic experimental data, even when applied in conjunction with rather modest basis sets. However, it is conjectured that the energetics of low-lying spin-states is beyond the capabilities of the DFT/B3LYP level, constituting a limitation to mechanistic studies of multinuclear oxomanganese complexes where until now the performance of DFT/B3LYP has raised little concern.

The performance of semilocal and hybrid density functionals in 3d transition-metal chemistry

We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic pi ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca(2), Mn(2), and Zn(2) are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground-state symmetry in the vast majority of cases and rarely fails qualitatively. This further confirms TPSS as a general purpose functional that works throughout the periodic table. We also give workstation timing comparisons for the 645-atom protein crambin.