Density functional theory study of small vanadium oxide clusters

Density functional theory is employed to study structure and stability of small neutral vanadium oxide clusters in the gas phase. BPW91/LANL2DZ level of theory is used to obtain structures of VOy (y=1-5), V2Oy (y=2-7), V3Oy (y=4-9), and V4Oy (y=7-12) clusters. Enthalpies of growth and fragmentation reactions of the lowest energy isomers of vanadium oxide molecules are also obtained to study the stability of neutral vanadium oxide species under oxygen saturated gas-phase conditions. Our results suggest that cyclic and cage-like structures are preferred for the lowest energy isomers of neutral vanadium oxide clusters, and oxygen-oxygen bonds are present for oxygen-rich clusters. Clusters with an odd number of vanadium atoms tend to have low spin ground states, while clusters with even number of vanadium atoms have a variety of spin multiplicities for their ground electronic state. VO2, V2O5, V3O7, and V4O10 are predicted to be the most stable neutral clusters under the oxygen saturated conditions. These results are in agreement with and complement previous gas-phase experimental studies of neutral vanadium oxide clusters.     

Oxidation pattern of small silicon oxide clusters: structures and stability of Si6On (n = 1-12)

We have performed systematic ab initio calculations to study the structures and stability of Si(6)O(n)() clusters (n = 1-12) in order to understand the oxidation process in silicon systems. Our calculation results show that oxidation pattern of the small silicon cluster, with continuous addition of O atoms, extends from one side to the entire Si cluster. Si atoms are found to be separated from the pure Si cluster one-by-one by insertion of oxygen into the Si-O bonds. From fragmentation energy analyses, it is found that the Si-rich clusters usually dissociate into a smaller pure Si clusters (Si(5), Si(4), Si(3), or Si(2)), plus oxide fragments such as SiO, Si(2)O(2), Si(3)O(3), Si(3)O(4), and Si(4)O(5). We have also studied the structures of the ionic Si(6)O(n)(+/-) (n = 1-12) clusters and found that most of ionic clusters have different lowest-energy structures in comparison with the neutral clusters. Our calculation results suggest that transformation Si(6)O(n)+(a) + O --> Si(6)O(n+1)+(a) should be easier.     

Si_2C_(m-2)(m=4～15)团簇的结构与稳定性

用密度泛函理论(DFT)的B3LYP/6-31G*方法, 对Si2Cm-2(m = 4~15)团簇的几何构型、电子结构、振动频率等性质进行了研究, 讨论了化学键的特征和热力学稳定性, 比较了Si2Cm-2团中环状和线状结构的差异。结果表明, m = 4~13的团簇为线状结构, m=14~15的团簇为环状结构。在线状结构中, 随着m增大, 自旋多重度出现1、3交替变化, 并且Si原子倾向于在C链端部成键;环状结构中, C原子形成环状, 2个Si原子处于椭圆环状构型的两端。m 为奇数的Si2Cm-2团簇比m为偶数的更为稳定。     

Long, directional interactions in cofacial silicon phthalocyanine oligomers

Single crystal structures have been determined for the three cofacial, oxygen-bridged, silicon phthalocyanine oligomers, [((CH(3))(3)SiO)(2)(CH(3))SiO](SiPcO)(2-4)[Si(CH(3))(OSi(CH(3))(3))(2)], and for the corresponding monomer. The data for the oligomers give structural parameters for a matching set of three cofacial, oxygen-bridged silicon phthalocyanine oligomers for the first time. The staggering angles between the six adjacent cofacial ring pairs in the three oligomers are not in a random distribution nor in a cluster at the intuitively expected angle of 45° but rather are in two clusters, one at an angle of 15° and the other at an angle of 41°. These two clusters lead to the conclusion that long, directional interactions (LDI) exist between the adjacent ring pairs. An understanding of these interactions is provided by atoms-in-molecules (AIM) and reduced-density-gradient (RDG) studies. A survey of the staggering angles in other single-atom-bridged, cofacial phthalocyanine oligomers provides further evidence for the existence of LDI between cofacial phthalocyanine ring pairs in single-atom-bridged phthalocyanine oligomers.     

Structures and bonding of B4O5 and B4O5− clusters: Emergence of boroxol ring and competition between rhombic B2O2 and hexagonal B3O3 cores

Boron oxide clusters are electron-deficient species with novel structures and bonding, in which the emergence of rhombic and boroxol rings is of interest. We report on computational prediction of the global-minimum structures for two boron oxide clusters: B₄O₅ and B₄O₅⁻. These structures differ distinctly, as established through global machine searches and electronic structure calculations at B3LYP and single-point CCSD(T) levels. While B₄O₅ neutral cluster has a rhombic B₂O₂ core, the B₄O₅⁻ anion features a boroxol B₃O₃ ring. One electron completely changes the potential landscapes. Bonding analyses show that the 4π electron-counting is crucial for a rhombic B O cluster, in contrast to π sextet for a boroxol ring, which underlies the competition between rhombic and boroxol rings in B₄O₅/B₄O₅⁻ clusters. A possible pathway for rhombic-to-hexagonal transformation is proposed based on intrinsic reaction coordinate calculations. Anion B₄O₅⁻ cluster, a new member of the inorganic benzene family, is among the smallest B O species with a free-standing boroxol ring, governed collectively by composition and electron-counting.     

Structure and stability of oxygen adsorption on Si(n) (n = 5-10) clusters

The structures, binding energies, and electronic properties of one oxygen atom (O) and two oxygen atoms (2O) adsorption on silicon clusters Si(n) with n ranging from 5 to 10 are studied systematically by ab initio calculations. Twelve stable structures are obtained, two of which are in agreement with those reported in previous literature and the others are new structures that have not been proposed before. Further investigations on the fragmentations of Si(n)O and Si(n)O2 (n = 5-10) clusters indicate that the pathways Si(n)O --> Si(n-1) + SiO and Si(n)O2 --> Si(n-2) + Si2O2 are most favorable from thermodynamic viewpoint. Among the studied silicon oxide clusters, Si8O, Si9O, Si5O2 and Si8O2 correspond to large adsorption energies of silicon clusters with respect to O or 2O, while Si8O, with the smallest dissociation energy, has a tendency to separate into Si7 + SiO. Using the recently developed quasi-atomic minimal-basis-orbital method, we have also calculated the unsaturated valences of the neutral Si(n) clusters. Our calculation results show that the Si atoms which have the largest unsaturated valences are more attractive to O atom. Placing O atom right around the Si atoms with the largest unsaturated valences usually leads to stable structures of the silicon oxide clusters.     

Multiple aromaticity and antiaromaticity in silicon clusters.

A series of silicon clusters containing four atoms but with different charge states (Si4(2+), Si4, Si4(2-), and NaSi4-) were studied by photoelectron spectroscopy and ab initio calculations. Structure evolution and chemical bonding in this series were interpreted in terms of aromaticity and antiaromaticity, which allowed the prediction of how structures of the four-atom silicon clusters change upon addition or removal of two electrons. It is shown that Si4(2+) is square-planar, analogous to the recently discovered aromatic Al4(2-) cluster. Upon of two electrons, neutral Si4 becomes sigma-antiaromatic and exhibits a rhombus distortion. Adding two more electrons to Si4 leads to two energetically close structures of Si4(2-): either a double antiaromatic parallelogram structure or an aromatic system with a butterfly distortion. Because of the electronic instability of doubly charged Si4(2-), a stabilizing cation (Na+) was used to produce Si4(2-) in the gas phase in the form of Na+[Si4(2-)], which was characterized experimentally by photoelectron spectroscopy. Multiple antiaromaticity in the parallelogram Na+[Si4(2-)] species is highly unusual.     

Molecular and electronic structures, Brönsted basicities, and Lewis acidities of group VIB transition metal oxide clusters

The structures and properties of transition metal oxide (TMO) clusters of the group VIB metals, (MO(3))(n) (M = Cr, Mo, W; n = 1-6), have been studied with density functional theory (DFT) methods. Geometry optimizations and frequency calculations were carried out at the local and nonlocal DFT levels with polarized valence double-zeta quality basis sets, and final energies were calculated at nonlocal DFT levels with polarized valence triple-zeta quality basis sets at the local and nonlocal DFT geometries. Effective core potentials were used to treat the transition metal atoms. Two types of clusters were investigated, the ring and the chain, with the ring being lower in energy. Large ring structures (n > 3) were shown to be fluxional in their out of plane deformations. Long chain structures (n > 3) of (CrO(3))(n) were predicted to be weakly bound complexes of the smaller clusters at the nonlocal DFT levels. For M(6)O(18), two additional isomers were also studied, the cage and the inverted cage. The relative stability of the different conformations of M(6)O(18) depends on the transition metal as well as the level of theory. Normalized and differential clustering energies of the ring structures were calculated and were shown to vary with respect to the cluster size. Br?nsted basicities and Lewis acidities based on a fluoride affinity scale were also calculated. The Br?nsted basicities as well as the Lewis acidities depend on the size of the cluster and the site to which the proton or the fluoride anion binds. These clusters are fairly weak Br?nsted bases with gas phase basicities comparable to those of H(2)O and NH(3). The clusters are, however, very strong Lewis acids and many of them are stronger than strong Lewis acids such as SbF(5). Br?nsted acidities of M(6)O(19)H(2) and M(6)O(18)FH were calculated for M = Mo and W and these compounds were shown to be very strong acids in the gas phase. The acid/base properties of these TMO clusters are expected to play important roles in their catalytic activities.     

Reactivity of sub 1 nm supported clusters: (TiO2)n clusters supported on rutile TiO2 (110)

Metal oxide clusters of sub-nm dimensions dispersed on a metal oxide support are an important class of catalytic materials for a number of key chemical reactions, showing enhanced reactivity over the corresponding bulk oxide. In this paper we present the results of a density functional theory study of small sub-nm TiO(2) clusters, Ti(2)O(4), Ti(3)O(6) and Ti(4)O(8) supported on the rutile (110) surface. We find that all three clusters adsorb strongly with adsorption energies ranging from -3 eV to -4.5 eV. The more stable adsorption structures show a larger number of new Ti-O bonds formed between the cluster and the surface. These new bonds increase the coordination of cluster Ti and O as well as surface oxygen, so that each has more neighbours. The electronic structure shows that the top of the valence band is made up of cluster derived states, while the conduction band is made up of Ti 3d states from the surface, resulting in a reduction of the effective band gap and spatial separation of electrons and holes after photon absorption, which shows their potential utility in photocatalysis. To examine reactivity, we study the formation of oxygen vacancies in the cluster-support system. The most stable oxygen vacancy sites on the cluster show formation energies that are significantly lower than in bulk TiO(2), demonstrating the usefulness of this composite system for redox catalysis.     

M4 @Si28 (M=Al,Ga): metal-encapsulated tetrahedral silicon fullerene

It is known that silicon fullerenes cannot maintain perfect cage structures like carbon fullerenes. Previous density-functional theory calculations have shown that even with encapsulated species, nearly all endohedral silicon fullerenes exhibit highly puckered cage structures in comparison with their carbon counterparts. In this work, we present theoretical evidences that the tetrahedral fullerene cage Si(28) can be fully stabilized by encapsulating a tetrahedral metallic cluster (Al(4) or Ga(4)). To our knowledge, this is the first predicted endohedral silicon fullerene that can retain perfectly the same cage structure (without puckering) as the carbon fullerene counterpart (T(d)-C(28) fullerene). Density-functional theory calculations also suggest that the two endohedral metallosilicon fullerenes T(d)-M(4)@Si(28) (M=Al and Ga) can be chemically stable because both clusters have a large highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap ( approximately 0.9 eV), strong spherical aromaticity (nucleus-independent chemical shift value of -36 and -44), and large binding and embedding energies.     

Silica nanoarchitectures with tailored pores based on the hybrid three- and four-membered rings

Inspired by the recent developments in the controlled synthesis of porous materials, we present herein the structural prediction of silica nanoarchitectures by using the three- (3MRs) and four-membered rings (4MRs), which are more frequently found in the nanometer-sized particles than in the bulk form, as building blocks. The proposed models include the active molecular rings, thin nanowires, hollow nanotubes, discrete fullerene-like cages, and porous zeolite-like three-dimensional networks. Their geometrical and electronic structures and properties were studied by performing density functional calculations. These silica nanostructures were proved, using molecular dynamics simulations, to possess intrinsic structural stabilities with highly symmetrical geometries and regular nanochannels. These atomically well-defined clusters, (SiO)(n), are chemically more reactive than those proposed earlier and are energetically more favorable for n > 20 in high-level density functional calculations over the corresponding two-membered ring (2MR) chains and rings as well as the pure 3MR networks. The nanoparticles and nanodevices based on them are expected to have potential technological applications that mainly make use of their characteristic geometrical structures (nanosized pores) and novel electronic properties.     

Probing the electronic structure of early transition-metal oxide clusters: polyhedral cages of (V2O5)n(-) (n = 2-4) and (M2O5(2)(-) (M = Nb, Ta)

Vanadium oxide clusters, (V2O5)n, have been predicted to possess interesting polyhedral cage structures, which may serve as ideal molecular models for oxide surfaces and catalysts. Here we examine the electronic properties of these oxide clusters via anion photoelectron spectroscopy for (V2O5)n(-) (n = 2-4), as well as for the 4d/5d species, Nb4O10(-) and Ta4O10(-). Well-resolved photoelectron spectra have been obtained at 193 and 157 nm and used to compare with density functional calculations. Very high electron affinities and large HOMO-LUMO gaps are observed for all the (V2O5)n clusters. The HOMO-LUMO gaps of (V2O5)n, all exceeding that of the band gap of the bulk oxide, are found to increase with cluster size from n = 2-4. For the M4O10 clusters, we find that the Nb/Ta species yield similar spectra, both possessing lower electron affinities and larger HOMO-LUMO gaps relative to V4O10. The structures of the anionic and neutral clusters are optimized; the calculated electron binding energies and excitation spectra for the global minimum cage structures are in good agreement with the experiment. Evidence is also observed for the predicted trend of electron delocalization versus localization in the (V2O5)n(-) clusters. Further insights are provided pertaining to the potential chemical reactivities of the oxide clusters and properties of the bulk oxides.     

Structural model of silica nanowire assembled from a highly stable (SiO2)8 unit

The ground-state structures of silica clusters (SiO2)n for n = 1-8 were studied by performing calculations at the B3LYP/6-311+G(d) level of density functional theory. The results indicate that the growth mode of a silica nanowire based on small silica clusters may change at different wire lengths. A linear chain might be assembled from the smallest clusters of rhombic two-membered ring (2MR) with n < or = 5, while the growth motif changes at n = 6 into a more compact form composed of three-membered-rings (3MRs). The 3MR-containing structures become energetically favorable configurations for even longer silica clusters. In particular, the closed molecular ring consisting of 3MRs at n = 8 (i.e., (SiO2)8) with a high symmetry shows extreme energetic stability and relatively high chemical reactivity and thus is considered to be an important building block to assemble into silica nanowires. The relative stability of so-assembled silica nanowires were evaluated and compared with the models of silica nanowires in the literature.     

多体展开势能函数在碳族元素原子簇研究中的应用(Ⅳ)──Si9─Si16原子簇的分子结构与稳定性

多体展开势能函数研究表明,Si₄-Si₁₆原子簇分子间的结构衍生关系为:依次增加一个二配位或三配位的表面原子,分子表面被四元蝶形环Si₄(D_2d)所覆盖;Si_n(n=5-16)结构中多含有Si₅(D_3h)、Si₆(D_2d)区域结构单元,笼状Si₁₀及Si₁₆的表面原子均为三配位或三配位以上,预计Si₅、Si₆、Si₁₀及Si₁₆是硅原子簇碎片化产物分布中丰度较高的序列;在这一范围内的分子结构呈与晶体硅结构(金刚石)无关的密堆积,最高配位数为5,在小于半球的立体角内形成六配位的硅中心,使该簇合物在能量上不稳定。     

Theoretical study on (Al2O3)n (n = 1-10 and 30) fullerenes and H2 adsorption properties

The structures and stabilities of (Al2O3)n (n = 1-10 and 30) clusters were studied by means of first principles calculations. The calculated results reveal that the global minima of small (Al2O3)n (n = 1-5) clusters are cage structures with high symmetries, in which Al and O atoms are three- and two-coordinated, respectively, and are linked to neighbors via single bonds. Beyond (Al2O3)5, we calculated both cage and cage-dimer structures for (Al2O3)n (n = 6-10), and the results show that, at this size range, cage-dimer structures are more stable than cage structures. Furthermore, an onion-like motif for (Al2O3)10 was studied, and it is interesting to find that, at this size, the onion structure is more favorable than cage and cage-dimer structures. For large clusters, a shell structure of Al60O90 is suggested. Electronic properties and calculations on hydrogen adsorption of these aluminum oxide structures are reported, and we discuss their possible use as hydrogen storage materials.     

Structure and stability of small boron and boron oxide clusters

To rationally design and explore a potential energy source based on the highly exothermic oxidation of boron, density functional theory (DFT) was used to characterize small boron clusters with 0-3 oxygen atoms and a total of up to ten atoms. The structures, vibrational frequencies, and stabilities were calculated for each of these clusters. A quantum molecular dynamics procedure was used to locate the global minimum for each species, which proved to be crucial given the unintuitive structure of many of the most stable isomers. Additionally, due to the plane-wave, periodic DFT code used in this study a straightforward comparison of these clusters to the bulk boron and B2O3 structures was possible despite the great structural and energetic differences between the two forms. Through evaluation of previous computational and experimental work, the relevant low-energy structures of all but one of the pure boron clusters can be assigned with great certainty. Nearly all of the boron oxide clusters are described here for the first time, but there are strong indications that the DFT procedure chosen is particularly well suited for the task. Insight into the trends in boron and boron oxide cluster stabilities, as well as the ultimate limits of growth for each, are also provided. The work reported herein provides crucial information towards understanding the oxidation of boron at a molecular level.     

硅氧团簇(SiO2)nO2H4的密度泛函理论研究

提出硅氧团簇(SiO2)nO2H4的两种新构型: 基于笼状结构和环状结构的构型, 并与链状构型相比较, 用密度泛函理论的B3LYP方法在6-31G(d)基组水平上计算了三种构型n=2～22(n取偶数)的几何结构、平均结合能、能隙以及能量的二次差分. 分析计算结果发现, 笼状构型不但在n=4和8处存在幻数团簇(实验上已经观察到), 而且预测在n=14处也存在类似的幻数团簇; 此外, 与(SiO2)n团簇不同的是, (SiO2)nO2H4团簇的环状构型的稳定性从n=4开始大于链状构型, 意味着水的加成对硅氧团簇的稳定性有着重要的影响.     

Stability of small Pdn (n=1-7) clusters on the basis of structural and electronic properties: a density functional approach

Density functional calculations within the generalized gradient approximation have been used to investigate the lowest energy electronic and geometric structures of neutral, cationic, and anionic Pd(n) (n=1-7) clusters in the gas phase. In this study, we have examined three different spin multiplicities (M=1, 3, and 5) for different possible structural isomers of each neutral cluster. The calculated lowest energy structures of the neutral clusters are found to have multiplicities, M=1 for Pd(1), Pd(3), Pd(5), Pd(6), and Pd(7), while M=3 for Pd(2) and Pd(4). We have also determined the lowest energy states of cationic and anionic Pd(n) (n=1-7) clusters, formed from the most stable neutral clusters, in three spin multiplicities (M=2, 4, and 6). Bond length, coordination number, binding energy, fragmentation energy, bond dissociation energy, ionization potential, electron affinity, chemical hardness, and electric dipole moment of the optimized clusters are compared with experimental and other theoretical results available in the literature. Based on these criteria, we predict the four-atom palladium cluster to be a magic-number cluster.