A density functional theory study of the Cu+·(CO)n (n = 1–3) complexes

Density functional theory calculations, with an effective core potential for the copper ion, and large polarized basis set functions have been used to construct the potential energy surface of the Cu⁺·(CO)_n (n = 1–3) complexes. A linear configuration is obtained for the global minimum of the Cu⁺·CO and Cu⁺·(CO)₂ complexes with a bond dissociation energy (BDE) of 35.9 and 40.0 kcal mol^-1, respectively. For the Cu⁺·(CO)₃ complex, a trigonal planar geometry is obtained for the global minimum with a BDE of 16.5 kcal mol^?1. C-coordinated copper ion complexes exhibit stronger binding energy than O-coordinated complexes as a result of C_lp → 4s σ-donation. The computed sequential BDEs of Cu⁺·(CO)_n (n = 1–4) complexes agree well with experimental findings, in which the electrostatic energy and σ-donation play an important role in the observed trend.     

A density functional theory investigation of CrSin (n=1–6) clusters

Clusters of the form CrSi_n (n=1–6) were investigated computationally using a density functional approach. In particular, geometry optimizations were carried out under the constraint of well-defined point group symmetries at the B3LYP level employing a pseudopotential method in conjunction with double zeta basis sets. In this article, the resulting total energies, Mulliken atomic net populations, overlap populations, fragmentation energies and geometries of CrSi_n (n=1–6) are presented and discussed, together with natural populations and natural electron configurations. In addition, we comment on the charge transfer within the clusters. From this analysis, the 3d orbital of the Cr atom in CrSi_n (n=1–6) cluster absorbs electrons. From this tendency, conclusions are drawn with respect to the electronic populations and the chemical bond between Si and Cr as well as Si and Si.     

The density functional theory study CO oxidation catalyzed by subnanometer AlAg n (n = 1–3) clusters

The efficiency of AlAg _n (n = 1–3) alloy clusters toward CO oxidation is demonstrated from first-principles theory. It is found that these subnanometer species transform into reaction complexes which catalyze CO oxidation through the Langmuir–Hinshelwood path. It is shown that mixing two different metals (Al and Ag) can have beneficial effects on the catalytic activity and the alloyed AlAg₃ cluster is proposed as the best effective nanocatalysts.     

Density functional calculation on a high energy density compound having the formula C2OH4?n (NO2) n

A series of ethylene oxide derivations, C₂OH_4?Cn (NO2) _n (x?=?1?C4), has been designed computationally. We calculated the heats of formation (HOFs), bond dissociation energy (BDE), and explosive performances (detonation velocity and detonation pressure) of the title compounds by using density functional theory with 6-311G** basis set. The results show that most of ethylene oxide derivations have positive HOFs values except I. All the calculated BDE and the bond dissociation energies without zero-point energy corrections (BDE⁰) are larger than 200?kJ?mol^?1, which indicate that all the ethylene oxide derivations have good thermal stabilities. The explosive performances of most of ethylene oxide derivations would rank up with cyclotrimethylenetrinitramine (RDX). The results have not only shown that these compounds may be used as high energy density compounds, but also provide some useful information for further syntheses.     

Gaussian binning of the vibrational distributions for the Cl + CH4(v(4/2) = 0, 1) → H + CH3Cl(n(1)n(2)n(3)n(4)n(5)n(6)) reactions

We test several binning techniques to obtain mode-specific final-state distributions for polyatomic reactions. Normal mode analysis is done after an exact transformation to the Eckart frame. Standard histogram binning (HB) and three different variants of the energy-based Gaussian binning (1GB) are employed to obtain the probabilities of the vibrational states. We consider the two major issues of the polyatomic quasiclassical product analysis, i.e., (1) rounding the classical action to the nearest integer can result in unphysical states and (2) the normal-mode analysis can break down for highly distorted geometries. We show that 1GB can handle issue 1 when the total vibrational energy is evaluated in the normal mode space using the harmonic approximation and both issues 1 and 2 can be solved when the total vibrational energy is calculated exactly in the Cartesian space. We found that anharmonicity in the quantized energy levels does not have a significant effect on the final-state distributions. Quasiclassical trajectory calculations are performed for the reactant ground-state and bending-excited Cl((2)P(3/2)) + CH(4)(v(4/2) = 0, 1) → H + CH(3)Cl reactions using an ab initio potential energy surface. The product analysis techniques are successfully applied to the CH(3)Cl product molecules and some qualitative features of the results are discussed.     

Is the uniform electron gas limit important for small Ag clusters? Assessment of different density functionals for Ag(n) (n < or = 4)

Twenty-three density functional theory (DFT) methods, including the second- and the third-generation functionals, are tested in conjunction with two basis sets (LANL2DZ and SDD) for studying the properties of neutral and ionic silver clusters. We find that DFT methods incorporating the uniform electron gas limit in the correlation functional, namely, those with Perdew's correlation functionals (PW91, PBE, P86, and TPSS), Becke's B95, and the Van Voorhis-Scuseria functional VSXC, generally perform better than the other group of functionals, e.g., those incorporating the LYP correlation functional and variations of the B97 functional. Strikingly, these two groups of functionals can produce qualitatively different results for the Ag3 and Ag4 clusters. The energetic properties and vibrational frequencies of Ag(n) are also evaluated by the different functionals. The present study shows that the choice of DFT methods for heavy metals may be critical. It is found that the exact-exchange-incorporated PBE functional (PBE1PBE) is among the best for predicting the range of properties.     

Structure and stability of hypervalent NLi n Na2 (n = 1–4) and related species at density functional theory level

A number of configurations of NLi _n Na₂ (n = 1–4) species were optimized using the B3LYP–density functional theory method; the 6-31G^* basis set was used in this calculation. In order to study all possible dissociation energies, some related species such as NLi₂Na, NLi _n (n = 1–4), Li _n (n = 1, 2) and Na _n (n = 1, 2) were also considered. Optimizations of these species were followed by fundamental frequency calculations at the same level. Global minima of these species were shown to adopt C _{2 v} (NLi₄Na₂, NLi₂Na₂), D _{3 h} (NLi₃Na₂) and C _s (NLiNa₂ and NLi₂Na) configurations. All possible dissociation energies were obtained. Received: 30 November 1998 / Accepted: 15 October 1999 / Published online: 14 March 2000     

(HAlNH)n(n=2～3)环状化合物的结构特征

用自洽场理论 (HF)和密度泛函理论 (DFT)的B3LYP方法 ,在 6 31G 的水平上对化合物(HAlNH) 2 和 (HAlNH) 3 的几何结构进行优化 ,并分别与环丁二烯C4 H4 和苯分子C6H6的结构和成键方式进行比较。以B3LYP STO 3G方法讨论其分子轨道波函数 (Ψ)。结果表明 :C4 H4 和 (HAlNH) 2 均为D2h对称 ,前者为长方形结构 ,形成两个孤立的π键 ;而后者为菱形结构 ,形成一个π44键。C6H6和 (HAlNH) 3分子点群分别为D6h和D3h,并均形成一个π66键。成键原子对分子轨道的贡献不同 ,其中C原子是完全等价的 ,而Al和N原子各不相同 ,N原子比Al的贡献要大得多     

Structures,stabilities, and electronic properties of Al n Au (n = 1–15) clusters: A density functional study

We present density functional calculations of Al _n Au clusters for n = 1–15. The growth pattern for Al _n Au (n = 1–7, 12, 14, 15) clusters is the Au atom occupying a peripheral position of Al _n clusters, and the growth pattern for Al _n Au (n = 8, 10 and 13) clusters is Au-substituted Al _n+1 clusters. It is found that the Au atom replaces the surface atom of an Al _n+1 cluster and occupies a peripheral position. In addition, the ground state structures of Al _n Au clusters are more stable than pure Aln clusters. It is found that the Al₁₃Au cluster exhibits high stability.     

Assessment of multicoefficient correlation methods, second-order Møller-Plesset perturbation theory, and density functional theory for H3O(+)(H2O)n (n = 1-5) and OH(-)(H2O)n (n = 1-4)

We have assessed the ability of 52 methods including 15 multicoefficient correlation methods (MCCMs), two complete basis set (CBS) methods, second-order M?ller-Plesset perturbation theory (MP2) with 5 basis sets, the popular B3LYP hybrid functional with 6 basis sets, and 24 combinations of local density functional and basis set to accurately reproduce reaction energies obtained at the Weizmann-1 level of theory for hydronium, hydroxide, and pure water clusters. The three best methods overall are BMC-CCSD, G3SX(MP3), and M06-L/aug-cc-pVTZ. If only microsolvated ion data is included, M06-L/aug-cc-pVTZ is the best method; it has errors only half as large as the other density functionals. The deviations between the three best performing methods are larger for the larger hydronium- and hydroxide-containing clusters, despite a decrease in the average reaction energy, making it impossible to determine which of the three methods is overall the best, so they might be ranked in order of increasing cost, with BMC-CCSD least expensive, followed by M06-L/aug-cc-pVTZ. However, the cost for M06-L will increase more slowly as cluster size increases. This study shows that the M06-L functional is very promising for condensed-phase simulations of the transport of hydronium and hydroxide ions in aqueous solution.     

Gaussian density-functional study for small neutral (Al n ), positive (Al n + ) and negative (Al n − ) aluminium clusters (n=2–5)

The structures and properties of Al _n , Al _n ⁺ , Al _n ^– (n=1,5) clusters have been investigated by using the Linear Combination of Gaussian Type Orbitals (LCGTO) method, considering Local (LSD) and Non Local (NLSD) Spin Density Approximations and employing a Model Core Potential (MCP) that allows the explicit treatment of 3s ² 3p ¹ valence electrons. For each system different geometrical structures and electronic states have been considered. For Al₃, Al ₃ ⁺ , Al ₃ ^– the most stable geometry proved to be the equilateral triangle (D _3h ). Al₄ and Al ₄ ⁺ prefer the rhombus (D _2h ) structure, while the corresponding anion prefers the square (D _4h ) one. The trapezoidal form (C _2v ) is the most stable isomer for Al₅, Al ₅ ⁺ and Al ₅ ^– clusters. The analysis of vibrational frequencies shows that these structures are minima in the potential energy surface. The binding energies (D _e), the adiabatic ionization potentials (IP) and electron affinities (EA), the chemical potentials or absolute hardnesses () and electronegativities () have been computed. Results are in good agreement with the available experimental data and the previous high level theoretical computations.     

Ion induced dipole clusters H(n)- (3 ≤ n-odd ≤ 13): density functional theory calculations of structure and energy

We investigate anew the possible equilibrium geometries of ion induced dipole clusters of hydrogen molecular ions, of molecular formula H(n)(-) (3 ≤ n-odd ≤ 13). Our previous publications [Sapse, A. M.; et al. Nature 1979, 278, 332; Rayez, J. C.; et al., J. Chem. Phys. 1981, 75, 5393] indicated these molecules would have a shallow minimum and adopt symmetrical geometries that accord with the valence shell electron pair repulsion (VSEPR) rules for geometries defined by electron pairs surrounding a central point of attraction. These earlier calculations were all based upon Hartree-Fock (HF) calculations with a fairly small basis of atomic functions, except for the H3(-) ion for which configuration interaction (CI) calculations were carried out. A related paper [Hirao, K.; et al., Chem. Phys. 1983, 80, 237] carried out similar calculations on the same clusters, finding geometries similar to our earlier calculations. However, although that paper argued that the stabilization energy of negative ion clusters H(n)(-) is small, vibration frequencies for the whole set of clusters was not reported, and so a definitive assertion of a true equilibrium was not present. In this paper we recalculate the energetics of the ion induced dipole clusters using density function theory (DFT) B3LYP method calculations in a basis of functions (6-311++G(d,p)). By calculating the vibration frequencies of the VSEPR geometries, we prove that in general they are not true minima because not all the resulting frequencies correspond to real values. By searching the energy surface of the B3LYP calculations, we find the true minimum geometries, which are surprising configurations and are perhaps counterintuitive. We calculate the total energy and binding energy of the new geometries. We also calculate the bond paths associated with the quantum theory of atoms in molecules (QTAIM). The B3LYP/6-311++G(d,p) results, for each molecule, deliver bond paths that radiate between each polarized H2 molecule and the polarizing H(-) ion.     

(BN)n(n≤12)团簇的结构及成键性质

利用遗传算法和Gastreich提出的经验势函数研究了(BN)_n(n≤12)团簇的可能稳定结构, 并对能量较低的异构体在HF/6-31G(d)水平进行优化, 得到了(BN)_n(n≤12)团簇的线状、蒲扇形、单环、双环、三环和笼状结构, 讨论了各种结构的特征及相对稳定性. 分析了BN团簇中原子的成键性质, 在单环结构中, N原子以sp2杂化成键, B原子以sp杂化成键, 而在节点处B原子以sp2杂化成键. (BN)6是唯一没有张力的单环结构.     

Heats of formation for Fe(CO)n(n=1–4)

The dissociation energies of Fe(CO)_n (n=2–4) are computed using correlation consistent basis sets and the CCSD(T) approach. The dissociation energies are extrapolated to the CBS limit and are corrected for core–valence (CV), scalar relativistic, spin–orbit, zero-point, and thermal effects. Our iron carbonyl bond strengths agree with experiment within the respective error bars. We use our dissociations energies at 298 K to obtain the heats of formation of Fe(CO)_n (n=1–4).     

Photoelectron imaging of tetrahydrofuran cluster anions (THF)(n) (-) (1≤n≤100)

Anionic tetrahydrofuran clusters (THF)(n) (-) (1≤n≤100) are studied with photoelectron imaging as gas-phase precursors for electrons solvated in THF. Photoelectron spectra of clusters up to n=5 show two peaks, one of which is attributed to a solvated open chain radical anion and the other to the closed THF ring. At n=6, the spectra change shape abruptly, which become more characteristic of (THF)(n) (-) clusters containing solvated electrons. From n=6-100, the vertical detachment energies (VDEs) of these solvated electron clusters increase from 1.96 to 2.71 eV, scaling linearly with n(-1/3). For fully deuterated (THF-d8)(n) (-) clusters, the apparent transition to a solvated electron cluster is delayed to n=11. Extrapolation of the VDEs to infinite cluster size yields a value of 3.10 eV for the bulk photoelectric threshold. The relatively large VDEs at onset and small stabilization with increasing cluster size compared to other solvated electron clusters may reflect the tendency of the bulk solvent to form preexisting voids that can readily solvate a free electron.     

A density functional study for the isomers of anionic sulfur clusters Sn− (n=3–20)

The signals of anionic sulfur clusters are intense in the mass spectrum of sulfur clusters generated in direct laser vaporization. We have acquired numerous isomers of sulfur clusters by means of the B3LYP DFT method. According to total energies, the most stable S_n⁻ (n=3–13) isomers are predicted. The helical S_n (n=14–20) structures are also calculated. Most of the anionic clusters are with chain configurations; the ring clusters with threefold atom(s) are higher in total energy. The most stable forms of isomers, from S₉⁻ to S₁₃⁻, show helical configurations that are completely different from those of the corresponding neutral and cationic clusters.     

Augmentation of the Generalised n×n Eigenvalue Equation to a Generalised (n+1)×(n+1) Eigenvalue Equation

Generalised n×n eigenvalue equation B| _i = _i S ^b | _i (i=1,...,n) where B and S ^b are n×n Hermitian matrices while S ^b is in addition positive definite is considered. This equation is augmented to a generalised (n+1)(n+1) eigenvalue equation H| _k = _k S| _k (k=1,...,n+1) where Hermitian matrices H and S represent matrices B and S ^b , respectively, augmented by one additional row and one additional column. It is shown how the eigenvalues _k and the eigenvectors | _k of the augmented eigenvalue equation can be expressed in terms of the eigenvalues _i and the eigenvectors | _i of the original eigenvalue equation. Operation count to obtain by this method all augmented eigenvalues and eigenvectors is of the order O(n ²). Unless matrices involved are of some special kind such as sparse matrices or alike, this operation count is one order of magnitude smaller than operation count required by other presently known methods. In many practical cases operation count to obtain a single selected eigenvalue and/or eigenvector by this method is of the order O(n). In the case of the generalised eigenvalue equation, all other methods usually require again O(n ³) operations, even if only a single eigenvalue and/or eigenvector is required. Thus in many cases of interest operation count to obtain a selected eigenvalue and/or eigenvector by this method is two orders of magnitude smaller than operation count required by other methods.     

(Cl2AlNH2)n和(H2AlNH2)n(n=1～5)簇结构及其热力学性质

用HF自洽场理论和密度泛函理论(DFT)的B3LYP方法,在6 31G水平上研究了低聚物(Cl2AlNH2)n和(H2AlNH2)n(n=1～5)簇的几何构型、电子结构和聚合反应热力学性质,比较了两个系列化合物中化学键的强度.结果表明,Cl2AlNH2和H2AlNH2分子为C2 (EC)平面型结构,其中Al－N为由一个σ键和一个键组成的双键.(Cl2AlNH2)n和(H2AlNH2)n(n=1～5)分子为Dnh对称,Al－N是典型的σ单键 .低聚物(Cl2AlNH2)n和(H2AlNH2)n的稳定性顺序分别为: 3 > 2 > 4> 5 > 1和8 > 7 > 9 > 11 > 6.     

羰基硼化合物(BCO)n(n==1～12)的理论研究

运用HF/3-21G方法和密度泛函理论(DFT)的B3LYP/6-31G*方法, 对羰基硼化合物(BCO)n(n=1~12)的各种可能结构进行了优化, 对在B3LYP/6-31G*水平上得到的几何构型、电子态、结合能、振动频率、核独立化学位移(NICS)、能量二次差分和热力学性质进行了理论研究, 得到了(BCO)n(n=1~12)结构的稳定性信息. 十二种基态结构都是端配位(μ1-CO), (BCO)n(n=1~3, 5, 6) 的基态是线型或平面结构, (BCO)n(n=4, 7~12)的基态是笼状结构; B—C平均键能呈现奇偶交替现象, 偶数的结构比奇数稳定; 能量二次差分得到同样的结论；羰基的振动频率与实验值非常吻合; 热力学性质的研究对实验具有重要的指导意义.