A density functional study of bare and hydrogenated platinum clusters

We perform density functional theory calculations using Gaussian atomic-orbital methods within the generalized gradient approximation for the exchange and correlation to study the interactions in the bare and hydrogenated platinum clusters. The minimum-energy structures, binding energies, relative stabilities, vibrational frequencies and the highest occupied and lowest unoccupied molecular-orbital gaps of Pt_nH_m (n = 1–5, m = 0–2) clusters are calculated and compared with previously studied pure platinum and hydrogenated platinum clusters. We investigate any magic behavior in hydrogenated platinum clusters and find that Pt₄H₂ is more stable than its neighboring sizes. The lowest energy structure of Pt₄ is found to be a distorted tetrahedron and that of Pt₅ is found to be a bridge site capped tetrahedron which is a new global minimum for Pt₅ cluster. The successive addition of H atoms to Pt_n clusters leads to an oscillatory change in the magnetic moment of Pt₃–Pt₅ clusters.     

A density functional study of the SERS spectra of pyridine adsorbed on silver clusters

The effect of binding pyridine to silver clusters has been studied by density functional calculations by adopting a hybrid functional. The calculations allow proposing an explanation for the different SERS spectra of the ligand observed on Ag colloids in the presence and in the absence of coadsorbed chloride anions. In the latter case, a better agreement is obtained modeling the system by adsorption of pyridine on a surface cluster.     

A density functional theory study of gold clusters supported on layered double hydroxides

The geometrical structure and electronic properties of a series of Au _N (N = 1–8) clusters supported on a Mg²⁺, Al³⁺-containing layered double hydroxides (MgAl–LDH) are investigated using density functional theory. The Au clusters are supported on two typical crystal faces of the LDH platelet, the basal {0001} and the lateral $ \{ 10\,\bar{1}\,0\} $ crystal face, respectively, corresponding to the top and edge site of monolayer MgAl–LDH lamella for the sake of simplicity. It is revealed that an increase in the charge transfer from the LDH lamella to the Au _N clusters at the edge site rather than clusters on the top surface, demonstrating a preferential adsorption for Au _N clusters at the edge of LDH lamella. Moreover, the calculated adsorption energy of the Au _N clusters on the LDH lamella increases with the cluster size, irrespective of the adsorption site. The investigation on the interaction between O₂ and Au _N clusters on the LDH lamella is further carried out for understanding the catalytic oxidation properties of the LDH-supported Au catalyst. The formation of reactive O₂ ^? species, a necessary prerequisite in catalytic oxidation of CO, by O₂ bridging two Au atoms of Au _N clusters indicates that the LDH-supported Au catalyst has the required characteristics of a chemically active gold catalyst in CO oxidation.     

Stability and structure of rare-gas ionic clusters using density functional methods: A study of helium clusters

Several modelings of exchange and correlation forces which can be carried out using density functional theory (DFT) methods have been analyzed to study their efficiency and reliability when evaluating possible competing structures of helium ionic clusters of increasing size. This study examines He_n⁺systems with n from 1 to 7 and compares the present calculations with earlier evaluations that used more conventional, and more computationally intensive, methods with configuration interaction (CI) approaches. The present results indicate that it is indeed possible to strike a fruitful balance between reduction of computational times and quality of the ensuing structural information. © 1996 John Wiley & Sons, Inc.     

Two hydrogen ligands on tetrairidium clusters: a relativistic density functional study

Structural and energetic properties of Ir(4)H(2) have been determined by applying a relativistic density functional method. As previously obtained for Ir(4)H, terminal coordination of H ligands is preferred, in contrast to some other transition metals. Square-planar Ir(4) isomers with an H binding energy of up to 318 kJ mol(-1) per atom were determined as the most stable structures of Ir(4)H(2). Isomers with a tetrahedral or butterfly structure of the metal framework exhibit average H atom binding energies of up to approximately 300 kJ mol(-1). For all three types of isomers, a surprisingly large number of stable minima was identified. Unexpectedly, structural as well as energetic properties of Ir(4)H(2) complexes are very similar to Ir(4)H. Thus binding of an H atom to Ir(4) is only slightly affected by the presence of a second H ligand. In all cases examined, the reaction H(2)+ Ir(4)--> H(2)Ir(4) was found to be exothermic with reaction energies of up to 170 kJ mol(-1).     

Catalytic hydrogenolysis of alkyl halides by sulfido-bridged molybdenum clusters: A density functional study

Density functional theory has been used to explore the mechanism of cleavage of H₂ at a sulfido-bridged molybdenum cluster, CpMo(μ-SH)(μ-S)(μ-S₂CH₂)MoCp. The addition occurs across a single Mo-S bond, and the disruption of the strong Mo-S π bonding in the ground state leads to a very high-lying transition state (+43 kcal mol⁻¹). Once formed, the adsorbed hydrogen migrates over the cluster via a series of hops from metal to sulphur, formally corresponding to a switch from hydridic to protic character. The low barrier (+15 kcal mol⁻¹) for migration leads to facile hydrogenolysis of coordinated substrates.     

A density functional study of YnAl (n=1-14) clusters

The geometries, stabilities, and electronic and magnetic properties of Y(n)Al (n=1-14) clusters have been systematically investigated by using density functional theory with generalized gradient approximation. The growth pattern for different sized Y(n)Al (n=1-14) clusters is Al-substituted Y(n+1) clusters and it keeps the similar frameworks of the most stable Y(n+1) clusters except for Y(9)Al cluster. The Al atom substituted the surface atom of the Y(n+1) clusters for n<9. Starting from n=9, the Al atom completely falls into the center of the Y-frame. The Al atom substituted the center atom of the Y(n+1) clusters to form the Al-encapsulated Y(n) geometries for n>9. The calculated results manifest that doping of the Al atom contributes to strengthen the stabilities of the yttrium framework. In addition, the relative stability of Y(12)Al is the strongest among all different sized Y(n)Al clusters, which might stem from its highly symmetric geometry. Mulliken population analysis shows that the charges always transfer from Y atoms to Al atom in all different sized clusters. Doping of the Al atom decreases the average magnetic moments of most Y(n) clusters. Especially, the magnetic moment is completely quenched after doping Al in the Y(13), which is ascribed to the disappearance of the ininerant 4d electron spin exchange effect. Finally, the frontier orbitals properties of Y(n)Al are also discussed.     

A density functional theory study on the water aggregation behaviour of fatty acid-based anionic surface active ionic liquids

Structural Chemistry - The hydrogen bond interactions between methyl-imidazolium cation (MIM+) and fatty acid anions (CmHnCOO–, where m?=?1–6; n-3–13) of ionic liquids...     

Au_n团簇催化水煤气变换反应机理的密度泛函理论研究

利用密度泛函理论(DFT)研究了Au10、Au13和Au20三类团簇的稳定性和对水煤气变换(WGSR)反应的催化活性,考察了各物质在Aun团簇上的吸附行为和微观反应机理。结果表明,三类Aun团簇的稳定性顺序为Au10Au13Au20,而Aun团簇中电子离域性及吸附能力大小趋势为Au13Au10Au20。在三类Aun团簇上,水煤气变换反应的控速步骤均为H2O的解离,但其反应机理路径有所不同。Au10团簇上为羧基机理,COOH*中间体直接解离;Au13团簇上为氧化还原机理,两个OH*发生歧化反应;Au20团簇上为羧基机理,COOH*和OH*发生歧化反应。通过对三类团簇上的最佳反应路径进行比较发现,Au13团簇在低温下具有较好的催化活性。     

A density functional theory study on pelargonidin

A complete conformational analysis on the isolated and polarizable continuum model (PCM) modeled aqueous solution cation, quinonoidal, and anion forms of pelargonidin, comprising the diverse tautomers of the latter forms, was carried out at the B3LYP/6-31++G(d,p) level. The results indicate that the most stable conformer of cationic and quinonoidal forms of pelargonidin are completely planar in the gas phase, whereas that of the anionic form is not planar. In contrast, PCM calculations show that the plane of the B ring is slightly rotated with regard to the AC bicycle in the most stable conformer of the cation and quinonoidal form. The most stable conformers of the cation, both in gas phase and aqueous solution, display anti and syn orientations for, respectively, C2-C3-O-H and C6-C5-O-H dihedral angles, whereas syn and anti orientation of hydroxyls at 7 and 4' positions are nearly isoenergetic. The most stable tautomer of quinonoidal pelargonidin is obtained by deprotonating hydroxyl at C5 in gas phase but at C7 according to PCM. Also, the most stable tautomer of the anion is different in gas phase (hydrogens are abstracted from hydroxyls at C5 and C4') and PCM simulation (C3 and C5). Tautomeric equilibria affect substantially the geometries of the AC-B backbone providing bond length variations that basically agree with the predictions of the resonance model. Most of the conformers obtained display an intramolecular hydrogen bond between O3 and H6'. Nevertheless, this interaction is not present in the most stable anions. Ionization potentials and O-H bond dissociation energies computed for the most stable conformers of cation, quinonoidal, and anion forms are consistent with an important antioxidant activity.     

Anion photoelectron spectroscopy and density functional study of small aluminum-vanadium oxide clusters

Small aluminum-vanadium oxide clusters, AlVO(y)(-) (y = 1-3) and Al(x)VO(2)(-) (x = 2, 3), were investigated with anion photoelectron spectroscopy and density functional calculations. The adiabatic detachment energies of AlVO(y)(-) were estimated to be 1.06 ± 0.05, 1.50 ± 0.08, and 2.83 ± 0.08 eV for y = 1, 2, and 3. Those of Al(2)VO(2)(-) and Al(3)VO(2)(-) were estimated to be 1.22 ± 0.08 and 1.25 ± 0.08 eV. Comparison of theoretical calculations with experimental measurement suggests that the most probable structure of AlVO(-) cluster is quasilinear with O atom in the middle. AlVO(2)(-) has an irregular chain structure of Al-O-V-O and a C(2v) cyclic structure very close in energy. The structure of AlVO(3)(-) cluster is evolved from the C(2v) cyclic AlVO(2)(-) structure by adding the third O atom to the V atom. Al(2)VO(2)(-) has a pair of nearly degenerate Al-O-V-O-Al chain structures that can be considered as cis and trans forms. Al(3)VO(2)(-) probably has two low-lying isomers each containing a four-membered ring. The structures of the corresponding neutral clusters are discussed.     

NO adsorption on MoS(x) clusters: a density functional theory study

The density functional theory (DFT) method has been used to investigate NO probe molecule adsorption on the stoichiometric (Mo(16)S(32)) and nonstoichiometric (Mo(16)S(34) and Mo(16)S(29)) clusters. The calculated adsorption energies indicate that the stoichiometric cluster has stronger NO affinity than the nonstoichiometric surfaces. It is also found that mononitrosyl adsorption is favored at low NO coverage, while dinitrosyl (germinal) and (NO)(2) dimer adsorption at high NO coverage are possible. Strong repulsive interaction has been found for the adsorbed dinitrosyl and (NO)(2) dimer species. In addition, the computed NO stretching frequencies for the mononitrosyl and dinitrosyl species agree well with the experimental data, while those of the dimer species are much lower than the suggested experimental data.     

Ionization potentials of tantalum-carbide clusters: an experimental and density functional theory study

We have used photoionization efficiency spectroscopy to determine the ionization potentials (IP) of the tantalum-carbide clusters, Ta3Cn (n = 1-3) and Ta4Cn (n = 1-4). The ionization potentials follow an overall reduction as the number of carbon atoms increases; however, the trend is not steady as expected from a simple electrostatic argument. Instead, an oscillatory behavior is observed such that clusters with an odd number of carbon atoms have higher IPs and clusters with an even number of carbon atoms have lower IPs, with the Ta4C4 cluster exhibiting the lowest IP. Excellent agreement is found with relative IPs calculated using density functional theory for the lowest energy structures, which are consistent with the development of a 2 x 2 x 2 face-centered nanocrystal. This work shows that IPs may be used as a reliable validation for the geometries of metal-carbide clusters calculated by theory. The variation in IP can also be interpreted qualitatively with application of a simple model based upon isolobal frontier orbitals.     

A density functional theory for studying ionization processes in water clusters

A generalized Kohn-Sham (GKS) approach to density functional theory (DFT), based on the Baer-Neuhauser-Livshits range-separated hybrid, combined with ab initio motivated range-parameter tuning is used to study properties of water dimer and pentamer cations. The water dimer is first used as a benchmark system to check the approach. The present brand of DFT localizes the positive charge (hole), stabilizing the proton transferred geometry in agreement with recent coupled-cluster calculations. Relative energies of various conformers of the water dimer cation compare well with previously published coupled cluster results. The GKS orbital energies are good approximations to the experimental ionization potentials of the system. Low-lying excitation energies calculated from time-dependent DFT based on the present method compare well with recently published high-level "equation of motion-coupled-cluster" calculations. The harmonic frequencies of the water dimer cation are in good agreement with experimental and wave function calculations where available. The method is applied to study the water pentamer cation. Three conformers are identified: two are Eigen type and one is a Zundel type. The structure and harmonic vibrational structure are analyzed. The ionization dynamics of a pentamer water cluster at 0 K shows a fast <50 fs transient for transferring a proton from one of the water molecules, releasing a hydroxyl radical and creating a protonated tetramer carrying the excess hole.     

Growth pattern and electronic properties of acetonitrile clusters: a density functional study

We report a systematic theoretical study on the growth pattern and electronic properties of acetonitrile clusters [(CH(3)CN)(n) (n = 1, 9, 12)] using density functional approach at the B3LYP6-31++G(d,p) level. Although we have considered a large number of configurations for each cluster, the stability of the lowest energy isomer was verified from the Hessian calculation. It is found that the lowest energy isomer of the dimer adopts an antiparallel configuration. For trimer and tetramer, cyclic ring structures were found to be favored over the dipole stabilized structure. In general, it is found that the intermolecular CH...N interactions play a significant role in the stabilization of the cyclic layered geometry of acetonitrile clusters. A critical comparison between trimer and tetramer clusters suggests that the three member cyclic ring is more stable than four member rings. The growth motif for larger clusters (n = 5-9, 12) follows a layered pattern consisting of three or four membered rings, which, in fact, is used as the building block. Based on the stability analysis, it is found that clusters with an even number of molecular entities are more stable than the odd clusters, except trimer and nonamer. The exceptional stability of these two clusters is attributed to the formation of trimembered cyclic rings, which have been found to form the building blocks for larger clusters.     

A density functional study on nitrogen-doped carbon clusters CnN3- (n=1-8)

Using molecular graphics software, we designed numerous models of CnN3- (n=1-8). Geometry optimization and calculation on vibration frequency were carried out by the B3LYP density functional method. After comparison of structure stability, we found that the structures of ground-state CN3- and C2N3- are bent chains with a nitrogen atom at either end, whereas when n=3-8, the ground-state clusters show three branches, each with a nitrogen atom located at the end. When n=5-8, the longest branch of CnN3- is polyacetylenelike. When n=5 or 7, the longest branch is connected to the central sp2 carbon in a nonlinear manner. The CnN3- (n=1-8) with an even number of carbon atoms are more stable than those with odd numbers, matching the peak pattern observed in laser-induced mass spectra of CnN3-. The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electron affinities, and incremental binding energies.     

On the accuracy of density functional theory for iron-sulfur clusters

A simple, yet powerful wave function manipulation method was introduced utilizing a generalized ionic fragment approach that allows for systematic mapping of the wave function space for multispin systems with antiferromagnetic coupling. The use of this method was demonstrated for developing ground state electronic wave function for [2Fe-2S] and [Mo-3Fe-4S] clusters. Using well-defined ionic wave functions for ferrous and ferric irons, sulfide, and thiolate fragments, the accuracy of various density functionals and basis sets including effective core potentials were evaluated on a [4Fe-4S] cluster by comparing the calculated geometric and electronic structures with crystallographic data and experimental atomic spin densities from X-ray absorption spectroscopy, respectively. We found that the most reasonable agreement for both geometry and atomic spin densities is obtained by a hybrid functional with 5% HF exchange and 95% density functional exchange supplemented with Perdew's 1986 correlation functional. The basis set seems to saturate only at the triple-zeta level with polarization and diffuse functions. Reasonably preoptimized structures can be obtained by employing computationally less expensive effective core potentials, such as the Stuttgart-Dresden potential with a triple-zeta valence basis set. The extension of the described calibration methodology to other biologically important and more complex iron-sulfur clusters, such as hydrogenase H-cluster and nitrogenase FeMo-co will follow.     

A density functional theory study of the dissociation of H2 on gold clusters: importance of fluxionality and ensemble effects

Density functional theory was employed to calculate the adsorption/dissociation of H2 on gold surfaces, Au(111) and Au(100), and on gold particles from 0.7 (Au14) to 1.2 nm (Au29). Flat surfaces of the bulk metal were not active towards H2, but a different effect was observed in gold nanoclusters, where the hydrogen was adsorbed through a dissociative pathway. Several parameters such as the coordination of the Au atoms, ensemble effects and fluxionality of the particle were analyzed to explain the observed activity. The effect of the employed functional was also studied. The flexibility of the structure, i.e., its adaptability towards the adsorbate, plays a key role in the bonding and dissociation of H2. The interaction with hydrogen leads to drastic changes in the structure of the Au nanoparticles. Furthermore, it appears that not only low coordinated Au atoms are needed because H2 adsorption/dissociation was only observed when a cooperation between several (4) active Au atoms was allowed.     

Adsorption of CO on oxygen preadsorbed neutral and charged gas phase Pd4 clusters: A density functional study

We present the results of a density functional calculation on adsorption of O₂, CO, and their coadsorption at various sites of neutral, cationic, and anionic Pd₄ clusters. For all the clusters, the dissociative adsorption of oxygen sitting on Pd bridge sites is found to be preferable. Both O₂ and CO binding energies are found to be higher for the anionic Pd₄ cluster followed by cationic and neutral cluster. However, binding energies of O₂ or CO in the coadsorption complexes follow the trend: anionic > neutral > cationic. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

A density functional study of uranyl monocarboxylates

We studied uranium(VI) monocarboxylate complexes by a relativistic density functional method using simple carboxylic acids as ligands, i.e. [UO2(OOCR)]+ (R = H, CH3, CH2CH3). These complexes exist in aqueous solution and, for R = CH3 and CH2CH3, may also be considered as models of uranyl complexated by humic substances. We investigated mono- and bidentate coordination modes. Short-range solvent effects were accounted for explicitly via aqua ligands of the first hydration shell and long-range electrostatic interactions were described via a polarizable continuum model. The calculated results for the uranyl U=O bond, the bond to aqua ligands, and the averaged uranium distances to equatorial oxygen atoms, U-Oeq, agreed quite well with EXAFS-derived interatomic distances. However, the uranyl-carboxylate bond was calculated to be notably shorter than the experimentally determined value. Experimental differences between mono- and bidentate coordination, obtained mainly from crystal structures, were qualitatively reproduced for the U-C distance but not for the average bond length, U-Oeq. We discuss these discrepancies between calculated and experimental results in some detail and suggest changes in the coordination number rather than variations of the coordination geometry as the main source of the experimentally observed variation of the U-Oeq distance.