C-N coupling on transition metal surfaces: a density functional theory study

We have investigated the formation of C-N bonds from individual atoms and single hydrogenated moieties on a series of transition metals. These reactions play a role in HCN formation at high oxygen coverage, also known as Andrussow oxidation, and they are fundamental to understand the ability of other materials to form part of alloys where Pt is the major component. Dehydrogenations take place quite easily under these high oxygen conditions and thus, the C+N, HC+N, and N+CH recombinations to form HCN or its isomer CNH might represent the rate-limiting steps for the reaction. For all the metals in the present study we have found that the activation energy for the reactions between H(x)C and NH(y) (x,y = 0,1) involved in C-N formation follow a linear relationship with the adsorption energy of the N atom. This is due to the common nature of all these transition states, where N-containing fragments get activated from three-fold hollow sites to bridge positions. The slopes of the linear dependence, though, depend on the valence of the N fragment, i.e., smaller slopes are found for NH moieties with respect to N ones.     

BFW: a density functional for transition metal clusters

Ionization potentials (IPs) or electron affinities (EAs) for transition metal clusters are an important property that can be used to identify and differentiate between clusters. Accurate calculation of these values is therefore vital. Previous attempts using a variety of DFT models have correctly predicted trends, but have relied on the use of scaling factors to compare to experimental IPs. In this paper, we introduce a new density functional (BFW) that is explicitly designed to yield accurate, absolute IPs for transition metal clusters. This paper presents the numerical results for a selection of transition metal clusters and their carbides, nitrides, and oxides for which experimental IPs are known. When tested on transition metal clusters, the BFW functional is found to be significantly more accurate than B3LYP and B3PW91.     

Effects of carbon on the stability and chemical performance of transition metal carbides: a density functional study

Density functional theory was employed to study the stabilities and chemical activities of transition metal carbides. Here we take the well-known Mo carbides and Ti carbides as an example. Different kinds of structures including the bulk surfaces [Mo(2)C(001), MoC(001), and TiC(001)] and metcars [Mo(8)C(12) and Ti(8)C(12)] are taken into consideration. Systematic studies show that by raising the C coordination number of the metal atoms in the carbides, in general the stability of the carbides increases (metcars are an exception since they include both high-coordinated and low-coordinated metal atoms.); at the same time, the chemical activities of the carbides decrease due to a downshift of the metal d-band center (ligand effect). Considering the better catalysts those that combine high stability and moderate chemical activity, our results suggest that the catalytic potential of Mo carbide systems should decrease in the following sequence: Mo(8)C(12)>Mo(2)C(001) or MoC(001)>pure Mo(110). In spite of having the largest C/Mo ratio, the metcar appears as the most attractive system. Our studies also indicate that the "magic" behavior of metcars is not unique for Mo carbides. Similar behavior is also observed for Ti carbides. This implies that nanoparticles like metcar species could exhibit better performances than the corresponding bulk metal carbides as catalysts.     

H2 adsorption on 3d transition metal clusters: a combined infrared spectroscopy and density functional study

The adsorption of H2 on a series of gas-phase transition metal (scandium, vanadium, iron, cobalt, and nickel) clusters containing up to 20 metal atoms is studied using IR-multiple photon dissociation spectroscopy complemented with density functional theory based calculations. Comparison of the experimental and calculated spectra gives information on hydrogen-bonding geometries. The adsorption of H2 is found to be exclusively dissociative on Sc(n)O+, V(n)+, Fe(n)+, and Co(n)+, and both atomic and molecularly chemisorbed hydrogen is present in Ni(n)H(m)+ complexes. It is shown that hydrogen adsorption geometries depend on the elemental composition as well as on the cluster size and that the adsorption sites are different for clusters and extended surfaces. In contrast to what is observed for extended metal surfaces, where hydrogen has a preference for high coordination sites, hydrogen can be both 2- or 3-fold coordinated to cationic metal clusters.     

A systematic study of CO oxidation on metals and metal oxides: density functional theory calculations

CO oxidation on Ru(0001), Rh(111), Pd(111), Os(0001), Ir(111), Pt(111), and their corresponding metal oxides is studied using density functional theory. It is found that (i) the reactivity of metal oxide is generally higher than that of the corresponding metal, and (ii) on both metals and metal oxides, the higher the chemisorption energy is in the initial state, the larger the reaction barrier. The barriers are further analyzed by decomposing them into electronic and geometric effects, and the higher reactivity of metal oxides is attributed mainly to the surface geometric effect. Moreover, the electronic effect on both metals and metal oxides follows the same pattern: the shorter the OC-O bond distance in the TS, the higher the barrier.     

Associative versus dissociative binding of CO to 4d transition metal trimers: A density functional study

Density functional calculations were performed to determine equilibrium geometrical structures, transition states and relative energies for M(3) clusters (M = Nb, Mo, Tc, Ru, Rh, Pd, Ag) reacting with CO, leading to proposed reaction pathways. For the Nb(3), Mo(3), and Tc(3) clusters, the lowest energy structure correlates to dissociated CO, with the C and O atoms bound on opposite sides of the metal triangle. For all other trimers, the lowest energy structures maintain the CO moiety. In the case of Pd(3) and Ag(3) the dissociated geometries lie higher in energy than the sum of the separated reactants. In most cases, several multiplicities were found to be similar in energy and for Mo(3)CO and Pd(3)CO singlet-triplet minimum energy crossing points were identified. In the case of Rh(3)CO, minimum energy crossing points for the doublet, quartet, and sextet reaction pathways were determined and compared. The electron densities of pertinent M(3)CO species were investigated using Natural Bond Order calculations. It was found that the effect of the metal trimer on the energy of the pure p-type pi* antibonding orbital of carbon monoxide directly correlates with the occurrence of CO dissociation.     

Vibrational corrections to geometries of transition metal complexes from density functional theory

Zero-point vibrational corrections are computed at the BP86/AE1 level for the set of 50 transition-metal/ligand bonds that have recently been proposed as testing ground for DFT methods, because of the availability of precise experimental gas-phase geometries (Bühl and Kabrede, J Chem Theory Comput 2006, 2, 1282). These corrections are indicated to be transferable to a large extent between various density-functional/basis-set combinations, so that they can be used to estimate zero-point averaged r0g distances from re values optimized at other theoretical levels. Applying this approach to a number of popular DFT levels does not, in general, improve their overall accuracy in terms of mean and standard deviations from experiment. The hybrid variant of the meta-functional TPSS is confirmed as promising choice for computing structures of transition-metal complexes.     

Charge state of metal atoms on oxide supports: a systematic study based on simulated infrared spectroscopy and density functional theory

A long standing question in the study of supported clusters of metal atoms in the properties of metal–oxide interfaces is the extent of metal–oxide charge transfer. However, the determination of this charge transfer is far from straight forward and a combination of different methods (both experimental and theoretical) is required. In this paper, we systematically study the charging of some adsorbed transition metal atoms on two widely used metal oxides surfaces [α-Al₂O₃ (0001) and rutile TiO₂ (110)]. Two procedures are combined to this end: the computed vibrational shift of the CO molecule, that is used as a probe, and the calculation of the atoms charges from a Bader analysis of the electron density of the systems under study. At difference from previous studies that directly compared the vibrational vawenumber of adsorbed CO with that of the gas phase molecule, we have validated the procedure by comparison of the computed CO stretching wavenumbers in isolated monocarbonyls (MCO) and their singly charged ions with experimental data for these species in rare gas matrices. It is found that the computational results correctly reproduce the experimental trend for the observed shift on the CO stretching mode but that care must be taken for negatively charged complexes as in this case there is a significative difference between the total charge of the MCO complex and the charge of the M atom. For the supported adatoms, our results show that while Cu and Ag atoms show a partial charge transfer to the Al₂O₃ surface, this is not the case for Au adatoms, that are basically neutral on the most stable adsorption site. Pd and Pt adatoms also show a significative amount of charge transfer to this surface. On the TiO₂ surface our results allow an interpretation of previous contradictory data by showing that the adsorption of the probe molecule may repolarize the Au adatoms, that are basically neutral when isolated, and show the presence of highly charged Au^δ+–CO complexes. The other two coinage metal atoms are found to significatively reduce the TiO₂ surface. The combined use of the shift on the vibrational frequency of the CO molecule and the computation of the Bader charges shows to be an useful tool for the study the charge state of adsorbed transition metal atoms and allow to rationalize the information coming from complementary tools.     

Spin quenching of transition metals deposited on MgO insulator and CdO semiconductor density functional calculations

We have analyzed spin quenching of first row transition metals deposited on (001) defect‐free and defect‐containing surfaces of MgO insulator and CdO semiconductor by means of density functional calculations and embedded cluster model. Clusters of moderate sizes were embedded in the simulated Coulomb fields that closely approximate the Madelung fields of the host surfaces, and relaxation of ions that surround the defect sites was taken into account. Spin states of metals deposited on the defect free surfaces were maintained as in the isolated metals except for Ti, V, and Co on MgO, and Ti, V, and Cr on CdO. On the defect containing surfaces, spin states were maintained too except for Fe on MgO, and V and Cr on CdO. The metal‐support interactions stabilize the low spin state of the adsorbed metal with respect to the isolated metal, but the effect was not in general enough to quench the spin. Spin polarization effects tend to preserve the spin states of the adsorbed metals relative to those of the isolated metals. Although charge transfer took place from the adsorbed metal to the insulator surface, it took place the other way round from the semiconductor surface to the adsorbed metal. The encountered variations in magnetic properties were attributed to the smaller band gap of the semiconductor, and the behavior of a single metal atom adsorbed on a particular surface was a result of a competition between Hund's rule for the adsorbed metal and the formation of a chemical bond at the interface. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Acetylene to vinylidene rearrangements on electron rich d6 metal centers: a density functional study

The acetylene to vinylidene isomerization on several Ru(II) d(6) metal fragments with different electron richness of the metal center has been investigated by means of density functional theory calculations. We considered the [(eta(5)-C(5)Me(5))Ru(dippe)](+), [(eta(5)-C(5)Me(5))Ru(dmpe)](+), [(eta(5)-C(5)H(5))Ru(PMe(3))(2)](+), [(eta(6)-C(6)Me(6))(PMe(3))ClRu](+), [(eta(5)-C(5)H(5))Ru(CO)(PPh(3))](+) and [eta(6)-C(6)H(6))(PMe(3))ClRu](+), species which are quite common in the chemistry of cationic Ru(II) complexes and span a wide range of electron-richness. For each of the considered fragments, the minima on the potential energy surfaces for the two possible isomerization mechanisms, i.e. through a direct 1,2-hydrogen shift or through a hydrido-alkynyl intermediate, have been localized. A linear correlation has been found between the C=C stretching frequencies of the vinylidene complexes, as an estimate of the electron richness, and the stability of the corresponding hydrido-alkynyl intermediates. For the most electron-rich among the considered fragments, [(Cp*)(dippe)Ru(HCCH)](+), the hydrido-alkynyl species has been found essentially isoenergetic with the alkyne complex (only 1.9 kcal mol(-1) higher), in agreement with the experimental evidence showing for this system an equilibrium between these two species. For the same [(Cp*)(dippe)Ru](+) fragment, a detailed analysis of the reaction profiles for the two possible acetylene rearrangement pathways has been performed. Our results show that once the eta(2)-C-H coordinated acetylene intermediate is accessed, the system can easily evolve towards a hydrido-alkynyl intermediate, this process being kinetically favored with respect to the direct 1,2-shift leading to the vinylidene product.     

Effect of surface pressure on the insulator to metal transition of a langmuir polyaniline monolayer

A remarkable change in the conductivity of a polyaniline (PAN) Langmuir monolayer in the conducting state, as a function of surface pressure, has been observed using scanning electrochemical microscopy (SECM). The film conductivity, as expressed by the SECM current response of a redox mediator, was measured in-situ in a Langmuir film balance. The conductivity of the film increases significantly with surface pressure, above a threshold value of ca. 20 mN m-1.     

Modeling of highly efficient drug delivery system induced by self-assembly of nanocarriers: A density functional study

Owing to the importance of drug delivery in cancer or other diseases’ therapy, the targeted drug delivery (TDD) system has been attracting enormous interest. Herein, we model the TDD system and design a novel rod-like nanocarrier by using the coarse grained model-based density functional theory, which combines a modified fundamental measure theory for the excluded-volume effects, Wertheim’s first-order thermodynamics perturbation theory for the chain connectivity and the mean field approximation for van der Waals attraction. For comparison, the monomer nanocarrier TDD system and the no nanocarrier one are also investigated. The results indicate that the drug delivery capacity of rod-like nanocarriers is about 62 times that of the no nanocarrier one, and about 6 times that of the monomer nanocarriers. The reason is that the rod-like nanocarriers would self-assemble into the smectic phase perpendicular to the membrane surface. It is the self-assembly of the rod-like nanocarriers that yields the driving force for the targeted delivery of drugs inside the cell membrane. By contrast, the conventional monomer nanocarrier drug delivery system lacks the driving force to deliver the drugs into the cell membrane. In short, the novel rod-like nanocarrier TDD system may improve the drug delivery efficiency. Although the model in this work is simple, it is expected that the system may provide a new perspective for cancer targeted therapy.     

Detection of nitrobenzene using transition metal doped C24: A DFT study

Structural Chemistry - First principles calculations have been used to verify the ability for sensing nitrobenzene by the host systems, pristine C24 and transition metal (TM) doped C24, viz.,...     

Characterization of methoxy adsorption on some transition metals: a first principles density functional theory study

Based on the gradient-density functional theory, calculation results of methoxy adsorption on Au(111), Ag(111), Cu(111), Pt(111), Pd(111), Ni(111), Rh(111), and Fe(100) surfaces are presented, and a consistent picture for some key physical properties determining the reactivity of metals appears. These eight metals belong to two groups: either with filled d electrons (group IB) or with unfilled but more than half filled d electrons (group VIII). The calculated adsorption energies are quite in agreement with the experimental data as well as the previous theoretical calculation results. Importantly, using the analysis of B. Hammer and J. K. Norskov, Nature (London) 376, 232 (1995) and in Chemisorption and Reactivity on Supported Clusters and Thin Films, edited by R. M. Lambert and G. Pacchioni (Kluwer Academic, Dordrecht, 1997), pp. 285-351, the binding energies have selectively been linearly correlated to the d-band center and to the size of the metal d-band orbital overlapping with the adsorbate (coupling matrix element) for these two groups of metals. And by analyzing the nature of the adsorption bonding, the possible reason of this difference is suggested.     

A screened hybrid density functional study on energetic complexes: Alkaline-earth metal carbohydrazide perchlorates

The molecular geometries, electronic structures and stabilities of a series of alkaline-earth metal carbohydrazide perchlorates were investigated using the Heyd–Scuseria–Ernzerhof (HSE) screened hybrid density functional. The results show that Be and Mg complexes have six-coordinated octahedron features, as previously reported for the transition metal complexes. However, Ca, Sr and Ba complexes have additional coordinated oxygen atoms from the perchlorate ion. Detailed NBO analyses indicate that the metal–ligand interactions are essentially ionic and play an important role in the stabilities of these energetic complexes. The donor–acceptor interactions result in a reduction of occupancies of σ_C=O and σ_N–H bond orbitals, and also their subsequent impact on bond length and bond order.     

The versatility of pentalene coordination to transition metals: a density functional theory investigation

DFT calculations with full geometry optimization have been carried out on a series of real and hypothetical compounds of the type [CpM(C8H6)], [(CO)3M(C8H6)], [M(C8H6)2], [(CpM)2(C8H6)], [[(CO)3M]2(C8H6)], and [M2(C8H6)2] (M = transition metal). The bonding in all the currently known compounds is rationalized, as well as in the (so far) hypothetical stable complexes. Depending on the electron count and the nature of the metal(s), eta2 (predicted), eta3, eta5, eta8, or intermediate coordination modes can be adopted. In the case of the mononuclear species, the most favored closed-shell electron counts are 18 and 16 metal valence electrons (MVE). In the case of the dinuclear species, an electron count of 34 MVEs is most favored. However, other electron counts can be stabilized, especially in the case of dinuclear complexes. Coordinated pentalene should most often be considered as formally being a dianion, but sometimes as a neutral ligand. In the former case it can behave as an aromatic species made of two equivalent fused rings, as a C5 aromatic ring connected to an allylic anion, or even as two allylic anions bridged by a C7=C8 double bond. In the latter case, it can behave as a bond-alternating cyclic polyene or as a C5 aromatic ring connected to an allylic cation.     

Time-dependent density functional theory study on optical properties of GaN doped with alkaline-earth atom

The absorption and the emission spectra of GaN and (Ga,AE)N (AE = Be, Mg, Ca) were calculated by using TDDFT method with cluster model. The calculation results show that Ga₂₆N₂₆H₅₀ cluster only has very weak absorption but has some emission bands in visible region. The band-to-band emission appears around 3.56 eV. If alkaline-earth atoms (Be, Mg and Ca) are doped into GaN, the blue emission band around 3.0 eV appears. Thereinto, Ga₂₅Ca₁N₂₆H₅₀ cluster shows a weaker yellow emission band and a stronger blue one (2.92 eV), which make (Ga,Ca)N material a promising candidate to emit blue light.     

Time-dependent density functional theory applied to ligand-field excitations and their circular dichroism in some transition metal complexes

Ligand-field transitions in [Co(en)₃]³⁺ and [Rh(en)₃]³⁺ as well as the low-energy part of the electronic spectrum of [Fe(phen)₃]²⁺ are investigated with time-dependent density functional theory (TDDFT). There is a strong functional dependence for [Co(en)₃]³⁺ and [Fe(phen)₃]²⁺. ΔSCF methods reproduce the ligand-field singlet excitation energies of [Co(en)₃]³⁺ and [Rh(en)₃]³⁺ very well. The case of [Co(en)₃]³⁺ is analyzed in some detail, in particular regarding the possibility of applying a charge-transfer (CT) correction [M.E. Casida, F. Gutierrez, J. Guan, F.-X. Gadea, D.R. Salahub, J.-P. Daudey, J. Chem. Phys. 113 (2000) 7062]. A simple CT correction would not be sufficient, but the magnitude of the charge transfer correction term in comparison with the calculated excitation energy appears to be indicative of self-interaction problems in the ground state electronic structure and in the calculated excitation energies. For the ligand-field transition of [Co(en)₃]³⁺ a hybrid functional with about 25% exact exchange performs well. Range separation/long range correction/Coulomb attenuation offers little improvement for the ligand-field transitions in [Co(en)₃]³⁺ because the occupied and unoccupied orbitals involved are in close spatial proximity.     

Coupling of photoactive transition metal complexes to a functional polymer matrix**

Conductive polymers represent a promising alternative to semiconducting oxide electrodes typically used in dye-sensitized cathodes as they more easily allow a tuning of the physicochemical properties. This can then also be very beneficial for using them in light-driven catalysis. In this computational study, we address the coupling of Ru-based photosensitizers to a polymer matrix by combining two different first-principles electronic structure approaches. We use a periodic density functional theory code to properly account for the delocalized nature of the electronic states in the polymer. These ground state investigations are complemented by time-dependent density functional theory simulations to assess the Franck-Condon photophysics of the present photoactive hybrid material based on a molecular model system. Our results are consistent with recent experimental observations and allow to elucidate the light-driven redox chemical processes – eventually leading to charge separation – in the present functional hybrid systems with potential application as photocathode materials.     

Substituent effects on benzdiyne: a density functional theory study

Density functional theory (DFT) studies were performed to investigate the effect of substituents on the properties of benzdiyne derivatives. Twelve substituted benzdiynes-C(6)X(2), where X = F, Cl, Br, Me, CF(3), CN, OH, NO(2), NH(2), OMe, NMe(2), and Ph-were considered along with the unsubstituted 1,4-benzdiyne. The structures, vibrational frequencies, and IR intensities of these benzdiynes were studied with a popular three-parameter hybrid density functional (B3LYP) combined with the split-valence 6-31G(d) basis set and Dunning's correlation-consistent polarized triple-zeta (cc-pVTZ) basis set. The relative stabilities of the substituted benzdiynes were studied with the help of reaction energies of isodesmic reactions, which showed that the electron-withdrawing groups destabilized the benzdiynes more than they did the corresponding benzenes, whereas the electron-donating groups stabilized the benzdiynes more than they did their benzene counterparts. Correlation analyses revealed that field/inductive effects played a more important role than did resonance effects. The changes in atomic charges and spin populations due to the substituents were also studied. The asymmetric nu(Ctbd1;C) stretching modes obtained were close to the 1500-cm(-)(1) mark. Reinvestigation of the experimental results supported these results; a weak IR band at 1486 cm(-)(1) was assigned to this asymmetric stretching mode in C(6)(CF(3))(2) F. Some other benzdiynes also had large IR intensity values for their asymmetric nu(Ctbd1;C) vibrational modes due to the coupling with other vibrational modes. Heats of formation for the substituted benzdiynes were obtained from the reaction energies calculated at the B3LYP/cc-pVTZ level of theory.